聶曉 丁玲玲

摘要 采用SWAT分布式水文模型，結(jié)合漢江上游流域1961—1990年實(shí)測(cè)氣象數(shù)據(jù)和CMIP5多模式RCP 2.6、RCP 4.5 和 RCP 8.5情景下的輸出數(shù)據(jù)集，在對(duì)漢江上游流域徑流模擬檢驗(yàn)的基礎(chǔ)上，分析未來(lái)氣候變化對(duì)漢江上游流域水資源的響應(yīng)。結(jié)果表明，SWAT模型能較好地模擬漢江上游流域徑流變化。2011—2100年不同氣候變化情景下，漢江上游流域多年平均降水、徑流與基準(zhǔn)期（1961—1990年）相比均呈增加趨勢(shì)，增加幅度從大到小依次為RCP 8.5、RCP 4.5、RCP 2.6，且降水量的增幅大于徑流的增幅。3種情景下的多年平均月降水、月徑流總體呈增加趨勢(shì)，且在冬季（12月至次年2月）枯水期增加幅度較為明顯。

關(guān)鍵詞 降水； 徑流； RCP氣候情景； SWAT模型； 漢江上游流域

中圖分類號(hào) X16 文獻(xiàn)標(biāo)識(shí)碼 A 文章編號(hào) 0517-6611（2017）04-0058-03

Response of Future Water Resources to Climate Change in the Upper Reaches of the Hanjiang Basin

NIE Xiao，DING Ling-ling （Institute of the Han River，Hubei University of Arts and Science，Xiangyang，Hubei 441053）

Abstract Based on the the measured meteorological data from 1961 to 1990 and the output data set under RCP scenarios （RCP 2.6，RCP 4.5 and RCP 8.5） from the phase 5 of the Coupled Model Intercomparison Project （CMIP5），the verified SWAT distributed hydrological model was used to analyze the impact of future climate change on the water resources of the upper reaches of the Hanjiang Basin.The results show that the SWAT model can well simulate the runoff variation of the upper reaches of the Hanjiang Basin.Compared with the baseline period （1961-1990），the average annual precipitation and runoff during 2011-2100 under different climate change scenarios all increase with the increase range of RCP 8.5>RCP 4.5>RCP 2.6，but the increment amplitude of precipitation will be more than that of runoff.The mean monthly precipitation and runoff all show an increasing tendency under three RCP scenarios，and the increase range is obvious in winter which is the dry period of the upper reaches of the Hanjiang Basin.

Key words Precipitation；Runoff；RCP scenarios；SWAT model；The upper reaches of the Hanjiang Basin

近年來(lái)，關(guān)于未來(lái)氣候變化對(duì)流域水資源影響的研究較多，且大多采用氣候情景驅(qū)動(dòng)水文模型的方法[1-3]。氣候情景的定義主要采用2種方法：一種是假定不同的降雨量和氣溫變化情景或組合，如假定年降水量變化±5%、±10%、±15%等，年平均氣溫變化±1、±2、±3 ℃等，以及將假定的不同降水量和氣溫變化進(jìn)行組合[4-5]；另一種是基于全球氣候模式（GCM）或區(qū)域氣候模式（REGCM）在各種溫室氣體排放情景下的氣候數(shù)據(jù)[6-9]，將其輸入水文模型中，進(jìn)行氣候變化影響的模擬量化。

CMIP5成果是聯(lián)合國(guó)政府間氣候變化專門委員會(huì)（IPCC）第五次評(píng)估報(bào)告的重要基礎(chǔ)，參與試驗(yàn)的全球多種氣候模型有了改進(jìn)和發(fā)展，并且形成了最新一代排放情景，為國(guó)際耦合模式的評(píng)估和后續(xù)發(fā)展提供了重要的平臺(tái)[10]。漢江上游流域?yàn)槟纤闭{(diào)中線工程水源區(qū)，未來(lái)水源區(qū)的水資源變化是流域水資源調(diào)配及跨流域調(diào)水可靠性的重要依據(jù)[11]。鑒于此，筆者以漢江上游流域?yàn)檠芯繉?duì)象，采用最新的CMIP5 RCP情景下的氣候數(shù)據(jù)驅(qū)動(dòng)SWAT分布式水文模型，模擬分析未來(lái)氣候變化對(duì)漢江上游流域水資源的影響，以期為流域水資源管理及跨流域調(diào)水提供科學(xué)依據(jù)。

1 數(shù)據(jù)與方法

1.1 研究區(qū)概況

漢江屬于長(zhǎng)江最大的支流，該研究流域?yàn)闈h江丹江口水庫(kù)以上的上游流域，位于106°15′～114°20′ E、30°10′～34 ° 20′ N，長(zhǎng)約925 km，占漢江總長(zhǎng)的59%，集水面積為95 200 km2。漢江上游河段位于秦嶺、大巴山之間，流域呈羽葉狀，干流走向?yàn)闁|西方向。漢江上游屬于北亞熱帶季風(fēng)氣候，多年平均氣溫14.6 ℃，多年平均降水量904 mm，降水年內(nèi)分布不均，多集中在7—9月，多年平均徑流量是367.8億m3。水面蒸發(fā)變化在700～1 000 mm，其分布特點(diǎn)是西南少，東北多。

1.2 數(shù)據(jù)資料

流域水文氣象數(shù)據(jù)來(lái)自流域站點(diǎn)的觀測(cè)數(shù)據(jù)。其中，氣象數(shù)據(jù)來(lái)自研究流域9個(gè)氣象站點(diǎn)1961—1990年（基準(zhǔn)期）的逐日氣溫、降水量數(shù)據(jù)，降水量數(shù)據(jù)還來(lái)自流域內(nèi)195個(gè)雨量站的逐日觀測(cè)數(shù)據(jù)。黃家港水文站是漢江上游的水文控制站，因此該研究歷史期水文數(shù)據(jù)采用黃家港水文控制站1981—1990年實(shí)測(cè)的逐月天然徑流量序列。

流域DEM資料來(lái)自美國(guó)的USGS網(wǎng)站，分辨率為500 m×500 m。土壤數(shù)據(jù)來(lái)源于中國(guó)科學(xué)院南京土壤研究所土壤數(shù)據(jù)庫(kù)，該土壤數(shù)據(jù)在輸入模型前，需參照模型標(biāo)準(zhǔn)進(jìn)行土壤類型的重新劃分，并賦予模型適用的代碼。土地利用/覆蓋數(shù)據(jù)來(lái)自國(guó)家測(cè)繪總局提供的30 m國(guó)家土地覆蓋分幅TIF圖，采用1990年漢江上游流域的土地利用數(shù)據(jù)，并根據(jù)SWAT模型中土地利用/覆蓋分類標(biāo)準(zhǔn)，將其分為12類。

漢江上游流域基準(zhǔn)期（1961—1990年）和未來(lái)時(shí)期（2011—2100年）RCP 2.6、RCP 4.5、RCP 8.5（RCP：典型濃度路徑）情景下的逐日平均溫度和降水量數(shù)據(jù)序列來(lái)源于北京師范大學(xué)全球變化與地球系統(tǒng)科學(xué)研究院的CMIP5多模式統(tǒng)計(jì)降尺度集合數(shù)據(jù)集。采用Delta方法[12]對(duì)獲得的數(shù)據(jù)序列進(jìn)行校正，以便消除氣候模式模擬的系統(tǒng)偏差。

1.3 研究方法

在我國(guó)的流域水文水資源模擬中，SWAT模型被證實(shí)能較好地模擬流域降水、徑流、蒸散發(fā)等水文過(guò)程[12-13]，因此筆者選用SWAT模型研究漢江上游流域水資源對(duì)未來(lái)氣候變化的響應(yīng)。模擬過(guò)程中，選擇SCS徑流曲線數(shù)法演算流域徑流過(guò)程，選擇Penman-Monteith 公式計(jì)算流域潛在蒸散發(fā)，采用變動(dòng)存儲(chǔ)系數(shù)法演算河道的匯流過(guò)程。采用自動(dòng)優(yōu)選與手工優(yōu)選相結(jié)合的方法對(duì)參數(shù)進(jìn)行率定。參數(shù)的率定將基于先調(diào)整上游再調(diào)整下游、先調(diào)整總水量和蒸發(fā)、地表徑流、地下徑流等各部分水量平衡再調(diào)整過(guò)程的原則[13]。

筆者采用試錯(cuò)法，基于模型的反復(fù)調(diào)試，生成的集水面積約93 215.6 km2，模型最終將漢江上游流域劃分為42個(gè)子流域。根據(jù)優(yōu)勢(shì)土地利用類型/土壤類型方法生成水文響應(yīng)單元（HRU），將占子流域面積5%以上的土地利用及占該類土地利用20%以上的土壤類型生成HRU，即設(shè)定5/20的閾值，按該閾值，將漢江上游流域劃分為241個(gè)HRU。

2 基于SWAT模型的流域徑流模擬

2.1 模型參數(shù)的校準(zhǔn)與驗(yàn)證

選用Nash-Sutcliffe模擬效率系數(shù)（Ens）和決定系數(shù)（R2）2個(gè)指標(biāo)評(píng)估模型的模擬效果，進(jìn)而評(píng)價(jià)模型在漢江上游流域的適用性。Ens的計(jì)算公式為：

Ens=1-ni=1（QEi-QOi）2

ni=1（QOi-QOi）2 （1）

式中，QOi為實(shí)測(cè)徑流，QEi為模擬徑流，QOi為實(shí)測(cè)徑流的平均值，n為模擬的時(shí)段數(shù)。通常取Ens>0.5作為模型在流域適用的評(píng)價(jià)標(biāo)準(zhǔn)[14]。

利用SWAT2005自動(dòng)敏感性分析工具確定敏感性等級(jí)最高的6個(gè)參數(shù)：SCS徑流曲線系數(shù)（CN2）、土壤有效含水量（SOL_AWC）、土壤蒸發(fā)補(bǔ)償系數(shù)（ESCO）、地下水再蒸發(fā)系數(shù)（GW_REVAP）、淺層地下水再蒸發(fā)或下滲的閾值深度（REVAPMN）、淺層地下水回歸流閾值深度（GWQMN）。因此，在參數(shù)率定中主要利用黃家港控制站實(shí)測(cè)徑流對(duì)以上6個(gè)參數(shù)進(jìn)行調(diào)整。其中，選取1981—1987年的實(shí)測(cè)逐月徑流數(shù)據(jù)對(duì)參數(shù)進(jìn)行率定，校準(zhǔn)模型；采用率定后的參數(shù)，選取1988—1990年的實(shí)測(cè)逐月徑流數(shù)據(jù)對(duì)模型進(jìn)行驗(yàn)證。

2.2 校準(zhǔn)與驗(yàn)證結(jié)果

校準(zhǔn)期（1981—1987年）和檢驗(yàn)期（1988—1990年）模型評(píng)價(jià)結(jié)果分別見(jiàn)圖1、2。結(jié)果顯示，校準(zhǔn)期和檢驗(yàn)期漢江上游流域徑流實(shí)測(cè)值與模擬值的變化趨勢(shì)一致，徑流過(guò)程吻合較好。在校準(zhǔn)期決定系數(shù)是0.92，模型的模擬效果較好，表征模型效率的Ens達(dá)0.90。對(duì)模型的相關(guān)參數(shù)進(jìn)行校準(zhǔn)后，對(duì)模型進(jìn)行檢驗(yàn)，結(jié)果表明，檢驗(yàn)期決定系數(shù)為0.85，表征模型效率的Ens達(dá)0.87，模型對(duì)徑流過(guò)程的模擬較好。模型的校準(zhǔn)與驗(yàn)證結(jié)果表明，在漢江上游流域的地理環(huán)境下，筆者構(gòu)建的分布式水文模型SWAT能較準(zhǔn)確地模擬流域的徑流過(guò)程，運(yùn)用該模型開(kāi)展?jié)h江上游流域水資源對(duì)未來(lái)氣候變化的響應(yīng)研究是可行的。

3 漢江上游流域未來(lái)水資源變化

3.1 未來(lái)降水變化

圖3顯示，未來(lái)時(shí)期（2011—2100年），漢江上游流域除春季月份（3—5月）外，其他月份在3種情景下的多年平均月降水量相對(duì)于基準(zhǔn)期均呈增加趨勢(shì)，其中，在冬季（12月至次年2月）枯水期增加幅度尤為明顯，RCP 2.6、RCP 4.5、RCP 8.5情景下平均增加幅度分別為16.36%、22.78%、28.18%。

由表1可知，其中，未來(lái)時(shí)期（2011—2100年）RCP 2.6、RCP 4.5、RCP 8.5情景下的多年平均降水量分別較基準(zhǔn)期增加4.88%、7.74%、12.73%。短期（2020s）除RCP 2.6情景下降水量較基準(zhǔn)期略有減少外，其他2個(gè)情景下的降水量均呈增加趨勢(shì)，其中RCP 8.5的增幅偏大；中期（2050s）和長(zhǎng)期（2080s）3個(gè)情景下的降水量相對(duì)基準(zhǔn)期均為增加，增幅從大到小依次為RCP 8.5、RCP 4.5、RCP 2.6。

3.2 未來(lái)徑流變化

圖4顯示，未來(lái)時(shí)期（2011—2100年）

多年平均月徑流量相對(duì)于基準(zhǔn)期總體呈增加趨勢(shì)，但不同月份之間的增加幅度不同。在RCP 2.6情景下，除了3、5、7月的月徑流量呈減小趨勢(shì)外，其他月份的月徑流量均呈增加趨勢(shì)，特別是在冬季（12月至次年2月）枯水期增幅較大。在RCP 4.5情景下，除了3—4、10月的月徑流量呈減小趨勢(shì)外，其他月份的月徑流量呈均呈增加趨勢(shì)，特別是在冬季（12月至次年2月）枯水期和8月的增幅較大。在RCP 8.5情景下，全年月徑流量均呈增加趨勢(shì)，特別是在11月至次年2月和7、9月增幅較大。

由表1可知，未來(lái)時(shí)期（2011—2100年）RCP 2.6、RCP 4.5、RCP 8.5情景下的多年平均徑流量分別較基準(zhǔn)期增加2.35%、4.74%、9.04%。短期（2020s）除RCP 8.5情景下徑流量較基準(zhǔn)期略有增加外，其他2個(gè)情景下的徑流量均呈減小趨勢(shì)，但減幅均不大。中期（2050s）和長(zhǎng)期（2080s）3個(gè)情景下的徑流量相對(duì)基準(zhǔn)期均為增加，增幅從大到小依次為RCP 8.5、RCP 4.5、RCP 2.6，且徑流量增幅小于降水量的增幅，這可能與氣候變暖導(dǎo)致蒸發(fā)量增大有關(guān)。

4 結(jié)論

該研究基于CMIP5多模式氣象輸出數(shù)據(jù)集，結(jié)合分布式水文模型SWAT，分析了3種氣候變化情景下未來(lái)時(shí)期（2011—2100年）漢江上游流域水資源相對(duì)基準(zhǔn)期（1961—1990年）的變化，得出以下結(jié)論：

（1） 水文模型的模擬結(jié)果表明，SWAT分布式水文模型對(duì)漢江上游流域的水資源模擬有較好地適用性。運(yùn)用該模型開(kāi)展未來(lái)水資源對(duì)氣候變化的響應(yīng)研究是可行的。

（2） 未來(lái)時(shí)期（2011—2100年）RCP 2.6、RCP 4.5、RCP 8.5情景下的多年平均降水量均呈增加趨勢(shì)，增幅依次為RCP 8.5、RCP 4.5、RCP 2.6，并且中期（2050s）和長(zhǎng)期（2080s）的增加幅度較為明顯。未來(lái)時(shí)期，多年平均月降水量在冬季（12月至次2月）枯水期增加幅度最為顯著。

（3）未來(lái)時(shí)期（2011—2100年）RCP 2.6、RCP 4.5、RCP 8.5情景下的多年平均徑流量均呈增加趨勢(shì)，增幅依次為RCP 8.5、RCP 4.5、RCP 2.6，但其增幅小于降水量的增幅。長(zhǎng)期（2080s）徑流量

相對(duì)基準(zhǔn)期的增加幅度較為顯著。并且，未來(lái)時(shí)期（2011—2100年），3種情景下的多年平均月徑流量在冬季（12月至次年2月）枯水期增幅較大，在洪水期增幅相對(duì)較小。

參考文獻(xiàn)

[1]CHRISTENSEN N S，WOOD A W，VOISIN N，et al.The effects of climate change on the hydrology and water resources of the Colorado River basin [J].Climatic change，2004，62（1）：337-363.

[2] 黃國(guó)如，武傳號(hào)，劉志雨，等.氣候變化情景下北江飛來(lái)峽水庫(kù)極端入庫(kù)洪水預(yù)估[J].水科學(xué)進(jìn)展，2015，26（1）：10-19.

[3] 張建云，王國(guó)慶.氣候變化對(duì)水文水資源影響研究[M].北京：科學(xué)出版社，2007.

[4] 賈仰文，高輝，牛存穩(wěn)，等.氣候變化對(duì)黃河源區(qū)徑流過(guò)程的影響[J].水利學(xué)報(bào)，2008，39（1）：52-58.

[5] 夏智宏，周月華，許紅梅.基于SWAT模型的漢江流域水資源對(duì)氣候變化的響應(yīng)[J].長(zhǎng)江流域資源與環(huán)境，2010，19（2）：158-163.

[6] HORTON P，SCHAEFLI B，MEZGHANI A，et al.Assessment of climate change impacts on Alpine discharge regimes with climate model uncertainty [J].Hydrol Process，2006，20（10）：2091-2109.

[7] CHIEW F H S，VAZE J，VINEY N R，et al.Estimation of runoff and the impact of climate change and development on runoff across the Murray-Darling Basin[C]//Proceedings of the water down under 2008.Adelaide，Australia：[s.n]，2008：1957-1968.

[8] ZHANG Z X，CHEN X，XU C Y，et al.Evaluating the non-stationary relationship between precipitation and streamflow in nine major basins of China during the past 50 years[J].Journal of hydrology，2011，409（1/2）：81-93.

[9] 胡珊珊，鄭紅星，劉昌明，等.氣候變化和人類活動(dòng)對(duì)白洋淀上游水源區(qū)徑流的影響[J].地理學(xué)報(bào)，2012，67（1）：62-70.

[10] 成愛(ài)芳，馮起，張健愷，等.未來(lái)氣候情景下氣候變化響應(yīng)過(guò)程研究綜述[J].地理科學(xué)，2015，35（1）：84-90.

[11] 張利平，秦琳琳，胡志芳，等.南水北調(diào)中線工程水源區(qū)水文循環(huán)過(guò)程對(duì)氣候變化的響應(yīng)[J].水利學(xué)報(bào)，2010，41（11）：1261-1271.

[12] 李峰平，章光新，董李勤.氣候變化對(duì)水循環(huán)與水資源的影響研究綜述[J].地理科學(xué)，2013，33（4）：457-464.

[13] 孫瑞，張雪琴.基于SWAT模型的流域徑流模擬研究進(jìn)展[J].水文，2010，30（3）：28-32，47.

[14] 張利平，胡志芳，秦琳琳，等.2050年前南水北調(diào)中線工程水源區(qū)地表徑流的變化趨勢(shì)[J].氣候變化研究進(jìn)展，2010，6（6）：391-397.