1957-2009年烏魯木齊河徑流與氣候變化的對應(yīng)關(guān)系

穆艾塔爾·賽地, 熱孜燕, 阿不都·沙拉木

(1.新疆水利水電科學(xué)研究院, 烏魯木齊 830049; 2.鞏留縣水利局, 新疆 鞏留 835400)

?

1957-2009年烏魯木齊河徑流與氣候變化的對應(yīng)關(guān)系

穆艾塔爾·賽地1, 熱孜燕2, 阿不都·沙拉木1

(1.新疆水利水電科學(xué)研究院, 烏魯木齊 830049; 2.鞏留縣水利局, 新疆 鞏留 835400)

以1957—2009年烏魯木齊河徑流量、氣溫和降水量數(shù)據(jù)為基礎(chǔ)，采用線性分析、Mann-Kendall檢驗(yàn)方法、相關(guān)性分析和灰色關(guān)聯(lián)度分析法等方法，研究烏魯木齊河年徑流量、氣溫和降水量的變化特征以及影響徑流量的氣象因子。結(jié)果表明：(1) 50多年間烏魯木齊河年徑流量變化呈增加趨勢，氣溫顯著升高而降水量緩慢減少。從徑流量的年際變化來看，烏魯木齊河年徑流量變化出現(xiàn)平水期(1957—1990年)、豐水期(1991—2000年)和平水期(2001—2009年)，總體呈現(xiàn)波動(dòng)增加趨勢；從年內(nèi)分配來看，烏魯木齊河年徑流量主要集中在夏季和秋季，而春季和冬季較少；從月分布來看，6月、7月、8月份徑流量最大，其他月份較少。(2) Mann-Kendall檢驗(yàn)表明，1957—2009年烏魯木齊河年徑流量、春季徑流量和冬季徑流量均顯著上升，而夏季和秋季徑流量上升趨勢不顯著。(3) 相關(guān)性分析和灰色關(guān)聯(lián)度分析表明，烏魯木齊河年徑流量變化主要受降水的影響。

烏魯木齊河; 年徑流量; Mann-Kendall檢驗(yàn); 相關(guān)性分析; 灰色關(guān)聯(lián)度分析

氣候變化直接影響水資源徑流的大小及空間分布，是徑流變化的主要驅(qū)動(dòng)因素[1-2]，而氣溫和降水變化是氣候變化的主要體現(xiàn)。氣溫變化可以引起流域潛在蒸發(fā)和實(shí)際蒸發(fā)變化，降水量年內(nèi)與年際變化在下墊面引起流域水分垂向和橫向再分布，由此會影響流域徑流情勢變化[3]。本文以烏魯木齊河上游出山口英雄橋水文站1957—2009年年徑流量、氣溫和降水量等數(shù)據(jù)為基礎(chǔ)，研究近50多年烏魯木齊河年徑流量變化趨勢及其與氣溫、降水等氣象因子的關(guān)系，以期為了解全球及區(qū)域氣候變化背景下的干旱區(qū)內(nèi)陸河流域水資源變化提供參考。

1 研究區(qū)概況

烏魯木齊河是新疆干旱區(qū)典型的內(nèi)陸河，位于天山山脈北坡中段，補(bǔ)給以冰雪融水、降雨及地下水混合補(bǔ)給為主，是新疆維吾爾自治區(qū)烏魯木齊市工、農(nóng)業(yè)生產(chǎn)和城市生活用水的主要水源[5]。烏魯木齊河多年平均徑流量為2.44億m3，年最大徑流量為3.44億m3。烏魯木齊河由烏魯木齊河水系組成，包括發(fā)源于南山的烏魯木齊河及水西溝等10條支流以及發(fā)源于東山的水磨溝、蘆草溝、鐵廠溝、白楊溝等15條小河和山泉。烏魯木齊河徑流量在出山口附近達(dá)到最大值，在中下游逐漸散失，并最終消失在米泉市西北部的古爾班通古特沙漠中。烏魯木齊河流域是新疆水利水電建設(shè)中水資源開發(fā)利用程度較高、工程效益發(fā)揮較好的地區(qū)，河流上游建有青年渠、知青渠、黃草梁渠以及大西溝水庫等；中游攔河筑壩建有烏拉泊水庫和紅雁池水庫；下游平原地區(qū)修建了猛進(jìn)水庫、八一水庫。流域水資源利用和開發(fā)程度高，水資源變率對城市生產(chǎn)、生活影響極大[4-6]。

2 數(shù)據(jù)與方法

2.1 數(shù) 據(jù)

利用烏魯木齊河出山口英雄橋水文站1957—2009年實(shí)測年徑流量、降水量和氣溫等數(shù)據(jù)，采用線性回歸分析、相關(guān)分析、Mann-Kendall檢驗(yàn)法和灰色關(guān)聯(lián)度分析等方法定量研究烏魯木齊河年徑流量、氣溫與降水量變化趨勢以及徑流量與氣溫、降水量變化的相關(guān)性，分析烏魯木齊河年徑流量變化的原因及與氣候因子的相關(guān)性。

2.2 灰色關(guān)聯(lián)度分析

灰色理論方法視不確定的信息為灰數(shù)，應(yīng)用灰色數(shù)學(xué)處理不確定量，可以最終使不確定量量化[7-8]。

假設(shè)有n個(gè)時(shí)間序列：

(1)

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

Mann-Kendall趨勢檢驗(yàn)法是一種非參數(shù)統(tǒng)計(jì)檢驗(yàn)方法，該方法不需要樣本遵從一定的分布，也不受少數(shù)異常值的干擾，計(jì)算比較簡單，是目前較為常用的氣候、氣象趨勢判斷方法[9-12]。

假設(shè)H0為時(shí)間序列數(shù)據(jù)(X1,…,Xn)，是n個(gè)獨(dú)立的、隨機(jī)變量同分布的樣本；備擇假設(shè)H1，是雙邊檢驗(yàn)，對于所有的k，j

(2)

統(tǒng)計(jì)變量S的取值由符號函數(shù)Sgn(Xj-Xk)決定，其中：

(3)

S服從正態(tài)分布，其均值為0，方差var(S)計(jì)算公式為：

(4)

當(dāng)n≥10時(shí)，Z收斂于標(biāo)準(zhǔn)正態(tài)分布，正態(tài)分布檢驗(yàn)統(tǒng)計(jì)量Z計(jì)算公式為：

(5)

Mann-Kendall趨勢檢驗(yàn)法可以進(jìn)行序列變化的突變檢驗(yàn)。對于時(shí)間序列x1,x2,…,xn定義統(tǒng)計(jì)量Tj，表明第j個(gè)樣本xj大于xi的累計(jì)數(shù)，其中j=1,2,…,n；i=1,2,…,j-1。

(6)

在原序列隨機(jī)獨(dú)立的假定下，定義統(tǒng)計(jì)變量U(t)為：

(7)

其中，E(Tj)，var(Tj)分別為累計(jì)統(tǒng)計(jì)量Tj的均值和方差，其計(jì)算為：

E(Tj)=[j(j-1)]/4

(8)

var(Tj)=[j(j-1)(2j+5)]/72

(9)

3 結(jié)果與分析

3.1 徑流量與氣象因子的變化

運(yùn)用距平分析法和時(shí)間序列法分析研究烏魯木齊河上游氣象要素的時(shí)程變化特征。由圖1可以看出，1957年以來烏魯木齊河年徑流量呈增加趨勢。其中：1957—1990年平均年徑流量為2.35億m3，接近于多年平均水平，總體呈波動(dòng)增加趨勢；1991—2000年為豐水期，徑流量大于多年平均徑流量，在1996年徑流量達(dá)到最大正距平值，為0.42；2001—2009年徑流量存在一定的波動(dòng)，在2001年徑流量達(dá)到最大負(fù)距平值，為-0.28，總體上接近于多年平均徑流量，相對于前一段而言是流域的平水期。如圖2所示，從徑流量的季節(jié)性分布來看，春(3—5月)、夏(6—8月)、秋(9—11月)和冬(12—次年2月)4個(gè)季節(jié)分別占年總徑流量的10.28%，69.37%，15.74%，4.73%。從徑流量的月平均值來看，6月、7月、8月份徑流量最大，分別占19.3%，28.6%，21.4%。

圖1 1957-2009年烏魯木齊河年徑流量、氣溫和降水量變化

運(yùn)用時(shí)間序列分析法分析烏魯木齊河上游1957—2009年氣溫和降水量的變化規(guī)律(圖1)。結(jié)果表明，烏魯木齊河上游出山口氣溫變化總體呈上升趨勢，線性傾向率為0.018℃/a，相當(dāng)于每10a平均氣溫升高0.18℃。根據(jù)1957—2009年的氣溫變化總體情況，將氣溫變化趨勢可以分成4個(gè)時(shí)段，分別為1957—1979年、1980—1983年、1984—2004年、2005—2009年。其中，1957—1979年和1984—2004年兩個(gè)時(shí)期氣溫均呈顯著上升趨勢，而1980—1983年和2005—2009年氣溫變化總體上呈上升趨勢，但個(gè)別年份氣溫發(fā)生明顯的下降(<-1℃)，其中在1979年、1982年、1983年，2005年出現(xiàn)氣溫低值(0～1℃)；從1957—2009年烏魯木齊河上游降水量歷年變化可以看出，降水量在整個(gè)研究時(shí)段總體上呈緩慢降低的趨勢，線性傾向率為1.36mm/a。

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

利用Mann-Kendall檢驗(yàn)法計(jì)算1957—2009年序列實(shí)測徑流量資料，在0.05的置信水平上獲得烏魯木齊河年徑流量和四季徑流量的變化趨勢特征，計(jì)算結(jié)果如表1所示。

計(jì)算結(jié)果表明，在0.05的置信水平上，烏魯木齊河上游年徑流量、春季和冬季徑流量的Z值分別為2.11，2.45，4.41，均大于檢驗(yàn)臨界值1.89，表明在1957—2009年烏魯木齊河上游年徑流量、春季徑流量和冬季徑流量均呈顯著的上升趨勢；與此對比，烏魯木齊河上游夏季和秋季徑流量的Z值分別為1.03，1.61，均小于檢驗(yàn)臨界值1.89，表明該時(shí)段徑流量的上升趨勢不顯著。

突變檢驗(yàn)分析表明(圖3)，在1957—2009年除1974年外，其他時(shí)間U(t)的統(tǒng)計(jì)變量取值均大于零，表明在研究期間烏魯木齊河上游年徑流量的變化呈顯著的上升趨勢。其中除了1963—1967年、1994—2005年U(t)的曲線超過臨界線外，在其他時(shí)間段U(t)和U(t)′的統(tǒng)計(jì)曲線均落在0.05置信水平上的檢驗(yàn)臨界線之間，表明烏魯木齊河上游年徑流量的變化趨勢不顯著。

計(jì)算可知，1986年以前烏魯木齊河上游年徑流量上升變化有較大的波動(dòng)，而1986年后其徑流量變化呈逐漸增長的趨勢。1957—2009年烏魯木齊河上游年徑流量的U(t)和U(t)′統(tǒng)計(jì)曲線落在兩條檢驗(yàn)臨界線之間的部分在1970年、1972年、1978年出現(xiàn)交點(diǎn)，代表其年徑流量變化的突變發(fā)生時(shí)間。

圖2 1957-2009年烏魯木齊河徑流量年內(nèi)變化特征

表1 烏魯木齊河上游徑流量變化趨勢檢驗(yàn)結(jié)果

圖3 烏魯木齊河上游出山口年徑流量Mann-Kendall檢驗(yàn)曲線

3.3 烏魯木齊河上游年徑流量與氣溫、降水的關(guān)系

通過徑流量和氣溫、降水量等氣象因子的相關(guān)分析可知，烏魯木齊河年徑流量與降水量的相關(guān)性較好(R2=0.857)，通過0.01置信水平上的顯著性檢驗(yàn)；而徑流量與氣溫的相關(guān)性不顯著(R2=0.457)，未通過0.01置信水平上的顯著性檢驗(yàn)，表明烏魯木齊河年徑流量變化受降水的影響較顯著，受氣溫的影響不顯著。為了進(jìn)一步驗(yàn)證以上相關(guān)分析結(jié)果，利用灰色關(guān)聯(lián)度分析方法分析烏魯木齊河上游地區(qū)年徑流量與降水和氣溫的關(guān)聯(lián)度。首先通過對1957—2009年烏魯木齊河上游年徑流量與氣溫、降水的原始數(shù)據(jù)進(jìn)行均值化轉(zhuǎn)換，采用如下公式：

(10)

原始數(shù)據(jù)標(biāo)準(zhǔn)化處理后，以年徑流量為參考序列X0，以氣象因子(氣溫和降水量)為比較序列X1，求得比較序列與參考序列的絕對差值，從而計(jì)算出烏魯木齊河流域上游地區(qū)年徑流量與氣象因子的關(guān)聯(lián)系數(shù)，比較其大小，確定年徑流量與氣象因子之間的相關(guān)性。比較時(shí)，以L01代表徑流量與降水量間的關(guān)聯(lián)系數(shù)，以L02代表徑流量與氣溫間的關(guān)聯(lián)系數(shù)。計(jì)算可得L01=0.851 7，L02=0.715 8；年徑流量與降水和氣溫的關(guān)聯(lián)度的大小順序?yàn)長01>L02，表明烏魯木齊河上游降水量對徑流量的影響比氣溫更顯著，該結(jié)論與相關(guān)性分析的結(jié)果一致。

表2 徑流量、氣溫和降水量的標(biāo)準(zhǔn)化值

注：Rt,Tt,Pt分別為徑流標(biāo)準(zhǔn)化值、氣溫標(biāo)準(zhǔn)化值、降水標(biāo)準(zhǔn)化值。

4 結(jié) 論

本文通過分析1957—2009年烏魯木齊河年徑流量變化及其與氣溫、降水量的關(guān)系得出如下結(jié)論：

(1) 從徑流量的年際變化來看，1957—1990年平均年徑流量接近于多年平均水平，1991—2000年徑流量大于多年平均徑流量，2001—2009年與多年平均徑流量基本一致，總體呈增加的趨勢，出現(xiàn)平水期(1957—1990年)、豐水期(1991—2000年)和平水期(2001—2009年)的變化過程。徑流量主要集中在夏季和秋季，春季和冬季徑流量比重較小。

(2) 烏魯木齊河上游徑流量的變化與降水量的變化趨勢一致；徑流量增長的同時(shí)氣溫顯著升高而降水量緩慢減少。

(3) Mann-Kendall檢驗(yàn)表明烏魯木齊河年徑流量、春季徑流量和冬季徑流量均呈顯著上升趨勢，而夏季和秋季徑流量上升趨勢不顯著。

(4) 相關(guān)分析表明，烏魯木齊河徑流量與降水量的相關(guān)性較好，而與氣溫的相關(guān)性不顯著。灰色關(guān)聯(lián)度法分析表明徑流量與降水量間關(guān)聯(lián)系數(shù)大于其與氣溫間關(guān)聯(lián)系數(shù)，表明降水量對徑流量影響更顯著。

[1] Dooge J C I. Sensitivity of runoff to climate change: A hortonian approach[J]. Bulletin of the American Meteorological Society,1992,73(12):2013-2024.

[2] Bergstrom S, Carlsson B, Gardelin M, et al. Climate change impacts on runoff in Sweden-assessments by global climate models, dynamical downscaling and hydrological modeling[J]. Climate Research,2001,16(2):101-112.

[3] 張士鋒,華東,孟秀敬,等.三江源氣候變化及其對徑流的驅(qū)動(dòng)分析[J].地理學(xué)報(bào),2011,66(1):13-24.

[4] 谷鑫,金兵.烏魯木齊河水文特性分析[J].黑龍江水利科技,2012,40(11):17-18.

[5] 姚海燕,張民.烏魯木齊河流域水管理信息系統(tǒng)[J].水文,2006,26(2):78-80,96.

[6] 穆艾塔爾·賽地,阿不都·沙拉木,崔春亮,等.新疆天山北坡山區(qū)流域水文特征分析[J].水文,2013,33(2):87-92.

[7] 秦昌波,鄭丙輝,秦延文,等.渤海灣天津段海岸帶水環(huán)境質(zhì)量灰色關(guān)聯(lián)度評價(jià)[J].環(huán)境科學(xué)研究,2006,19(6):94-99.

[8] 李香云,王玉杰.不同植被類型對地表徑流影響的灰色關(guān)聯(lián)度分析[J].水土保持通報(bào),2007,27(2):83-86.

[9] 曹輝,黃強(qiáng),暢建霞,等.黑河徑流時(shí)空演變規(guī)律分析[J].水資源與水工程學(xué)報(bào),2008,19(5):69-72.

[10] 符淙斌,王強(qiáng).氣候突變的定義和檢測方法[J].大氣科學(xué),1992,16(4):482-493.

[11] 陳忠升,陳亞寧,徐長春.近50 a來塔里木河干流年徑流量變化趨勢及預(yù)測[J].干旱區(qū)地理,2011,34(1):43-51.

[12] Gerstengarbe F W, Werner P C. Estimation of the beginning and end of recurrent events within a climate regime[J]. Climate Research,1999,11(2):97-107.

Corresponding Relationship Between the Runoff and Climate Change in the Upper Urumqi River Basin During the Period from 1957 to 2009

Muattar Saydi1, Reziyan2, Shalamu Abudu1

(1.XinjiangResearchInstituteofWaterResources,Urumqi830049,China;2.WaterConservancyBureauofGongliu,Gongliu,Xinjiang835400,China)

We analyzed the annual and interannual variability of runoff discharge and the related influential factors in Urumqi River basin during the hydrological years from 1957 to 2009 using linear analysis, Mann-Kendall test, correlation and grey correlation analysis methods based on the annual runoff, temperature, and precipitation data from the hydrological station located at the river outlet. Results showed that: (1) annual runoff had the increasing trend during the study period, and the annual average temperature increased significantly while the annual precipitation decreased slightly; the runoff had relatively greater value and was sufficient during the period from 1991 to 2000, and had average value during the hydrological years from 1957 to 1990 and from 2001 to 2009; in terms of seasonal distribution, the runoff in Urumqi River basin mainly concentrated in summer and fall seasons while runoff was relatively small in winter and spring; the runoff was the greatest in June, July and August within a year; (2) Mann-Kendall analysis showed that the annual, spring and winter runoff had significant increasing trend while the increasing trend for runoff in summer and fall seasons was insignificant; (3) the correlation and gray correlation analysis showed that the increased runoff in Urumqi River basin was mainly caused by precipitation during the study period.

Urumqi River; annual runoff; Mann-Kendall test; correlation analysis; grey correlation analysis

2014-09-23

2014-10-22

新疆維吾爾自治區(qū)科技計(jì)劃項(xiàng)目(20126013);國家自然科學(xué)資助項(xiàng)目(51069017);新疆大學(xué)綠洲生態(tài)重點(diǎn)實(shí)驗(yàn)室項(xiàng)目(XJDX0201-2013-04);新疆維吾爾自治區(qū)公益性科研院所基本科研業(yè)務(wù)經(jīng)費(fèi)自助項(xiàng)目(KY2014044)

穆艾塔爾·賽地(1984—),女,新疆策勒縣人,工程師,主要從事水文水資源研究。E-mail:895798730@qq.com

阿不都·沙拉木(1968—),男,新疆鞏留人,研究員,主要從事水文水資源研究。E-mail:shalamu3@gmail.com

P339

1005-3409(2015)05-0289-05