楊 柳， 陳興偉， 許有鵬， 劉梅冰， 胡春生

(1.南京大學(xué) 地理與海洋科學(xué)學(xué)院， 江蘇 南京210093; 2.福建師范大學(xué) 地理科學(xué)學(xué)院， 福建 福州 350007)

東南沿海地區(qū)小流域土地利用/覆被變化的水文效應(yīng)

楊 柳1,2， 陳興偉2， 許有鵬1， 劉梅冰2， 胡春生1

(1.南京大學(xué) 地理與海洋科學(xué)學(xué)院， 江蘇 南京210093; 2.福建師范大學(xué) 地理科學(xué)學(xué)院， 福建 福州 350007)

資助項目：福建省教育廳項目“流域施肥的時空分配對山美水庫N污染特性的影響”(JB12038); 福建科技廳省屬公益類科研院所青年科研項目(2014R1034-4); 國家自然科學(xué)基金重點項目(40730635)

第一作者：楊柳(1986—),女(漢族),湖南省永州市人,博士研究生,研究方向為水文水資源環(huán)境。E-mail:palapalayl2006@126.com。

摘要：[目的] 研究我國東南沿海小流域土地利用/覆被變化的水文效應(yīng)，為該區(qū)水土資源合理配置和可持續(xù)發(fā)展提供科學(xué)依據(jù)。[方法] 以東南沿海福建省晉江山美水庫集水區(qū)為例，運用SWAT水文模型。[結(jié)果] (1) 研究區(qū)土地利用變化表現(xiàn)為林地轉(zhuǎn)向園地和建設(shè)用地、耕地轉(zhuǎn)向建設(shè)用地的特征。 (2) 固定其他條件，將構(gòu)建的模型2006年土地利用數(shù)據(jù)替換為1995年，對比分析發(fā)現(xiàn)2000—2010年全流域和典型子流域的年徑流深增加，洪水期平均徑流增加，枯水期徑流基本減少。 (3) 從全流域空間角度看，徑流量均增大，特別是研究區(qū)西部地區(qū)和東南部地區(qū)；同時徑流量變率增大，穩(wěn)性減小。 (4) 園地、林地、建設(shè)用地和水田是驅(qū)動區(qū)內(nèi)水文變化的主要因子。[結(jié)論] 研究區(qū)土地利用與覆被變化導(dǎo)致“洪水量更大，枯水量更小”的徑流變化趨勢。河川徑流的變化受各地類變化的組合疊加效應(yīng)影響，這種疊加效應(yīng)還須進一步研究。

關(guān)鍵詞：SWAT模型； 日徑流過程模擬； 土地利用/覆被變化； 水文效應(yīng)； 東南沿海區(qū)； 小流域

土地利用/覆被變化(LUCC)對徑流過程影響是當(dāng)前人們關(guān)注的熱點問題。SWAT模型因為其輸入輸出方便，數(shù)據(jù)獲取相對容易，可模擬長時間序列變化等優(yōu)勢[1]被廣泛應(yīng)用于LUCC水文效應(yīng)研究中。在SWAT日徑流模擬研究中，李慧[2]等在瑪納斯流域得出該模型適用于冰雪融水區(qū)的日徑流模擬；李麗嬌[3]、黃清華[4]等日徑流模擬研究表明雨季徑流量年際變化趨勢與降雨量的年際變化一致，豐水年份的模型效率較高，其次為平水年和枯水年；李占玲[5]在黑河運用WASMOD模型與SWAT模型進行日徑流模擬，比較兩種模型的優(yōu)劣等等。由此可見，目前SWAT在日徑流模擬的分析研究涉及多方面，但還鮮有在日尺度模擬基礎(chǔ)上進行LUCC水文響應(yīng)的研究。另外，從空間上看，已有日徑流模擬的研究多集中于北方地區(qū)的黃河、海河流域以及南方地區(qū)長江流域，在東南沿海諸河流域的研究較為少見[6]。而我國東南沿海地區(qū)以中小流域為主，多獨流入海，以東南沿海小流域為研究對象，具有一定的典型性[7]。本研究應(yīng)用SWAT模型，構(gòu)建山美水庫集水區(qū)日徑流分布式水文模型，分析研究區(qū)土地利用面積、地類轉(zhuǎn)移情況，繼而研究LUCC下時間(年徑流、洪水期日徑流和枯水期日徑流)和空間(徑流系數(shù)等值線)尺度上的水文響應(yīng)情況。

1研究區(qū)概況

山美水庫集水區(qū)地理位置在25°32′—25°9′N，118°3′—118°30′E之間，位于福建晉江東溪中游，集水面積1 023 km2。區(qū)內(nèi)多山地和丘陵，海拔約在40～1 360 m之間。氣候?qū)賮啛釒ШＱ笮约撅L(fēng)氣候，年均溫約為20.9 ℃；年日照時數(shù)約為1 800～2 200 h，年降雨量約為1 600 mm[8]。區(qū)內(nèi)降水量年際變化大，少雨年份降水量不及多雨年份1/2。土壤類型有紅壤、赤紅壤、黃壤、水稻土和石質(zhì)土等，其中硅鋁質(zhì)紅壤、滲育水稻土、鋁硅質(zhì)紅壤占78.2%，為區(qū)內(nèi)主要土壤類型。

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

2.1　數(shù)據(jù)來源

模型輸入數(shù)據(jù)包括空間數(shù)據(jù)和水文氣象數(shù)據(jù)?？臻g數(shù)據(jù)有土壤、DEM和土地利用數(shù)據(jù)。其中土壤資料來源于福建省1∶50萬數(shù)字化土壤數(shù)據(jù)；DEM來自國際科學(xué)數(shù)據(jù)中心的SRTM數(shù)據(jù)(分辨率為30 m×30 m)；土地利用數(shù)據(jù)通過對研究區(qū)2006年和1995年TM遙感影像進行目視解譯獲得。水文數(shù)據(jù)采用2000—2010年研究區(qū)的逐日入庫流量進行模型率定和驗證，氣象數(shù)據(jù)采用永春和德化2個氣象站2000—2010年的逐日最高、最低氣溫，風(fēng)速，相對濕度以及16個雨量站的逐日降水?dāng)?shù)據(jù)。

2.2　研究方法

2.2.1SWAT模型SWAT(soil and water assessment tool)模型是美國農(nóng)業(yè)部農(nóng)業(yè)研究所(USDA-ARS)[9]開發(fā)的運行以日為時間步長，可進行連續(xù)長時段模擬計算的分布式流域水文模型，目前已經(jīng)在不同地區(qū)得到廣泛應(yīng)用，本研究選用該模型對研究區(qū)的日徑流過程進行模擬。模型用2000—2005年的逐日入庫流量進行率定，2006—2010年日流量數(shù)據(jù)進行驗證，尋求晉江東溪流域山美水庫集水區(qū)徑流模型的最優(yōu)參數(shù)，構(gòu)建該地區(qū)日尺度徑流模擬的SWAT模型[6]。地表徑流計算選用徑流曲線數(shù)方法(SCS runoff curve number)，河道匯流過程采用馬斯京根法(Muskingum)進行演算。子流域劃分基于數(shù)字高程模型，采用“burn-in”算法。水文響應(yīng)單元劃分采用優(yōu)勢地面覆蓋/優(yōu)勢土壤類型方法，閾值設(shè)為5/10。此外，為評價模型的模擬精度，選用Nash—Suttcliffe效率系數(shù)[10](Ens)和確定性系數(shù)(R2)進行表征。

2.2.2土地利用數(shù)據(jù)解譯與變化分析首先對研究區(qū)1995年和2006年同時段TM影像進行幾何校正，再結(jié)合研究區(qū)實際與模型的需要，將土地利用類型劃分為8類，包括水田(RICE)、旱地(WWHT)、林地(FRST)、園地(ORCD)、草地(RNGE)、未利用地(SWRN)、水域(WATR)和建設(shè)用地(URMD)。為提高解譯精度，通過兩種方法進行驗證：一是在研究區(qū)中選取52個典型樣點，野外利用手持GPS定位，實地確定目標(biāo)位置的土地利用類型，通過反復(fù)修正建立起適合研究區(qū)的解譯標(biāo)志。二是按行政單位將整個永春縣2006年LUCC情況解譯，與權(quán)威部門土地利用各地類面積數(shù)據(jù)進行比對，調(diào)整解譯標(biāo)志，修正解譯結(jié)果。解譯完成后用ArcGIS 9.3對數(shù)據(jù)進行拓?fù)錂z查，用Arcview空間分析模塊分析研究區(qū)1995—2006年土地利用變化情況。

2.2.3日徑流時空響應(yīng)分析構(gòu)建研究區(qū)SWAT日徑流模型，固定其它條件，僅將土地利用數(shù)據(jù)替換為1995年數(shù)據(jù)，得出水文響應(yīng)的模型用以探討土地利用/覆被變化后年徑流和日徑流響應(yīng)情況。日徑流變化主要從洪水徑流變化和枯水徑流變化兩方面進行分析，其中，洪水徑流變化選用最大1 d和連續(xù)5 d最大的徑流深進行研究，選取后者主要是由于本流域洪水過程一般在5 d之內(nèi)，而枯水徑流變化通過常用的最小1 d和連續(xù)7 d最小的徑流深進行分析。同時考慮到整個流域LUCC的疊加效應(yīng)和復(fù)雜性，可能會對分析結(jié)果造成影響，為進一步研究土地利用/覆被變化強烈區(qū)域的水文效應(yīng)，研究在37個子流域中遴選出林地轉(zhuǎn)化園地的面積最大子流域(第26號子流域，記為S26)作為典型子流域，進一步分析研究區(qū)典型日徑流變化。

空間響應(yīng)部分的分析，首先計算出37個子流域各地類面積凈變化量與年徑流深變化量，再運用ArcGIS 9.3地統(tǒng)計分析模塊，采用普通克里格方法，對研究區(qū)兩期土地利用下的37個子流域徑流系數(shù)進行插值，得出全流域的徑流系數(shù)等值線，結(jié)果須通過該方法的精度評價檢驗。

3結(jié)果與分析

3.1　土地利用/覆被變化

研究區(qū)林地面積最大，占總面積的65%以上；園地和水田為主要地類，分別占8%以上，其余地類面積較小(表1)。由表1可知，1995—2006年，各地類變化程度表現(xiàn)為林地>園地>建設(shè)用地>旱地>水田>未利用地>水域>草地，其中林地減少面積占總減少面積的92.31%，而園地增加面積則占總增加面積的73.43%。此外，建設(shè)用地的增加是另外一個重要特征，其增加面積占總增加面積的13.83%。

表1　研究區(qū)不同土地利用類型面積及其變化

為獲取研究區(qū)土地利用變化除數(shù)量以外的隱含信息，通過轉(zhuǎn)移矩陣和累計百分比分析研究區(qū)各地類空間轉(zhuǎn)移情況，發(fā)現(xiàn)研究區(qū)最顯著的轉(zhuǎn)變?yōu)榱值叵驁@地的轉(zhuǎn)化，其次是耕地、園地和林地之間的相互轉(zhuǎn)化，而第3層次地類相互轉(zhuǎn)化的變化率在3%～1%之間，變化率合計達15.055%，共有7種地類之間的相互轉(zhuǎn)化；第4層次的地類相互轉(zhuǎn)化有27種之多，但轉(zhuǎn)化量很小，變化面積總和不超過整個變化面積的5%。

3.2　徑流模擬

SWAT模型對山美水庫集水區(qū)流域日尺度的徑流過程模擬效果良好，其中在率定期內(nèi)Ens和R2分別為0.85和0.86；驗證期內(nèi)Ens和R2分別為0.86和0.87。從整個徑流模擬期內(nèi)的流量過程線來看，模擬效率評價結(jié)果也較好，Ens和R2均達到了0.85，模擬與實測流量過程線吻合良好，區(qū)內(nèi)徑流對于降水的響應(yīng)強烈(圖1)。分析表明，基于SWAT構(gòu)建的日尺度徑流模型在研究區(qū)的適應(yīng)性較好，為進一步分析土地利用/覆被變化水文響應(yīng)奠定了堅實基礎(chǔ)。

3.3　土地利用/覆被變化的水文時空響應(yīng)

3.3.1時間響應(yīng)分析研究區(qū)土地利用變化的水文響應(yīng)從時間尺度(表2)來看，年徑流深增加，日徑流變化中，最大1 d和連續(xù)5 d最大徑流均明顯增加，增幅達3%以上，即洪水期的日徑流增加；同時值得注意的是最小1d枯水徑流出現(xiàn)了最大幅度的增加，變化率達到8.18%，只有連續(xù)7 d最小徑流變化量表現(xiàn)為減小，但變化率非常小，僅為0.21%，幾乎可忽略不計。

可見，全流域尺度下，洪水期和枯水期日徑流深均出現(xiàn)了較大幅度的增加。典型子流域的分析表明，S26洪水徑流平均增加5.54%，枯水徑流平均減少12.84%，呈現(xiàn)出“洪水量更大，枯水量更小”的兩極化變化趨勢，與全流域尺度下的日徑流變化存在一定的差異。這種差異是由于全流域各地類之間的相互轉(zhuǎn)化較子流域的變化更復(fù)雜，水文響應(yīng)中正負(fù)效應(yīng)相互抵消情況嚴(yán)重。另外，降水中心的空間變化對徑流的影響[11]也是不容忽視的因素。

圖1　研究區(qū)日徑流實測值與模擬值

研究區(qū)域項目年徑流深洪水期(4—10月)最大1d連續(xù)5d最大枯水期(11月至翌年3月)最小1d連續(xù)7d最小全流域1995年1436.8150.57295.980.739.542006年1470.0852.23301.690.799.52變化量33.271.665.710.06-0.02變化率/%2.323.294.228.18-0.21典型子流域 (S26)1995年1104.9484.2035.180.060.072006年1133.9989.7136.770.050.06變化量29.055.511.60-0.01-0.01變化率/%2.636.554.53-12.22-13.45

當(dāng)研究區(qū)地表植被郁閉度降低，地表覆被和冠層最大截流量減少，地表植被最大截流量的減少，使得降水轉(zhuǎn)為徑流的量增加，年徑流深增加，日徑流深在洪水期增加，枯水期基本減少。

事實上，在其他流域的研究[12-14]得出相似結(jié)論，建設(shè)用地增加導(dǎo)致不透水面積增加，造成入滲減少，洪峰流量增大，徑流系數(shù)增大。葉寶瑩[15]以泥鰍河和科洛河流域為例的分析得出，林地、耕地面積的減少導(dǎo)致地表徑流增加。有關(guān)梭磨河[16]的研究也表明地表覆被和冠層最大截流量的增加，導(dǎo)致徑流深減少，洪水期和枯水期徑流深均減少，且枯水期減少幅度小。

3.3.2空間響應(yīng)分析選用普通克里格[17]插值，并通過該方法的4個參數(shù)對插值結(jié)果進行精度評價。評價結(jié)果(表3)顯示，插值結(jié)果符合研究區(qū)徑流系數(shù)的克里格最優(yōu)參數(shù)模型標(biāo)準(zhǔn)。且研究結(jié)果與林文嬌[18]在山美庫區(qū)以上集水區(qū)的插值數(shù)值上相近；等值線布局的不同，與研究時段和插值方法不同有關(guān)。研究區(qū)空間上，土地利用變化下的水文響應(yīng)結(jié)果如圖2所示。

根據(jù)圖2可知，兩期土地利用下徑流系數(shù)的差異體現(xiàn)為2006年土地利用覆被下徑流系數(shù)大于1995年，表明2006年土地利用/覆被變化后，山美水庫集水區(qū)降水量中有更多水量變成徑流，而其余耗損于植物截留、填洼、入滲和蒸發(fā)的水量相應(yīng)減少，這恰好驗證了研究區(qū)林地減少、建設(shè)用地增加的現(xiàn)狀。此外，從等值線分布來看，2006年徑流系數(shù)等值線比1995年等值線更密集，特別是西部地區(qū)和東南部地區(qū)，這說明研究區(qū)土地利用/覆被變化后，徑流系數(shù)變率增大、穩(wěn)定性下降。

表3　普通克里格插值結(jié)果精度評價

圖2　研究區(qū)徑流系數(shù)等值線圖

3.3.3徑流變化影響因子分析徑流的變化從來都不是單一地類影響所致。因各地類變化并不相互獨立，單個變量包含信息較片面，而相關(guān)分析只能區(qū)分出各地類對研究區(qū)徑流變化的正負(fù)效應(yīng)，不能確定其影響程度。所以，為進一步識別研究區(qū)徑流量變化的主要影響因子，對各地類變化做主成分分析，結(jié)合主成分載荷矩陣，確定各地類變化對區(qū)內(nèi)徑流變化的影響程度，進而識別出主要影響因子。

主成分分析結(jié)果表明，4個主成分可以解釋81.97%的地類變化信息，第1到第4主成分的貢獻率依次為31.02%，21.20%，16.89和12.86%。由載荷矩陣(表4)可知，園地、林地、建設(shè)用地和水田是影響區(qū)內(nèi)水文變化的主要因素，草地雖然與主成分具有較高的相關(guān)系數(shù)，但其總面積很小，因此影響也相對較小。

表4　主成分載荷矩陣

4結(jié) 論

(1) LUCC分析得出，在全流域園地、建設(shè)用地增加，林地減少，并以林地轉(zhuǎn)為園地為主，即研究區(qū)的地表覆被和冠層最大截留量減少。這樣的變化導(dǎo)致時間水文響應(yīng)上，年徑流深無論全流域或子流域均增加；全流域在日尺度的洪水期和枯水期日徑流深一致增加；對典型子流域(S26)的研究發(fā)現(xiàn)，洪水期日徑流深增加，而枯水期研究區(qū)土壤下滲量增加，導(dǎo)致徑流深減少，枯水徑流減少。體現(xiàn)子流域洪水與枯水徑流深變化的不一致性，呈現(xiàn)出“洪水量更大，枯水量更小”的特點。

(2) 空間水文響應(yīng)分析得出，同一地點2006年土地利用下的徑流系數(shù)大于1995年土地利用下的徑流系數(shù)，且前者較后者等值線更為密集。這與研究區(qū)研究期內(nèi)林地減少，建設(shè)用地的增加密不可分，表明降水量中耗損的水量減少，形成的徑流更多；研究區(qū)土地利用/覆被變化后，徑流系數(shù)變率大，不似從前穩(wěn)定。

(3) 研究期內(nèi)，研究區(qū)徑流時空變化的主要影響因子為園地、林地、建設(shè)用地和水田的變化。土地利用/覆被變化與流域徑流量的關(guān)系，絕非某一種土地利用類型的變化能闡述清楚，而是由各種地類變化的組合疊加效應(yīng)作用于河川徑流。這種疊加效應(yīng)還須進一步研究并加以揭示。

[參考文獻]

[1]Neitsch S L, Arnold J G, Kiniry J R, et al. Soil and water assessment tool theoretical documentation version 2009[R]. Texas： Texas Water Resources Institute, 2011.

[2]李慧,雷曉云,包安明,等.基于SWAT模型的山區(qū)日徑流模擬在瑪納斯河流域的應(yīng)用[J].干旱區(qū)研究,2010,27(5):686-690.

[3]李麗嬌,薛麗娟,張奇.基于SWAT的西苕溪流域降雨—徑流關(guān)系及水量平衡分析[J].水土保持通報,2008,28(5):81-85.

[4]黃清華,張萬昌.SWAT模型參數(shù)敏感性分析及應(yīng)用[J].干旱區(qū)地理,2010(1):8-15.

[5]李占玲,徐宗學(xué).黑河流域上游山區(qū)徑流模擬及模型評估[J].北京師范大學(xué)學(xué)報:自然科學(xué)版,2010(3):344-349.

[6]楊柳,劉梅冰,陳瑩,等.山美水庫集水區(qū)徑流模擬的日尺度SWAT模型[J].南水北調(diào)與水利科技,2013,11(S1):3-5.

[7]許有鵬,葛小平,張立峰,等.東南沿海中小流域平原區(qū)洪水淹沒模擬[J].地理研究,2005,24(1):38-45.

[8]福建師范大學(xué)地理系《福建自然地理》編寫組.福建自然地理[M].福州:福建人民出版社,1987:64-131.

[9]賴格英,吳敦銀,鐘業(yè)喜,等.SWAT模型的開發(fā)與應(yīng)用進展[J].河海大學(xué)學(xué)報:自然科學(xué)版,2012,40(3):243-251.

[10]Nash J E, Sutcliffe J V. River flow forecasting through conceptual models(part I)：A discussion of principles[J]. Journal of hydrology, 1970, 10(3): 282-290.

[11]林木生,陳興偉,陳瑩.晉江西溪流域洪水與暴雨時空分布特征的相關(guān)分析[J].資源科學(xué),2011,33(12):2226-2231.

[12]Bhaduri B L. A geographic information system-based model of the long-term impact of land use change on nonpoint-source pollution at a watershed scale[D]. West Lafayette： Purdue University, 1998.

[13]陳瑩,許有鵬,陳興偉.長江三角洲地區(qū)中小流域未來城鎮(zhèn)化的水文效應(yīng)[J].資源科學(xué),2011(1):64-69.

[14]史曉亮,李穎,嚴(yán)登華,等.流域土地利用/覆被變化對水文過程的影響研究進展[J].水土保持研究,2013,20(4):301-308.

[15]葉寶瑩,張養(yǎng)貞,張樹文,等.嫩江流域土地覆被變化對徑流量的影響分析[J].水土保持通報,2003,23(2):15-18.

[16]鄧慧平,李秀彬.流域土地覆被變化水文效應(yīng)的模擬：以長江上游源頭區(qū)梭磨河為例[J].地理學(xué)報,2003,58(1):53-62.

[17]湯國安.楊昕.ArcGIS地理信息系統(tǒng)空間分析實驗教程[M].北京：科學(xué)出版社,2006:401-408.

[18]林文嬌.晉江東溪流域農(nóng)業(yè)非點源污染模擬分析[D].福建 福州：福建師范大學(xué),2009.

Hydrological Responses to Land Use and Cover Change in a Small Watershed of Southeast China Coast Area

YANG Liu1,2, CHEN Xingwei2, XU Youpeng2, HU Chunsheng1

(1.GeographicandOceanographicSciences,NanjingUniversity,Nanjing,Jiangsu210093,China; 2.SchoolofGeographicSciences,FujianNormalUniversity,Fuzhou,Fujian350007,China)

Abstract：[Objective] Studying the hydrological responses to land use and cover change of Southeast China coast in a small watershed, in order to provide references for the rational allocation and sustainable development of water and soil resources. [Methods] The confluence catchment of Jinjiang Shanmei Reservoir in Fujian Province of China, a place located at the southeast coast of China, was selected to illustrate the effects of land use and cover change(LUCC) on hydrological processes by SWAT model. [Results] (1) From 1995 to 2006, the forest land was mainly converted to garden land and construction land, while farm land was mainly transformed to construction land; (2) In the model, if other conditions were fixed, only data of land use in 2006 was replaced with the corresponding input variable in 1995, the runoff could be outputted as contrast: from 2000 to 2010, the annual runoffs in the whole watershed and also in the tributary drainage area consequently increased; and the averaged runoff during flooding increased, while the runoff during rainless period decreased; (3) The runoff volume in the whole watershed increased with an increased variation as contrast with the previous corresponding value; (4) Garden land, forest land, construction land and paddy field were the main driving factors that caused the regional hydrological changes. [Conclusion] Land use and cover change in the study area were found to result in a runoff trend that are “the flood volume is larger, the runoff yield is smaller”. The changes of river runoff was resulted from the combined effects of all types of land use change. This combination effects need to be further studied.

Keywords:SWAT model; daily runoff process simulation; land use and cover change; hydrological responses; Southeast China coastal area; small watershed

文獻標(biāo)識碼：A

文章編號：1000-288X(2015)02-0070-06

中圖分類號：

收稿日期：2014-01-21修回日期：2014-04-12