洪波　海米提?依米提　馬蓉等

摘要以開都河徑流出山口以上部分作為研究區(qū)域，基于Mann-Kendall以及雙累積曲線等方法對(duì)1960～2007年氣候變化和人類活動(dòng)對(duì)徑流影響的貢獻(xiàn)率進(jìn)行分離。結(jié)果表明，開都河流域氣候趨向暖濕方向發(fā)展，且徑流增加；以MK檢驗(yàn)得到的徑流量突變點(diǎn)，將研究時(shí)段劃為1960～1996和1997～2007年2個(gè)階段，后半段與前半段相比降水增加18.98 mm，溫度升高0.92 ℃，且徑流增加29.78 m3/s；氣候變化導(dǎo)致開都河徑流增加，但人類活動(dòng)則導(dǎo)致徑流減少，其貢獻(xiàn)率分別為122%和-22%；氣候變化是導(dǎo)致開都河出山口徑流量增加的主要驅(qū)動(dòng)力。

關(guān)鍵詞氣候變化；人類活動(dòng)；開都河；徑流；影響

中圖分類號(hào)S181文獻(xiàn)標(biāo)識(shí)碼A文章編號(hào)0517-6611（2015）16-201-04

Effects of Climate Change and Human Activities on Path of Kaidu River

HONG Bo1，2， HAIMITI Yimiti2，3*， MA Rong1，2 et al

（1. School of Resource and Environment Science， Xinjiang University， Urumqi， Xinjiang 830046； 2.Oasis Ecology Key Laboratory of Ministry of Education， Urumqi， Xinjiang 830046； 3. Key Laboratory of Lake Environment and Resource in Arid Region， Urumqi， Xinjiang 830054）

AbstractTaking water source area of Kaidu River as the study area， based on MannKendall and double cumulative curve methods， the contribution rate of climate change and human activities on path during 1960-2007 was separated. The results showed that in Kaidu River Basin climate became warm and wet， and the runoff increased； according to runoff mutation point， the research period was divided into two stages of 1960-1996 and 1997-2007； with the first periods comparison， the precipitation of the last period increased 18.98 mm， temperature increased 0.92 ℃， and the runoff increased 29.78 m3/s.Human activities bring about the runoff cut down， but climate change generate the runoff increase， the contribution rate was respectively 122% and -22%. Climate change is the main driving force causing runoff increase of Kaidu River.

Key wordsClimate change； Human activities；Kaidu River； Path； Influence

在全球氣候變暖的當(dāng)下，我國(guó)西北地區(qū)氣候發(fā)生了很大變化，逐漸由暖干向暖濕型轉(zhuǎn)變[1-3]。西北干旱區(qū)氣溫對(duì)氣候變暖的敏感性明顯強(qiáng)于東部季風(fēng)區(qū)，是氣候變化的敏感地帶[4]。對(duì)于西北干旱區(qū)而言，地表徑流對(duì)氣候變化的響應(yīng)極為敏感，由于空氣濕度較低，輕微的降水和溫度變化均會(huì)引起河流徑流量的明顯改變；與此同時(shí)人類活動(dòng)對(duì)水資源的重新分配同樣對(duì)水資源的變化有重大影響。我國(guó)最大的內(nèi)陸淡水湖是位于新疆庫(kù)爾勒市的博斯騰湖，湖區(qū)蒸發(fā)量大且降水量少。開都河是博斯騰湖的最大入湖河流，近年來(lái)，有的學(xué)者對(duì)開都河流域中上游的氣溫和降水變化進(jìn)行了一些研究[5-10]。但氣候變化以及人類活動(dòng)對(duì)開都河水源區(qū)流量影響的程度未見報(bào)道。筆者在此利用MannKendall和雙累積曲線等方法分離了1960～2007年氣候變化和人類活動(dòng)對(duì)開都河徑流影響的貢獻(xiàn)率。

1資料與方法

1.1研究區(qū)概況

開都河處于新疆維吾爾自治區(qū)的巴州（82°58′～86°55′E、41°47′～43°21′ N），全長(zhǎng)610 km，是雪融水和雨水補(bǔ)給為主的河流[11]。流域內(nèi)的冰川總數(shù)為722條，面積可達(dá)445 km2，主要分布于開都河水源地4 000～4 500 m地區(qū)。積雪、冰川和地下水補(bǔ)給在河流水源中占相當(dāng)大的比重，冰川融水補(bǔ)給占15%左右[12]。開都河水量年內(nèi)變化較為分散，呈單峰弧線，高水量多發(fā)生在高溫多雨的夏季，約占年總量的50%，而在寒冷且降水稀少的冬季，主要補(bǔ)給源是地下水回流。該研究的研究區(qū)域?yàn)榇笊娇谒恼疽陨蠀^(qū)域，區(qū)域面積約為18 311 km2（圖1）。

圖1研究區(qū)概況

1.2數(shù)據(jù)處理與研究方法

1.2.1氣象觀測(cè)數(shù)據(jù)。

開都河流域的氣象觀測(cè)數(shù)據(jù)來(lái)源于新疆氣象局，包括了巴音布魯克、和碩、和靜、焉耆以及大山口氣象站1960～2007年的逐日氣溫、降水量等氣象要素，其中氣溫和降水?dāng)?shù)據(jù)均求取算術(shù)平均值作為流域的平均值。

1.2.2水文氣象要素的MK突變檢驗(yàn)。采用MK檢驗(yàn)法以及距平百分比來(lái)研究水文氣象要素變化特征。距平百分比能夠輕易獲得各水文氣象要素的變化趨勢(shì)，并可對(duì)各要素間的相互關(guān)系進(jìn)行簡(jiǎn)要分析；MK檢驗(yàn)可以發(fā)現(xiàn)水文氣象要素的變化趨勢(shì)和突變年份，為雙累積曲線法提供臨界年份。在此利用統(tǒng)計(jì)軟件DPS對(duì)大山口水文站1960～2007年的觀測(cè)徑流數(shù)據(jù)進(jìn)行MK檢驗(yàn)。檢驗(yàn)的置信水平設(shè)為0.05。

1.2.3氣候變化和人類活動(dòng)對(duì)徑流影響的分離。

采用雙累積曲線法將氣候變化和人類活動(dòng)對(duì)開都河水量的影響進(jìn)行分離[13-14]。根據(jù)MK檢驗(yàn)獲得的突變分析結(jié)果，將實(shí)測(cè)時(shí)段分為突變點(diǎn)前段和突變點(diǎn)后段兩段（簡(jiǎn)稱為標(biāo)準(zhǔn)期和變化期）。

將標(biāo)準(zhǔn)期的累積徑流量R1以及累積降水量P1通過線性回歸得到以下方程：

R1=kP1+b（1）

將（1）中的方程式帶入到變化期（即包含氣候變化和人類活動(dòng)的影響成分的部分），以變化期的累積降水量P2得到累積徑流量R2，此時(shí)計(jì)算得到的累計(jì)徑流量R2可認(rèn)為是與標(biāo)準(zhǔn)期下墊面相同情況下的人類活動(dòng)較小時(shí)期的徑流。

用計(jì)算得到的R2可以得出年平均徑流，再分別算出標(biāo)準(zhǔn)期和變化期的實(shí)測(cè)和演算得到的年徑流量的均值，并對(duì)這些均值求差，即

δ1=R2實(shí)測(cè)-R1實(shí)測(cè)（2）

δ2=R2計(jì)算-R2實(shí)測(cè)（3）

式中，R1實(shí)測(cè)、R2實(shí)測(cè)分別表示標(biāo)準(zhǔn)期和變化期的實(shí)際測(cè)量平均徑流量，R2計(jì)算表示變化期由計(jì)算得到的平均徑流量，δ1表示氣候變化和人類活動(dòng)協(xié)同下的徑流變化，在（3）中消除了降水因素變化的影響，因此，δ2表示由非降水要素（如人類活動(dòng)）導(dǎo)致的徑流變化。計(jì)算其比值：

α=δ2/δ1（4）

α值則表示非降水要素變化在總變化中所占的比重。

通過（4）式獲取的非降水變化要素的影響大小，就可知降水要素對(duì)徑流的影響，前提是要判斷非降水和降水要素對(duì)徑流變化影響的趨勢(shì)。將這2種變化要素的綜合影響當(dāng)作是1，其中非降水要素對(duì)徑流的影響看作α，當(dāng)這2種要素導(dǎo)致徑流向相同方向變化時(shí)，即同時(shí)增加或同時(shí)減小，那么降水對(duì)徑流的影響為（1-|α|）。當(dāng)2種要素使徑流的大小向不同方向變化時(shí)，還需要看徑流量是向大變還是向小變化，當(dāng)非降水要素對(duì)徑流的影響趨勢(shì)增大（或減?。覐搅髁客瑫r(shí)增大（或減?。r(shí)，降水對(duì)徑流的影響為（|α|-1）；當(dāng)非降水要素對(duì)徑流的影響趨勢(shì)增大（或減?。?，而徑流卻減小（或增大）時(shí)，降水要素對(duì)徑流的影響是（1+|α|）。

2結(jié)果與分析

2.1水文氣象要素的變化分析

1960～2007年開都河流域平均降水、徑流和溫度分別是121.19 mm、110.09 m3/s、6.06 ℃，從極值比、標(biāo)準(zhǔn)偏差和變異系數(shù)（表1）來(lái)看，三者均有較

表11960～2007年開都河流域水文氣象要素的特征值

要素均值最大值最小值極值比標(biāo)準(zhǔn)偏差變異系數(shù)

降水∥mm121.19173.4467.382.625.760.21

徑流∥m3/s110.09180.2877.822.320.290.18

氣溫∥℃6.067.214.711.50.610.10

為緩和的變化。徑流的變異程度和降水的變異程度差不多，說(shuō)明研究區(qū)的徑流變化受降水因素影響程度較高。

從年際變化的線性趨勢(shì)（圖2）來(lái)看，降水和溫度均有升高的趨勢(shì)，說(shuō)明在未來(lái)一段時(shí)間內(nèi)開都河流域的氣候會(huì)向暖濕方向發(fā)展。降水的變化增加了河流水源的補(bǔ)給，而溫度的增高會(huì)使水源區(qū)的冰雪融水增加，因此，開都河流域的這種氣候變化導(dǎo)致了河流徑流的增加。

圖21960～2007年開都河流域徑流（a）、降水（b）和溫度（c）線性變化趨勢(shì)

2.2水文氣象要素的距平百分比變化

分析1960～2007年研究區(qū)水文氣象要素的距平百分比（圖3）可以看出，開都河流域的溫度變化相對(duì)平穩(wěn)，最大和最小變幅均為多年平均值的20%左右。降水變化較大，20世紀(jì)60～80年代中期以負(fù)距平為主，說(shuō)明這一時(shí)期開都河流域降水低于多年平均降水，可能會(huì)出現(xiàn)旱災(zāi)；80年代中期～21世紀(jì)初，距平值以正值為主，說(shuō)明這一時(shí)期的降水量高于多年平均水平，水量充足；正負(fù)距平的極值分別為1992年的43.12%和1985年的-44.40%，說(shuō)明1985～1992年降水變化明顯，可能是徑流量發(fā)生突變的主要原因之一。徑流變化規(guī)律與降水變化規(guī)律基本一致，一定程度上其變幅小于降水，

圖31960～2007年開都河流域年徑流、降水和氣溫的距平變化

正負(fù)距平的極值分別為2002年的63.75%和1986年的-29.31%?？傮w上，降水量和徑流量的變化方向是一致的，表明降水對(duì)徑流有決定性的影響。一般而言降水變化會(huì)被放大到徑流變化中[15]，但開都河流域降水的變化與徑流的變化相差不大，有些時(shí)段降水距平百分率的絕對(duì)值大于徑流，說(shuō)明非降水要素如水利工程建設(shè)、農(nóng)業(yè)用水等擾動(dòng)了降水與徑流的關(guān)系，且緩和了降水與徑流之間的關(guān)系，如部分降水低于平均值的年份（2001、2004和2006年）徑流卻高于平均值。

2.3水文氣象要素的突變分析

從MK突變性檢驗(yàn)結(jié)果（圖4）可看出，徑流自1985年以來(lái)不斷減少，在2000年左右開始增加，并在1997年發(fā)生突變性增加；在1979年以來(lái)，開都河流域的降水在總體上呈上升趨勢(shì)，并在1988年發(fā)生突變性升高；而溫度自1986年以來(lái)一直有升高趨勢(shì)，且在1998年發(fā)生升高的突變；從這3個(gè)要素的變化趨勢(shì)來(lái)看，徑流增加與降水增加和溫度升高的趨勢(shì)基本同步，說(shuō)明1960～2007年氣候變化必然會(huì)對(duì)開都河徑流產(chǎn)生重要影響。

圖41960～2007年開都河流域徑流（a）、降水（b）和氣溫（c）的M-K突變分析

以1997年為分界，將1960～1996和1997～2007年2個(gè)階段的水文氣象要素進(jìn)行對(duì)比分析，降水、徑流和氣溫分別從116.84 mm、103.26 m3/s和5.85 ℃變化至135.82 mm、133.05 m3/s和6.77 ℃，其中降水量增加了16.3%、徑流量增加28.8%，降水量的輕微變化將導(dǎo)致徑流相應(yīng)的變化，這與相關(guān)結(jié)論一致[15]；溫度上升0.92 ℃，這會(huì)使蒸發(fā)增加而導(dǎo)致徑流量減少，但有研究指出溫度每提高1 ℃可導(dǎo)致5%～6%空氣濕度的增加[16]。那么，以研究區(qū)的多年平均降水量121.19 mm作為實(shí)際蒸散量，且以溫度升高1 ℃來(lái)計(jì)算的話，最多只可能減少7 mm左右的水量，而降水增幅（18.98 mm）和徑流增幅（29.79 mm）的差值為-10.81 mm，可見單一的溫度變化是不可能使徑流變化為目前的狀態(tài)。溫度變化導(dǎo)致的冰雪融水以及人類活動(dòng)必然在其中扮演重要角色。

43卷16期

洪 波等氣候變化和人類活動(dòng)對(duì)開都河徑流的影響

2.4氣候變化和人類活動(dòng)的水文效應(yīng)評(píng)價(jià)

根據(jù)MK突變性檢驗(yàn)得到的結(jié)論，將1997年定為臨界年份，把水文要素劃分為1960～1996和1997～2007年兩段。根據(jù)雙累積曲線，由公式（1）建立起來(lái)的1960～1996年累積降水量與累積徑流量的線性關(guān)系（圖5）相關(guān)性較好（R2= 0.996），說(shuō)明可以用于將1997～2007年徑流還原至1960～1996年期間的受人類活動(dòng)影響水平的狀態(tài)。根據(jù)該回歸關(guān)系，可以計(jì)算得到1997～2007年的累積徑流量（圖6），計(jì)算值大致等同于臨界年份之前的受人類活動(dòng)影響較小時(shí)的狀態(tài)值，可見，該時(shí)段得到的計(jì)算徑流值大于實(shí)際測(cè)量值。然后，按照前面提到的方法，對(duì)1997～2007年的徑流進(jìn)行反推、

圖51960～1996年的累積降水量與累積徑流量的線性關(guān)系

圖61960～2007年徑流的計(jì)算值和實(shí)測(cè)值變化

還原并根據(jù)公式（2）和（3）計(jì)算2個(gè)差值，結(jié)果顯示（表2），通過還原得到的前段年徑流量均值為106.31 m3/s，對(duì)比相應(yīng)的實(shí)測(cè)值103.26 m3/s，誤差僅3%，說(shuō)明該方法的精確性較好。后段的實(shí)測(cè)徑流均值為133.05 m3/s，計(jì)算均值為126.60 m3/s，因而人類活動(dòng)的影響將會(huì)導(dǎo)致徑流量減小6.55 m3/s，但由表2可知實(shí)測(cè)的徑流共減小29.78 m3/s，因此，根據(jù)公式（4）可得人類活動(dòng)對(duì)徑流的最大貢獻(xiàn)率為-22%，同時(shí)從表2可以看出氣候變化對(duì)徑流的貢獻(xiàn)率為122%。

表2氣候變化和人類活動(dòng)對(duì)開都河徑流的影響

時(shí)段實(shí)測(cè)徑流m3/s計(jì)算徑流m3/s總變化量m3/s人類活動(dòng)影響的貢獻(xiàn)率∥%降雨影響的貢獻(xiàn)率%

1960～1996103.26106.31 ---

1997～2007133.05 126.60 29.78 -22122

3結(jié)論

1960～2007年開都河水源區(qū)氣候向暖濕方向發(fā)展，以1997年為界，降水增加了18.98 mm，溫度升高0.92 ℃，徑流在1997年開始突變性增加，截至2007年平均增加29.78 m3/s。氣候變化使水量增加而人類活動(dòng)致使水量減少，其貢獻(xiàn)率分別為122%和-22%。盡管氣候變化對(duì)開都河水源區(qū)徑流的影響具有主導(dǎo)性，但人類活動(dòng)是徑流減少的重要影響因素，可見，開都河水源區(qū)在進(jìn)行生態(tài)建設(shè)時(shí)應(yīng)充分考慮其不利影響。該研究采用雙累積曲線法建立受人類活動(dòng)影響較小時(shí)期的降雨徑流模型后，將受人類活動(dòng)影響較大時(shí)期的徑流還原至人類活動(dòng)影響較小的前一階段水平，進(jìn)而可分別求出氣候變化和人類活動(dòng)對(duì)水量變化的影響大小。一般情況

下，采用該方法可能會(huì)導(dǎo)致人類活動(dòng)貢獻(xiàn)率的最終結(jié)果變高，而氣候變化的貢獻(xiàn)率偏低。這是因?yàn)樵摲椒ǚ蛛x的是非

降水要素和降水要素的水文效應(yīng)，影響徑流變化的決定性因

素是降水，但溫度的升高將導(dǎo)致蒸發(fā)增加而徑流減少，該方法卻沒有考慮到溫度升高將導(dǎo)致的效應(yīng)。但由于該研究區(qū)存在冰川和積雪，溫度升高可增加冰雪融水對(duì)徑流的補(bǔ)給，這將縮小氣候變化對(duì)徑流的影響?？偟膩?lái)說(shuō)，這種方法獲得的評(píng)估結(jié)果是基本準(zhǔn)確的，具有重要的參考價(jià)值?？梢姡撗芯恐袣夂蜃兓蛩貙?duì)徑流的影響占主導(dǎo)地位，這主要是由于開都河水源區(qū)海拔高、地理環(huán)境復(fù)雜，極大程度上限制了人類活動(dòng)導(dǎo)致的。

參考文獻(xiàn)

[1] 竇睿音，延軍平，王鵬濤.全球變化背景下甘肅近半個(gè)世紀(jì)氣溫時(shí)空變化特征[J].干旱區(qū)研究，2015，32（1）：73-79.

[2] 郝振純，王加虎，李麗，等.氣候變化對(duì)黃河源區(qū)水資源的影響[J].冰川凍土，2006，28（1）：1-7.

[3] 陳亞寧，李稚，范煜婷，等.西北干旱區(qū)氣候變化對(duì)水文水資源影響研究進(jìn)展[J].地理學(xué)報(bào)，2014，69（9）：1295-1304.

[4] 張強(qiáng)，鄧振鏞，趙映東，等.全球氣候變化對(duì)我國(guó)西北地區(qū)農(nóng)業(yè)的影響[J].生態(tài)學(xué)報(bào)，2008，28（3）：1210-1218.

[5] 薛杰，李蘭海，李雪梅，等.開都河流域降水與徑流年內(nèi)分配特征及其變化的同步性分析[J].干旱區(qū)資源與環(huán)境，2014，28（12）：99-104.

[6] 趙直，徐晗.新疆開都河流域近50a徑流量年際年內(nèi)變化及其對(duì)氣候變化的響應(yīng)分析[J].干旱區(qū)資源與環(huán)境，2014，28（10）：151-156.

[7] 李新國(guó)，阿斯耶姆·圖爾迪，劉彬，等.近50年開都河流域下游降水和氣溫變化特征分析[J].干旱地區(qū)農(nóng)業(yè)研究，2013，31（3）：252-259.

[8] 王維霞，王秀君，姜逢清，等.開都河流域上下游過去50a氣溫降水變化特征分析[J].干旱區(qū)地理，2012，35（5）：746-753.

[9] 黃粵，陳曦，包安明，等.開都河流域山區(qū)徑流模擬及降雨輸入的不確定性分析[J].冰川凍土，2010，32（3）：567-572.

[10] 姚亞楠，徐俊榮，郭玲鵬，等.氣溫與降水形態(tài)關(guān)系研究——以開都河流域?yàn)槔齕J].干旱區(qū)研究，2015，32（1）：94-101.

[11] 王莉.基于SWAT模型的氣候變化對(duì)水資源響應(yīng)研究——以開都河流域?yàn)槔齕D].南京：南京信息工程大學(xué)，2012：9-10.

[12] 黃粵，陳曦，包安明.干旱區(qū)資料稀缺流域日徑流過程模擬[J].水科學(xué)進(jìn)展，2009，20（3）：332-336.

[13] 穆興民，張秀勤，高鵬，等.雙累積曲線方法理論及在水文氣象領(lǐng)域應(yīng)用中應(yīng)注意的問題[J].水文，2010，30（4）：47-51.

[14] 冉大川，劉斌.雙累積曲線計(jì)算水土保持減水減沙效益方法探討[J].人民黃河，1996，18（6）：24-25.

[15] CHIEW F H S，MCMAHON T A.Modelling the impacts of climate change on Australian streamflow[J].Hydrological Processes，2002，16（6）：1235-1245.

[16] 陳亞寧，徐宗學(xué).全球氣候變化對(duì)新疆塔里木河流域水資源的可能性影響[J].中國(guó)科學(xué)D輯，2004，34（11）：1047-1053.