夏智宏，劉 敏，王 苗，王 凱，秦鵬程

(武漢區(qū)域氣候中心，武漢430074)

隨著全球氣候變化和人類活動影響的加劇，水資源短缺和環(huán)境惡化已經(jīng)成為世界性的問題，是當前人類面臨的嚴峻挑戰(zhàn)之一［1］.對洪湖流域而言，近50年來，年平均氣溫呈極顯著的上升趨勢;年降水量呈緩慢增加趨勢，但時空分布不均;春、夏、冬季和年蒸發(fā)量均呈上升趨勢;極端旱澇事件交替轉(zhuǎn)換周期縮至4年左右［2］.人類活動同樣顯著影響著洪湖流域水資源與生態(tài)環(huán)境，1950s中期到1970s末大規(guī)模的水利建設(shè)使流域水系河網(wǎng)發(fā)生重大變化，1980s末圍網(wǎng)養(yǎng)殖致使流域水質(zhì)逐漸惡化，進入1990s城鎮(zhèn)化快速發(fā)展使流域不透水層范圍不斷擴大.氣候變化和人類活動的雙重影響已經(jīng)導(dǎo)致洪湖流域湖泊面積縮小，旱澇頻率明顯升高，枯水期提前，時間延長，水位波幅增大，水質(zhì)富營養(yǎng)化，進而威脅到洪湖流域生態(tài)安全.如1998和2010年的洪澇以及2011年的冬春季連旱，使洪湖流域及其生態(tài)系統(tǒng)遭受極大危害.因此客觀定量地評價氣候變化和人類活動對洪湖流域水資源的影響有助于未來水資源的合理開發(fā)和可持續(xù)利用，對于適應(yīng)氣候變化的對策研究也具有重要意義.有關(guān)氣候變化和人類活動對洪湖流域影響的研究，大多是從流域水資源、水環(huán)境、河湖環(huán)境演變、土地覆被、氣象災(zāi)害對氣候變化和人類活動的響應(yīng)以及生態(tài)環(huán)境綜合評價等方面著手［3-10］，缺乏氣候變化和人類活動對洪湖流域水資源影響貢獻率定量分析的研究成果.國內(nèi)外有關(guān)氣候變化和人類活動對流域水資源影響貢獻的區(qū)分量化研究，采取的主要技術(shù)手段有長序列資料對比分析法、試驗對比分析法及流域水文模擬等［11］.由于水文模擬法具有能彌補缺乏長序列實測徑流資料和能保證評價期與基準期自然徑流成因一致等特點［12］，因此基于統(tǒng)計和物理的水文模擬成為目前區(qū)分量化氣候變化和人類活動對河川徑流影響貢獻的主要方法，并廣泛應(yīng)用于我國的北方地區(qū)，但在長江中下游地區(qū)應(yīng)用較少［13-20］.本文選擇長江中下游具有代表性的洪湖流域為研究對象，采用基于SWAT模型的徑流影響分割分析法，定量辨識1990s洪湖流域城鎮(zhèn)化迅速發(fā)展以來，氣候變化和人類活動對洪湖流域地表徑流的影響程度，為保護洪湖流域水資源、促進其生態(tài)可持續(xù)發(fā)展提供科學(xué)依據(jù).

1 研究區(qū)概況

洪湖流域位于亞熱帶中緯度邊緣，流域面積11152 km2(圖1)，為典型的北亞熱帶濕潤季風氣候區(qū)，光熱充足，雨熱同季，無霜期長，適宜多種農(nóng)作物生長，有“魚米之鄉(xiāng)”的美稱.常年年均氣溫17.0℃，年均降水量1278.3 mm，無霜期250 d以上，夏半年(4-10月)降水量占全年降水總量的77%［21-22］.降水形成的地表徑流與地下水以及長江與漢水帶來的過境客水是洪湖流域水資源的主要構(gòu)成，年均地表徑流量達到36.40×108m3.

圖1 洪湖流域地理位置Fig.1 Location of Lake Honghu basin

洪湖是洪湖流域最大的淡水湖，也是中國7大淡水湖之一，湖面面積348.2 km2，位于長江中游北岸，長湖、三湖、白鷺湖的下游，通過四湖總干渠承接這3個湖泊與區(qū)間來水，是長江和漢水支流東荊河之間的洼地壅塞湖，最高水位24.58 ～27.18 m，最低水位22.87 ～23.92 m，多年平均水位24.43 m，年水位變幅為 1 ～4 m［23］.洪湖具有調(diào)蓄、供水、灌溉、濕地生物棲息地、養(yǎng)殖、航運和旅游等多種功能，在長江中下游淺水湖泊濕地中具有代表性.

2 資料收集與處理

2.1 地理空間數(shù)據(jù)

2.1.1 地形數(shù)據(jù) 洪湖流域地形數(shù)據(jù)來源于國家測繪與地理信息局1∶5萬基本比例尺地形圖，用于提取水系、匯流區(qū)域、地形特征參數(shù)等信息.

2.1.2 土地利用/覆被數(shù)據(jù) 洪湖流域土地利用/覆被數(shù)據(jù)來源于中國科學(xué)院武漢地球物理與測量研究所.包括比例尺為1∶25萬的1980年和2005年土地利用/覆被數(shù)據(jù)、比例尺為1∶5萬的1996年土地利用/覆被數(shù)據(jù).本研究中，分別以1980、1996、2005年的土地利用/覆被數(shù)據(jù)表征1980s、1990s、2000s的土地利用/覆被情況.

2.1.3 土壤數(shù)據(jù) 洪湖流域土壤數(shù)據(jù)來源于中國科學(xué)院南京土壤研究所1995年編制的1∶100萬中國土壤數(shù)據(jù)庫.考慮到其土壤物理特性與SWAT模型采用的美制標準的不同(如粒徑等級)，在充分采納前人的經(jīng)驗和知識的前提下［24-25］，將土壤屬性轉(zhuǎn)換至美制標準.

2.2 氣象數(shù)據(jù)

洪湖流域氣象數(shù)據(jù)來源于湖北省氣象信息保障中心，包括石首、監(jiān)利、洪湖等11個氣象臺站1961-2011年逐日平均氣溫、最高/最低氣溫、降水、相對濕度、風速等要素，用于獲取SWAT天氣發(fā)生器并驅(qū)動SWAT模型進行徑流模擬.

2.3 水文數(shù)據(jù)

水文數(shù)據(jù)來源于湖北省水文局，包括習(xí)家口水文站1964-1974年、2008-2011年逐日流量和福田寺水文站2008-2010年逐日流量(圖1).考慮到長湖對洪湖流域中下游的調(diào)蓄作用，利用習(xí)家口水文站流量數(shù)據(jù)計算長湖的調(diào)蓄參數(shù)，以反映人為調(diào)蓄對洪湖濕地水資源模擬的影響.福田寺水文站流量數(shù)據(jù)則用于率定模型的參數(shù)和評價模型在洪湖流域的適用性.

3 研究方法

根據(jù)研究區(qū)域資料情況，本文采用基于SWAT模型的水文模擬分割分析法對氣候變化與人類活動對洪湖流域地表徑流影響貢獻率進行辨識.

3.1 SWAT模型徑流模擬計算方法

SWAT模型的產(chǎn)流模擬不同于一般的統(tǒng)計模型，它包括氣溫控制降水形態(tài)、冰雪融化、蒸散發(fā)、地下徑流補給等物理過程，并最終導(dǎo)致產(chǎn)流量的變化，因此具有很強的物理機制，適用于模擬不同土地利用及多種土地管理措施對水文的長期影響，其水文計算基于如下水量平衡方程:

式中，swt為土壤最終含水量(mm)，sw0為土壤前期含水量(mm)，t為時間步長(d)，Rday為第 i日降水量(mm)，Qsurf為第i日的地表徑流(mm)，Ea為第i日的蒸發(fā)量(mm)，Wseep為第i日土壤剖面地層的滲透量和側(cè)流量(mm)，Qgw為第 i日的基流量(mm)［26-27］.

SWAT模型適用性評價采用Nash模型效率系數(shù)Ens［28］和相關(guān)系數(shù)R2來評價［29］:

式中，Ens為 Nash-Sutcliffe 效率系數(shù)，R2為相關(guān)系數(shù)平方，Qobs，i為天然流量，Qsim，i為模擬流量，Qobs為多年平均實測流量為多年平均模擬流量，n為模擬流量序列長度.

3.2 徑流影響分割分析方法

本文將人類活動和氣候變化視為影響該洪湖濕地地表徑流變化的兩個相互獨立因子，其中人類活動以土地利用/覆被變化表征.選用SWAT模型，應(yīng)用實測徑流數(shù)據(jù)率定模型參數(shù)，獲取基準期的模擬徑流深WB.保持基準期土地利用/覆被類型不變，輸入評價期的氣候要素驅(qū)動SWAT模型，獲得僅考慮氣候變化影響的評價期模擬徑流深WHN;同時輸入評價期氣候要素和土地利用/覆被資料驅(qū)動SWAT模型，獲取考慮氣候變化和人類活動共同影響的評價期模擬徑流深WHR，采用分割分析方法(如公式(4)～(8))分別得到評價期氣候變化和人類活動對洪湖流域徑流影響的貢獻率ηH和ηC.

式中，ΔWT為評價期的人類活動和氣候變化共同作用對徑流深的影響量，ΔWH為評價期的人類活動對徑流深的影響量，ΔWC為評價期的氣候變化對徑流深的影響量.

4 結(jié)果與分析

4.1 水文模型適用性評價

應(yīng)用AVSWAT2005自帶的基于LH-OAT方法的敏感性分析工具［30］，選取27個模型輸入?yún)?shù)進行參數(shù)敏感性分析試驗，設(shè)置參數(shù)值采樣間隔為10，進行了280次比較判斷，確定了敏感性等級最高的5個參數(shù)，分別為 CN2、SOL_AWC、ESCO、GW_REVAP 和 GW_DELAY.

利用2000s土地利用/覆被資料，選用2008-2009年福田寺水文站徑流實測數(shù)據(jù)，根據(jù)“先調(diào)整水量平衡，再調(diào)整過程;先調(diào)整地表徑流，再調(diào)整土壤水、蒸散發(fā)和地下徑流”的原則對模型參數(shù)進行逐月水量平衡校正，模型率定后的模擬值和實測值擬合很好(圖2)，率定結(jié)果顯示R2=0.94，Ens=0.86.選用2010年福田寺水文站徑流實測數(shù)據(jù)對模型進行驗證，逐月模擬結(jié)果顯示R2=0.94，Ens=0.88，表明SWAT模型適用于洪湖濕地水資源模擬.參數(shù)率定結(jié)果如表1所示.

圖2 模型率定期與驗證期實測與模擬徑流量對比Fig.2 Comparing runoff of calibration period with verification period

表1 SWAT模型率定的主要參數(shù)Tab.1 Main calibration parameters in SWAT model

4.2 氣候變化和人類活動對徑流影響辨識

考慮到1990s以來洪湖流域城鎮(zhèn)化進程加快，人類活動以城鎮(zhèn)化帶來的土地覆被類型變化為主要表現(xiàn)形式，因此本文用洪湖流域土地利用/覆被變化表征人類活動的影響.以1980s作為基準期，采用基于SWAT模型的分割分析法分階段辨識1990-2011年氣候變化和人類活動對洪湖流域徑流深影響的貢獻率.結(jié)果表明，受環(huán)境變化的影響，相對基準期而言，1990s徑流深呈增加趨勢，2000s呈減少趨勢，減少量是1990s增加量的1.7倍.不同時期人類活動和氣候變化對徑流深影響的貢獻率不同，1990s氣候變化對徑流深的影響貢獻率高于人類活動，受其影響的徑流深增加量占總增加量的52.35%;2000s人類活動影響的貢獻率高于氣候變化，受其影響的徑流深減少量占總減少量的55.62%.2000年以后氣候因素和人類活動在對徑流深的絕對影響量均高于1990s.總體上看，1990-2011年，人類活動是地表徑流減少的主要原因，其影響的貢獻率占徑流減少量的63.72%，氣候變化影響的貢獻率占36.28%(表2).

表2 1990s和2000s人類活動與氣候變化對洪湖流域地表徑流變化的影響*Tab.2 The impact of climate change and human activity on runoff in 1990s and 2000s in Lake Honghu basin

4.3 土地利用/覆蓋變化對徑流影響分析

圖3 1980s、1990s和2000s洪湖流域土地利用空間分布Fig.3 Spatial distribution of landuse in 1980s，1990s and 2000s in Lake Honghu basin

分別提取洪湖流域1980s、1990s和2000s土地利用/覆被數(shù)據(jù)(圖3)，研究土地利用類型變化與徑流變化的關(guān)系.結(jié)果表明(表3)，1990s相對1980s，農(nóng)村居民地、林地、水體、城鎮(zhèn)建設(shè)用地面積分別增長146.53%、16.44%、17.42% 和 6.33%，灘地、沼澤地、耕地、草地分別減少了 25.79%、25.29%、8.86% 和8.15%，徑流量增加了3.86%.分析其原因，一方面是由于降水量的增加，另一方面，由于人口增長加速、經(jīng)濟擴張式發(fā)展和城市化進程加快，大量耕地和部分林地、草地、水體、沼澤地轉(zhuǎn)變?yōu)槌擎?zhèn)建設(shè)用地和農(nóng)村居民用地，導(dǎo)致洪湖流域不透水層面積增加，抵消了由林地和水體面積增加引起的徑流減少的幅度.2000s相對1990s，草地、水體、灘地、城鎮(zhèn)建設(shè)用地和沼澤地面積分別增加了 40.27%、22.21%、11.06%、7.65%和6.28%，農(nóng)村居民點、林地和耕地面積分別減少了58.25%、13.13%和0.98%，徑流量減少8.92%.分析其原因，一方面是由于降水量減少，另一方面是由于草地、水體、沼澤地覆蓋率的增加使土壤貯水能力增強、攔截作用加大，盡管城鎮(zhèn)建設(shè)用地面積繼續(xù)擴張，但是流域前期的水土保持措施的蓄水功能效益得到發(fā)揮，抵消了城鎮(zhèn)化引起的徑流增加.

總體上看，近20年來洪湖流域水土保持工作得到加強，洪湖流域的水體、草地面積大幅增加，林地的面積得到有效恢復(fù)與保護，盡管城市建設(shè)用地持續(xù)增加，但徑流量仍然下降了5.41%，可見水土保持措施發(fā)揮了較好的徑流調(diào)節(jié)和保水效益.

表3 洪湖流域不同年代土地利用與徑流變化(%)Tab.3 The rate of runoff and landuse change in different decades in Lake Honghu basin

5 結(jié)論與討論

1)SWAT模型在洪湖流域徑流模擬驗證期的效率系數(shù)和相關(guān)系數(shù)平方分別達到0.88和0.94，表明該模型適用于洪湖流域徑流模擬，可用于辨識氣候變化和人類活動對洪湖流域徑流的影響.

2)受氣候變化和人類活動共同影響，相對于基準期，1990-2011年洪湖流域年均地表徑流深減少了11.63 mm.其中1990s徑流深呈增加趨勢;而2000s呈減少趨勢，減少量是1990s增加量的1.7倍.

3)近20年來，人類活動是洪湖流域地表徑流減少的主要原因，其影響的貢獻率占徑流減少量的63.72%，氣候變化影響貢獻率占36.28%.但不同階段人類活動與氣候變化對洪湖流域地表徑流影響的程度不同，1990s氣候變化對洪湖流域地表徑流影響的貢獻率高于人類活動，2000s則相反.

4)相對于基準期，近20年來洪湖濕地的水體、草地面積大幅增加，林地的面積得到有效恢復(fù)與保護.盡管城市建設(shè)用地持續(xù)增加，但年均徑流量仍然下降了5.41%，表明水土保持措施發(fā)揮了較好的徑流調(diào)節(jié)和保水效益.

5)廣義的人類活動包括人口增長、城市化水平發(fā)展、植被增減變化、硬覆蓋變化、水利設(shè)施對截留的影響、流域內(nèi)水體小循環(huán)等，本文為研究城鎮(zhèn)化快速發(fā)展階段的人類活動對洪湖流域水資源的影響，以土地覆被類型的變化來表征人類活動，在后續(xù)研究中將進一步區(qū)分不同人類活動對徑流變化的影響.

［1］李子君，李秀彬.近45年來降水變化和人類活動對潮河流域年徑流量的影響.地理科學(xué)，2008，28(6):809-813.

［2］武漢區(qū)域氣候中心.中國氣象局氣候變化專項(CCSF2012-05)“氣候變化對長江中游典型湖泊濕地的影響”，2012.

［3］蔡述明，王學(xué)雷.江漢平原四湖地區(qū)生態(tài)環(huán)境綜合評價.長江流域資源與環(huán)境，1993，2(4):356-364.

［4］李昌峰，張鴻輝.建國以來人類活動對湖北省四湖地區(qū)水環(huán)境的影響研究.國土與自然資源研究，2003，(4):68-70.

［5］李昌峰，高俊峰，張鴻輝.近50年來人類活動對四湖地區(qū)河湖環(huán)境演變的影響.地域研究與開發(fā)，2004，23(5):120-124.

［6］何 浩，劉章勇.自然與人為干擾對四湖地區(qū)澇漬生態(tài)系統(tǒng)的影響.長江大學(xué)學(xué)報:自科版，2006，3(3):179-182.

［7］桂 峰，于 革，賴格英.洪湖流域自然農(nóng)耕條件下營養(yǎng)鹽沉積輸移演化模擬研究.沉積學(xué)報，2006，24(3):334-338.

［8］桂 峰，于 革.20世紀80年代以來洪湖流域營養(yǎng)物質(zhì)輸移模擬研究.長江流域資源與環(huán)境，2007，16(22):179-182.

［9］吳啟俠，朱建強，耿顯波.江漢平原四湖流域降水特征分析.中國農(nóng)業(yè)氣象，2008，29(2):146-150.

［10］金衛(wèi)斌，熊勤學(xué)，薛 蓮.基于MODIS-EVI時序數(shù)據(jù)的江漢平原四湖地區(qū)土地覆蓋動態(tài)分析.湖北農(nóng)業(yè)科學(xué)，2011，50(11):2220-2224.

［11］張建云，王國慶.氣候變化對水文水資源影響研究.北京:科學(xué)出版社，2007.

［12］張愛靜，王本德，曹明亮.氣候變化和人類活動對徑流影響貢獻的研究概況.東北水利水電，2012，(1):6-9.

［13］劉吉峰，霍世青，李世杰等.SWAT模型在青海湖布哈河流域徑流變化成因分析中的應(yīng)用.河海大學(xué)學(xué)報:自然科學(xué)版，2007，35(2):159-163.

［14］胡宏昌，王根緒，李 志等.渭河徑流對LUCC和氣候波動的響應(yīng)研究.人民黃河，2008，30(7):25-26.

［15］王國慶，張建云，劉九夫等.氣候變化和人類活動對河川徑流影響的定量分析.中國水利，2008，(2):55-58.

［16］黃玉霞，王寶鑒，張 強等.氣候變化和人類活動對石羊河流域水資影響評價.高原氣象，2008，27(4):867-872.

［17］葉許春，張 奇，劉 健等.氣候變化和人類活動對鄱陽湖流域徑流變化的影響研究.冰川凍土，2009，31(5):836-842.

［18］江善虎，任立良，雍 斌等.氣候變化和人類活動對老哈河流域徑流的影響.水資源保護，2010，26(6):1-4.

［19］朱穎潔，郭純青，黃夏坤.氣候變化和人類活動影響下西江梧州降水徑流演變研究.水文，2010，30(3):50-55.

［20］王 俊，陳亞寧，陳忠升.氣候變化與人類活動對博斯騰湖入湖徑流影響的定量分析.新疆農(nóng)業(yè)科學(xué)，2012，49(3):581-587.

［21］項國榮.四湖地區(qū)濕地農(nóng)業(yè)持續(xù)發(fā)展研究.北京:科學(xué)出版社，1997:271.

［22］王蘇民，竇鴻身.中國湖泊志.北京:科學(xué)出版社，1998:580.

［23］王學(xué)雷，寧龍梅，肖 銳.洪湖濕地恢復(fù)中的生態(tài)水位控制與江湖聯(lián)系研究.濕地科學(xué)，2008，6(2):316-320.

［24］龔子同，張甘霖，陳志誠等.以中國土壤系統(tǒng)分類為基礎(chǔ)的土壤參比.土壤通報，2002，33(1):1-5.

［25］王慶云，徐能海.湖北省土系概要.武漢:湖北科技出版社，2001.

［26］楊宏偉，許崇育.東江流域典型子流域土地利用/覆被變化對地表徑流影響.湖泊科學(xué)，2011，23(6):991-996.

［27］孫 瑞，張雪芹.基于SWAT模型的流域徑流模擬研究進展.水文，2010，30(3):29-32，47.

［28］Xie XH，Cui YL.Development and test of SWAT for modeling hydrological processes in irrigation districts with paddy rice.Journal of Hydrology，2011，396(1/2):61-71.

［29］夏智宏，周月華，許紅梅.基于SWAT模型的漢江流域水資源對氣候變化的響應(yīng).長江流域資源與環(huán)境，2010，19(2):158-163.

［30］Mulungu DMM，Munishi SE.River catchment parameterization using SWAT model.Physics and Chemistry of the Earth，2007，32(15/16/17/18):1032-1039.