胡 倩，王軍霞，劉世強(qiáng)，吳嘉鈴，唐仲華*，成建梅

(1.中國(guó)地質(zhì)大學(xué)(武漢)環(huán)境學(xué)院，湖北 武漢 430078;2.中國(guó)地質(zhì)大學(xué)(武漢)數(shù)學(xué)與物理學(xué)院，湖北 武漢 430074;3.山東正元建設(shè)工程有限責(zé)任公司，山東 臨沂 276006)

洞庭湖是我國(guó)第二大淡水湖[1]。洞庭湖作為長(zhǎng)江流域重要的調(diào)蓄型湖泊，保護(hù)洞庭湖湖泊水資源并掌握洞庭湖水文水資源狀況，對(duì)于洞庭湖生態(tài)保護(hù)和水資源合理調(diào)配具有十分重要的意義。水資源短缺、水污染加劇、地下水超采和用水效率低下正在不斷加劇我國(guó)水資源供需矛盾[2]。隨著全球氣候變暖以及人類(lèi)活動(dòng)的干擾，洞庭湖的面積和來(lái)水量不斷減小，其生態(tài)環(huán)境和生態(tài)資源都受到了嚴(yán)重影響，湖區(qū)內(nèi)水環(huán)境也發(fā)生了變化，這些已引起學(xué)者們的廣泛關(guān)注[3-7]。面對(duì)氣候變化背景下，學(xué)者們對(duì)氣候變化環(huán)境下洞庭湖水資源評(píng)價(jià)進(jìn)行了大量的研究，如徐志[8]分析了洞庭湖水資源現(xiàn)狀，并探討了氣候變化對(duì)洞庭湖水資源的影響；周慧[9]基于氣象、水文觀測(cè)資料和氣候模式CMIP5模擬數(shù)據(jù)，采用統(tǒng)計(jì)學(xué)分析方法和水文學(xué)模擬方法，分析了洞庭湖流域未來(lái)集雨范圍內(nèi)以及過(guò)境水資源的變化情況。在區(qū)域水資源評(píng)價(jià)方法方面，學(xué)者們多采用水均衡法、解析法和數(shù)值模型法來(lái)對(duì)地下水資源進(jìn)行評(píng)價(jià)[10-12]。但現(xiàn)有研究利用水均衡法計(jì)算洞庭湖平原水資源量以及氣候變化對(duì)水資源量影響方面的研究較少。

鑒于此，本文通過(guò)調(diào)查和收集洞庭湖平原地形、氣象水文、土壤、土地利用等資料，利用ArcSWAT建立洞庭湖平原分布式水文模型，在評(píng)價(jià)該模型對(duì)洞庭湖平原建模適用性的基礎(chǔ)上，對(duì)洞庭湖平原水資源量進(jìn)行計(jì)算，并通過(guò)假定氣候情景法來(lái)模擬預(yù)測(cè)未來(lái)氣候變化對(duì)洞庭湖平原水資源量的影響，對(duì)洞庭湖平原生態(tài)經(jīng)濟(jì)發(fā)展具有重要的意義。

1 研究區(qū)概況

洞庭湖平原位于湖南省北部，主要由長(zhǎng)江通過(guò)松滋、太平、藕池、調(diào)弦四口輸入的泥沙和洞庭湖水系湘江、資水、沅江、澧水等帶來(lái)的泥沙沖積而成，北由東面岳陽(yáng)城陵磯注入長(zhǎng)江[13]。研究區(qū)范圍西起石門(mén)—常德市，東抵湘陰，南起常德—益陽(yáng)，北邊是江漢平原，地理坐標(biāo)為東經(jīng)111°14′～113°15′、北緯28°15′～30°45′，面積為2.14×104km2。洞庭湖平原內(nèi)地形以丘陵、平原為主，該區(qū)域?qū)儆趤啛釒Ъ撅L(fēng)性濕潤(rùn)氣候，年平均氣溫在15～17℃之間。研究區(qū)地形特征及水系分布，見(jiàn)圖1。

圖1 洞庭湖平原地形特征及水系分布圖

氣候變化會(huì)導(dǎo)致水文循環(huán)的變化，加劇水資源分配不均，對(duì)流域水資源產(chǎn)生一定的影響[14-17]。為了分析研究區(qū)內(nèi)水資源量與氣候要素特征之間的關(guān)系，本文以洞庭湖平原城陵磯水文站為代表，從月尺度上分析了2010—2013年年平均氣溫、降雨量與流量之間的變化特征(見(jiàn)圖2)，結(jié)果表明：研究區(qū)降雨量和流量6月份達(dá)到最大，平均氣溫在8月份達(dá)到最大，降雨量和流量以及氣溫在月尺度內(nèi)存在明顯的相關(guān)性關(guān)系。

圖2 2010—2013年洞庭湖平原城陵磯水文站年均氣候要素變化特征圖

2 數(shù)據(jù)來(lái)源與研究方法

2.1 數(shù)據(jù)來(lái)源與處理

建立研究區(qū)的SWAT模型需要空間數(shù)據(jù)和屬性數(shù)據(jù)兩類(lèi)數(shù)據(jù)。其中，空間數(shù)據(jù)包括洞庭湖平原的數(shù)字高程模型(DEM)、土地利用類(lèi)型、土壤類(lèi)型數(shù)據(jù)、河流水系圖，以及水文站、氣象站點(diǎn)位信息；屬性數(shù)據(jù)包括土壤物理化學(xué)屬性信息、植被覆蓋等調(diào)查統(tǒng)計(jì)數(shù)據(jù)，以及水文站監(jiān)測(cè)點(diǎn)、氣象站觀測(cè)點(diǎn)數(shù)據(jù)等。本文主要使用數(shù)據(jù)的精度，見(jiàn)表1。

表1 SWAT模型主要數(shù)據(jù)的精度

(1) 數(shù)字高程模型(DEM)數(shù)據(jù)。研究區(qū)DEM數(shù)據(jù)由中國(guó)科學(xué)院地理空間數(shù)據(jù)云下載的分辨率為30 m×30 m的ASTER GDEM數(shù)據(jù)，投影方式為墨卡托投影，在建立數(shù)字高程模型時(shí)要對(duì)加載的柵格數(shù)據(jù)進(jìn)行裁剪等預(yù)處理，DEM數(shù)據(jù)可以提取數(shù)字河網(wǎng)、劃分子流域以及水文響應(yīng)單元(HRU)。

(2) 土地利用數(shù)據(jù)。土地利用數(shù)據(jù)來(lái)源于中國(guó)科學(xué)院地理空間數(shù)據(jù)云的Landsat 4-5 TM衛(wèi)星遙感影像，影像精度為30 m×30 m。通過(guò)對(duì)遙感影像進(jìn)行圖像鑲嵌、裁剪、遙感解譯，得到研究區(qū)2013年土地利用數(shù)據(jù)。將洞庭湖平原劃分為7種土地利用類(lèi)型，分別為林地、旱地、水田、水體、建設(shè)用地、裸地和草地。

(3) 土壤類(lèi)型數(shù)據(jù)。土壤類(lèi)型數(shù)據(jù)來(lái)源于1∶100萬(wàn)的世界和諧土壤數(shù)據(jù)庫(kù)(HWSD)，將土壤類(lèi)型數(shù)據(jù)進(jìn)行裁剪并重分類(lèi)，得到洞庭湖流域內(nèi)分布有13種土壤類(lèi)型。

(4) 土壤屬性和植被覆蓋數(shù)據(jù)。植被覆蓋和土壤屬性參數(shù)在SWAT自帶數(shù)據(jù)庫(kù)中獲得，其中土壤屬性數(shù)據(jù)庫(kù)中的物理屬性大部分由HWSD土壤數(shù)據(jù)庫(kù)獲得，一部分利用SPAW軟件計(jì)算得出。

(5) 氣象水文數(shù)據(jù)。本文收集到的氣象數(shù)據(jù)來(lái)自中國(guó)氣象數(shù)據(jù)共享服務(wù)中心，包括13個(gè)氣象站2010—2013年的日最高氣溫、日最低氣溫、日相對(duì)濕度、日降雨量、日照時(shí)數(shù)和日均風(fēng)速實(shí)測(cè)數(shù)據(jù)。氣象站點(diǎn)包括松滋、公安、石首、澧縣、安鄉(xiāng)、南縣、岳陽(yáng)、益陽(yáng)、常德、漢壽、華容、沅江和湘陰。本文將收集到的石龜山(澧水控制站)、南縣(藕池河北支控制站)、南咀(南洞庭湖入湖控制站)、城陵磯(洞庭湖出口控制站)4個(gè)站點(diǎn)2010—2013年的日均徑流量觀測(cè)數(shù)據(jù)，用于后期模型的率定及驗(yàn)證。

2.2 研究方法

本文采用SWAT模型構(gòu)建洞庭湖流域分布式水文模型，對(duì)模型進(jìn)行率定后，運(yùn)用水均衡法計(jì)算研究區(qū)水資源量。

2.2.1 SWAT模型基本原理

根據(jù)水循環(huán)原理，SWAT模型將隔水層以上的非飽和帶以及淺層含水層視為一個(gè)整體，其水量平衡方程為

(1)

式中：SWt、SW0分別為時(shí)段末和時(shí)段初土壤含水量(mm)；t為時(shí)間步長(zhǎng)(d)；Rday,i為第i天的降雨量(mm)；Qsurf,i為第i天的地表徑流量(mm)；Ea,i為第i天的蒸散發(fā)量(mm)；Qlat,i為第i天的側(cè)向壤中流流量(mm)；Qgw,i為第i天的地下水回歸流量(mm)。

(1) 土壤水均衡方程式為

ΔS=PREC-SURQ-PERC-LATQ-ET

(2)

式中：ΔS為土壤水儲(chǔ)存量變化量(mm)；PREC為降雨量(mm)；SURQ為地表徑流量(mm)；PERC為土壤對(duì)地下水的補(bǔ)給量(mm)；LATQ為側(cè)向流量(mm)；ET為實(shí)際蒸散發(fā)量(mm)。

(2) 河道徑流水均衡方程式為

WYLD=SURQ+GWQ+LATQ

(3)

式中：WYLD為河道徑流總量(mm)；GWQ為基流量(mm)。

(3) 地下水水均衡方程式為

ΔS=PERC-GWQ-REAVP-DARCHG

(4)

式中：ΔS為淺層含水層儲(chǔ)存量變化量(mm)；REAVP為淺層含水層流向包氣帶的流量(mm)；DARCHG為深層含水層的滲漏量(mm)。

上述公式(1)至(4)中的各項(xiàng)公式可以參考文獻(xiàn)[18]。

2.2.2 評(píng)價(jià)模型適用性分析方法

本文選用Nash-Sutcliffe確定性系數(shù)Ens和相關(guān)系數(shù)R2來(lái)評(píng)價(jià)模型的適用性，具體計(jì)算公式如下：

(5)

(6)

Ens和R2因子評(píng)價(jià)模型的可信度分布區(qū)間，見(jiàn)表2。

表2 Ens和R2因子評(píng)價(jià)模型的可信度分布區(qū)間

2.2.3 模型中水資源量計(jì)算方法

根據(jù)子流域輸出文件分別將所有子流域地表水各項(xiàng)水資源量分別乘以對(duì)應(yīng)子流域的面積，并將所有子流域各項(xiàng)計(jì)算結(jié)果進(jìn)行相加，可求得整個(gè)研究區(qū)地表水資源量；地下水補(bǔ)給資源量只需計(jì)算每個(gè)子流域土壤對(duì)地下水的補(bǔ)給量分別乘以對(duì)應(yīng)子流域面積，再將所有子流域各項(xiàng)計(jì)算結(jié)果相加即可。

3 結(jié)果分析

3.1 SWAT模型模擬過(guò)程及適用性評(píng)價(jià)

在建模中，由于缺乏洞庭湖流域模擬期內(nèi)的水位數(shù)據(jù)，根據(jù)劉萬(wàn)千等[19]提到洞庭湖流域年均水位未見(jiàn)趨勢(shì)性變化，湖區(qū)面積變化不大且對(duì)整個(gè)區(qū)域的影響較小，因此本文在模擬中將湖區(qū)概化成河道，未考慮洞庭湖流域水位和面積對(duì)徑流的調(diào)節(jié)作用。

加載DEM數(shù)據(jù)并設(shè)置掩膜區(qū)域，添加burn-in對(duì)河流進(jìn)行修正，設(shè)置集水面積閾值為5 000 hm2，添加水文站站點(diǎn)信息，設(shè)置城陵磯水文站點(diǎn)為模型總出水口，自動(dòng)生成模型邊界，最終將研究區(qū)共劃分為43個(gè)子流域(見(jiàn)圖3)。加載土壤類(lèi)型和土地利用類(lèi)型數(shù)據(jù)，將坡度分為兩級(jí)(0%～20%和大于20%)，對(duì)土壤、土地利用類(lèi)型和坡度進(jìn)行重分類(lèi)；設(shè)置土地利用類(lèi)型比例閾值為5%、土壤類(lèi)型比例閾值為10%、坡度分類(lèi)閾值為10%，進(jìn)行空間疊加，劃分得到444個(gè)水文響應(yīng)單元。

圖3 研究區(qū)邊界生成及子流域劃分圖

3.1.1 參數(shù)敏感性分析

為了降低參數(shù)率定的復(fù)雜程度，通過(guò)全局參數(shù)敏感性分析來(lái)挑選對(duì)模擬結(jié)果影響較大的參數(shù)。利用輔助軟件SWAT-CUP2012進(jìn)行全局參數(shù)敏感性分析，選取8個(gè)敏感性較大的參數(shù)進(jìn)行率定，通過(guò)參數(shù)率定后所得的模型滿足精度要求，最終得到率定參數(shù)的取值，見(jiàn)表3。

表3 率定參數(shù)的取值

3.1.2 模型率定及驗(yàn)證

將率定后的參數(shù)代入到構(gòu)建好的SWAT模型中，并將2010—2012年設(shè)為率定期，2012—2013年設(shè)為驗(yàn)證期，利用石龜山、南縣、南咀、城陵磯水文站2010—2013年實(shí)測(cè)流量數(shù)據(jù)對(duì)率定后的徑流模擬值與水文站觀測(cè)值進(jìn)行對(duì)比，得出典型水文站點(diǎn)率定期和驗(yàn)證期流量對(duì)比曲線(見(jiàn)圖4)，并計(jì)算了模擬精度(見(jiàn)表4)?？紤]到城陵磯水文站的徑流過(guò)程還受到洞庭湖上游、長(zhǎng)江三口來(lái)水的影響，所以在徑流模擬的過(guò)程中，將長(zhǎng)江三口水文站(彌陀寺站、新江口站、管家鋪站、沙道觀站、康家崗站)和湘、資、沅、澧四水水文站(湘潭水文站、桃江水文站、桃源水文站、石門(mén)水文站)的流量作為境外輸入量進(jìn)行水資源量計(jì)算。

由表4和圖4可知：各水文站點(diǎn)率定期和驗(yàn)證期的相關(guān)系數(shù)R2均大于0.7，確定性系數(shù)Ens均大于0.65，達(dá)到模型模擬精度且在誤差允許范圍內(nèi)，整體擬合效果良好，說(shuō)明SWAT模型適用于洞庭平原徑流模擬；石龜山和南咀水文站點(diǎn)的擬合效果較好，南縣水文站點(diǎn)的擬合效果較差，但是整體上枯水期的擬合效果要優(yōu)于豐水期。

圖4 典型水文站點(diǎn)率定期和驗(yàn)證期流量的對(duì)比曲線

表4 典型水文站點(diǎn)率定期和驗(yàn)證期擬合結(jié)果的精度評(píng)價(jià)

3.2 水資源量計(jì)算

3.2.1 降雨量計(jì)算

通過(guò)將模型校準(zhǔn)率定后的參數(shù)代入到SWAT模型中再次進(jìn)行模擬，得出模型運(yùn)行結(jié)果，并依據(jù)土壤水均衡方程[見(jiàn)式(2)]，計(jì)算2010—2013年洞庭湖平原土壤水均衡各項(xiàng)水資源量，其計(jì)算結(jié)果見(jiàn)表5。

表5 2010—2013年洞庭湖平原土壤水均衡各項(xiàng)水資源量計(jì)算結(jié)果(單位：×108 m3)

由表5可知：洞庭湖平原年均降雨量為1 186.24 mm，約為254.27×108m3，土壤水呈負(fù)均衡；在模擬期間2011年為降雨量最小的年份，2012年為降雨量最大的年份，2012年與2011年降雨量之差為148.41×108m3，相比2011年,2012年地表徑流量增加了37.97×108m3，土壤水對(duì)地下水的補(bǔ)給量增加了85.07×108m3，兩者增加的量分別占2011年地表徑流量和土壤對(duì)地下水補(bǔ)給量的133.65%和210.83%，說(shuō)明洞庭湖平原的大氣降雨主要以蒸散發(fā)的形式排泄，而且在降雨量變化幅度較大時(shí)流域內(nèi)土壤對(duì)地下水的補(bǔ)給比地表徑流受降雨的影響也更大。

3.2.2 地表水資源量計(jì)算

研究區(qū)內(nèi)地表水總資源量包括大氣降雨經(jīng)過(guò)地表產(chǎn)流流入河流的流量和土壤水側(cè)向流入河流的流量以及區(qū)域總基流量，2010—2013年洞庭湖平原河道徑流水均衡各項(xiàng)水資源量的計(jì)算結(jié)果，見(jiàn)表6。

表6 2010—2013年洞庭湖平原河道總徑流水均衡各項(xiàng)水資源量計(jì)算結(jié)果(×108 m3)

由表6可知：洞庭湖平原地表水資源量主要為地表徑流量和基流量，地表水資源總量為120.84×108m3，其中地表徑流量為49.69×108m3，占地表水資源量的41.12%，基流量為69.03×108m3，占地表水資源量的57.13%，而側(cè)向流量為2.12×108m3，占比僅為地表水資源量的1.75%。由此可見(jiàn)，流域內(nèi)基流對(duì)河道徑流的補(bǔ)給量大于地表徑流對(duì)河道徑流的補(bǔ)給量。

3.2.3 地下水資源量計(jì)算

利用公式(4)計(jì)算研究區(qū)地下水資源量，得到2010—2013年洞庭湖平原地下水水均衡各項(xiàng)水資源量，見(jiàn)表7。

表7 2010—2013年洞庭湖平原地下水水均衡各項(xiàng)水資源量計(jì)算結(jié)果(×108 m3)

由表7可知：土壤對(duì)地下水的補(bǔ)給是研究區(qū)地下水補(bǔ)給的主要來(lái)源，模擬期內(nèi)年均地下水資源總量為78.64×108m3，其中年均基流量為69.09×108m3，占地下水資源量的87.77%，而淺層含水層流向包氣帶的流量和深層含水層的滲漏量總共占地下水資源量的7.82%。土壤對(duì)地下水的補(bǔ)給量與基流量成正相關(guān)關(guān)系，在整個(gè)模擬期內(nèi)2012年的地下水總補(bǔ)給量和基流量達(dá)到最大；從整個(gè)模擬期來(lái)看，區(qū)域內(nèi)年均淺層含水層儲(chǔ)存量變化量為3.45×108m3，占地下水總補(bǔ)給量的4.41%，說(shuō)明地下水均衡呈正均衡。

3.3 氣候變化對(duì)水資源量的影響分析

根據(jù)聯(lián)合國(guó)政府間氣候變化專(zhuān)門(mén)委員會(huì)(IPCC)第五次評(píng)估報(bào)告，并參考相關(guān)文獻(xiàn)[20-23]，預(yù)測(cè)并設(shè)置了2030年氣候情景：降雨量在原來(lái)降雨量(P)基礎(chǔ)上增加或減少10%，氣溫在原來(lái)氣溫(T)基礎(chǔ)上升高或降低1℃，降雨量和氣溫變化兩兩組合(不考慮單一變化)共5種氣候情景組合模式(包括基礎(chǔ)情景C3)，具體氣溫和降雨量情景組合見(jiàn)表8。

表8 氣候情景設(shè)置表

輸入5種氣候情景模式對(duì)應(yīng)的氣溫和降雨量數(shù)據(jù)，保持其他輸入不變，建立模型并進(jìn)行徑流模擬，得到水均衡各項(xiàng)的年均流量變化(見(jiàn)圖5)，并計(jì)算不同氣候情景模式下的土壤水、地表水和地下水資源量，其計(jì)算結(jié)果見(jiàn)表9至表12。

圖5 不同氣候情景模式下洞庭湖平原水均衡各項(xiàng)年均流量變化

由圖5和表9至表12可知：研究區(qū)地表徑流量、基流量、河道徑流總量與氣溫之間呈負(fù)相關(guān)關(guān)系，與降雨量之間呈正相關(guān)關(guān)系，而側(cè)向流量、淺層含水層流向包氣帶的流量和深層含水層滲漏量的變化均較小。不同降雨量和氣溫的組合情景下對(duì)研究區(qū)水資源量的影響均不相同：相對(duì)基礎(chǔ)情景C3,氣候情景C5下均衡期內(nèi)年均降雨量增加最多，增加了30.42×108m3，而C1減少最多；相對(duì)基礎(chǔ)情景C3,氣候情景C2下地表水資源量增加最多，增加了23.77×108m3，地下水資源量變化最大，增加了4.20×108m3，而氣候情景C4下地表水資源量和地下水資源量減少最多，分別減少了19.18×108m3、2.70×108m3。

表9 氣候變化對(duì)洞庭湖平原土壤水均衡各項(xiàng)水資源量的影響(×108 m3)

表10 氣候變化對(duì)洞庭湖平原河道徑流水均衡各項(xiàng)水資源量的影響(×108 m3)

表11 氣候變化對(duì)洞庭湖平原地下水均衡各項(xiàng)水資源量的影響(×108 m3)

表12 氣候變化下洞庭湖平原水資源量變化量(×108 m3)

氣溫對(duì)水均衡各項(xiàng)水資源量的影響比降雨量小，水均衡各項(xiàng)水資源量隨著降雨量的增加而增大，隨著降雨量的減少而減小。當(dāng)降雨量均增加或者減少，氣溫減小或增加1℃時(shí)，水均衡各項(xiàng)水資源量的變化幅度小，例如氣候情景C1與C4、C2與C5；當(dāng)氣溫均增加或減少，降雨量增加或減少10%時(shí)，水均衡各項(xiàng)水資源量的變化幅度大，例如氣候情景C1與C2、C4與C5。因此，影響研究區(qū)水均衡對(duì)氣候變化響應(yīng)的主要因素為降雨量的變化。

4 結(jié) 論

本文利用ArcSWAT構(gòu)建洞庭湖平原分布式水文模型，通過(guò)對(duì)模型進(jìn)行率定和驗(yàn)證，對(duì)研究區(qū)進(jìn)行了水均衡分析，計(jì)算了水資源量，并分析了氣候變化對(duì)水資源量的影響，得出以下主要結(jié)論：

(1) 根據(jù)模型率定和驗(yàn)證結(jié)果證明，在誤差允許范圍內(nèi)該模型較適用于洞庭湖平原；根據(jù)參數(shù)敏感性分析結(jié)果表明，對(duì)洞庭湖平原徑流過(guò)程影響較大的參數(shù)主要為CN2(SCS徑流曲線數(shù))、ALPHA_BF(基流衰減因子)、ESCO(土壤補(bǔ)償蒸發(fā)因子)、SOL_AWC(土壤有效含水量)，其余的參數(shù)相對(duì)來(lái)說(shuō)影響較小。

(2) 從土壤水、地表水(河道徑流)、地下水三個(gè)方面分析了洞庭湖平原4年的年平均水均衡特征并計(jì)算了年均水資源量，結(jié)果表明：研究區(qū)年均降雨量為254.27×108m3，大氣降雨以蒸散發(fā)的形式排泄，而且在降雨量變化幅度較大時(shí)研究區(qū)內(nèi)土壤對(duì)地下水的補(bǔ)給比地表徑流受降雨的影響更大；區(qū)內(nèi)年均地表水資源量為120.84×108m3，其中流域內(nèi)基流量對(duì)河道徑流的補(bǔ)給量大于地表徑流對(duì)河道的補(bǔ)給量；土壤對(duì)地下水的補(bǔ)給是研究區(qū)地下水補(bǔ)給的主要來(lái)源，模擬期內(nèi)年均土壤對(duì)地下水的補(bǔ)給量為78.64×108m3，主要排泄方式是地下水對(duì)河流的補(bǔ)給，地下水呈正均衡。

(3) 在5種氣候變化情景下，地表徑流量、基流量與降雨量成正相關(guān)關(guān)系，與氣溫成負(fù)相關(guān)關(guān)系。氣溫對(duì)地表徑流量、河道徑流總量、基流量等其他水均衡各項(xiàng)水資源量的影響與降雨量相關(guān)，即隨降雨量的增加而變大，隨降雨量的減少而減小，因此影響研究區(qū)年徑流量變化的主要因素為降雨量。

(4) 不同氣候變化情景下降雨量、地表水資源量和地下水資源量的變化明顯，當(dāng)降雨量減少10%(年均降雨量減少118.62 mm)、溫度升高1℃時(shí)，年均地表水資源量減少19.18×108m3，地下水資源量減少2.7×108m3；當(dāng)降雨量增加10%(年均降雨量增加118.62 mm)、溫度下降1℃時(shí)，年均地表水資源量增加23.77×108m3，地下水資源量增加4.2×108m3。

本文雖然模擬了洞庭湖平原的徑流過(guò)程并計(jì)算了水資源量，但是未考慮洞庭湖水位和面積對(duì)徑流的調(diào)節(jié)作用，這將是今后研究的方向。