胡文峰，陳玲玲，姚俊強(qiáng)，何 清

（1.阜陽(yáng)師范大學(xué)，安徽 阜陽(yáng)236037；2.中國(guó)氣象局烏魯木齊沙漠氣象研究所，新疆 烏魯木齊830013；3.中亞大氣科學(xué)研究中心，新疆 烏魯木齊830013）

氣候變化及其對(duì)生態(tài)環(huán)境的影響已經(jīng)引起科學(xué)家和各國(guó)政府的關(guān)注與重視[1]。氣溫、徑流和降水是水文循環(huán)過(guò)程研究中的主要核心內(nèi)容，而地表徑流變化又是水文循環(huán)過(guò)程中最敏感的因素。河川徑流作為水循環(huán)的關(guān)鍵環(huán)節(jié)，是水資源科學(xué)管理、綜合開(kāi)發(fā)利用、優(yōu)化調(diào)度的重要依據(jù)[2，3]。

不少學(xué)者對(duì)新疆地區(qū)近50 a 的水文與氣候變化開(kāi)展了相關(guān)研究[4-6]，研究表明新疆天山西部地區(qū)在20 世紀(jì)80—90 年代氣候轉(zhuǎn)暖，造成冰川融化加速，部分地區(qū)降水變多，從而導(dǎo)致徑流增加的趨勢(shì)，1987—2000 年平均年徑流總量比1956—1986 年增加約7%。但對(duì)于天山北坡年徑流量在3 億m3以下的中小河流特點(diǎn)研究不夠深入。天山北坡的精河是艾比湖流域的一部分，是新疆北部生態(tài)環(huán)境最為敏感的地區(qū)之一，也是天山西部氣候變化的指示器。本文以天山北坡精河為研究對(duì)象，利用近50 a 來(lái)精河氣象、水文資料，研究精河流域徑流變化對(duì)氣候變化的響應(yīng)，揭示精河水文過(guò)程與氣候變化相互作用的耦合機(jī)制，為其生態(tài)系統(tǒng)保護(hù)和可持續(xù)發(fā)展提供理論依據(jù)與科學(xué)指導(dǎo)，為流域水資源開(kāi)發(fā)利用及制定“三條紅線”提供參考依據(jù)。

1 研究區(qū)概況

精河流域位于新疆準(zhǔn)噶爾盆地西南邊緣，婆羅科努山北麓（圖1），其主要發(fā)源于婆羅科努山北坡，向北注入艾比湖，河流集水面積1419 km2，流域面積2150 km2。水系呈扇形，河網(wǎng)密度為0.091，坡降為25‰～8.3‰，全長(zhǎng)114 km，徑流由冰雪融水、降水和地下水組成，其中降水補(bǔ)給、地下水補(bǔ)給、冰川和永久積雪融水補(bǔ)給分別占總補(bǔ)給量的70%、20%、10%；流域內(nèi)有129 條冰川，冰川面積91 km2，冰川總儲(chǔ)量54.598 億m3，冰雪年融水量為0.96 億m3，占河川徑流量的20.6%。據(jù)精河山口水文站近50 a 的水文資料統(tǒng)計(jì)，精河流域多年平均降水量156 mm，多年平均年徑流量4.74 億m3。春旱限制了灌溉面積擴(kuò)大，導(dǎo)致在夏季有相對(duì)豐富水資源[7]。精河水資源總量為12.9 億m3，水能資源理論儲(chǔ)藏量為16.27 萬(wàn)kW[8]。

圖1 精河流域水系分布

2 資料與方法

精河流域水文數(shù)據(jù)來(lái)自精河山口水文站，時(shí)間序列從1957—2012 年，氣象數(shù)據(jù)來(lái)自精河山口水文站和精河氣象站，包括逐月溫度和降水量資料；代表四季的尺度分別是春季為3—5 月，夏季為6—8 月，秋季為9—11 月，冬季為12 月—次年2 月。精河的出山口控制站為山口水文站，位于精河縣托里鄉(xiāng)山口村，建立于1956 年，其控制斷面在1966 年上遷了550 m，坐標(biāo)為85°55′E、44°22′N，海拔619.2 m。

內(nèi)陸河徑流一般產(chǎn)流位于人類活動(dòng)干擾較少的山區(qū)，而平原盆地產(chǎn)生徑流較少，故河流出山口的徑流變化能直觀反映氣候要素決定的天然徑流的變化[9]。精河流域的水文站位于出山口，站點(diǎn)控制斷面以上，其水資源受人類干擾較少，因此，該水文站點(diǎn)的水文資料能夠反映精河流域的天然徑流變化趨勢(shì)。該站點(diǎn)水文資料有50 a，其變化可以代表長(zhǎng)時(shí)間系列的變化趨勢(shì)。采用相關(guān)系數(shù)、小波分析等[10-12]方法進(jìn)行分析。

3 結(jié)果與分析

3.1 徑流年內(nèi)變化

精河流域徑流在年內(nèi)分配不均勻（表1），其中6—9 月為徑流較大的4 個(gè)月，是主要產(chǎn)流時(shí)期，占全年徑流量的74%，其中7、8 月的徑流量更是占全年徑流量的50%以上，7 月徑流量最大，其多年平均徑流量是最小月（2 月）徑流量的12.7 倍。因此，在6—9 月期間出現(xiàn)洪水的概率也最多，其中在1999 年7 月實(shí)測(cè)的最大洪峰流量達(dá)到390 m3/s，占當(dāng)年徑流量的3.63%。此外，12 月—次年4 月為徑流較小時(shí)期，在長(zhǎng)達(dá)6 個(gè)月的多年平均總徑流量?jī)H占其全年徑流量的17.46%，枯水期時(shí)間長(zhǎng)、徑流量小。

表1 精河流域1957—2012 年月平均徑流量變化

從季節(jié)變化來(lái)看，精河流域徑流量的季節(jié)變化顯著，集中于夏季，冬季徑流量極小。春汛在精河流域出現(xiàn)情況不多，春季徑流量?jī)H占全年徑流量的10.5%；夏季是豐水期，徑流量占全年徑流量的63.9%；冬季是枯水期，僅占全年徑流量的6.8%，主要是由于在冬季降水是以降雪為主，難以形成徑流，其徑流主要依靠淺層地下水補(bǔ)給，徑流量很?。磺锛緩搅髁扛哂诙簝杉?，占全年徑流量的18.7%。

3.2 徑流的年際及年代際變化

從精河流域1957—2012 年徑流量距平變化（圖2）可以看出，精河流域徑流量在20 世紀(jì)50 年代末、60 年代初為其豐水期，呈現(xiàn)增加趨勢(shì)；在60 年代后期呈下降趨勢(shì)，開(kāi)始進(jìn)入枯水期，說(shuō)明在此期間，精河徑流量呈減少趨勢(shì)；這種減少趨勢(shì)持續(xù)到80 年代初開(kāi)始回升，再進(jìn)入豐水期；期間徑流量波動(dòng)變化，在80 年代末期至90 年代前中期呈現(xiàn)緩慢下降趨勢(shì)，到90 年代后期呈上升趨勢(shì)，說(shuō)明精河流域開(kāi)始進(jìn)入又一個(gè)豐水期，進(jìn)入21 世紀(jì)以后，河川徑流波動(dòng)較大，豐枯相間變化。

圖2 精河流域1957—2012 年徑流量距平變化

20 世紀(jì)60 年代之后，精河徑流量較多年平均值偏高1%～3%，其中60 年代偏高1.04%（表2），70年代徑流量較多年平均值偏少3.18%，是近50 多年來(lái)徑流量最小的10 a，與60 年代相比，10 a 內(nèi)減少了0.2×108m3的來(lái)水量。80 年代以來(lái)徑流增加趨勢(shì)顯著，在此期間，1987 年和1988 年為最大連年豐水年，1988 年年徑流量達(dá)到6.05×108m3；90 年代持續(xù)增加，并趨于穩(wěn)定，高于多年平均徑流量的2.28%，但在1992 年出現(xiàn)了最大枯水年，年徑流量?jī)H為3.68×108m3，比多年平均值偏少22%。21 世紀(jì)以來(lái)精河徑流量有下降的趨勢(shì)，21 世紀(jì)最初10 a 是僅次于20 世紀(jì)70 年代的第二枯水年代，徑流量比20世紀(jì)90 年代減少了0.16×108m3的來(lái)水量。從20 世紀(jì)80 年代開(kāi)始精河流域徑流量增加，持續(xù)至90 年代，但21 世紀(jì)以來(lái)有減小的趨勢(shì)；進(jìn)一步分析發(fā)現(xiàn)自1981 年開(kāi)始徑流量呈現(xiàn)明顯增加趨勢(shì)，1981—2005 年平均年徑流總量比1957—1980 年增加了3.24%。也有研究[13]表明，天山北坡河流大多徑流量在20 世紀(jì)80 年代以后呈增加趨勢(shì)，增加幅度為10%～20%。

表2 精河徑流量年代際變化

3.3 徑流的多時(shí)間尺度分析

復(fù)值Morlet 小波可用于時(shí)間序列分析，小波系數(shù)實(shí)部能夠反映徑流在不同時(shí)間尺度的周期變化及在時(shí)間域中的分布，其變化趨勢(shì)與信號(hào)的起伏是一致的，中心值的大小可以反映波動(dòng)的振蕩強(qiáng)度，正值代表豐水期，負(fù)值表示枯水期。而Morlet 小波系數(shù)模表示不同時(shí)間尺度變化周期所對(duì)應(yīng)的能量密度在時(shí)間域上的分布，反映了對(duì)應(yīng)時(shí)間尺度的周期性強(qiáng)弱，模值的大小與周期顯著性的強(qiáng)弱成正比。

圖3 為精河流域年均徑流量序列的復(fù)值Morlet小波變換實(shí)部及對(duì)應(yīng)的模值，能夠明顯地反映出周期振蕩及時(shí)間尺度特征，圖4 是小波方差圖，可以看出主周期。從圖3、圖4 結(jié)合圖5 可以看出，精河流域年徑流量序列在21 a 和13 a 左右的振蕩周期最為明顯，其次是32 a 和9 a；而其中的21 a 和13 a時(shí)間尺度上的振蕩是全時(shí)域的。在21 a 時(shí)間尺度上存在著4 個(gè)明顯的豐水中心和4 個(gè)枯水中心，豐水中心分別是1957、1971、1985 年和1999 年，枯水中心分別是1964、1978、1993 年和2006 年，流域徑流變化的平均周期為14 a 左右，大約經(jīng)歷了4 個(gè)周期的豐—枯變化；13 a 尺度上存在6 個(gè)豐水中心和6個(gè)枯水中心，其中豐水中心為1964、1972、1981、1989、1998 年和2007 年，枯水中心為1959、1968、1977、1985、1994 年和2003 年，流域徑流變化的平均周期為8 a 左右，大約經(jīng)歷了6 個(gè)周期的豐—枯轉(zhuǎn)換期。9 a 時(shí)間尺度上的振蕩在20 世紀(jì)50 年代后期和80 年代之后時(shí)段內(nèi)相對(duì)較強(qiáng)，其他時(shí)段較弱，存在著8 個(gè)豐水中心和8 個(gè)枯水中心，平均周期為6 a 左右，大約經(jīng)歷了9 個(gè)周期的豐—枯轉(zhuǎn)換期。此外，精河流域年徑流量還存在準(zhǔn)32 a 的振蕩周期。

3.4 徑流變化對(duì)氣候變化的響應(yīng)

圖3 精河流域年均徑流量的小波分析

圖4 精河流域年徑流量變化小波方差圖

圖5 精河流域年徑流量變化主要時(shí)間尺度小波實(shí)部過(guò)程線

河流出山徑流的年際變化與山區(qū)氣候變化有密切的關(guān)系。以精河山口水文站的氣溫和降水時(shí)間序列資料作為山區(qū)的主要?dú)夂蛞蜃?，以出山天然徑流量?lái)研究徑流變化對(duì)氣候變化的響應(yīng)。氣溫和降水是基本的氣候因子，降水直接影響徑流，是產(chǎn)生徑流的直接原因，而氣溫通過(guò)影響蒸散發(fā)間接影響徑流，兩者共同影響著流域徑流過(guò)程。

1957—2012 年精河流域山區(qū)的氣溫、降水與出山徑流的相關(guān)分析（圖6、圖7）可知，氣溫與徑流的相關(guān)系數(shù)為0.13，呈不顯著的相關(guān)關(guān)系。而降水與出山徑流呈正相關(guān)關(guān)系，降水曲線的波峰、波谷走向與徑流量的波動(dòng)變化基本吻合，相關(guān)系數(shù)為0.40，通過(guò)了P＜0.01 的顯著性水平檢驗(yàn)，表明兩者之間有著極顯著的正相關(guān)，山區(qū)降水越多，出山徑流量越大；反之，降水越少?gòu)搅髁恳簿拖鄳?yīng)地越小。這也在精河流域得到體現(xiàn)，在1958、1964、1969、1978、1981、1988、1990、1993、1998、2000、2005 年等降水較多，對(duì)應(yīng)的出山徑流也越大；而1961、1967、1976、1985、1989、1992、1994、2003 年等典型年降水較少，則相應(yīng)年份的出山徑流量較小[14]。

圖6 精河流域山區(qū)氣溫、降水和徑流的變化

河川徑流對(duì)氣候變化的響應(yīng)過(guò)程雖然十分復(fù)雜的，但通過(guò)建立模型，近似地去模擬原型仍然是一種可行的途徑[15-16]?；谫Y料的代表性和完整性，對(duì)精河山口水文站1957—2012 年近50 a 的徑流、氣溫和降水資料進(jìn)行分析，并利用1957—2000 年的徑流、氣溫和降水資料初步建立徑流與降水、氣溫的多變量自回歸CAR 模型：

圖7 精河流域山區(qū)氣溫、降水和徑流的相關(guān)性

式中，t 為時(shí)間序列編號(hào)，t＞1；R 為年徑流量，P 為年降水量，T 為年均氣溫。圖8 為1957—2000 年精河流域的實(shí)測(cè)徑流量和上式模擬徑流量的比較，發(fā)現(xiàn)兩者趨勢(shì)基本一致，該模型通過(guò)了F-test 檢驗(yàn)和殘差檢驗(yàn)，說(shuō)明建立的多變量自回歸CAR 模型是合理的。

圖8 1957—2000 年精河流域的實(shí)測(cè)徑流量和上式模擬徑流量的比較

為了進(jìn)一步驗(yàn)證模型的適用性，采用精河出山口水文站2001—2012 年年徑流量觀測(cè)資料。圖9 為年徑流量的模擬值和觀測(cè)值的相對(duì)誤差圖和觀測(cè)值與模擬值的擬合圖?？梢钥闯瞿陱搅髁康哪M值和觀測(cè)值的平均絕對(duì)誤差為0.84×108m3，平均相對(duì)誤差為6.54%，均方根誤差為0.039?？梢钥闯鲇^測(cè)值和模擬值的擬合效果很好，整體趨勢(shì)吻合，偏離程度很小。用CAR 模型模擬河流年徑流量誤差在可接受范圍內(nèi)，可以利用該模型對(duì)精河流域年徑流量進(jìn)行模擬。

4 結(jié)論

（1）精河徑流在年內(nèi)分配不均，季節(jié)變化明顯，夏季集中，枯水期長(zhǎng)且枯季徑流量小，6—9 月為豐水期，占全年徑流量的74%。從20 世紀(jì)80 年代開(kāi)始，河川徑流量增加，持續(xù)至90 年代，21 世紀(jì)以來(lái)有減小的趨勢(shì)。1981 年開(kāi)始徑流量出現(xiàn)了明顯增加趨勢(shì)，1981—2005 年平均年徑流總量比1957—1980年增加了3.24%。

（2）精河流域年徑流量序列在21 a 和13 a 左右的振蕩周期最為明顯，其次是32 a 和9 a，21 a 和13 a 時(shí)間尺度上的振蕩是全時(shí)域的，其中，在21 a時(shí)間尺度上存在4 個(gè)明顯的豐水中心和4 個(gè)枯水中心，13 a 尺度上存在6 個(gè)豐水中心和6 個(gè)枯水中心。

（3）利用精河流域徑流、氣溫和降水資料，建立了徑流與降水和氣溫的多變量自回歸CAR 模型，模擬的年徑流量平均相對(duì)誤差為6.54%，均方根誤差為0.039，模擬效果較好。

由于氣候和下墊面條件是影響流域水循環(huán)的兩個(gè)重要方面，在不同的流域二者對(duì)流域徑流影響作用存在差異，在本文中僅僅只考慮了氣候因素中降水和氣溫對(duì)精河徑流的影響，而沒(méi)有考慮潛在蒸散發(fā)量、日照時(shí)數(shù)、近地表平均風(fēng)速等因素對(duì)徑流的影響，下墊面條件的變化在本文中也讓沒(méi)有考慮，這也將是下一步研究的內(nèi)容。

圖9 2001—2012 年年徑流量的模擬值和觀測(cè)值的誤差分析