劉樹鋒 陳記臣 關(guān)帥 王會(huì)

摘要以楊溪河流域?yàn)槔?，利用未來土地利用變化情景模擬（GeoSOS-FLUS）模型對(duì)未來土地利用情況進(jìn)行模擬，并通過InVEST模型對(duì)水源涵養(yǎng)功能進(jìn)行分析。結(jié)果表明：整體上楊溪河流域內(nèi)未來各類土地利用類型更趨于集中，耕地面積大量減少并趨于消失，林地面積大量增加并普遍大于其他土地類型;楊溪河流域在現(xiàn)狀和未來各情景下均表現(xiàn)為上游和下游的產(chǎn)水量較高;人類活動(dòng)會(huì)對(duì)楊溪河流域產(chǎn)水量產(chǎn)生影響，注重全球化的變化和環(huán)境保護(hù)可保持流域良好的水源涵養(yǎng)功能。水源涵養(yǎng)功能在流域的開發(fā)利用中有重要作用，需要保護(hù)林地并協(xié)調(diào)經(jīng)濟(jì)、環(huán)境和資源等各方面之間的關(guān)系，以維持楊溪河流域等我國南方中小流域較強(qiáng)的水源涵養(yǎng)功能。

關(guān)鍵詞InVEST模型;未來土地利用變化;水源涵養(yǎng);楊溪河流域

中圖分類號(hào)X171.1文獻(xiàn)標(biāo)識(shí)碼A文章編號(hào)0517-6611（2020）15-0067-04

doi：10.3969/j.issn.0517-6611.2020.15.020

開放科學(xué)（資源服務(wù)）標(biāo)識(shí)碼（OSID）：

Effect of Future Land Use Change on Water Conservation Function Based on InVEST Model—Taking Yangxi River Basin as an Example

LIU Shufeng1，2，3， CHEN Jichen1，2，3， GUAN Shuai1，2，3et al

（1.Guangdong Research Institute of Water Resources and Hydropower， Guangzhou， Guangdong 510635;2.State and Local Joint Engineering Laboratory of Estuarine Hydraulic Technology， Guangzhou， Guangdong 510635;3.Guangdong Key Laboratory of Hydrodynamic Research， Guangzhou， Guangdong 510635）

AbstractTaking Yangxi River Basin as an example， this paper used GeoSOSFLUS model for simulating land use change in the future， and used InVEST model to analyse the water conservation function. The results showed that： As a whole， the land use types in Yangxi River Basin would be more concentrated in the future， and the area of farmland would be greatly reduced and tend to disappear. The area of forest land would increase greatly and would be generally larger than other types of land. Yangxi River Basin had higher water yield in the upstream and downstream areas in current and all future scenarios. Human activities would affect the water yield of Yangxi River Basin. Paying attention to the change of globalization and environmental protection could maintain a good water conservation function of the basin. Water conservation function would play an important role in the development and utilization of the basin. It is necessary to protect forest land and coordinate the relationship among economy， environment and resources to maintain the strong water conservation function of Yangxi River Basin and other small and mediumsized river basins in the south of China.

Key wordsInVEST model;Future land use change;Water conservation;Yangxi River Basin

基金項(xiàng)目廣東省水利科技創(chuàng)新項(xiàng)目（2017-12）。

作者簡介劉樹鋒（1978—），男，廣東韶關(guān)人，高級(jí)工程師，博士，從事水文水資源研究。

收稿日期2019-12-26

在全球水資源需求不斷增加的情況下，對(duì)不同土地利用條件下的產(chǎn)水量評(píng)估和水源涵養(yǎng)功能評(píng)價(jià)研究具有巨大意義。目前常用的產(chǎn)水量評(píng)估模型有MIKE模型[1]和TOPMODEL模型[2]等，但這些模型所需數(shù)據(jù)量較大，難以適用于中小流域。InVEST（integrated valuation of ecosystem services and tradeoffs）模型[3]被引入我國后，因其具有很強(qiáng)的空間分析能力，在產(chǎn)水量評(píng)估和水源涵養(yǎng)功能研究等方面取得了良好的模擬效果。已有不少學(xué)者在我國黃土高原[4]、橫斷山區(qū)[5]和漢江上流[6]等地區(qū)運(yùn)用InVEST模型對(duì)當(dāng)?shù)禺a(chǎn)水量或水源涵養(yǎng)的空間格局等進(jìn)行了分析與評(píng)價(jià)。該研究利用InVEST模型的水源涵養(yǎng)模塊（Water Yield），對(duì)韶關(guān)市乳源瑤族自治縣楊溪河流域現(xiàn)狀和未來不同土地利用情景進(jìn)行水源涵養(yǎng)功能分析，以期為我國南方地區(qū)中小流域的開發(fā)利用提供參考依據(jù)。

1區(qū)域概況

楊溪河屬于珠江流域北江水系武江的一級(jí)支流，總集雨面積498 km2，河長61.29 km，坡降1.19%，總落差1 637 m。該流域發(fā)源于廣東省韶關(guān)市乳源瑤族自治縣的老鵬頂山（海拔高程1 737 m，位于南嶺國家級(jí)自然保護(hù)區(qū)），自西向東流經(jīng)大橋鎮(zhèn)、必背鎮(zhèn)和桂頭鎮(zhèn)，在楊溪口匯入武江。楊溪河流域?qū)儆谥衼啛釒Ъ撅L(fēng)氣候，雨量豐富，流域多年平均降水量1 700 mm，上游降水量比中下游大;流域多年平均蒸發(fā)量1 089 mm。

2研究方法

2.1InVEST模型算法

InVEST模型Water Yield模塊的基礎(chǔ)是Budyko水熱耦合平衡，并考慮年均降水量和實(shí)際蒸散發(fā)，運(yùn)用水量平衡法結(jié)合氣候、地形、根系深度和土壤深度等計(jì)算得出各柵格的水源涵養(yǎng)量[7-8]。模型的計(jì)算原理如公式（1）～（4）：

式中，Y?xj為第j種土地利用類型柵格x的年產(chǎn)水量;AET?xj為第j種土地利用類型柵格x的年實(shí)際蒸散發(fā)量;P?x為柵格x的年降雨量;ω?x為自然氣候-土壤性質(zhì)的非物理參數(shù);R?xj為Budyko干燥指數(shù)，即潛在蒸散與降水量的比值;Z為季節(jié)系數(shù)，是與降雨季節(jié)分布有關(guān)的浮點(diǎn)型數(shù)據(jù)（1～30）;AWC?x為柵格x的年平均植物可利用含水量，由土壤深度和理化性質(zhì)決定;K?xj為第j種土地利用類型柵格x的植被蒸散發(fā)系數(shù);ET?0為年平均參考作物蒸散發(fā)量。

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

該研究對(duì)各參數(shù)和數(shù)據(jù)作了區(qū)域處理以增強(qiáng)模型適用性。地形數(shù)據(jù)來源于地理空間數(shù)據(jù)云（http：∥www.gscloud.cn），空間分辨率為30 m×30 m;年平均降水量（P?x）和年實(shí)際蒸散量（AET?xj）等氣候數(shù)據(jù)根據(jù)《廣東省水文圖集》獲得多年平均值，并進(jìn)行克里金插值而得;土壤數(shù)據(jù)來源于世界土壤數(shù)據(jù)庫（harmonized world soil database，HWSD）和南京土壤所1∶100萬的土壤資料;2010年的土地利用數(shù)據(jù)來源于中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心（http：∥www.resdc.cn），空間分辨率為30 m×30 m;年平均植物可利用含水量（AWC?x）等模型參數(shù)參考周文佐[9]的研究結(jié)果，使用SPAW Hydrology模型中的Open Soil Water Characteristics進(jìn)行計(jì)算獲得。

2.3未來土地利用變化情景模擬模型

未來土地利用數(shù)據(jù)采用未來土地利用變化情景模擬模型（GeoSOS-FLUS）進(jìn)行模擬而得，該模型是中山大學(xué)LI等[10]在2017年開發(fā)的多類土地利用變化情景模擬軟件。GeoSOS-FLUS模型的原理源自元胞自動(dòng)機(jī)（CA），并在傳統(tǒng)元胞自動(dòng)機(jī)的基礎(chǔ)上作了較大改進(jìn)。首先，GeoSOS-FLUS模型采用神經(jīng)網(wǎng)絡(luò)算法（ANN）從一期土地利用數(shù)據(jù)與包含人類活動(dòng)與自然效應(yīng)的多種驅(qū)動(dòng)力因子中獲取各種用地類型在研究范圍內(nèi)的適宜性概率。其次，GeoSOS-FLUS模型采用從一期土地利用分布數(shù)據(jù)中采樣的方式，能較好地避免誤差傳遞的發(fā)生。另外，在土地變化模擬過程中，GeoSOS-FLUS模型提出一種基于輪盤賭選擇的自適應(yīng)慣性競(jìng)爭(zhēng)機(jī)制，該機(jī)制能有效處理多種土地利用類型在自然作用與人類活動(dòng)共同影響下發(fā)生相互轉(zhuǎn)化時(shí)的不確定性與復(fù)雜性，使得該模型具有較高的模擬精度并能獲得與現(xiàn)實(shí)土地利用分布相似的結(jié)果[11]。

3結(jié)果與分析

3.1未來土地利用變化模擬

LI等[10-11]選用2010年土地利用數(shù)據(jù)對(duì)未來土地利用變化進(jìn)行模擬，均發(fā)現(xiàn)GeoSOS-FLUS模型的模擬結(jié)果具有較高的準(zhǔn)確性。由于未來?xiàng)l件的不確定性，該研究選擇基于聯(lián)合國政府間氣候變化專門委員會(huì)（IPCC）所發(fā)布的4個(gè)SRES情景來進(jìn)行未來土地利用變化的預(yù)測(cè)。在不同的人口與經(jīng)濟(jì)增長及資源使用等情況下對(duì)未來土地利用變化進(jìn)行模擬，分為A?1、A?2、B?1和B?2 4種不同情景，每一個(gè)情景代表不同水平的人口增長、經(jīng)濟(jì)發(fā)展和其他社會(huì)經(jīng)濟(jì)、環(huán)境變量，其中A強(qiáng)調(diào)經(jīng)濟(jì)，B注重環(huán)境發(fā)展;1更注重全球化的變化，2則更注重局部發(fā)展。

根據(jù)我國土地資源分類系統(tǒng)，土地利用類型共分為兩級(jí)，該研究參照一級(jí)分類標(biāo)準(zhǔn)，將土地利用類型分為耕地、林地、草地、水域和居民用地5類。該研究利用區(qū)域2010年的土地利用類型，考慮不同的影響因素，模擬得出不同情景下2050年和2100年的土地利用變化（圖1）。楊溪河流域現(xiàn)狀（2010年）與未來（2050和2100年）不同情景下的土地利用類型統(tǒng)計(jì)情況見表1。模擬結(jié)果顯示，整體上流域內(nèi)各類土地利用類型更趨于集中，耕地面積大量減少并趨于消失，林地面積大量增加并遠(yuǎn)大于其他土地利用類型。其中，從2050年到2100年，A?1和B?1情景下草地面積均大幅減少，林地面積均大幅增加，到2100年該2種情景下林地面積占比均在96%以上，因此全球化的變化將使該流域內(nèi)的草地面積大幅減少、林地面積大幅增加。但B?1更追求綠色發(fā)展，其林地面積最大，為482.81 km2，占流域總面積的96.95%。到2100年，A?2情景下的水域面積增加最多，說明在注重經(jīng)濟(jì)和局部發(fā)展的情景下，該流域內(nèi)小水電等水利工程的建設(shè)將增多并導(dǎo)致其水域面積增加。B?2情景表示在注重環(huán)境發(fā)展的同時(shí)，亦保持人口增長處于中間水平[11]，因此該情景最顯著的一個(gè)特點(diǎn)是從2050到2100年其居民用地顯著增多。

3.2產(chǎn)水量和水源涵養(yǎng)量分析

基于InVEST模型對(duì)楊溪河流域的水源涵養(yǎng)功能進(jìn)行評(píng)價(jià)，并利用ArcGIS平臺(tái)對(duì)楊溪河流域水源涵養(yǎng)的空間分布和產(chǎn)水量情況進(jìn)行處理。該研究根據(jù)現(xiàn)狀條件中楊溪河干流上橫溪、釣魚臺(tái)和銀溪3宗水電站的壩址位置，將楊溪河流域劃分為Ⅰ、Ⅱ、Ⅲ和Ⅳ4個(gè)區(qū)域，這4個(gè)區(qū)域的面積分別為281.78、183.49、14.00和18.74 km2。圖2為根據(jù)楊溪河流域在2010年以及模擬得到的2050和2100年不同情景的土地利用下的產(chǎn)水量。楊溪河流域不同情景下Ⅰ、Ⅱ、Ⅲ和Ⅳ4個(gè)區(qū)域的產(chǎn)水量和水源涵養(yǎng)量情況見表2。根據(jù)圖2和表2可知，楊溪河流域2010年的產(chǎn)水量平均值為707.24 mm，水源涵養(yǎng)量為35 220.55萬m3;該流域2050年A?1、A?2、B?1和B?2情景下的產(chǎn)水量平均值分別為741.43、692.17、747.95和708.42 mm，水源涵養(yǎng)量分別為36 923.00萬、34 470.02萬、37 247.96萬和35 279.56萬m3;該流域2100年A?1、A?2、B?1和B?2情景下的產(chǎn)水量平均值分別為734.25、628.35、738.63和629.08 mm，水源涵養(yǎng)量分別為36 565.68萬、31 292.05萬、36 783.62萬和31 328.31萬m3。從各區(qū)的產(chǎn)水量來看，2010年、2050年B?1以及2100年A?2和B?2情景下各區(qū)的產(chǎn)水量大小情況為：Ⅳ區(qū)>Ⅰ區(qū)>Ⅱ區(qū)>Ⅲ區(qū);2050年A?1、A?2和B?2以及2100年A?1和B?1情景下各區(qū)的產(chǎn)水量大小情況為：Ⅰ區(qū)>Ⅳ區(qū)>Ⅱ區(qū)>Ⅲ區(qū)。因此，在現(xiàn)狀和未來4種情景下，該流域均表現(xiàn)為上游和下游的產(chǎn)水量較高。

4討論

土地利用類型的變化通過改變下墊面狀況進(jìn)而影響楊溪河流域的水源供給。下墊面是水文循環(huán)的重要場(chǎng)所，下墊面狀況的改變會(huì)影響陸面的下滲量、實(shí)際蒸散發(fā)量以及土壤理化性質(zhì)等。

楊溪河流域內(nèi)土地利用變化主要是人為因素造成的，人類的社會(huì)經(jīng)濟(jì)活動(dòng)和對(duì)大自然的改造會(huì)改變流域的土地利用類型，進(jìn)而對(duì)其產(chǎn)水量和水源涵養(yǎng)量造成影響。該研究發(fā)現(xiàn)，楊溪河流域2010年的產(chǎn)水量和水源涵養(yǎng)量與2050年B?2情景較接近，高于2050年A?2情景，低于2050年A?1和B?1情景。同時(shí)，該流域各情景下2050年的產(chǎn)水量和水源涵養(yǎng)量均比各情景下2100年的大，說明人類活動(dòng)會(huì)對(duì)流域產(chǎn)水量產(chǎn)生影響。該流域2050年和2100年各情景下的產(chǎn)水量和水源涵養(yǎng)量均表現(xiàn)為：B?1>A?1>B?2>A?2。因此注重全球化的變化，楊溪河流域能夠更好地保持其水源涵養(yǎng)功能;同時(shí)，注重環(huán)境保護(hù)，亦可保持流域良好的水源涵養(yǎng)功能。龔文明等[12-13]研究表明，林地覆蓋地區(qū)的水源涵養(yǎng)基礎(chǔ)條件優(yōu)越，能夠?qū)涤赀M(jìn)行有效蓄滲，其水源涵養(yǎng)功能較好;而居民用地因?yàn)槭瓜聣|面硬質(zhì)化，導(dǎo)致降雨消耗率偏高，其水源涵養(yǎng)功能較弱。

因此，未來人類活動(dòng)將對(duì)楊溪河流域的土地利用情況造成改變，為了維持楊溪河流域較強(qiáng)的水源涵養(yǎng)功能，一方面，流域相關(guān)部門需加強(qiáng)對(duì)林地的保護(hù)，嚴(yán)禁亂砍濫伐森林和亂占濫用林地等行為，另一方面則需要科學(xué)調(diào)控居民用地的擴(kuò)張速度和規(guī)模，協(xié)調(diào)好經(jīng)濟(jì)、環(huán)境和資源等各方面的關(guān)系。

5結(jié)語

該研究先利用GeoSOS-FLUS模型對(duì)楊溪河流域未來的土地利用變化進(jìn)行模擬，再利用InVEST模型的Water Yield模塊對(duì)楊溪河流域2010年、不同情景下2050年和2100年的土地利用情況進(jìn)行水源涵養(yǎng)功能評(píng)價(jià)。主要結(jié)論如下：

（1）GeoSOS-FLUS模型模擬結(jié)果顯示，整體上楊溪河流域內(nèi)未來各類土地利用類型更趨于集中，耕地面積大量減少并趨于消失，林地面積大量增加并普遍大于其他類型土地。

（2）楊溪河流域2010年各區(qū)的產(chǎn)水量情況為：Ⅳ區(qū)>Ⅰ區(qū)>Ⅱ區(qū)>Ⅲ區(qū);2050年和2100年各區(qū)的產(chǎn)水量情況均為：Ⅰ區(qū)、Ⅳ區(qū)>Ⅱ區(qū)>Ⅲ區(qū)。在現(xiàn)狀和未來4種情景下，該流域均表現(xiàn)為上游和下游的產(chǎn)水量較高。

（3）楊溪河流域各情景下2050年的產(chǎn)水量和水源涵養(yǎng)量均比各情景下2100年的大，說明人類活動(dòng)會(huì)對(duì)流域產(chǎn)水量產(chǎn)生影響。2050年和2100年各情景下的產(chǎn)水量和水源涵養(yǎng)量情況均為：B?1>A?1>B?2>A?2。注重全球化的變化和環(huán)境保護(hù)，可保持流域良好的水源涵養(yǎng)功能。

參考文獻(xiàn)

[1] 王盛萍，張志強(qiáng)，SUN G，等.基于MIKESHE模型的潮河流域土地利用與降水變化對(duì)水文的影響評(píng)價(jià)[J].生態(tài)與農(nóng)村環(huán)境學(xué)報(bào)，2012，28（3）：320-325.

[2] 王潤，劉洪斌，武偉.TOPMODEL模型研究進(jìn)展與熱點(diǎn)[J].水土保持研究，2005，12（1）：47-48，169.

[3] TALLIS H，RICKETTS T.InVEST 1.0 beta users guide：Integrated valuation of ecosystem services and tradeoffs[Z].San Francisco：Stanford，2010.

[4] 包玉斌，李婷，柳輝，等.基于InVEST模型的陜北黃土高原水源涵養(yǎng)功能時(shí)空變化[J].地理研究，2016，35（4）：664-676.

[5] 竇苗，孫建國，陳海鵬.基于InVEST模型的橫斷山區(qū)產(chǎn)水量模擬[J].安徽農(nóng)業(yè)科學(xué)，2017，45（21）：54-58.

[6] 王紀(jì)偉，劉康，甕耐義.基于In-VEST模型的漢江上游森林生態(tài)系統(tǒng)水源涵養(yǎng)服務(wù)功能研究[J].水土保持通報(bào)，2014，34（5）：213-217.

[7] 徐建寧，孫建國，陳海鵬.基于InVEST模型的生境質(zhì)量評(píng)估：以橫斷山小江流域?yàn)槔齕J].安徽農(nóng)業(yè)科學(xué)，2016，44（15）：105-108，134.

[8] 陳海鵬，竇苗.基于InVEST模型的云南小江流域生態(tài)系統(tǒng)碳儲(chǔ)量及空間格局研究[J].安徽農(nóng)業(yè)科學(xué)，2017，45（12）：51-54.

[9] 周文佐.基于GIS的我國主要土壤類型土壤有效含水量研究[D].南京：南京農(nóng)業(yè)大學(xué)，2003.

[10] LI X，CHEN G Z，LIU X P，et al.A new global landuse and landcover change product at a 1km resolution for 2010 to 2100 based on humanenvironment interactions[J].Annals of the American association of geographers，2017，107（5）：1040-1059.

[11] 劉祖發(fā)，姚寒梅，陳曉越，等.土地利用變化對(duì)茂名市降水入滲補(bǔ)給量的影響[J].生態(tài)環(huán)境學(xué)報(bào)，2017，26（12）：2102-2111.

[12] 龔文明.不同林分類型凋落物及土壤水源涵養(yǎng)功能差異分析[J].安徽農(nóng)業(yè)科學(xué)，2013，41（15）：6763-6766.

[13] 陳駿宇，劉鋼，白楊.基于InVEST模型的太湖流域水源涵養(yǎng)服務(wù)價(jià)值評(píng)估[J].水利經(jīng)濟(jì)，2016，34（2）：25-29.