通信作者:薛聯(lián)青 教授。E-mail:lqxue@hhu.edu.cn

氣候變化條件下洪澤湖以上流域水資源演變趨勢(shì)

王苓如1,薛聯(lián)青1,2,王思琪1,張榮3

(1.河海大學(xué)水文水資源與水利工程科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,江蘇 南京210098;

2.石河子大學(xué)水利建筑工程學(xué)院,新疆 石河子832003; 3.泰安市水文局,山東 泰安271000)

摘要：基于IPCC對(duì)全球和中國(guó)的氣候變化趨勢(shì),利用1990—2011年氣象資料,采用增量情景設(shè)置方法,分析氣候變化情景下洪澤湖以上流域水資源的演變趨勢(shì)。結(jié)果表明:該流域水資源量對(duì)降雨變化有較強(qiáng)的敏感性,實(shí)際蒸散發(fā)對(duì)溫度變化的敏感性較強(qiáng)。與基準(zhǔn)期相比,在氣溫同等條件、降水增加情景下,流域水資源量呈增加趨勢(shì);在降水同等條件、氣溫升高情景下,流域的實(shí)際蒸發(fā)會(huì)增加,導(dǎo)致水資源量呈減少的趨勢(shì)。徑流年內(nèi)分配受降水變化影響較大;隨著降水增加,徑流年內(nèi)分配更集中,加大年內(nèi)徑流分配差異,可能加大流域湖泊調(diào)蓄壓力。

關(guān)鍵詞：氣候變化;水資源演變;SWAT模型;洪澤湖以上流域

基金項(xiàng)目:國(guó)家自然科學(xué)基金(41371052, U1203282);江蘇省“青藍(lán)工程”及兵團(tuán)創(chuàng)新團(tuán)隊(duì)項(xiàng)目

作者簡(jiǎn)介:王苓如(1989—),女,碩士研究生,研究方向?yàn)樗h(huán)境保護(hù)及水資源演變。E-mail:wanglingru891228@126.com

中圖分類(lèi)號(hào)：TV211文獻(xiàn)標(biāo)志碼：A

收稿日期：(2014-09-28編輯：徐娟)

Evolutionary tendency of water resources in the basin above Hongze Lake

under the condition of climate change

WANG Lingru1, XUE Lianqing1,2,WANG Siqi1,ZHANG Rong3

(1.CollegeofHydrologyandWaterReseources,HohaiUniversity,Nanjing210098,China;

2.ShiheziUniversit,Shihezi832003,China;

3.Tai'anCityBureauofHydrology,Tai’an271000,China)

Abstract：Climate change has caused re-allocation of water resources in the spatial and temporal pattern, leading to more frequent occurrence of extreme hydrological events. Based on the predictions to global and China climate change by the Intergovernmental Panel on Climate Change (IPCC), the authors of this paper analyzed the evolutionary tendency of water resources in the basin above Hongze Lake under the condition of climate change by the method of setting incremental scene and using the climatological data of the years of 1990 to 2011. The results show that the catchment runoff was sensitive to precipitation changes; the evaporation was relatively sensitive to temperature changes. Compared with the baseline period, in the scenario of increasing precipitation and the same temperature, the water resources in the basin will be in the tendency of increasing; in the scenario of the same precipitation and temperature rise, the actual evaporation in the basin will increase, causing the tendency of water resources decreasing; the annual distribution of runoff is greatly influenced by the change of precipitation; runoff annual distribution becomes more concentrated with the increase of precipitation, which may result in the increasing the difference of runoff annual distribution, and increasing the regulating pressure of Hongze Lake.

Key words： climate change; evolutionary of water resources; SWAT model; basin above Hongze Lake

氣候變化引起了水文循環(huán)的變化,導(dǎo)致水資源時(shí)空分布上的重新分配及水資源總量的改變,引發(fā)了部分流域極端氣候及水文事件頻率和強(qiáng)度的增加,加劇了洪旱災(zāi)害頻發(fā)的風(fēng)險(xiǎn)[1-2]。全球氣候變化既影響水文與生態(tài)系統(tǒng)、經(jīng)濟(jì)與生活,甚至影響國(guó)家、地區(qū)的可持續(xù)發(fā)展。因此,全球氣候變化對(duì)水資源的影響已引起國(guó)際普遍關(guān)注[3-4]。水文學(xué)者一般利用氣候模式輸出的氣溫、降水等結(jié)果或是氣候情景作為輸入資料,再通過(guò)水文模型來(lái)分析流域水資源變化趨勢(shì)。有許多學(xué)者運(yùn)用SWAT模型來(lái)研究氣候變化的水文效應(yīng):Cruise等[5]研究不同氣候變化條件下,美國(guó)東南部現(xiàn)在及未來(lái)水資源的變化情況;Stone等[6]結(jié)合密蘇里河流域的歷史數(shù)據(jù),模擬CO2濃度增加兩倍情況下的水文變化;李道峰等[7]探討研究了SWAT模型在黃河河源區(qū)的適用性,假定了5種土地利用變化情景和24種氣候變化組合情景,研究區(qū)域水資源的變化趨勢(shì)。

流域內(nèi)未來(lái)水資源變化將直接影響洪澤湖的調(diào)蓄功能,是水資源配置與管理運(yùn)行的重要依據(jù)。筆者以洪澤湖以上流域?yàn)檠芯繀^(qū),假設(shè)其下墊面狀況在水文響應(yīng)預(yù)測(cè)期內(nèi)保持不變;根據(jù)IPCC(the Intergovernmental Panel on Climate Change)第4次報(bào)告(the fourth assessment report,AR4),全球大氣二氧化碳的含量增加1倍,氣溫增幅為1～4℃,降水變幅為20%,設(shè)置不同氣候變化情景[8-9];運(yùn)用SWAT(soil and water assessment tool)分布式水文模型模擬分析該流域水資源演變情況。

1氣候、水文變化概況

淮河流域位于東經(jīng)112°～121°,北緯31°～36°,處于北亞熱帶向暖溫帶過(guò)渡區(qū),屬大陸性季風(fēng)氣候。淮河干流至洪澤湖三河閘(中渡)出口為洪澤湖以上流域,匯水面積為15.8萬(wàn)km2。淮河干流為洪澤湖入湖水量的主要來(lái)源,約占87%。流域內(nèi)降水時(shí)空分布不均,極易產(chǎn)生水旱災(zāi)害。洪澤湖以上流域年平均氣溫在11~16℃,大部分地區(qū)氣溫呈上升趨勢(shì);年降水量總體呈上升趨勢(shì),增速緩慢,四季降水有增有減,春季和秋季降水量減少,而夏季和冬季降水量增加。1969—2011年流域徑流經(jīng)M-K檢驗(yàn)表明:該流域徑流總體呈顯著減少趨勢(shì);年內(nèi)徑流在春夏季節(jié)減少明顯,在秋冬季節(jié)略呈上升趨勢(shì)。

淮河流域多年平均地表水資源量595億m3。近幾十年來(lái),人均地表水資源量發(fā)生銳減,因缺水造成的農(nóng)業(yè)旱災(zāi)日益嚴(yán)重。加之水資源時(shí)空分配極為不均,導(dǎo)致水資源短缺問(wèn)題更加突出[10-11]。洪澤湖當(dāng)?shù)睾^(qū)徑流逐漸減少,區(qū)域降水和淮河干流水資源時(shí)空分配不均,導(dǎo)致湖區(qū)地表徑流年際、年內(nèi)分配不均,給洪澤湖調(diào)蓄帶來(lái)較大壓力[12]。

2資料與方法

2.1資料說(shuō)明

以洪澤湖以上流域?yàn)檠芯繀^(qū),基于數(shù)字高程(圖1)、土地利用(圖2)和土壤類(lèi)型(圖3)建立SWAT模型,輸入19個(gè)氣象站點(diǎn)(圖4)1990—2011年氣象數(shù)據(jù)資料,模擬流域基準(zhǔn)期的逐月徑流,采用淮河干流控制站蚌埠水文站1990—2011年天然徑流量進(jìn)行率定和驗(yàn)證。在率定最佳月徑流模擬結(jié)果參數(shù)的基礎(chǔ)上,分析25種氣候變化組合情景下流域水資源演變趨勢(shì)。

圖1　洪澤湖以上流域DEM

圖2　研究區(qū)土地利用分布

圖3　研究區(qū)土壤分布

圖4　研究區(qū)水文氣象站點(diǎn)分布

其中研究區(qū)土地利用圖參照SWAT模型中美國(guó)土地利用參數(shù)庫(kù),將其重分類(lèi)分為10類(lèi)(詳見(jiàn)GB_T2010—2007《土地利用現(xiàn)狀分類(lèi)》);將土壤分布圖按照模型要求及數(shù)據(jù)說(shuō)明整理(詳見(jiàn)世界土壤數(shù)據(jù)庫(kù)中國(guó)土壤數(shù)據(jù)集),建立SWAT模型土壤數(shù)據(jù)庫(kù)。

由于模型中天氣發(fā)生器為月尺度,為了便于進(jìn)行氣候變化情景模擬,筆者將日降水?dāng)?shù)據(jù)、最高和最低氣溫?cái)?shù)據(jù)、濕度數(shù)據(jù)以及輻射數(shù)據(jù)通過(guò)EGNmaker 4.1軟件進(jìn)行時(shí)間上的降尺度處理,得到月降水?dāng)?shù)據(jù)及相關(guān)氣象參數(shù)。

2.2水文要素敏感性定義

水文要素敏感性分析是在定量研究氣候變化情景下(如降水變化0,±10%,±20%;氣溫變化0℃,+1℃,+2℃,+3℃,+4℃),通過(guò)水文要素的變化率來(lái)表征不同氣候情景相對(duì)于基準(zhǔn)期氣候情景(即ΔT=0,ΔP=0) 的水資源變化趨勢(shì)。其中水文要素的變化率采用式(1)計(jì)算:

(1)

氣候變化幅度相同的情況下,水文要素的變化率越大表明其敏感性越強(qiáng)。敏感性可以揭示流域不同水文要素對(duì)氣候相應(yīng)的機(jī)理和差異。

3SWAT模型的適用性

3.1模型概述

SWAT模型在流域水系基礎(chǔ)上根據(jù)閾值劃分出子流域,在各個(gè)子流域內(nèi)運(yùn)用傳統(tǒng)概念模型推求凈雨、計(jì)算匯流,得出各子流域出口斷面徑流。該模型主要用于水文方面的模擬,具有很強(qiáng)的物理基礎(chǔ),適用于多土地利用和土壤類(lèi)型的大型流域。至今,SWAT模型被廣泛應(yīng)用于模擬流域水文過(guò)程對(duì)氣候變化、下墊面條件變化的響應(yīng),評(píng)價(jià)人類(lèi)活動(dòng)對(duì)流域生態(tài)環(huán)境的影響,為流域水資源的規(guī)劃管理提供依據(jù)[13]。

3.2模型模擬方法的選擇

在模型模擬時(shí)選擇“日降雨數(shù)據(jù)SCS徑流曲線(xiàn)數(shù)方法以日為時(shí)間單位”方法進(jìn)行徑流演算。選擇“Priestly-Taylor”方法計(jì)算流域的潛在蒸發(fā)量,相關(guān)研究表明該方法適用于濕潤(rùn)地區(qū)的潛在蒸發(fā)量計(jì)算[14]。降雨量分布模擬方法選擇偏正態(tài)分布,河道演算采用馬斯京根法。模型輸出選擇月時(shí)間尺度,選取相關(guān)指標(biāo)對(duì)參數(shù)率定進(jìn)行評(píng)價(jià)。

3.3模型參數(shù)的率定

本次采用SWAT-CUP 2012軟件中SUFI-2優(yōu)化算法進(jìn)行參數(shù)率定。SWAT模型中的參數(shù)眾多,由于本文側(cè)重流域的徑流模擬,故參數(shù)選擇主要是徑流參數(shù)。在參數(shù)率定時(shí)遵循“從上游至下游;先調(diào)整水量平衡,再調(diào)整徑流過(guò)程;先調(diào)整徑流,后調(diào)整蒸發(fā)”的原則。由于沒(méi)有流域?qū)崪y(cè)蒸散發(fā)資料,故本次研究中僅進(jìn)行徑流校準(zhǔn)?？紤]到模型中參數(shù)數(shù)量較多且模型結(jié)構(gòu)復(fù)雜,故本次參數(shù)率定將手動(dòng)優(yōu)化與自動(dòng)優(yōu)化相結(jié)合,根據(jù)流域參數(shù)敏感性,來(lái)率定和驗(yàn)證模型。

3.4模擬結(jié)果評(píng)價(jià)

本文選用決定系數(shù)R2、Nash效率系數(shù)NNS和相對(duì)誤差r這3個(gè)指標(biāo)來(lái)評(píng)價(jià)該模型在洪澤湖以上流域的適用性。

(3)

(4)

(5)

運(yùn)用淮河控制站—蚌埠(吳家渡)水文站1990—2003年月徑流資料進(jìn)行模型參數(shù)率定,2004—2011年月徑流資料用來(lái)模型檢驗(yàn)(表1，圖5)。結(jié)果表明:分布式水文模型SWAT在洪澤湖以上游流域適用于徑流模擬,為探討氣候變化對(duì)流域水資源影響奠定了基礎(chǔ)。

表1　1990—2011年模擬驗(yàn)證期月徑流評(píng)價(jià)

圖5　1990—2011年模擬驗(yàn)證期月徑流模擬值與實(shí)測(cè)值比較

4研究區(qū)水資源對(duì)未來(lái)氣候變化的響應(yīng)分析

根據(jù)IPCC提出的氣候變化預(yù)測(cè)結(jié)果和研究區(qū)可能出現(xiàn)的氣候變化條件,假定25 種氣候變化組合情景(表2),即在研究區(qū)22年平均日氣溫和降水的基礎(chǔ)上,將溫度增加0℃、1℃、2℃、3℃、4℃和降水變化0、±10%、±20%進(jìn)行兩兩組合。假定在水文響應(yīng)期內(nèi)下墊面條件不變,分別對(duì)25種氣候變化情景進(jìn)行月徑流模擬,得到每種情景下的月平均徑流量。

表2　氣候變化情景類(lèi)型

經(jīng)計(jì)算得到每種氣候情景下研究區(qū)年水資源量和蒸發(fā)量相對(duì)于基準(zhǔn)期模擬值的變化率(表3、表4),其中蚌埠站變化率表示淮河干流水資源變化,三河閘(中渡)變化率表示研究區(qū)水資源整體變化。

表3　不同氣候變化情景下水資源變化率

a. 不同氣候變化情景下研究區(qū)水資源量變化差異明顯。水資源量減少最多的氣候情景是降水減少20%、氣溫增加4℃的情景方案(S21);增加最多的氣候情景是降水增加20%、氣溫不變的情景方案(S5)。

b. 不同氣候變化情景下研究區(qū)實(shí)際蒸散發(fā)變化沒(méi)有水資源量變化明顯。實(shí)際蒸散發(fā)對(duì)氣溫變化的敏感性強(qiáng),其減少最多的氣候情景是氣溫不變、降水減少20%的情景方案(S1);增加最多的氣候情景是氣溫增加4℃、降水增加20%的情景方案(S25)。

c. 降水增加或氣溫降低都會(huì)導(dǎo)致研究區(qū)水資源量的增加。流域水資源量變化與降水變化呈正相關(guān),與氣溫變化呈負(fù)相關(guān);其對(duì)降雨變化的敏感性強(qiáng),對(duì)氣溫變化的敏感性相對(duì)較弱。說(shuō)明洪澤湖以上流域未來(lái)徑流變化主要受到降水的影響,受氣溫影響相對(duì)較小。

d. 淮河干流水資源對(duì)氣候變化響應(yīng)的敏感性強(qiáng)于流域水資源。即在同一氣候變化情景下,洪澤湖來(lái)水量的增加或減少程度比下泄水量的增加或減少程度要大。在氣候變化較大的情景(如S4或S22)下,可能會(huì)導(dǎo)致洪澤湖發(fā)生洪澇或干旱事件,需要提前做好防護(hù)措施。

表4　不同氣候變化情景下研究區(qū)實(shí)際蒸散發(fā)變化率

降水、氣溫的變化除了導(dǎo)致流域年際徑流的變化,也將將影響其年內(nèi)徑流分布。由圖6可知:

圖6　氣溫不變、降水變化情景下研究區(qū)徑流年內(nèi)分布

a. 在相同氣溫、降水變化情景下,隨著降水增加,年內(nèi)各月徑流與降水變化趨勢(shì)一致,但月增幅各不相同,在秋季增加最多。流域徑流在夏季隨著降水增加,變幅十分明顯:在1-6月除降水增加20%情況下變化明顯外,其余降水變化情景下變化趨勢(shì)一致。

b. 在相同氣溫、降水變化情景下,淮河干流徑流受降水變化影響比流域出口徑流變化大;淮河干流徑流在夏秋兩季增加明顯;研究區(qū)出口徑流只在降水增加20%的氣候變化情景下增幅明顯。在氣溫相同,降水變化情景下,降水增加會(huì)導(dǎo)致流域徑流分配更集中,年內(nèi)各月徑流分配差異加大;洪澤湖在夏秋兩季來(lái)水量可能增加明顯,下泄量變化較小,加大了洪澤湖調(diào)蓄壓力。

圖7　降水不變、氣溫變化情景下研究區(qū)徑流年內(nèi)分布

由圖7可知：

a. 在降水相同、氣溫變化情景下,隨著溫度升高,各月徑流均呈減小趨勢(shì);其中在氣溫增加2℃的情景下,流域徑流在秋季減少最多。原因可能是當(dāng)溫度升高超過(guò)一定水平時(shí),蒸發(fā)量增大到一定程度后,空氣中水氣含量增加,云層增厚,導(dǎo)致蒸發(fā)能力降低。研究區(qū)在氣溫增加2℃時(shí),實(shí)際蒸發(fā)較大,使得徑流值減小。

b. 氣溫對(duì)徑流的影響主要是通過(guò)蒸發(fā)體現(xiàn),其對(duì)年內(nèi)徑流分配的影響程度小于降水變化。其中冬季徑流減少程度最小,造成該現(xiàn)象的原因很復(fù)雜,可能是由于在降水相同條件下,氣溫上升導(dǎo)致蒸發(fā)隨著溫度升高而增加,造成徑流減少。

4總結(jié)與討論

a. 運(yùn)用SWAT模型模擬洪澤湖以上流域月徑流過(guò)程,模擬精度達(dá)到了評(píng)價(jià)標(biāo)準(zhǔn),表明該模型適用于洪澤湖以上流域徑流模擬,為分析未來(lái)氣候變化下淮河流域流域水資源演變規(guī)律奠定基礎(chǔ)。

b. 不同氣候組合情景下,流域徑流變化比實(shí)際蒸散發(fā)變化明顯。降水對(duì)徑流的影響要明顯大于溫度對(duì)徑流的影響。溫度對(duì)流域?qū)嶋H蒸散發(fā)的影響大于降水變化對(duì)流域蒸散發(fā)的影響。隨著降水條件的變化,流域?qū)嶋H蒸散發(fā)變化差異較大,其中降水變化±20%時(shí)對(duì)蒸散發(fā)影響最大。

c. 由于沒(méi)有流域的月蒸散發(fā)實(shí)測(cè)資料,故在模型模擬研究中沒(méi)有用實(shí)測(cè)月蒸散發(fā)數(shù)據(jù)對(duì)模擬結(jié)果進(jìn)行校準(zhǔn)和驗(yàn)證。希望在下一步的研究中可以對(duì)此作出改善。

d. 降水相同、氣溫變化情景下，流域徑流年內(nèi)分布在氣溫增加2℃時(shí)減小明顯,具體原因還有待進(jìn)一步驗(yàn)證。

e. 運(yùn)用SWAT模型模擬未來(lái)氣候變化下流域月徑流時(shí),沒(méi)有考慮未來(lái)下墊面變化可能對(duì)流域水資源的影響,使氣候情景模擬的預(yù)測(cè)結(jié)果具有一定的不確定性。

參考文獻(xiàn)：

[1] 鄧慧平,吳正方,唐來(lái)華.氣候變化對(duì)水文和水資源影響研究綜述[J].地理學(xué)報(bào),1996,51(增刊):161-170.(DENG Huiping,WU Zhengfang,TANG Laihua.A review of the studies of the impacts of climate change on hydrology and water resources[J].Journal of Acta Geographica Sinica,1996,51(sup):161-170.(in Chinese))

[2] 夏軍,劉春蓁,任國(guó)玉.氣候變化對(duì)我國(guó)水資源影響研究面臨的機(jī)遇與挑戰(zhàn)[J].地球科學(xué)進(jìn)展,2011,26(1):1-12.(XIA Jun,LIU Chunzhen,REN Guoyu.Opprotunity and challenge of the climate change impact on the water resource of China[J]. Journal of Advances in Earth Science,2011,26(1):1-12.(in Chinese))

[3] 袁飛,謝正輝,任立良,等.氣候變化對(duì)海河流域水文特性的影響[J].水利學(xué)報(bào),2005,36(3):274-279.(YUAN Fei,XIE Zhenghui,REN Liliang,et al.Hydrological variation in Haihe River Basin due to climate change[J]. Journal of Hydraulic Engineering,2005,36(3):274-279.(in Chinese))

[4] 張利平,秦琳琳,胡志芳,等.南水北調(diào)中線(xiàn)工程水源區(qū)水文循環(huán)過(guò)程對(duì)氣候變化的響應(yīng)[J].水利學(xué)報(bào),2010,41(11):1261-1271.(ZHANG Liping,QIN Linlin,HU Zhifang,et al.Simulated hydrologic respinses to climate change of water source area in the Middle Route of South-to-North Water Transfer Project[J]. Journal of Hydraulic Engineering,2010,41(11):1261-1271.(in Chinese))

[5] CRUISE J F,LIMAYE A S. Assessment of impacts of climate change on water quality in the southeastern United States [J]. Journal of the American Water Resources Association,1999,35(6):1539-1550.

[6] STONE FELT M C,FONTAME T A,HOTCHKISS R H. Impacts of climate change on Missouri river basin water yield [J]. Journal of the American Water Resources Association,2001,37:1119-1129.

[7] 李道峰,田英,劉昌明.黃河河源區(qū)變化環(huán)境下分布式水文模擬[J].地理學(xué)報(bào),2004,59(4):565-573.(LI Daofeng,TIAN Ying,LIU Changming.Distributed hydrological simulation of the source region of the Yellow River under environment changes[J]. Journal of Acta Geographica Sinica,2004,59(4):565-573.)

[8] IPCC. IPCC forth report of assessment of climate change [R].Cambridge,United Kingdom and New York:Cambridge University Press,2007.

[9] 夏智宏,周月華,許紅梅.基于SWAT模型的漢江流域水資源對(duì)氣候變化的響應(yīng)[J]. 長(zhǎng)江流域資源與環(huán)境,2010,19(2):158-163.(XIA Zhihong,ZHOU Yuehua,XU Hongmei.Water resources responses to climate changes in Hanjiang River Basin based on SWAT model[J]. Journal of Resources and Environment in the Yangtze Basin,2010,19(2):158-163.(in Chinese))

[10] 費(fèi)永法,張炎齋.淮河流域水資源問(wèn)題及對(duì)策[J].水利規(guī)劃與設(shè)計(jì),2005(1):23-25.(FEI Yongfa,ZHANG Yanzhai.Questions and suggestions on water resouces in Huaihe River Basin[J]. Journal of Water Resources Planning and Design,2005(1):23-25.(in Chinese))

[11] 胡瑞,左其亭,竇明,等. 淮河流域水資源承載力研究[C]//四川大學(xué).第六屆中國(guó)水論壇:河流開(kāi)發(fā)、保護(hù)與水資源可持續(xù)利用.北京:中國(guó)水利水電出版社,2008:347-352.

[12] 楚恩國(guó).洪澤湖水資源現(xiàn)狀分析及對(duì)策[J].中國(guó)水利,2007(23):33-35.(CHU Enguo. Anaysis of present situation of water resources in Hongze Lake and its counter measures[J]. Journal of China Water Resources,2007(23):33-35.(in Chinese))

[13] 肖軍倉(cāng),周文斌,羅定貴,等.非點(diǎn)源污染模型-SWAT用戶(hù)應(yīng)用指南[M]. 北京:地質(zhì)出版社,2009.

[14] 彭世彰,徐俊增.參考作物蒸發(fā)蒸騰量計(jì)算方法的應(yīng)用比較[J].灌溉排水學(xué)報(bào),2004,23(6):5-9.(PENG Shizhang,XU Junzeng. Comparsion of reference crop evapotranspiration computing methods[J]. Journal of Irrigation and Drainage,2004,23(6):5-9.(in Chinese))