侯 蕾，彭文啟，董 飛，劉經(jīng)強(qiáng)，陳渠昌，黃愛(ài)平，陳學(xué)凱

(1.中國(guó)水利水電科學(xué)研究院，流域水循環(huán)模擬與調(diào)控國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京100038；2.山東農(nóng)業(yè)大學(xué)水利土木工程學(xué)院，山東 泰安 271000)

0 引 言

IPCC第五次評(píng)估報(bào)告：1982-2012年間，全球的地表平均溫度升高了0.85 ℃，到21世紀(jì)末，全球地表溫度變化可能超過(guò)1.5 ℃[1]。氣候變化對(duì)陸面水文循環(huán)產(chǎn)生了顯著影響[2,3-6]。研究全球氣候變暖背景下河川徑流的變化已成為各界普遍關(guān)注的熱點(diǎn)問(wèn)題[7, 8]。河川徑流是重要的水資源形式，受氣候變化和下墊面變化的影響較為明顯，在時(shí)間和空間上多尺度研究氣候變化對(duì)水資源的影響程度，利用長(zhǎng)序列氣象水文數(shù)據(jù)分析流域水文氣象要素的演變特征，是實(shí)現(xiàn)水資源開(kāi)發(fā)利用及管理的重要基礎(chǔ)性工作[9]。永定河流域是我國(guó)北方干旱半干旱區(qū)頗受關(guān)注的水資源及生態(tài)環(huán)境脆弱地區(qū)。近年來(lái)，永定河流域年徑流量出現(xiàn)了銳減趨勢(shì)[10-12]，多數(shù)河段出現(xiàn)斷流現(xiàn)象，尤其是2000年以后官?gòu)d水庫(kù)年均徑流量?jī)H為0.9 億m3，較1955-1981年豐水期，衰減了93%，較1955-1981年平水期衰減了60%～70%。對(duì)于永定河流域的水資源問(wèn)題已有學(xué)者和專家進(jìn)行了大量的研究[13-17]，目前為止，永定河流域的研究多集中于徑流衰減成因的識(shí)別[18-20]，對(duì)于永定河上游流域及主要支流的氣候水文要素演變特征的系統(tǒng)研究較少，并且時(shí)間序列較短。本文系統(tǒng)分析了流域氣溫和降水的空間分布格局，揭示了永定河上游流域及主要支流的氣溫、降水、徑流的變化特征，分析了流域徑流變化的主要影響因素，為流域水資源保護(hù)、水資源的開(kāi)發(fā)利用及管理提供了技術(shù)支撐。

1 流域概況及研究方法

1.1 流域概況

永定河流域位于112°～117°45′E，39°～41°20′N之間，是海河流域的七大支流之一，流域面積約為4.70 萬(wàn)km2，約為海河流域面積的14.7%[21]。流域地跨內(nèi)蒙古、山西、河北、北京、天津等5個(gè)省(自治區(qū)、直轄市)，主要有桑干河、洋河兩大支流，在河北懷來(lái)縣匯合后稱永定河，在官?gòu)d水庫(kù)納入媯水河，經(jīng)官?gòu)d水庫(kù)流入官?gòu)d山峽于三家店進(jìn)入平原。本研究的區(qū)域在永定河上游流域，研究區(qū)面積約為4.3 萬(wàn)km2，位于112°～116°22′E，38°55′～41°15′N，官?gòu)d水庫(kù)以上稱為永定河上游，約為永定河流域面積的93%，見(jiàn)圖1。永定河上游流域跨燕山山脈、內(nèi)蒙古高原和華北平原，高程有很大差異，山區(qū)面積約占流域面積的80%，以山區(qū)為主。

圖1 永定河上游流域位置

1.2 研究方法

基于1960-2012年永定河上游流域官?gòu)d水庫(kù)、響水堡、石匣里三個(gè)水文站(來(lái)源于《海河流域水文年鑒》)的年徑流量以及流域范圍及周邊15個(gè)國(guó)家氣象站的逐日最高氣溫、最低氣溫、平均氣溫、降水等氣象數(shù)據(jù)，采用Mann-Kendall 非參數(shù)檢驗(yàn)[22,23]及線性回歸法對(duì)流域水文氣象要素進(jìn)行了趨勢(shì)性分析；采用累積距平法[24]對(duì)流域水文氣象要素進(jìn)行了突變性分析；采用Morlet小波分析了流域水文氣象要素的周期性演變特征。文中官?gòu)d水庫(kù)水文站控制流域簡(jiǎn)稱永定河上游流域，響水堡水文站控制流域簡(jiǎn)稱洋河流域，石匣里水文站控制流域簡(jiǎn)稱桑干河流域。

2 結(jié)果與討論

2.1 降水、氣溫的空間分布特征

采用趨勢(shì)面分析法，利用1960-2012年永定河上游流域及周圍15個(gè)氣象臺(tái)站的統(tǒng)計(jì)數(shù)據(jù)來(lái)繪制流域多年的平均降水和多年平均氣溫的等值線圖(圖2、圖3)。在空間上，降水量地域性差異較大，年降水量自東南向西北遞減；受燕山山脈迎風(fēng)坡多雨帶的影響，在桑干河上游以及媯水河流域的多年平均降水量較大，其中，東南部五臺(tái)山氣象站的降水量達(dá)到744 mm；受大陸性半干旱氣候的影響，洋河的上游流域多年平均降水量相對(duì)較低，其中，化德氣象站的多年平均降水量約為314 mm。多年平均氣溫地域性差異也較大，受流域地形和高程的影響，流域多年平均氣溫以北京氣象站為中心向外延伸逐漸遞減。受北京城市熱島效應(yīng)[25]的影響，官?gòu)d水庫(kù)及媯水河區(qū)域氣溫偏高；張北地區(qū)受北寒帶干旱氣候以及東南地區(qū)地形的影響，洋河和桑干河流域的上游地區(qū)氣溫較低。

圖2 永定河上游流域年降水等值線圖

圖3 永定河上游流域年平均氣溫空間分布圖

2.2 降水、氣溫、徑流的年內(nèi)及年代際變化特征

根據(jù)1960-2012年永定河上游流域、桑干河流域、洋河流域各月的降水統(tǒng)計(jì)分析(見(jiàn)圖4)，結(jié)果表明，流域降水量年內(nèi)分配極不均勻，6-9月降水量占全年降水的75%左右；12月份和1月份的降水量最小，各流域月降水量不到全年的0.5%，均小于4 mm。通過(guò)對(duì)降水量年代距平的變化特征分析(見(jiàn)圖5)。永定河河上游流域在20世紀(jì)60、70年代是豐水期，尤其60年代，降水量為452 mm，比多年平均降水值增加了6.4%，80年代、2000-2012年為枯水期，80年代的年降水量比多年平均降水量減少了15.89 mm，2000-2012年降水量較多年平均值減少了21.79 mm，90年代月降水量與多年平均的月降水變化趨勢(shì)較為一致；桑干河流域在60、70年代也是豐水期，60年代的降水量比多年平均降水增加了9.5%，80年代以后為枯水期，尤其是90年代出現(xiàn)了銳減現(xiàn)象，該時(shí)期降水量較多年平均值減少4.9%；洋河流域在60年代、70年代以及90年代是豐水期，其中70年代、90年代的降水量分別比多年平均值增加了6.3%、5.8%，在80年代、2000-2012年間，流域降水偏枯，尤其是80年代，該時(shí)期降水量比多年平均值減少了7.7%。

圖5 永定河上游流域、桑干河流域和洋河流域降水距平值年代變化圖

通過(guò)分析不同年代的月降水量的變化特征(圖4)，90年代流域月降水最大值出現(xiàn)在7月份，其余各年代降水量的分布不均衡，最大值主要出現(xiàn)在7月份和8月份；60、70年代各流域的降水偏豐，非汛期的降水量基本與多年平均降水量一致，汛期的月降水量均高于多年平均降水量；從2000-2012開(kāi)始各流域月降水的分配情況發(fā)生了明顯的變化，其中非汛期的降水量均大于多年平均降水量，而汛期的月降水量均小于多年平均值。

圖4 永定河上游流域、桑干河流域和洋河流域各年代降水月均值變化圖

通過(guò)分析1960-2012年各流域月平均氣溫的變化特征(圖6)。各流域四季氣溫變化分明，1月份各流域的氣溫最低，分別為-10.51、-10.96、-12.03 ℃；7月份各流域的氣溫最高，分別為21.71、20.90、21.25 ℃；秋季以后氣溫逐漸下降。通過(guò)分析不同年代各月平均氣溫的變化特征(圖6)，20世紀(jì)60、70和80年代，各月的平均氣溫均小于多年平均氣溫；90年代，各月平均氣溫與多年平均值一致；2000年以后，各月平均氣溫均大于多年平均值。通過(guò)對(duì)各流域氣溫年代距平的變化特征分析(圖7)，60、70、80年代，各流域的氣溫均降低；90年代，氣溫開(kāi)始逐漸升高；2000年以后，各流域的氣溫劇增。

圖6 永定河上游流域、桑干河流域和洋河流域各年代氣溫月均值變化圖

圖7 永定河上游流域、桑干河流域和洋河流域氣溫距平值年代變化圖

通過(guò)分析各流域不同年代各月的徑流量變化特征(圖8)。研究表明，各流域的在7、8月份(夏汛)和3、4月份(春汛)徑流量較大，總體上與華北地區(qū)的氣候變化趨勢(shì)一致。20世紀(jì)70年代，各流域的月徑流均大于相應(yīng)的月平均值；80、90年代，各流域的月徑流量與相應(yīng)的月均值變化一致；2000年以后，各流域的月徑流偏枯，均小于相應(yīng)的月均值；與不同年代的月降水、月氣溫的變化進(jìn)行對(duì)比分析，70年代各流域月徑流變化特征與月降水和月氣溫的變化一致，80年代以后，各流域月徑流變化特征與月降水和月氣溫的變化特征差異性較大；結(jié)果表明，70年代，氣象條件對(duì)流域徑流量的影響較為顯著；80年代以后，人類活動(dòng)成為了流域徑流變化的主導(dǎo)因素。

圖8 永定河上游流域、桑干河流域和洋河流域各年代徑流月均值變化對(duì)比圖

通過(guò)對(duì)各流域徑流量年代距平的變化特征分析(圖9)。結(jié)果表明，在20世紀(jì)60、70年代，各流域處于豐水期，尤其是60年代各流域水量豐沛，永定河上游、桑干河、洋河流域的徑流量分別比相應(yīng)的多年平均徑流量增加了119.39%、147.74%、72.25%；80年代以后，各流域的徑流量處于偏枯期，尤其是2000年以后各流域徑流量出現(xiàn)了銳減，永定河上游、桑干河、洋河流域的徑流量分別比相應(yīng)的多年平均值減少了84.7%、79.12%、75.15%。通過(guò)與降水、氣溫的年代際變化對(duì)比分析，在60、70年代，永定河上游和桑干河流域流域的徑流量變化的決定性因素是降水量，該時(shí)期降水降水較多，因此徑流量也較大；洋河流域60年代的降水量小于70年代，而60年代的徑流量明顯高于70年代，同時(shí)氣溫不斷升高，氣候干暖化趨勢(shì)顯著，因此徑流量逐漸減少；80年代以后，各流域的徑流量出現(xiàn)了大幅度的減少；90年代，洋河流域的降水量明顯高于多年平均值，而徑流量比多年平均值減少了2.89 億m3，減少率為23.14%，研究表明，90年代洋河流域的氣溫升高、大規(guī)模的生態(tài)建設(shè)以及水資源的過(guò)度開(kāi)發(fā)導(dǎo)致了該時(shí)期徑流量的銳減。

圖9 各流域徑流量距平值年代際變化對(duì)比圖

2.3 降水、氣溫、徑流的年際變化特征

2.3.1 趨勢(shì)性分析

永定河上游流域、桑干河流域和洋河流域的多年平均降水量分別為452、470、365 mm。根據(jù)擬合方程的斜率可知，永定河上游流域、桑干河流域、洋河流域的降水量減少量分別為8.21、10.37、6.92 mm/10 a。如圖10所示。通過(guò)Mann-Kendall 非參數(shù)檢驗(yàn)，各流域降水量統(tǒng)計(jì)值均為負(fù)值，分別為-1.18、-1.08、-1.1，|ZMK|<1.96,表明在1960-2012年間各流域的降水量均呈下降趨勢(shì)，而下降趨勢(shì)不顯著(β<0.05)。對(duì)1960-2012年15個(gè)氣象站點(diǎn)的年降水量進(jìn)行M-K趨勢(shì)性檢驗(yàn)分析(圖11)，結(jié)果表明，在β=0.05的顯著性水平，其中有11個(gè)氣象站點(diǎn)的年降水量呈減少下降趨勢(shì)，有4個(gè)氣象站點(diǎn)的年降水量呈增加趨勢(shì)，除五臺(tái)山站以外，其余站點(diǎn)的年降水量變化趨勢(shì)不顯著。

圖10 永定河上游流域、桑干河流域和洋河流域降水量年際變化

圖11 永定河上游流域年降水量M-k檢驗(yàn)結(jié)果空間分布

永定河上游流域、桑干河流域、洋河流域多年平均氣溫分別為6.79、6.22、5.74 ℃，根據(jù)擬合方程的斜率可知，氣溫分別增加了1.85、1.98、1.78 ℃，氣溫增加率分別為0.36、0.38、0.34 ℃/10 a。如圖12所示。通過(guò)Mann-Kendall 非參數(shù)檢驗(yàn)，各流域平均氣溫統(tǒng)計(jì)值均為正值，分別為4.86、4.36、4.91，表明在1960-2012年間各流域的氣溫均呈升高趨勢(shì)，而升高趨勢(shì)顯著(β>0.05顯著水平)。對(duì)1960-2012年15個(gè)氣象站點(diǎn)的年平均氣溫進(jìn)行M-K趨勢(shì)性檢驗(yàn)分析(圖13)，在β=0.05的顯著性水平，除河曲氣象站外，其余14個(gè)氣象站的年平均氣溫均呈顯著升高趨勢(shì)?？傮w來(lái)看，永定河上游流域年平均氣溫呈現(xiàn)顯著升高趨勢(shì)，與線性回歸分析的結(jié)果一致。

圖12 永定河上游流域、桑干河流域和洋河流域氣溫年際變化

圖13 永定河上游流域年氣溫M-k檢驗(yàn)結(jié)果空間分布

永定河上游流域、桑干河流域和洋河流域多年平均徑流量分別為6.02、3.18、2.89 億m3。通過(guò)線性回歸分析(圖14)，各流域的年徑流量均呈減少趨勢(shì)，徑流量的減少量分別為2.64、1.53、1.02 億m3/10 a。通過(guò)Mann-Kendall 非參數(shù)檢驗(yàn)，各流域年徑流量統(tǒng)計(jì)值均為負(fù)值，分別為-7.6、-7.77、-7.17，|ZMK|>2.58表明在1960-2012年間各流域的年平均徑流量均呈下降趨勢(shì)，而下降趨勢(shì)通過(guò)了可信度99%的顯著性檢驗(yàn)。

圖14 永定河上游流域、桑干河流域和洋河流域徑流量年際變化

2.3.2 突變性分析

通過(guò)累積距平分析，結(jié)果表明，總體上各流域的年降水量在1960-1980年呈增加趨勢(shì)，2000-2012年呈下降趨勢(shì)，豐、枯變化沒(méi)有規(guī)律性，且波動(dòng)較大；各流域的年降水量沒(méi)有顯著的突變性(圖15)?？傮w上各流域的年平均氣溫在1960-1990年呈下降趨勢(shì)，2000-2012年呈升高趨勢(shì)；永定河上游流域有兩個(gè)年平均氣溫突變點(diǎn)，分別為1988年、1993年，桑干河流域、洋河流域各有一個(gè)年平均氣溫突變點(diǎn)，分別為1996年、1986年(圖16)?？傮w上各流域的年徑流量在1960-1980年呈增加趨勢(shì)，1990-2012年呈下降趨勢(shì)，永定河上游流域、桑干河流域和洋河流域年徑流量的突變點(diǎn)分別為1981年、1982年和1983年。

圖15 永定河上游流域、桑干河流域和洋河流域降水量累積距平圖

圖16 永定河上游流域、桑干河流域和洋河流域年平均氣溫累積距平

圖17 永定河上游流域、桑干河流域和洋河流域年徑流量累積距平圖

2.4 流域徑流變化的主要影響因素

已有研究表明，20世紀(jì)80年代以后，永定河上游流域人類活動(dòng)對(duì)徑流的影響占主導(dǎo)作用，人類活動(dòng)和氣候變化對(duì)流域徑流的影響的貢獻(xiàn)占比約為7∶3，人類活動(dòng)多樣性決定了人類活動(dòng)對(duì)水文循環(huán)影響的復(fù)雜性[26]。20世紀(jì)80年代，永定河上游被列為國(guó)家水土保持重點(diǎn)治理區(qū)域，已有研究表明，80年代以后，流域耕地面積減少，林地增加，導(dǎo)致地表產(chǎn)流減少，而土地利用的變化對(duì)地表徑流的影響貢獻(xiàn)較小。20世紀(jì)50年代開(kāi)始，山西省(永定河上游)建設(shè)了大量的水利工程，大中型水庫(kù)的總庫(kù)容達(dá)到了52.8 億m3，蓄水容積占多年平均地表水資源量的75.33%，因此上游來(lái)水逐年減少，永定河出現(xiàn)了連續(xù)斷流現(xiàn)象[27,28]。另外，20世紀(jì)80年代以后，隨著經(jīng)濟(jì)社會(huì)的發(fā)展，該區(qū)域人口逐漸增加，因此工農(nóng)業(yè)、生活用水量也隨之增加，而地下水是該地區(qū)的主要供水水源，地下水超采嚴(yán)重，在該區(qū)域形成了地下漏斗，地下水水位下降嚴(yán)重，改變了流域的水文過(guò)程，使地表徑流減少[29]。

3 結(jié) 論

文中系統(tǒng)的分析了永定河上游流域及主要支流的水文氣象要素的變化特征以及主要的影響因素，為流域水資源的管理和開(kāi)發(fā)利用提供了技術(shù)支撐。形成以下主要結(jié)論。

(1)流域的多年平均氣溫和多年平均降水呈現(xiàn)地域性差異，總體上，流域的多年平均氣溫以從官?gòu)d水庫(kù)向河源地區(qū)逐漸遞減，降水量由東南向西北逐漸減少。

(2)從趨勢(shì)性角度分析，各流域的年平均氣溫均呈顯著升高趨勢(shì)，年徑流量均呈顯著下降趨勢(shì)，而年降水量的變化趨勢(shì)不明顯；從突變性的角度分析，各流域的年平均氣溫及徑流量均有突變點(diǎn)，而降水量無(wú)顯著突變點(diǎn)。

(3)流域氣溫四季分明，年內(nèi)降水量分配極不均衡，在20世紀(jì)60、70年代，徑流量變化的主導(dǎo)因素是氣象要素；而80年代以后，流域的徑流量變化的主要影響因素是由人類活動(dòng)引起的。

(4)影響永定河流域水資源變化的主要因素是水利工程攔蓄及高強(qiáng)度的社會(huì)取用水，如果恢復(fù)永定河自然的河流水文過(guò)程，一方面，應(yīng)通過(guò)恢復(fù)河流的自然屬性等措施，逐漸實(shí)現(xiàn)洪水的資源化，緩解流域水資源短缺問(wèn)題，另一方面，應(yīng)從轉(zhuǎn)變用水結(jié)構(gòu)，大力節(jié)水入手。

□