李鴻儒，盧　遠(yuǎn)，何　文

(廣西師范學(xué)院 地理科學(xué)與規(guī)劃學(xué)院 北部灣環(huán)境演變與資源利用省部共建教育部重點(diǎn)實(shí)驗(yàn)室，廣西 南寧 530001)

?

欽江流域土地利用/覆被變化對(duì)徑流產(chǎn)沙的影響

李鴻儒，盧遠(yuǎn)，何文

(廣西師范學(xué)院 地理科學(xué)與規(guī)劃學(xué)院 北部灣環(huán)境演變與資源利用省部共建教育部重點(diǎn)實(shí)驗(yàn)室，廣西 南寧 530001)

SWAT；土地利用類型；水文效應(yīng)；欽江流域

以欽江流域?yàn)檠芯繉?duì)象，在GIS軟件的支持下，構(gòu)建研究區(qū)SWAT分布式水文模型，以徑流和泥沙實(shí)測(cè)數(shù)據(jù)對(duì)模型進(jìn)行校準(zhǔn)和驗(yàn)證，結(jié)果顯示相關(guān)系數(shù)R2高于0.9，Nash-Sutcliffe系數(shù)高于0.8，模型模擬精度較高。在固定天氣因子等不變的情況下，基于1990、2008年兩期土地利用數(shù)據(jù)，模擬研究區(qū)1968—2008年徑流產(chǎn)沙過(guò)程，并分析土地利用/覆被變化對(duì)水文效應(yīng)的影響，結(jié)果表明：1990—2008年土地利用變化主要表現(xiàn)為林地、水田面積減少，旱地、草地、水域和建設(shè)用地面積增加；與1990年土地利用情景相比，2008年土地利用情景下徑流和輸沙總量均增加，土地利用變化對(duì)產(chǎn)流影響較小，對(duì)產(chǎn)沙影響較顯著；由土地利用變化引起的月均徑流產(chǎn)沙變化量為汛期增幅大于非汛期；土地利用變化對(duì)枯水年徑流產(chǎn)沙過(guò)程影響最大，平水年次之，豐水年影響最小。

氣候變化和人類活動(dòng)是影響流域水文過(guò)程的兩個(gè)重要因素，了解自然變化和人類活動(dòng)對(duì)水文行為乃至周邊環(huán)境的影響是人類社會(huì)發(fā)展面臨的共同問(wèn)題，其中土地利用/覆被的變化直接反映了人類活動(dòng)的影響程度[1]。隨著經(jīng)濟(jì)社會(huì)的快速發(fā)展，受城市化、人口增長(zhǎng)、政府政策等影響，人類對(duì)自然資源的利用和改造必然會(huì)促使土地利用結(jié)構(gòu)發(fā)生改變，不當(dāng)?shù)耐恋乩梅绞綄?dǎo)致流域水量失衡，從而引發(fā)水土流失、土地退化等生態(tài)環(huán)境問(wèn)題，因此開(kāi)展土地利用/覆被變化的水文響應(yīng)研究對(duì)流域水資源合理開(kāi)發(fā)與管理、保證水資源可持續(xù)利用具有重要意義[2]。SWAT模型是美國(guó)農(nóng)業(yè)部農(nóng)業(yè)研究中心1993年為美國(guó)水文模型項(xiàng)目開(kāi)發(fā)的大、中尺度流域環(huán)境模擬模型，能夠利用GIS和RS提供的空間數(shù)據(jù)信息模擬產(chǎn)水產(chǎn)沙、非點(diǎn)源污染等過(guò)程。我國(guó)一些學(xué)者自2003年開(kāi)始運(yùn)用該模型對(duì)不同尺度流域展開(kāi)了大量研究，取得許多有意義的成果[3-12]?？紤]到廣西境內(nèi)流域利用分布式水文模型的相關(guān)研究成果較少，本研究以欽江流域?yàn)檠芯繉?duì)象，構(gòu)建流域分布式水文模型，在固定天氣因子、DEM、土壤類型等因素不變的情況下，模擬真實(shí)土地利用情景下的徑流產(chǎn)沙過(guò)程，并分析不同土地利用狀況對(duì)水文效應(yīng)的影響，以期為欽江流域水土資源與生態(tài)環(huán)境的可持續(xù)性綜合管理提供決策依據(jù)和參考。

1　研究區(qū)概況

欽江流域位于廣西南部，地處E108°34′～109°29′、N21°53′～22°36′，欽江從東北向西南穿過(guò)，流域面積2 241.69 km2，主干流全長(zhǎng)173.54 km。流域呈狹長(zhǎng)狀，地勢(shì)東北高、西南低，地形呈舟狀，除流域出口西南部外，其他三面高、中間低，東北部及四周多為低山丘陵地區(qū)。該流域?qū)賮啛釒Ъ撅L(fēng)氣候區(qū)，春夏季多暴雨，受熱帶氣旋影響，雨量豐沛，多年平均降水量在1 600 mm左右，降水年內(nèi)分配不均，多集中在4—9月，降水量由下游河口向上游河源遞減。

2　水文模型構(gòu)建

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

研究區(qū)DEM數(shù)據(jù)來(lái)源于國(guó)家1 ∶5萬(wàn)地形圖，經(jīng)ArcGIS數(shù)字化后重采樣生成30 m×30 m Grid數(shù)據(jù)；1990和2008年兩期土地利用數(shù)據(jù)通過(guò)TM影像人工解譯得到，同上將數(shù)據(jù)轉(zhuǎn)化成Grid格式，網(wǎng)格大小為30 m×30 m，并對(duì)土地利用類型進(jìn)行重分類，以符合SWAT模型模擬要求；土壤類型圖來(lái)源于中國(guó)1 ∶100萬(wàn)土壤數(shù)據(jù)庫(kù)，坐標(biāo)統(tǒng)一為Albers Krasovsky_1940。

土壤屬性數(shù)據(jù)主要來(lái)源于《廣西土種志》，部分通過(guò)SPAW軟件計(jì)算得到。氣象數(shù)據(jù)采用陸屋、小江、坡郎平、鎮(zhèn)龍、邕寧5個(gè)站點(diǎn)1961—2008年長(zhǎng)序列數(shù)據(jù)，包括日降雨量、相對(duì)濕度、最高最低氣溫、平均風(fēng)速和太陽(yáng)輻射，其中太陽(yáng)輻射由日照時(shí)數(shù)計(jì)算得到。水文數(shù)據(jù)來(lái)源于《珠江流域水文資料》，由于獲取和采集的難度較大，因此只收集到陸屋站水文數(shù)據(jù)，時(shí)間范圍為2003—2005、2001和2008 年。

2.2模型構(gòu)建

采用ArcGIS 10.1環(huán)境下的SWAT-2012，首先將研究區(qū)DEM數(shù)據(jù)輸入ArcSWAT以驅(qū)動(dòng)模型運(yùn)行，設(shè)置集水面積閾值為3 500 hm2，輸入水文站點(diǎn)和流域總出口，自動(dòng)將研究區(qū)劃分為44個(gè)子流域；其次疊加不同精度的土地利用、土壤類型和地形圖再將子流域劃分為355個(gè)水文響應(yīng)單元；最后將氣象數(shù)據(jù)輸入天氣發(fā)生器，寫入數(shù)據(jù)則完成模型的初步構(gòu)建。

2.3參數(shù)率定與驗(yàn)證

參數(shù)率定的目的在于盡可能地縮小觀測(cè)值與模擬值之間的差距，以提高模型的精準(zhǔn)度。利用SWAT-CUP軟件提供的SUFI2算法對(duì)影響水文過(guò)程的參數(shù)進(jìn)行敏感性分析，其中CN2(徑流曲線系數(shù))、SOL_K(土壤飽和水力傳導(dǎo)系數(shù))、ESCO(植物吸收補(bǔ)償因子)和SOL_AWC(土壤可利用有效水)對(duì)產(chǎn)流過(guò)程影響較大，USLE_C(植被覆蓋因子)、USLE_P(水土保持措施因子)和HRU_SLP(平均坡度)等對(duì)輸沙過(guò)程影響較大。在此基礎(chǔ)上，對(duì)參數(shù)進(jìn)行有針對(duì)性的微調(diào)，使用相關(guān)系數(shù)R2和Nash-Sutcliffe效率系數(shù)NS對(duì)模型參數(shù)率定的合理性進(jìn)行評(píng)價(jià)，計(jì)算公式為

(1)

(2)

式中：Q1i為模擬值；Q2i為實(shí)測(cè)值；Qavg為模擬均值；Pavg為實(shí)測(cè)均值；n為實(shí)測(cè)數(shù)據(jù)數(shù)量。

R2取值介于0～1，越接近于1，表明擬合效果越好；NS值越接近于1，表明模型的精準(zhǔn)度越高?；谑占降乃馁Y料，以2003—2005年為率定期，2001年和2008年為驗(yàn)證期，將校準(zhǔn)期調(diào)整好的參數(shù)代回SWAT模型進(jìn)行驗(yàn)證，評(píng)價(jià)結(jié)果見(jiàn)表1和圖1、2。由圖、表可知，陸屋站徑流、泥沙率定期和驗(yàn)證期的NS系數(shù)都在0.80以上，且R2都大于0.90， 產(chǎn)流產(chǎn)沙變化趨勢(shì)一致，說(shuō)明模擬精度較高，水文分布式模型適用于欽江流域，可對(duì)流域缺失的水文資料進(jìn)行模擬。

表1　陸屋站產(chǎn)流產(chǎn)沙量月尺度模擬評(píng)價(jià)

圖1　陸屋站逐月日均流量模擬與實(shí)測(cè)情況

圖2　陸屋站逐月日均輸沙量模擬與實(shí)測(cè)情況

3　結(jié)果分析

3.1土地利用類型及其動(dòng)態(tài)變化

欽江流域主要土地利用類型有水田、旱地、林地、草地、水域、建設(shè)用地等6種，以水田、旱地、林地和草地為主，1990、2008年合計(jì)分別占到流域總面積的92.62%、92.40%，尤其林地占比例最高，說(shuō)明流域植被覆蓋狀況較好。利用ArcGIS和Excel軟件，定量分析各土地利用類型在1990—2008年期間的變化情況，結(jié)果見(jiàn)表2。由表2可知，林地、水田面積減少，旱地、草地、水域、建設(shè)用地面積增加；水田、旱地、林地、草地不同程度地轉(zhuǎn)變?yōu)榻ㄔO(shè)用地，并出現(xiàn)水域干涸轉(zhuǎn)變?yōu)樗?、旱地和林地的情況。將退耕還林、植樹(shù)造林、減少建設(shè)用地等視為正向轉(zhuǎn)移，反之為負(fù)向轉(zhuǎn)移，通過(guò)統(tǒng)計(jì)可知?dú)J江流域28年間正向轉(zhuǎn)移面積為13.66 km2，負(fù)向轉(zhuǎn)移面積為38.81 km2，總體來(lái)看流域1990—2008年植被覆蓋度呈下降趨勢(shì)。

表2　欽江流域1990—2008年土地利用轉(zhuǎn)移矩陣　km2

3.2不同土地覆被對(duì)徑流產(chǎn)沙的影響

將2008年土地利用數(shù)據(jù)代入已校準(zhǔn)好的SWAT模型，固定氣候因子等不變，對(duì)欽江流域1968—2008年產(chǎn)流產(chǎn)沙量進(jìn)行模擬。與1990年土地利用情景相比，2008年徑流和輸沙量均增加，其中：徑流增加量較少，為1.19萬(wàn)m3，變化率為0.17%；輸沙量增加559.89萬(wàn)t，變化率高達(dá)4.01%。通過(guò)對(duì)比分析不同土地利用類型及真實(shí)情景下的產(chǎn)流產(chǎn)沙總量，發(fā)現(xiàn)欽江流域1990—2008年水文效應(yīng)改變與土地利用類型變化有關(guān)，植被覆蓋度下降引起流域產(chǎn)流產(chǎn)沙量增加，說(shuō)明土地利用結(jié)構(gòu)改變是影響欽江流域水文效應(yīng)的重要原因。

3.3不同土地覆被對(duì)汛期和非汛期徑流產(chǎn)沙的影響

為分析研究區(qū)土地利用變化對(duì)汛期(4—9月)和非汛期產(chǎn)流產(chǎn)沙過(guò)程的影響，計(jì)算了1968—2008年各月產(chǎn)流產(chǎn)沙平均增量(表3)。由表3可知，與1990年相比，2008年土地利用情景下流域汛期和非汛期各月平均產(chǎn)流產(chǎn)沙變化量趨勢(shì)基本一致，即均有所增加，且汛期增幅大于非汛期。汛期產(chǎn)流產(chǎn)沙量平均增量分別為243.91 m3、10.04萬(wàn)t，非汛期分別為46.29 m3、3.62萬(wàn)t，月平均輸沙增量較明顯，而徑流增量甚微。由此可見(jiàn)，土地利用覆被條件的變化是引起流域汛期產(chǎn)流產(chǎn)沙量增加的重要原因。由于植被具有巨大的保持水土功能，特別是枯枝落葉層阻礙消耗了降雨濺蝕的動(dòng)能，減少了面層徑流的擾動(dòng)和挾沙能力，因此恢復(fù)流域植被能起到滯洪和減少洪峰流量的作用，從而有助于減少汛期洪水、泥石流等災(zāi)害的發(fā)生。

表3　1990、2008年土地利用情景下各月產(chǎn)流產(chǎn)沙平均增量

3.4不同土地覆被對(duì)不同降水年型徑流產(chǎn)沙的影響

依據(jù)欽江流域陸屋站1961—2008年的降雨資料，用保證率進(jìn)行統(tǒng)計(jì)分析確定不同降水年型，以分析不同土地利用情景下各典型水文年份對(duì)流域產(chǎn)流產(chǎn)沙的影響(表4)。豐水年(降雨頻率10%)、平水年(降雨頻率50%)、枯水年(降雨頻率90%)分別選取與之相接近的1970、1978和2005年作為典型水文年份，降雨量分別為2 162.9、1 652.0、1 101.1 mm，平水年降雨量與48年的平均降雨量1 632.72 mm相近。從徑流量來(lái)看，2005年變化率最大，增幅為0.19%，1978和1970年次之，增幅分別為0.18%和0.12%。從輸沙量來(lái)看，2005年變化率最大，1978和1970年次之，趨勢(shì)與徑流量一致，增幅分別為6.01%、4.20%和3.34%。由此推斷，土地覆被變化對(duì)枯水年影響最大，平水年次之，豐水年影響最小。

表4　1990、2008年土地利用情景下不同降水年型產(chǎn)流產(chǎn)沙量

4　結(jié)　論

本研究基于SWAT分布式水文模型，利用DEM、土地利用和土壤屬性數(shù)據(jù)、氣象水文數(shù)據(jù)，構(gòu)建欽江流域水文分布式模型，在固定天氣因子等因素不變的情況下，對(duì)1990、2008年兩期不同土地利用情景下的產(chǎn)流產(chǎn)沙過(guò)程進(jìn)行模擬，結(jié)論如下：

(1)基于欽江流域陸屋站的實(shí)測(cè)徑流輸沙數(shù)據(jù)，采用R2、NS系數(shù)對(duì)模型的可用性進(jìn)行評(píng)價(jià)，陸屋站徑流、泥沙率定期和驗(yàn)證期的NS系數(shù)都在0.8以上，R2都大于0.9，模擬精度較好，徑流和輸沙變化趨勢(shì)趨于一致，表明模型校準(zhǔn)較成功，模擬精度滿足研究的需要，可對(duì)流域缺失的徑流泥沙數(shù)據(jù)進(jìn)行模擬和分析。固定天氣因子等不變的情況下，在校準(zhǔn)好的模型中輸入2008年土地利用數(shù)據(jù)，模擬流域1968—2008年徑流、輸沙量。

(2)對(duì)1990、2008年兩期土地利用數(shù)據(jù)進(jìn)行對(duì)比分析可知，主要地類水田、旱地、林地和草地占到流域總面積的92%以上，其中林地分別占到40.83%和39.84%，流域植被覆蓋度整體較高。相對(duì)于1990年，2008年土地利用變化表現(xiàn)為林地、水田面積減少，旱地、草地、水域、建設(shè)用地面積增加，流域植被覆蓋度總體呈下降趨勢(shì)。

(3)根據(jù)1968—2008年的模擬結(jié)果，與1990年相比，2008年土地利用情景下徑流和輸沙總量均增加，其中：徑流增加量較少，為1.19萬(wàn)m3，變化率為0.17%；輸沙量增加559.89萬(wàn)t，變化率高達(dá)4.01%。2008年土地利用情景下的流域汛期和非汛期各月平均產(chǎn)流產(chǎn)沙變化量趨勢(shì)基本一致，即均有所增加，且汛期增幅大于非汛期，土地覆被變化對(duì)枯水年影響最大，平水年次之，豐水年影響最小。

(4)本研究在固定天氣因子等不變的情況下，定量分析了欽江流域土地利用/覆被變化對(duì)徑流產(chǎn)沙的影響，未能考慮天氣變化、人類活動(dòng)等因素影響，具有一定的局限性。今后進(jìn)一步分析流域的水文效應(yīng)，還需綜合考慮天氣、河道外取用水、水土保持措施等因素的影響。

[1] 歐春平,夏軍,王中根,等.土地利用/覆被變化對(duì)SWAT模型水循環(huán)模擬結(jié)果的影響研究——以海河流域?yàn)槔齕J].水力發(fā)電學(xué)報(bào),2009,28(4):124-129.

[2] 王學(xué),張祖陸,寧吉才.基于SWAT模型的白馬河流域土地利用變化的徑流響應(yīng)[J].生態(tài)學(xué)雜志,2013,32(1):186-194.

[3] 王中根,劉昌明,黃友波.SWAT模型的原理、結(jié)構(gòu)及應(yīng)用研究[J].地理科學(xué)進(jìn)展,2003,22 (1):79-86.

[4] 郝芳華,陳利群,劉昌明,等.土地利用變化對(duì)產(chǎn)流和產(chǎn)沙的影響分析[J].水土保持學(xué)報(bào),2004,18(3):6-8.

[5] 吳軍,張萬(wàn)昌. SWAT徑流模擬及其對(duì)流域內(nèi)地形參數(shù)變化的響應(yīng)研究[J].水土保持通報(bào),2007,27(3):52-58.

[6] 龐靖鵬,劉昌明,徐宗學(xué).密云水庫(kù)流域土地利用變化對(duì)產(chǎn)流和產(chǎn)沙的影響[J].北京師范大學(xué)學(xué)報(bào):自然科學(xué)版,2010,46(3):290-299.

[7] 段超宇,張生,李錦榮,等.基于SWAT模型的內(nèi)蒙古錫林河流域降水—徑流特征及不同水文年徑流模擬研究[J].水土保持研究,2014,15(10):292-297.

[8] 袁媛,雷曉輝,蔣云鐘,等. 基于SWAT模型的西江流域徑流模擬研究[J].中國(guó)農(nóng)村水利水電，2015,15(3):14-21.

[9] 依熱下提·卡米力,李建柱,馮平.潘家口水庫(kù)流域土地利用變化的水文響應(yīng)分析[J].水資源與水工程學(xué)報(bào),2015,26(2):47-51.

[10] 王林,陳興偉.退化山地生態(tài)系統(tǒng)植被恢復(fù)水文效應(yīng)的SWAT模擬[J].山地學(xué)報(bào),2008,26(1):71-75.

[11] 曾赟,魏琳.土地利用和氣候變化對(duì)四川省紫色丘陵區(qū)徑流泥沙變化的影響研究[J].水土保持通報(bào),2013,33(3):1-6.

[12] 高楊,陳興偉,張紅梅.基于子流域土地覆被變化的產(chǎn)流產(chǎn)沙效應(yīng)模擬[J].水土保持通報(bào),2012,32(4):75-79.

(責(zé)任編輯李楊楊)

2015-07-15

廣西自然科學(xué)基金項(xiàng)目(2015GXNSFAA139234)

P332；S157

A

1000-0941(2016)04-0040-04

李鴻儒(1991—)，女(壯族)，廣西河池市人，碩士研究生，研究方向?yàn)榱饔蛩淖兓?；通信作者盧遠(yuǎn)(1971—)，男，廣西橫縣人，副研究員，博士，從事GIS應(yīng)用、環(huán)境遙感研究。