魏瀟娜，龍愛(ài)華,，尹振良，於嘉聞

(1.中國(guó)水利水電科學(xué)研究院流域水循環(huán)模擬與調(diào)控國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京 100038；2.石河子大學(xué)水利建筑工程學(xué)院，新疆 石河子 832003；3.中國(guó)科學(xué)院西北生態(tài)環(huán)境資源研究院內(nèi)陸河流域生態(tài)水文重點(diǎn)實(shí)驗(yàn)室，甘肅 蘭州 730000)

我國(guó)西北地區(qū)的內(nèi)陸河大多發(fā)源于高山冰川[1]，冰川融水是西北干旱區(qū)流域重要的徑流來(lái)源，在西北地區(qū)水資源的開(kāi)發(fā)利用中占有重要地位[2]。冰川作為“固體水庫(kù)”，對(duì)維持區(qū)域生態(tài)穩(wěn)定和調(diào)節(jié)河川徑流具有重要意義[3-4]，已有研究表明西北干旱區(qū)冰川融水量占全國(guó)冰川融水量的39%，占3省區(qū)(甘肅、新疆、青海)地表徑流的14%[5]。近50年來(lái)，新疆以氣溫明顯高于世界和全國(guó)的平均增幅變暖[3, 6]，其冰川整體上呈加速退縮之勢(shì)[3]，大大改變了區(qū)域的徑流組成、水資源量[7]等。氣候變化背景下，以典型流域?yàn)槔?，?duì)冰雪補(bǔ)給型河流出山徑流、冰川徑流的變化以及冰川徑流對(duì)出山徑流的貢獻(xiàn)開(kāi)展評(píng)估分析對(duì)揭示流域整體變化具有重要作用。塔里木河流域是我國(guó)最大的內(nèi)陸河流域，其河川徑流主要以冰川積雪融水補(bǔ)給為主，冰川融水徑流量約占總徑流的40%[7]。和田河是塔里木河四大主要源流之一，也是唯一一條完整發(fā)育著穿越塔克拉瑪干沙漠綠色生態(tài)走廊的河流，在變暖變濕趨勢(shì)下，流域生態(tài)環(huán)境變得異常脆弱[8]。因此開(kāi)展氣候變化背景下和田河流域冰川徑流變化研究，對(duì)保障流域社會(huì)經(jīng)濟(jì)發(fā)展與生態(tài)安全顯得十分迫切[9-11]。但是，長(zhǎng)期以來(lái)，學(xué)界對(duì)和田河流域冰川徑流變化的研究相對(duì)薄弱。

目前冰川徑流的研究方法大致分為直接觀測(cè)法、冰川物質(zhì)平衡法、水量平衡方程法、水化學(xué)示蹤法、水文模型法等5種，水文模型法是最常用的方法[12-13]。趙求東等[14]利用VIC(variable infiltration capacity)分布式水文模型估算得出天山南坡木扎提河流域冰川徑流占總流量的66.6%；李洪源等[15]利用包含冰雪消融模塊的寒區(qū)分布式水文模型——SPHY(spatial processes in hydrology)模型對(duì)疏勒河流域的徑流過(guò)程進(jìn)行定量模擬，得出冰川徑流占總徑流的年均比例為30.5%。SWAT(soil and water assessment tool)模型能夠充分利用遙感數(shù)據(jù)提供的空間信息，模擬復(fù)雜流域條件下氣候變化、土地利用變化和不同管理措施等對(duì)流域水量、水質(zhì)的影響[16-18]。Adnan等[16]在中國(guó)青藏高原南部南柯湖冰川和非冰川區(qū)子流域進(jìn)行的SWAT模型適應(yīng)性評(píng)價(jià)表明，常規(guī)的SWAT模型在有冰川的流域模擬效果并不理想，將SWAT模型與冰川模塊相結(jié)合才能更好地模擬冰川地區(qū)的流域徑流。近年來(lái)，有學(xué)者試圖在SWAT模型中嵌入冰川模塊，從而提高SWAT模型在高冰川覆蓋率地區(qū)的適用性。如Luo等[17]提出了冰川算法并將其嵌入SWAT模型中，在瑪納斯河流域的應(yīng)用中驗(yàn)證了該方法的可靠性；Yin等[18]利用嵌入冰川模塊的SWAT模型定量估算了葉爾羌河流域冰川徑流對(duì)出山徑流的貢獻(xiàn)，認(rèn)為嵌入冰川模塊的SWAT模型在該流域有較好的適用性。

本文針對(duì)和田河流域上游冰川徑流構(gòu)建了嵌入冰川模塊的SWAT模型來(lái)模擬分析冰川徑流變化規(guī)律及其對(duì)氣候變化的響應(yīng)，以期為塔里木河水資源利用與規(guī)劃提供參考。

1 研究區(qū)概況與數(shù)據(jù)來(lái)源

1.1 研究區(qū)概況

和田河流域位于新疆塔里木盆地的南部，流域面積約7.95萬(wàn)km2[19]，多年平均出山徑流量為44.90億m3。流域上游由玉龍喀什河(簡(jiǎn)稱玉河)與喀拉喀什河(簡(jiǎn)稱喀河)兩支流組成(圖1)，出山口分別設(shè)有同古孜洛克與烏魯瓦提水文站。主河流玉河發(fā)源于昆侖山，河長(zhǎng)504 km，多年(1967—2017年)平均出山徑流量為22.64億m3?？影l(fā)源于喀喇昆侖山，河長(zhǎng)808 km，多年(1967—2017年)平均出山徑流量為22.26億m3。兩支流于闊什拉什(88°33′E，38°05′N(xiāo))匯合后始稱和田河，并進(jìn)入塔克拉瑪干沙漠，和田河長(zhǎng)319 km[19]，最終經(jīng)由肖塔水文站匯入塔里木河。和田河流域氣候干旱少雨，寒冷干燥，常年處于極端干旱狀態(tài)[20]。和田河及其支流屬于冰雪融水及降水混合補(bǔ)給河流，冰川融水補(bǔ)給量占年平均徑流的40%以上[7]，夏季(6—8月)徑流占年總量的70%以上。據(jù)第一次冰川編目，和田河流域上游的冰川面積為5 687.35 km2，冰川覆蓋率為15.84%，其中玉河流域與喀河流域冰川覆蓋率分別為22.72%和11.13%。

圖1 研究區(qū)位置與水文站、氣象站以及冰川分布Fig.1 Location of study area and distributions of hydrological and meteorological stations as well as glacier

1.2 數(shù)據(jù)來(lái)源

a.DEM數(shù)據(jù)。在地理空間數(shù)據(jù)云平臺(tái)(http:∥www.gscloud.cn/)獲取了90 m分辨率的DEM數(shù)據(jù)，并對(duì)數(shù)據(jù)進(jìn)行了投影變換處理，將兩個(gè)出山口水文站作為流域出口控制斷面得到和田河上游的DEM數(shù)據(jù)。通過(guò)DEM數(shù)據(jù)進(jìn)行水流方向、集水面積的計(jì)算以及河網(wǎng)與子流域的劃分[21-22]。

b.土壤數(shù)據(jù)。土壤類型與土壤的物理屬性數(shù)據(jù)來(lái)自中國(guó)科學(xué)院南京土壤研究所發(fā)布的1∶1 000 000土壤數(shù)據(jù)集，其中土壤容重(SOL_BD)、飽和導(dǎo)水率(SOL_K)和田間持水量(SOL_AWC)等參數(shù)通過(guò)SPAW(soil-plant-atmosphere-water)模型中的土壤水文特性模塊計(jì)算得到。根據(jù)研究區(qū)情況，共劃分了10種土壤類型。

c.土地利用數(shù)據(jù)。土地利用數(shù)據(jù)來(lái)自地理空間數(shù)據(jù)云(http:∥www.gscloud.cn/)提供的2015年全國(guó)土地利用數(shù)據(jù)，分辨率為30 m × 30 m。研究區(qū)的數(shù)據(jù)經(jīng)過(guò)SWAT重分類為12種數(shù)據(jù)。

d.冰川編目數(shù)據(jù)。冰川數(shù)據(jù)來(lái)自國(guó)家冰川凍土沙漠科學(xué)數(shù)據(jù)中心(http:∥www.ncdc.ac.cn/)中國(guó)第一次與第二次冰川編目數(shù)據(jù)，第一次冰川編目數(shù)據(jù)主要依據(jù)1987—2004年航測(cè)地形圖制作而成，第二次冰川編目數(shù)據(jù)主要由2006—2013年Landsat TM/ETM+和ASTER遙感影像提取而成。本研究的模型率定期(1967—1992年)采用第一次冰川編目數(shù)據(jù)，驗(yàn)證期(1993—2017年)采用第二次冰川編目數(shù)據(jù)。

e.水文氣象數(shù)據(jù)。本文采用中國(guó)地面降水日值0.5°×0.5°格點(diǎn)數(shù)據(jù)集(V2.0)、中國(guó)地面氣溫日值0.5°×0.5°格點(diǎn)數(shù)據(jù)集(V2.0)，時(shí)間序列為1961—2017年，其中包括了日最高、平均、最低氣溫與日降水?dāng)?shù)據(jù)，提取20個(gè)格點(diǎn)作為虛擬氣象站(圖1)，格點(diǎn)數(shù)據(jù)經(jīng)過(guò)校正后作為模型的輸入數(shù)據(jù)。水文數(shù)據(jù)來(lái)自和田河流域管理局，包括烏魯瓦提與同古孜洛克水文站1961—2017年的逐月徑流數(shù)據(jù)。

2 研究方法

2.1 冰川融化算法

采用修正的溫度指數(shù)法對(duì)冰川融化過(guò)程進(jìn)行模擬[18，23]。根據(jù)消融因子與氣溫之間的線性關(guān)系計(jì)算冰川融化量，基于太陽(yáng)輻射和地形兩個(gè)因素的影響對(duì)空間異質(zhì)性進(jìn)行修正[23]。計(jì)算公式為

(1)

式中：M為日消融量，mm；FM為溫度消融因子，mm/(℃·d)；ri為輻射消融因子，m2·mm/(W·d·℃)；Ip為潛在太陽(yáng)總輻射，W/m2；T為日平均氣溫，℃；TMT為冰川消融閾值溫度，℃。

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

Mann-Kendall(M-K)檢驗(yàn)法作為非參數(shù)檢驗(yàn)方法，能夠很好地揭示時(shí)間序列的變化趨勢(shì)[24-25]。已知某時(shí)間序列數(shù)據(jù)，可計(jì)算得到兩組數(shù)據(jù)序列UF和UB。若UF>0，序列呈上升趨勢(shì)；若UF<0，序列呈下降趨勢(shì)。若UF和UB曲線交點(diǎn)位于顯著性水平線之間，交點(diǎn)對(duì)應(yīng)時(shí)間為突變點(diǎn)，該方法廣泛應(yīng)用于水文和氣象時(shí)間序列的趨勢(shì)變化分析[26-28]。

2.3 模型評(píng)價(jià)指標(biāo)

選用納什效率系數(shù)(ENS)、均方根誤差與實(shí)測(cè)值標(biāo)準(zhǔn)差的比值(RRS)及百分比誤差(PBIAS)對(duì)模型進(jìn)行評(píng)價(jià)[18]。ENS越接近1，表明模擬值對(duì)實(shí)測(cè)值的擬合精度越高，一般認(rèn)為月尺度上ENS>0.6時(shí)模擬結(jié)果較可靠[29]。當(dāng)0.75≤ENS≤1、0≤RRS≤0.5、|PBIAS|<10%時(shí)，認(rèn)為模型模擬效果為優(yōu)[24]。

3 結(jié)果與分析

3.1 模型評(píng)價(jià)

分別選取烏魯瓦提水文站、同古孜洛克水文站逐月徑流數(shù)據(jù)，進(jìn)行喀河流域與玉河流域的率定與驗(yàn)證，選取1961—1966年為模型預(yù)熱期，1967—1992年為模型率定期，1993—2017年為模型驗(yàn)證期。

3.1.1模型參數(shù)率定

采用SUFI2算法進(jìn)行參數(shù)敏感性分析[30-31]，選取最敏感的參數(shù)進(jìn)行率定，最終根據(jù)參數(shù)的物理意義和該流域的水文過(guò)程進(jìn)行參數(shù)的確定與調(diào)整，參數(shù)率定結(jié)果見(jiàn)表1。

表1 參數(shù)率定結(jié)果

3.1.2模型評(píng)價(jià)結(jié)果

將率定期的最優(yōu)參數(shù)代入SWAT模型和嵌入冰川模塊的SWAT模型進(jìn)行驗(yàn)證，模型評(píng)價(jià)結(jié)果如表2所示，率定期與驗(yàn)證期月徑流實(shí)測(cè)值與模擬值對(duì)比結(jié)果如圖2所示。結(jié)果表明，嵌入冰川模塊后，率定期與驗(yàn)證期ENS均在0.75以上，PBIAS均在10%以內(nèi)，RRS均在0.5以內(nèi)；模擬的出山徑流與實(shí)測(cè)徑流變化過(guò)程基本一致，其中玉河流域率定期峰值處模擬值偏低，主要原因是氣溫偏低導(dǎo)致冰川融水模擬值偏低，進(jìn)而導(dǎo)致峰值模擬值偏小。以上均表明嵌入冰川模塊的SWAT模型在和田河流域適用性較好。

圖2 1967—2017年月出山徑流量模擬結(jié)果

表2 模型評(píng)價(jià)結(jié)果

3.2 氣候與水文過(guò)程

3.2.1氣候特征

1961—2017年，玉河與喀河流域年平均氣溫均呈顯著上升趨勢(shì)(p<0.01)，上升速率分別為0.028 ℃/a和0.021 ℃/a；玉河流域降水量呈顯著增加趨勢(shì)，喀河流域降水量呈不顯著增加趨勢(shì)，降水量增長(zhǎng)速率分別為0.997 mm/a和0.646 mm/a(圖3)。就季節(jié)變化而言，玉河、喀河流域均在秋季升溫最明顯，夏季降水量增加最快，但玉河流域氣溫、降水量的增長(zhǎng)速率普遍高于喀河流域(表3)。

表3 和田河流域上游平均氣溫、降水量增長(zhǎng)速率Table 3 Average temperature and precipitation increase rates in the upper Hotan River Basin

3.2.2水文特征

1967年以來(lái)，玉河流域出山徑流和冰川徑流模擬值均呈顯著增加趨勢(shì)(p<0.01)，增長(zhǎng)速率分別為0.248億m3/a和0.077億m3/a；喀河流域出山徑流和冰川徑流模擬值均呈不顯著增加趨勢(shì)，增長(zhǎng)速率分別為0.061億m3/a和0.035億m3/a。玉河流域冰川徑流量的增長(zhǎng)速率遠(yuǎn)小于出山徑流量的增長(zhǎng)速率，冰川徑流增加對(duì)出山徑流增加的貢獻(xiàn)率約為31.1%，冰川徑流對(duì)出山徑流的貢獻(xiàn)(簡(jiǎn)稱冰川徑流貢獻(xiàn)率)呈顯著下降趨勢(shì)；而喀河流域冰川徑流量的增長(zhǎng)速率超過(guò)出山徑流量增長(zhǎng)速率的50%，冰川徑流增加對(duì)出山徑流增加的貢獻(xiàn)率約為57.38%，冰川徑流貢獻(xiàn)率呈不顯著的上升趨勢(shì)(圖4)。

圖4 年徑流量年際變化趨勢(shì)Fig.4 Annual runoff trends

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

分析表明，玉河、喀河流域冰川徑流量與氣溫的相關(guān)系數(shù)分別為0.70和0.43，均通過(guò)0.01的顯著性水平檢驗(yàn)；玉河、喀河流域冰川徑流量與降水量的相關(guān)系數(shù)分別為0.42和0.32，分別通過(guò)0.01和0.05的顯著性水平檢驗(yàn)。由此可知，玉河、喀河流域冰川徑流對(duì)氣溫的變化均更為敏感，而玉河流域的冰川覆蓋率約為喀河流域的2倍，使得玉河流域冰川徑流對(duì)氣溫變化的響應(yīng)更敏感，一定程度上解釋了玉河、喀河流域氣溫均呈顯著上升趨勢(shì)時(shí)，冰川徑流卻表現(xiàn)出顯著性不同的變化趨勢(shì)。冰川徑流量年代變化趨勢(shì)如圖5所示(圖中1970 s表示20世紀(jì)70年代)，從圖5可以看出，玉河流域冰川徑流量與氣溫變化規(guī)律基本一致；喀河流域冰川徑流量變化與氣溫變化趨勢(shì)之間關(guān)系并不明顯，但與降水量的變化規(guī)律相反，與趙求東等[14]得到的降水量增加對(duì)冰川消融有抑制作用的結(jié)論基本一致。由此可知，在氣溫、降水的共同作用下，氣溫上升加速了冰川的消融，但冰川萎縮與降水量增加均對(duì)冰川消融產(chǎn)生了抑制作用。

圖5 降水量、冰川徑流量、氣溫年代際變化趨勢(shì)Fig.5 Decadal trends of precipitation, glacier runoff, and air temperature

冰川徑流主要發(fā)生在4—10月，對(duì)冰川徑流量與季節(jié)性氣候的相關(guān)性分析(圖6)表明，玉河流域冰川徑流量與夏季氣溫的相關(guān)系數(shù)為0.76，高于與秋季氣溫的相關(guān)系數(shù)(0.64)，與夏季降水量的相關(guān)系數(shù)為0.38；喀河流域冰川徑流量與夏季氣溫的相關(guān)相關(guān)系數(shù)為0.67，與秋季氣溫的相關(guān)系數(shù)為0.40。氣溫是冰川消融最關(guān)鍵的因子，由于夏季氣溫基本在0℃以上，因此夏季氣溫對(duì)冰川徑流產(chǎn)生了較大的影響。同時(shí)，玉河、喀河流域海拔較高，秋季氣溫基本在0℃以下，秋季降水通過(guò)降雪的形式降落到地表，促進(jìn)冰川的積累，使得玉河、喀河流域的冰川徑流受秋季氣溫的影響也較大。

圖6 冰川徑流與氣象因子相關(guān)系數(shù)Fig.6 Correlation coefficient between glacier runoff and meteorological factors

4 討 論

4.1 模型的適用性

從模型評(píng)價(jià)結(jié)果可以看出，未嵌入冰川模塊的SWAT模型模擬結(jié)果與實(shí)測(cè)徑流相比，誤差較大，嵌入冰川模塊的SWAT模型在玉河與喀河流域的模擬效果較好。玉河流域多年(1967—2017年)平均年徑流量模擬值為22.63億m3，其中冰川徑流量為11.02億m3，冰川徑流貢獻(xiàn)率為48.70%(圖4)，與陳亞寧等[32]估算的塔里木河流域的冰川融水貢獻(xiàn)率(50%)較為接近。喀河流域多年(1967—2017年)平均年徑流量模擬值為20.92億m3，冰川徑流量為9.51億m3，冰川徑流貢獻(xiàn)率為45.46%(圖4)，與楊針娘[33]估算的(45.7%)較為接近。因此，改進(jìn)的SWAT模型不僅可以對(duì)高冰川覆蓋率流域的水文過(guò)程進(jìn)行較精確的模擬，也可為高寒山區(qū)水文過(guò)程研究提供參考。

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

在氣候暖濕化背景下，和田河流域的冰川徑流量呈現(xiàn)出與塔里木盆地其他冰川流域類似的增加趨勢(shì)[34-35]。但由于玉河流域與喀河流域氣溫、降水的增加速率不同，雖然冰川徑流量與出山徑流量均呈增加趨勢(shì)，但冰川徑流占出山徑流的比例呈現(xiàn)不同的變化趨勢(shì)，表明和田河流域出山徑流組分發(fā)生了變化。其中，玉河流域冰川徑流貢獻(xiàn)率下降的原因之一可能是冰川徑流量的上升速率變緩，冰川徑流量可能已到達(dá)拐點(diǎn)附近。而喀河流域的冰川徑流貢獻(xiàn)率的上升趨勢(shì)表明，冰川徑流量可能還處在加速上升的階段，但需進(jìn)一步確證研究。

近幾十年來(lái)，受全球變暖的影響，冰川普遍萎縮。據(jù)第一期、第二期冰川編目數(shù)據(jù)，和田河上游冰川面積減少了16.72%；1993—2017年和田河兩支流實(shí)測(cè)出山徑流序列的年際變差系數(shù)相較于1967—1992年增加了21%，表明冰川萎縮降低了冰川的多年調(diào)節(jié)作用[2]，冰川的調(diào)豐補(bǔ)枯作用正在發(fā)生著顯著變化；同期，徑流年內(nèi)分配不均勻系數(shù)減少了6%，這是由于氣溫上升引起冰川融化提前，一定程度上緩解了過(guò)去徑流年內(nèi)分配不均的問(wèn)題。

玉河與喀河同屬和田河流域的兩條支流，但各自徑流對(duì)氣候變化的響應(yīng)卻存在一定差異，主要原因可能是兩條支流發(fā)源于不同的山系?？錾脚c昆侖山氣候環(huán)境不同，其冰川結(jié)構(gòu)(大小、分布)、動(dòng)態(tài)變化規(guī)律也存在一定差異。加之玉河流域?qū)崪y(cè)徑流年內(nèi)分配不均勻系數(shù)(0.53)高于喀河流域(0.45)，從水資源高效利用與生態(tài)保護(hù)的角度，應(yīng)該采取不同的人工調(diào)控措施對(duì)兩支流的徑流進(jìn)行調(diào)節(jié)[36-39]。

5 結(jié) 論

a.利用嵌入冰川模塊的SWAT模型模擬和田河上游玉河流域與喀河流域的冰川徑流，率定期與驗(yàn)證期納什效率系數(shù)均在0.75以上，百分比誤差均在10%以內(nèi)，均方根誤差與實(shí)測(cè)值標(biāo)準(zhǔn)差的比值均在0.5以內(nèi)，表明嵌入冰川模塊的SWAT模型在玉河流域與喀河流域適用性較好。

b.1967年以來(lái)，玉河流域氣溫、降水與冰川徑流均呈顯著增加趨勢(shì)，冰川徑流貢獻(xiàn)率呈顯著下降趨勢(shì)，冰川徑流增加對(duì)出山徑流增加的貢獻(xiàn)率約為31.1%，冰川徑流對(duì)出山徑流的影響呈下降的趨勢(shì)；喀河流域氣溫呈顯著上升趨勢(shì)，降水量呈不顯著增加趨勢(shì)，冰川徑流、冰川徑流貢獻(xiàn)率均呈不顯著上升趨勢(shì)，冰川徑流增加對(duì)出山徑流增加的貢獻(xiàn)率約為57.38%，冰川徑流對(duì)出山徑流的影響呈上升的趨勢(shì)。

c.1993—2017年和田河兩支流實(shí)測(cè)出山徑流序列的年際變差系數(shù)相較于1967—1992年增加了21%，同期年內(nèi)徑流分配不均勻系數(shù)減少了6%，冰川面積萎縮降低了冰川的多年調(diào)節(jié)作用，一定程度上緩解了徑流年內(nèi)分配不均的問(wèn)題。