劉軍

(水利部新疆維吾爾自治區(qū)水利水電勘測(cè)設(shè)計(jì)研究院，烏魯木齊 830000)

近年來(lái)，全球氣候變化不僅加速了水分循環(huán)，還加劇了極端水文事件的發(fā)生，從而導(dǎo)致了水資源在時(shí)空上的重新分布和水資源數(shù)量的改變，進(jìn)而影響到部分流域綠洲生態(tài)環(huán)境和經(jīng)濟(jì)社會(huì)的發(fā)展[1]。流域徑流的產(chǎn)生源于一個(gè)復(fù)雜的水文物理過(guò)程，在此期間影響流域徑流變化的因素有很多，例如氣象條件、地理環(huán)境以及流域特性等。正是由于這些因素的干擾使得流域徑流的變化往往呈現(xiàn)出多種不確定、非線(xiàn)性、非平穩(wěn)的波動(dòng)性變化[2－3]。

地表水資源作為大氣環(huán)流系統(tǒng)和水文循環(huán)中不可或缺的要素，既受氣候變化的影響，同時(shí)又反作用于氣候變化。近年來(lái)諸多學(xué)者對(duì)流域徑流的分析從定性到定量，由簡(jiǎn)單的徑流演變規(guī)律與趨勢(shì)分析[4－6]到氣候影響下的徑流趨勢(shì)分析[7－11]，最后再到氣候變化和人類(lèi)活動(dòng)對(duì)于流域徑流變化影響的定量分析[12－16]。如劉靜等[4]利用多年實(shí)測(cè)徑流資料以塔河流域三條源流為研究對(duì)象，運(yùn)用多種方法研究多時(shí)間尺度下徑流的變化特征，探討環(huán)境變化對(duì)徑流趨勢(shì)變化的影響機(jī)理。柏玲[8]綜合運(yùn)用多方法、多尺度的建模技術(shù)對(duì)開(kāi)都河與阿克蘇河源流區(qū)出山口徑流和相關(guān)氣候因子對(duì)研究區(qū)氣候—徑流過(guò)程及其機(jī)理進(jìn)行對(duì)比分析，得出兩河徑流對(duì)氣溫和降水敏感性的差異可能與徑流的組成比例有關(guān)，不同時(shí)段和尺度的差異是由于氣溫和降水兩個(gè)主要?dú)夂蛞匾虏煌瑫r(shí)段補(bǔ)給徑流來(lái)源的變化，進(jìn)一步證明了氣候變化對(duì)河道徑流的影響。Zhang 等[13]發(fā)現(xiàn)在氣候變化和人類(lèi)活動(dòng)的雙重影響下，河流系統(tǒng)的水文過(guò)程處在不斷的變化過(guò)程中。王曉穎等[14]通過(guò)分析氣候變化和人類(lèi)活動(dòng)對(duì)白河流域徑流變化的貢獻(xiàn)率得出研究區(qū)2000 年以后氣候變化對(duì)徑流變化的影響呈上升趨勢(shì)。

位于干旱區(qū)的典型河流——庫(kù)山河源區(qū)以山地、盆地相間的地貌格局為特點(diǎn)，對(duì)氣候變化尤為敏感。簡(jiǎn)單地說(shuō)，氣候變化不僅加劇了流域內(nèi)水循環(huán)過(guò)程，相反，徑流的變化也在一定程度上使得研究區(qū)的氣候更加復(fù)雜多變，同時(shí)這種關(guān)系也增加了該地區(qū)水資源的不確定性?？梢哉f(shuō)，氣溫和降水變化都可以對(duì)庫(kù)山河流域地表徑流產(chǎn)生影響[8]。因此，本文利用不同的演變分析方法開(kāi)展庫(kù)山河流域在氣候影響下徑流長(zhǎng)期演變規(guī)律的研究，厘清徑流量的增減趨勢(shì)，并在此基礎(chǔ)上揭示多年氣候變化對(duì)徑流的影響。研究結(jié)論對(duì)于研究區(qū)水資源合理配置及可持續(xù)利用、流域水資源規(guī)劃以及促進(jìn)研究區(qū)生態(tài)可持續(xù)發(fā)展具有十分重要的現(xiàn)實(shí)意義。

1 材料及研究方法

1.1 研究區(qū)概況

庫(kù)山河是喀什噶爾河水系的第三大河。流域涉及克孜勒蘇柯?tīng)柨俗巫灾沃莸陌⒖颂湛h和喀什地區(qū)的英吉沙縣、疏勒縣。流域東與塔里木盆地銜接，南與葉爾羌河支流塔什庫(kù)爾干河及卡龍溝相鄰，西以慕土塔格山和公格爾山與蓋孜河分水，北與蓋孜河流域毗連，為喀什噶爾河水系的源流之一。庫(kù)山河流域面積為10924km2，干流發(fā)源于海拔7546m 的慕土塔格山和海拔7719m 的公格爾山的高山冰川，徑流以冰雪融水補(bǔ)給為主、降雨補(bǔ)給為輔，枯季為地下水和泉水補(bǔ)給，河流流經(jīng)英吉沙、阿克陶、疏勒等縣境內(nèi)，最后消失于布克拉沙漠中。近年來(lái)，由于下游河道整治，灌溉引水以及沙漠區(qū)的蒸發(fā)滲漏，庫(kù)山河地表水已難匯入喀什噶爾河干流，逐漸成為獨(dú)立的水系。

庫(kù)山河流域氣候呈明顯的垂直地帶性分布，山區(qū)一年只有冷、暖兩季，夏季短且涼爽，冬季長(zhǎng)而嚴(yán)寒，年平均氣溫10℃左右。海拔5000m 以上地區(qū)降水較多，年降水量達(dá)500mm～700mm，永久性冰川發(fā)育，形成較大的固體水庫(kù)，從而成為流域徑流的主要補(bǔ)給來(lái)源；低山區(qū)和平原區(qū)降水稀少，年降水量小于70mm，而蒸發(fā)量高達(dá)2700mm。

1.2 研究方法

本文通過(guò)獲取并收集庫(kù)山河流域1956 年設(shè)置的沙曼水文站近60 年的降水和水文資料以及近50 年的氣溫資料，運(yùn)用趨勢(shì)分析法、Mann－Kendall 趨勢(shì)檢驗(yàn)法及集合經(jīng)驗(yàn)?zāi)B(tài)分解法(EEMD，Ensemble Empirical Mode Decomposition)分析庫(kù)山河流域在多年氣候影響下的徑流演變規(guī)律及增減趨勢(shì)。

1.2.1 Mann－Kendall 趨勢(shì)檢驗(yàn)法

Mann－Kendall[17](以下簡(jiǎn)稱(chēng)M－K)方法所需樣本不需要遵循特定的分布規(guī)律，也不受數(shù)據(jù)異常值的影響，同時(shí)還適用于多種類(lèi)型及順序變量的分析。因此，作為一種非參數(shù)統(tǒng)計(jì)檢驗(yàn)方法，由于其計(jì)算簡(jiǎn)單、適用性強(qiáng)，M－K 趨勢(shì)檢驗(yàn)法被廣泛應(yīng)用于時(shí)間序列的變化趨勢(shì)分析。同時(shí)，該方法還可以通過(guò)計(jì)算統(tǒng)計(jì)量UF和UB以應(yīng)用于檢驗(yàn)時(shí)間序列數(shù)據(jù)是否發(fā)生了突變[18]。

對(duì)于時(shí)間序列Y，M－K 趨勢(shì)檢驗(yàn)的統(tǒng)計(jì)量如下:

式中，n為樣本數(shù)據(jù)的個(gè)數(shù)；xi為時(shí)間序列的第i個(gè)數(shù)值；sgn 為符號(hào)函數(shù)。另外，Mann(1945)和Kendall(1975)證明，當(dāng)n≥8 時(shí)，統(tǒng)計(jì)量P大致服從均值為0 的正態(tài)分布。方差為:

式中，ti為第i組數(shù)據(jù)點(diǎn)的數(shù)目。標(biāo)準(zhǔn)化統(tǒng)計(jì)量用下式計(jì)算:

式中，Zc服從標(biāo)準(zhǔn)正態(tài)分布，β用于衡量趨勢(shì)大小:

式中，10 時(shí)，數(shù)據(jù)呈現(xiàn)“上升趨勢(shì)”，反之則呈現(xiàn)“下降趨勢(shì)”。

1.2.2 集合經(jīng)驗(yàn)?zāi)B(tài)分解法(EEMD)

集合經(jīng)驗(yàn)?zāi)B(tài)分解(EEMD)方法[19]是在經(jīng)驗(yàn)?zāi)B(tài)分解(Empirical Mode Decomposition，EMD)的基礎(chǔ)上發(fā)展起來(lái)的，其首先對(duì)時(shí)間序列的數(shù)據(jù)進(jìn)行局部平穩(wěn)化處理，然后通過(guò)希爾伯特變換獲得時(shí)頻譜圖，從而獲得一些具有物理意義的頻率分量。其分解表達(dá)式為:

式中，Ci(t)為信號(hào)在i固有尺度上的本征模函數(shù)(IMFs)；rn(t)為信號(hào)分解后的趨勢(shì)項(xiàng)(RES)。

集合經(jīng)驗(yàn)?zāi)B(tài)分解法的過(guò)程可概括為:(1)在原數(shù)據(jù)的基礎(chǔ)上加入白噪聲序列；(2)對(duì)處理后的數(shù)據(jù)進(jìn)行經(jīng)驗(yàn)?zāi)B(tài)分解；(3)反復(fù)上述步驟得到不同的本征模態(tài)函數(shù)(IMFs)；(4)將每次分解得到的IMFs 進(jìn)行集合平均，用于抵消加入的白噪聲使其作為最終的分解結(jié)果。

2 氣候、徑流變化特征分析

已有研究表明氣候—水文系統(tǒng)是一種非線(xiàn)性的復(fù)雜系統(tǒng)[1，4，11]，而用來(lái)表征區(qū)域內(nèi)氣候—水文變化的因子如降雨量、溫度及徑流量等在長(zhǎng)期的時(shí)間序列上的變化也表現(xiàn)為一種非線(xiàn)性、非平穩(wěn)的復(fù)雜變化過(guò)程，同時(shí)還伴隨有多種尺度或周期性的振蕩?；诖?，本文以庫(kù)山河流域?yàn)槔?，采用多種方法聯(lián)合對(duì)比分析，從不同年際、年代際尺度分析庫(kù)山河徑流的變化特征。

2.1 趨勢(shì)性分析

通過(guò)收集到的1956—2018 年庫(kù)山河沙曼水文站年降雨量、徑流量及1962—2010 年的氣溫資料(部分缺失的數(shù)據(jù)運(yùn)用調(diào)整插補(bǔ)方法獲取)，運(yùn)用趨勢(shì)分析法及M－K 分析法對(duì)氣候要素和徑流量的變化趨勢(shì)進(jìn)行分析，結(jié)果如圖1 所示。庫(kù)山河流域降雨量、溫度及徑流量均呈增加趨勢(shì)，但可以明顯看出其增加趨勢(shì)并不是顯著的線(xiàn)性關(guān)系，這與劉靜[4]、金慶日[11]等得出的結(jié)論一致。其中線(xiàn)性回歸結(jié)果表明，研究流域內(nèi)降雨量、溫度和徑流量年增加率分別為1.02、0.99 和0.01，說(shuō)明流域內(nèi)徑流量波動(dòng)較大，其增加趨勢(shì)相對(duì)于降雨量及溫度并不顯著，降雨量和溫度的增量要略大于徑流量增量。

圖1 庫(kù)山河流域氣候、徑流量年際變化特征

距平分析發(fā)現(xiàn)，流域內(nèi)降雨量整體呈上升趨勢(shì)，溫度和徑流量卻呈現(xiàn)先下降后上升的趨勢(shì)，溫度及徑流量分別自1996 年和2000 年以后出現(xiàn)了較為明顯的上升趨勢(shì)。M－K 趨勢(shì)檢驗(yàn)結(jié)果如表1 所示，降雨量和溫度的Zc值分別為2.46和4.44，均≥2.32(P＝0.01)，徑流量Zc值為1.72≥1.64(P＝0.05)。以上表明降雨量、溫度通過(guò)了置信度為0.99 的顯著性檢驗(yàn)，徑流量通過(guò)了置信度為0.95 的顯著性檢驗(yàn)。根據(jù)衡量趨勢(shì)變化幅度的β值結(jié)果，降雨量、溫度和徑流量分別為1.20、0.28 和0.05，均大于0，表明近年來(lái)庫(kù)山河流域年平均降雨量增加趨勢(shì)顯著，溫度及徑流量上升趨勢(shì)雖不顯著，但不可否認(rèn)的是流域降雨量、溫度及徑流量均呈增加趨勢(shì)，與趨勢(shì)分析法結(jié)果一致，說(shuō)明該結(jié)論較為合理。

表1 庫(kù)山河流域氣候、徑流趨勢(shì)分析結(jié)果

2.2 突變性分析

突變分析是我們辨識(shí)流域徑流波動(dòng)變化的主要方法之一，對(duì)流域內(nèi)水文變化的性質(zhì)預(yù)測(cè)具有重要意義[2]。運(yùn)用M－K 突變分析在顯著性水平α＝0.05 的情況下將統(tǒng)計(jì)量UF和UB序列數(shù)據(jù)相結(jié)合做出圖2。結(jié)果表明，庫(kù)山河降雨量、溫度及徑流量UF和UB兩個(gè)統(tǒng)計(jì)量分別在1996 年、1981 年及2002 年相繼出現(xiàn)突變，其突變時(shí)間與上述距平分析時(shí)間較為接近，具有一致性。

圖2 庫(kù)山河流域氣候、徑流量突變檢驗(yàn)結(jié)果

根據(jù)圖中UF的變化趨勢(shì)可以發(fā)現(xiàn)，研究區(qū)降雨量從1956 年開(kāi)始到1962 年前，其值基本都小于0，說(shuō)明在此期間庫(kù)山河降雨量呈下降趨勢(shì)[18]；從1963 年至2018 年，UF的統(tǒng)計(jì)值都大于0，表明庫(kù)山河在該時(shí)段降雨表現(xiàn)為較明顯的增加趨勢(shì)。以上結(jié)果與上述分析結(jié)果相一致。突變前、后降雨量的平均值分別為140.48mm 和170.98mm，相對(duì)變化為21.71%。氣溫從1962 年開(kāi)始到1969 年前，其值均小于0，說(shuō)明在此期間庫(kù)山河流域溫度呈下降趨勢(shì)；從1970 年至2010 年，UF的統(tǒng)計(jì)值都大于0，表明庫(kù)山河流域在該時(shí)段氣溫呈現(xiàn)增加趨勢(shì)。突變前、后氣溫的平均值分別為8.59℃和9.37℃，相對(duì)變化為9.08%。徑流量從1960 年開(kāi)始到2002 年前，其值基本都小于0，說(shuō)明在此期間庫(kù)山河徑流波動(dòng)呈下降趨勢(shì)；從2003 年至2018 年，UF的統(tǒng)計(jì)值都大于0，表明庫(kù)山河在該時(shí)段徑流表現(xiàn)為增加的趨勢(shì)。突變前、后徑流量平均值分別為6.38 億m3/a 和7.03 億m3/a，相對(duì)變化為10.19%。該結(jié)果也反映出流域內(nèi)降雨量變化明顯強(qiáng)于溫度及徑流量。而流域內(nèi)降雨量的增量明顯大于徑流量的增量，這可能是受全球氣候變暖的影響，大氣降水和冰川融水逐漸增多所致。該結(jié)論雖與上述趨勢(shì)性分析得出的流域氣溫上升趨勢(shì)不顯著的結(jié)果略有差異，但不可否認(rèn)的是庫(kù)山河流域內(nèi)降雨量、溫度及徑流量都呈增加趨勢(shì)。

2.3 波動(dòng)性分析

對(duì)于距平后的降雨、溫度及徑流時(shí)間序列，運(yùn)用集合經(jīng)驗(yàn)?zāi)B(tài)分解方法進(jìn)行分解，均得到了4 個(gè)IMF 分量和一個(gè)趨勢(shì)項(xiàng)RES。為驗(yàn)證數(shù)據(jù)分解后的完備性和可靠性，將EEMD 分解得到的各分量和趨勢(shì)項(xiàng)疊加合成一個(gè)降雨、溫度及徑流的重構(gòu)序列與原始數(shù)據(jù)進(jìn)行擬合，結(jié)果如圖3 所示。從中發(fā)現(xiàn)各變量相對(duì)誤差均在0.1%之內(nèi)，表明EEMD 方法在此次分解中具有很好的完備性，能較好地反映原始數(shù)據(jù)在時(shí)間序列下的多尺度變化特征。

圖3 各參數(shù)EEMD 重構(gòu)序列與原始數(shù)據(jù)擬合圖

各分量結(jié)果依次反映了庫(kù)山河流域多年氣候要素及徑流量從高頻到低頻不同時(shí)間尺度下的波動(dòng)特征及流域降雨量、溫度及徑流量復(fù)雜的多時(shí)間尺度性[5]。結(jié)合獲取的趨勢(shì)項(xiàng)RES，其表明庫(kù)山河流域降雨量、溫度及徑流量隨時(shí)間的變化整體呈增加趨勢(shì)，究其原因可能與區(qū)域內(nèi)降水增加及氣候變暖等因素有關(guān)，這與前文所得結(jié)論一致。同時(shí)，表2 和圖4 中各分量顯示研究區(qū)降雨量、溫度及徑流量在年際尺度上具有2 年—3 年和5 年—7 年的準(zhǔn)周期性變化，在年代際尺度上具有11 年—17 年和27 年—29 年的準(zhǔn)周期性變化。其中降雨量、溫度及徑流量在IMF1 分量上方差貢獻(xiàn)率最大，分別為46.84%、44.48%和38.56%，說(shuō)明高頻振蕩在波動(dòng)頻率上對(duì)流域降雨及徑流總體變化特征影響程度最大，該分量降雨、溫度及徑流分別在1977—1990 年、1974—1988 年以及1997—2002 年振幅較小，其余時(shí)段振幅相對(duì)較高；同時(shí)，降雨量、溫度及徑流量中IMF3 分量均對(duì)各變量整體波動(dòng)影響最小，其方差貢獻(xiàn)率分別為2.58%、6.60%和3.03%。以上分析表明，庫(kù)山河流域內(nèi)年際振蕩(IMF1)在流域降水及徑流長(zhǎng)期變化中占據(jù)主導(dǎo)地位，可能是因?yàn)閹?kù)山河徑流受河源冰川固體水庫(kù)的調(diào)節(jié)作用，徑流量年際變化穩(wěn)定。這與柏玲[8]在天山南坡典型流域阿克蘇河和開(kāi)都河研究結(jié)果一致。另外結(jié)合各分量與原始數(shù)據(jù)的相關(guān)性發(fā)現(xiàn)，不論是年際還是年代際變化均與實(shí)測(cè)數(shù)據(jù)存在顯著的正相關(guān)關(guān)系，說(shuō)明流域內(nèi)短期的年際變化與長(zhǎng)期的年代際變化趨勢(shì)均有一定的研究意義。

表2 庫(kù)山河流域氣候、徑流集合經(jīng)驗(yàn)?zāi)B(tài)分解結(jié)果表

圖4 庫(kù)山河流域氣候、徑流集合經(jīng)驗(yàn)?zāi)B(tài)分解圖

2.4 氣候變化對(duì)徑流的影響

為定量分析庫(kù)山河流域氣候變化對(duì)徑流的影響，本文利用沙曼水文站實(shí)測(cè)資料采用相關(guān)性分析及通徑分析方法對(duì)研究流域進(jìn)行研究。但由于溫度資料的不足，本次在分析時(shí)選取1962—2010 年氣候要素及徑流數(shù)據(jù)進(jìn)行分析。

2.4.1 相關(guān)性分析

流域源流區(qū)氣候變化與徑流的相關(guān)關(guān)系結(jié)果如表3 所示。不難看出，流域年徑流量與年降水量及溫度之間存在著較為顯著的相關(guān)關(guān)系。分析認(rèn)為流域徑流量受到補(bǔ)給來(lái)源的影響較大，而庫(kù)山河流域補(bǔ)給來(lái)源主要為上游山區(qū)的冰雪融水，而這又恰恰受到氣候的影響。例如在氣候干燥的年份，降雨少，氣溫高，冰川融水多；在氣候較為濕潤(rùn)的年份，降雨多，氣溫較低，冰川融水也相對(duì)較少，這種關(guān)系不僅使得流域內(nèi)年際徑流相對(duì)穩(wěn)定，同時(shí)也間接證明了氣候變化會(huì)影響徑流從而起到調(diào)節(jié)作用，這與陳忠升[6]在西北干旱區(qū)河川徑流變化及歸因定量辨識(shí)中的分析結(jié)論一致。

表3 庫(kù)山河流域氣候變化與徑流之間的關(guān)系

2.4.2 通徑分析

使用通徑分析時(shí)通常不僅要求數(shù)據(jù)之間存在顯著的相關(guān)關(guān)系，同時(shí)還要求因變量服從正態(tài)分布[20]。上述分析表明流域內(nèi)徑流量與降雨量及溫度存在明顯的相關(guān)關(guān)系。另外，本次采用Kolmogorov－Smirnov 和Shapiro－Wilk 正態(tài)性檢驗(yàn)方法分析流域徑流量是否滿(mǎn)足正態(tài)性分布，結(jié)果如表4 所示。其顯著性分別為0.200 及0.967，均大于0.05，說(shuō)明庫(kù)山河徑流量在95%的置信區(qū)間下服從正態(tài)分布。因此可在流域內(nèi)徑流量與氣候的年際尺度上對(duì)其進(jìn)行通徑分析。

表4 徑流量正態(tài)性檢驗(yàn)統(tǒng)計(jì)表

溫度對(duì)徑流的間接通徑系數(shù)為溫度與降雨量的相關(guān)系數(shù)乘以降雨量對(duì)徑流的直接通徑系數(shù)，同理降雨量對(duì)徑流的間接通徑系數(shù)也與此相對(duì)應(yīng)，據(jù)此得出通徑分析結(jié)果如表5 所示。降雨量及溫度對(duì)徑流的通徑系數(shù)分別為0.295 及0.654，表明降雨及溫度對(duì)徑流的產(chǎn)生均具有驅(qū)動(dòng)作用，但溫度是徑流產(chǎn)生的主要驅(qū)動(dòng)因子，這與冰雪融水為庫(kù)山河流域主要補(bǔ)給來(lái)源相符。

表5 庫(kù)山河降雨量和氣溫與徑流量的通徑分析

3 討論

上述研究發(fā)現(xiàn)，庫(kù)山河流域地表徑流量受氣候因素的調(diào)節(jié)作用年際波動(dòng)較為穩(wěn)定，同時(shí)徑流量呈上升趨勢(shì)。但流域供水量調(diào)查統(tǒng)計(jì)數(shù)據(jù)顯示，庫(kù)山河流域地表水開(kāi)發(fā)利用率達(dá)到85.8%，地下水開(kāi)發(fā)利用率達(dá)到138.6%，水資源總量開(kāi)發(fā)利用率達(dá)129.8%。由此可見(jiàn)，喀什噶爾河流域水資源總量開(kāi)發(fā)利用程度最高。這也使得庫(kù)山河流域仍面臨著如下問(wèn)題:(1)水資源開(kāi)發(fā)利用量超出控制指標(biāo)要求，經(jīng)濟(jì)發(fā)展和生態(tài)修復(fù)建設(shè)存在爭(zhēng)水矛盾；(2)現(xiàn)狀水資源利用超過(guò)流域水資源承載能力；(3)流域用水結(jié)構(gòu)不合理，水利用效益低；(4)流域農(nóng)業(yè)用水效率低下。針對(duì)上述問(wèn)題，結(jié)合庫(kù)山河流域地表水資源演變規(guī)律及發(fā)展趨勢(shì)，在進(jìn)行區(qū)域水資源分配及流域規(guī)劃時(shí)應(yīng)本著節(jié)水優(yōu)先、高效利用，量水而行、協(xié)調(diào)發(fā)展，合理控制、增加調(diào)控，統(tǒng)籌協(xié)調(diào)、合理配置等原則，在滿(mǎn)足流域各行業(yè)用水需求的前提下，保障流域供水安全，更好地落實(shí)流域水資源配置方案，確保流域工農(nóng)業(yè)用水和下游生態(tài)需水。

本文僅探討了多年氣候變化對(duì)庫(kù)山河流域徑流變化趨勢(shì)及演變規(guī)律的影響，但庫(kù)山河大部分河道處于天然狀態(tài)，流域內(nèi)徑流動(dòng)態(tài)還受到現(xiàn)狀水利工程如渠首工程、堤防護(hù)岸工程、水庫(kù)工程及下游灌溉等人類(lèi)活動(dòng)的影響。因此，僅討論多年氣候條件變化對(duì)流域的動(dòng)態(tài)影響略顯不足。為有效地保護(hù)及利用庫(kù)山河流域地表水資源，下一步應(yīng)著重考慮氣候變化和人類(lèi)活動(dòng)疊加后對(duì)徑流變化的綜合影響。

4 結(jié)論

(1)基于多年實(shí)測(cè)數(shù)據(jù)，采用多種方法綜合分析研究區(qū)受氣候變暖的影響，研究表明降雨量、氣溫與徑流量均呈明顯的增加趨勢(shì)，但并不是呈線(xiàn)性增加而是存在著非線(xiàn)性、非平穩(wěn)的復(fù)雜波動(dòng)。

(2)通過(guò)EEMD 方法的波動(dòng)性分析發(fā)現(xiàn)，流域內(nèi)徑流受氣候變化的影響年際振蕩在其長(zhǎng)期變化中占據(jù)主導(dǎo)地位。同時(shí)，結(jié)合相關(guān)分析發(fā)現(xiàn)不論是年際還是年代際變化均與實(shí)測(cè)數(shù)據(jù)存在較顯著的正相關(guān)關(guān)系，說(shuō)明流域內(nèi)短期的年際變化與長(zhǎng)期的年代際變化趨勢(shì)均有一定的研究意義。

(4)研究流域內(nèi)氣候變化對(duì)流域徑流起調(diào)節(jié)作用，同時(shí)降雨及溫度均對(duì)徑流的產(chǎn)生具有驅(qū)動(dòng)作用。但溫度是徑流產(chǎn)生的主要驅(qū)動(dòng)因子，這與冰雪融水為庫(kù)山河流域主要補(bǔ)給來(lái)源相符。