基于InVEST模型的疏勒河上游產(chǎn)水量時空變化特征分析

魏培潔,吳明輝,賈映蘭,高雅月,徐浩杰,劉章文,陳生云1,,6,*

1 中國科學(xué)院西北生態(tài)環(huán)境資源研究院 冰凍圈科學(xué)國家重點實驗室 疏勒河源冰凍圈與生態(tài)環(huán)境綜合監(jiān)測研究站,蘭州 730000 2 中國科學(xué)院大學(xué),北京 100049 3 祁連山國家公園國家長期科研基地,西寧 810000 4 蘭州大學(xué)草地農(nóng)業(yè)科技學(xué)院草地農(nóng)業(yè)生態(tài)系統(tǒng)國家重點實驗室,蘭州 730020 5 中國科學(xué)院西北生態(tài)環(huán)境資源研究院 黑河上游生態(tài)水文試驗研究站,蘭州 730000 6 青海師范大學(xué)高原科學(xué)與可持續(xù)發(fā)展研究院 地理科學(xué)學(xué)院,西寧 810008

生態(tài)系統(tǒng)服務(wù)是指人類從各種生態(tài)系統(tǒng)中獲取的所有惠益,由生態(tài)系統(tǒng)的支持功能、供給功能、文化功能和調(diào)節(jié)功能及其相互作用形成[1—2]。生態(tài)系統(tǒng)功能是生態(tài)系統(tǒng)服務(wù)的基礎(chǔ)和物質(zhì)保障,其對人類生存和社會發(fā)展至關(guān)重要,特別是與水資源相關(guān)的功能[3]。隨著全球氣候變化、地表水資源短缺和水環(huán)境問題加劇,空間定量化、可視化評估區(qū)域產(chǎn)水量已成為當(dāng)前研究的熱點議題之一[4—8]。產(chǎn)水量是流域特別是干旱地區(qū)重要的調(diào)節(jié)功能,它對工業(yè)、農(nóng)業(yè)、漁業(yè)、人類消費和娛樂活動至關(guān)重要[9—10]。一方面,產(chǎn)水量是人類生存和區(qū)域經(jīng)濟可持續(xù)發(fā)展的重要物質(zhì)基礎(chǔ)；另一方面,產(chǎn)水量與區(qū)域自然地理條件和人類活動密切相關(guān)[10]。因此,研究區(qū)域產(chǎn)水量變化的驅(qū)動因素變得十分必要,其可為區(qū)域水資源的有效管理和可持續(xù)發(fā)展提供科學(xué)依據(jù)[10—11]。

自1970年以來,遙感技術(shù)和水文模型發(fā)展迅速,越來越多的學(xué)者嘗試通過模型模擬方法來定量化、可視化和精細(xì)化地評估和分析區(qū)域產(chǎn)水量,如MIKESHE(MIKE System Hydrological European)模型、TOPMODEL模型(TOPgraphy based hydrological Model)、SWAT模型(Soil and Water Assessment Tool)和InVEST模型等[12—13]。InVEST是一款由美國斯坦福大學(xué)、世界自然基金會、大自然保護(hù)協(xié)會和其他相關(guān)機構(gòu)聯(lián)合開發(fā)用于生態(tài)系統(tǒng)服務(wù)功能評估的模型方法[14]。InVEST模型產(chǎn)水量模塊基于水量平衡原理,充分考慮不同土地利用類型下土壤滲透性的空間差異及地形、地表粗糙程度等因素對徑流影響,以柵格為單元定量估算不同土地利用類型的產(chǎn)水量[12—13]。其模型由于參數(shù)調(diào)整靈活、評估結(jié)果空間表達(dá)性強,且能夠?qū)崿F(xiàn)時空間和多目標(biāo)權(quán)衡[15],被國外許多國家和地區(qū)廣泛應(yīng)用。如：Leh等[16]、Marquèsa等[17]和Redhead等[18]利用InVEST模型分別對西非科特迪瓦、西班牙東北部弗蘭科利河流域和英國多條河流的產(chǎn)水量進(jìn)行了評估,并取得了很好的應(yīng)用效果。此外,中國學(xué)者利用InVEST模型的產(chǎn)水量模塊分別對北京山區(qū)[19]和貴水流域[20]、陜西省商洛市五大水系流域[21]、三江并流區(qū)(云南、四川和西藏的交界處)[22]以及高原湖泊[23]等區(qū)域的產(chǎn)水量及其時空分布特征進(jìn)行了科學(xué)評價,模型在各研究區(qū)均取得了很好的應(yīng)用效果。同時,InVEST模型已成功用于青藏高原地區(qū)。其中,潘濤等[7]和呂樂婷等[24]分別評估了三江源區(qū)1981—2010、1980—2016年產(chǎn)水量的時間變化特征與空間分布狀況。王雅等[25]和趙亞茹等[10]分別對位于祁連山區(qū)中部的黑河流域、東部的石羊河流域的產(chǎn)水量空間分異特征進(jìn)行了分析,這些研究對西北地區(qū)水資源規(guī)劃和管理都具有重要意義。

目前,氣候變化和土地利用變化是影響產(chǎn)水量的關(guān)鍵因素[26]。其中,氣候變化在宏觀尺度上可以通過改變流域的降水量和蒸散量來影響產(chǎn)水量[27—28],土地利用可通過改變流域下墊面間接影響產(chǎn)水量[29—30]。一些研究已探究了氣候變化和土地變化對流域產(chǎn)水量的影響。例如：Pessacg等[9]利用降水?dāng)?shù)據(jù)集評估了楚布特河流域的產(chǎn)水量,表明降水量與產(chǎn)水量之間存在正相關(guān)關(guān)系；Stone等[31]重點研究了氣候變化對密蘇里河流域產(chǎn)水量的影響,證實氣候變化對產(chǎn)水量影響顯著；Zhang等[32]和Li等[33]通過分析土地利用變化對產(chǎn)水量影響,得出不同土地利用方式的改變對產(chǎn)水量的影響不同。目前,全球氣候變化對中國水資源的影響日益加劇[34]。青藏高原作為對全球環(huán)境變化特別是氣候變化響應(yīng)最敏感的地區(qū)之一,近年來受氣候變暖影響,50年來氣溫顯著升高0.3 ℃/(10a),部分區(qū)域降水波動較大[35],并且氣候暖濕化促進(jìn)了植被生長,NDVI顯著增加[36]。然而,青藏高原地區(qū)氣候變化和土地利用變化對產(chǎn)水量的相關(guān)研究仍然缺乏,尤其植被根系深度和NDVI對產(chǎn)水量的影響更是鮮有研究。

水在維持流域生態(tài)系統(tǒng)平衡和生態(tài)承載力可持續(xù)性方面扮演著極其重要的角色[37]。干旱區(qū)水資源空間分布不均衡,降水多集中于上游山區(qū),上游的出山口徑流量基本代表著流域的總水資源[38]。疏勒河流域地處青藏高原東北緣祁連山西段,是我國河西走廊內(nèi)陸干旱區(qū)三大內(nèi)陸河流域之一,也是河西農(nóng)牧民群眾的“生命線”和“天然水塔區(qū)”[39]。然而,疏勒河流域產(chǎn)水量相關(guān)的研究報道甚少。開展上游產(chǎn)水量研究,對整個流域水資源的動態(tài)評估、生態(tài)平衡的維系以及氣候調(diào)節(jié)具有十分重要的意義[10,15]。鑒于此,基于InVEST模型產(chǎn)水量模塊,本研究對疏勒河上游高寒山區(qū)產(chǎn)水量時空變化特征進(jìn)行了分析,并進(jìn)一步探討其對環(huán)境要素變化的響應(yīng),旨在為流域水資源的有效管理和可持續(xù)發(fā)展提供科學(xué)參考。

1 研究區(qū)概況

疏勒河流域地處我國西北干旱區(qū)腹地,是河西走廊僅次于黑河的第二大內(nèi)陸河[40—41],其水源補給主要源于冰川和凍土融水以及山區(qū)降水[37,42]。疏勒河上游(96°37′12″—98°59′24″ E,38°13′12″—39°52′12″ N)是指以昌馬堡水文站(96°51′ E,39°49′ N)為出水控制點的上部區(qū)域,為疏勒河干流的形成區(qū)[42],流域面積約為10973.9 km2,海拔2031—5763 m(圖1),多年平均徑流量為13.81×108m3。地形上主要由疏勒南山、托來南山和疏勒河谷地組成,山地地形陡峭,河谷地形相對低緩[38]。氣候類型為典型的大陸性氣候,按陸地表層干濕指數(shù)區(qū)劃屬于半干旱區(qū)—干旱區(qū)邊緣地帶,呈現(xiàn)干寒氣候條件下的草原—荒漠景觀[38—40]。土壤類型主要為高山寒漠土、高山草甸草原土、栗鈣土、淡栗鈣土和山地灰鈣土[40],土壤質(zhì)地主要以壤土為主,砂土和粘土占比均較小。疏勒河上游是整個流域的水源涵養(yǎng)區(qū)和產(chǎn)流區(qū),為中下游山前平原和綠洲地區(qū)賴以生存和發(fā)展提供了水源保障。

圖1 研究區(qū)位置和水文站點分布圖Fig.1 Location of study area and hydrological station

2 研究方法與數(shù)據(jù)處理

2.1 產(chǎn)水量模型

InVEST產(chǎn)水量模塊基于水量平衡原理,用各柵格單元的降水量減去實際蒸散后的水量,進(jìn)而得到該柵格產(chǎn)水量。該模型假設(shè)每個柵格單元所有的產(chǎn)水量通過地下徑流或地表徑流的方式到達(dá)流域出水口[14]。模型主要算法如下[7—8,43]：

(1)

式中,Yxj為第j土地利用類型柵格x的產(chǎn)水量(mm),AETxj為第j土地利用類型柵格x的年實際蒸散發(fā)量(mm),Px為柵格x的多年平均降水量(mm)。

(2)

式中,Rxj為第j土地利用類型柵格x的Budyko干燥指數(shù),ωx表示植被有效含水量與年均降水量的比值：

(3)

(4)

AWCx=Min(DS,Dr)×PAWCx

(5)

式中,AWCx為柵格x的植被有效含水量(mm)(由有效土壤深度和土壤質(zhì)地決定)；Z為季節(jié)參數(shù)(即Zhang系數(shù))用于表征降水的季節(jié)性特征,kxj為第j土地利用類型柵格x的植被蒸散系數(shù),ET0x表示柵格x的參考蒸散量(mm);Min為minimum的縮寫,DS為土壤深度(mm),Dr為根系深度(mm),PAWCx為柵格單元x的植被可利用水分(mm)。

2.2 參數(shù)來源及分析

2.2.1模型參數(shù)

產(chǎn)水量模塊所需輸入?yún)?shù)有氣象、土地利用、土壤深度、植被可利用水、生物物理參數(shù)表、流域以及Zhang系數(shù)。另外,相關(guān)分析涉及NDVI。本研究所有數(shù)據(jù)的空間分辨率均重采樣到500 m,且坐標(biāo)系均使用世界大地坐標(biāo)系統(tǒng)(WGS1984)。

(1)降水量和參考蒸散發(fā)

利用中國氣象科學(xué)數(shù)據(jù)中心(http://data.cma.cn/user/toLogin.html)收集周邊氣象站2001—2019年的降水量和氣溫數(shù)據(jù)。用ANUSPLIN模型進(jìn)行空間插值,得到降水量和氣溫的空間柵格數(shù)據(jù)。本文選用Modified—Hargreaves法計算參考蒸散發(fā),與其它方法相比,該方法只需最高氣溫、最低氣溫和太陽大氣頂層輻射便可計算[45]。公式如下：

(6)

式中,Tmax和Tmin分別為月最高氣溫和月最低氣溫(℃),Ra是月太陽大氣頂層輻射(MJ m-2d-1),C、E和T是常數(shù),其值參考胡慶芳等[46]推薦青藏高原C、E和T值。

(2)土地利用

土地利用數(shù)據(jù)源于美國宇航局地球科學(xué)數(shù)據(jù)系統(tǒng)(https://search.earthdata.nasa.gov/search)中的MCD12Q1數(shù)據(jù)產(chǎn)品。它是是根據(jù)一年的Terra和Aqua觀測所得的數(shù)據(jù)經(jīng)過處理,描述土地覆蓋的類型。MCD12Q1采用了五種不同的土地覆蓋分類方案,信息提取的主要技術(shù)是監(jiān)督?jīng)Q策樹分類。本文采用國際地圈—生物圈計劃(IGBP)的分類方法。研究區(qū)土地利用類型主要包括草地、永久性濕地、農(nóng)田、建設(shè)用地、農(nóng)作物—自然植被交錯帶、永久性冰雪和裸地。

(3)土壤數(shù)據(jù)和PAWC

基于蘭州市寒區(qū)旱區(qū)科學(xué)數(shù)據(jù)中心下載的第二次土壤調(diào)查數(shù)據(jù)庫生成土壤深度和土壤質(zhì)地柵格圖。PAWC是指土壤土層中為植物生長提供的水量所占比例,為0—1的小數(shù)。周文佐等[47]提出了一種根據(jù)物理和化學(xué)性質(zhì)計算PAWC的方法,具體算法如下;

PAWC=54.509-0.132Sand-0.003(Sand)2-0.055Silt-0.006(Silt)2-0.738Clay+0.007(Clay)2-2.668C+

0.501(C)2

(7)

式中,Sand為土壤砂粒含量(%),Silt為土壤粉粒含量(%),Clay為土壤粘粒含量(%)；C為土壤有機質(zhì)含量(%)。

(4)流域和生物物理參數(shù)表

數(shù)字高程模型(DEM)由資源環(huán)境科學(xué)與數(shù)據(jù)中心(http://www.resdc.cn/)提供的雷達(dá)地形測繪SRTM(90 m)數(shù)據(jù)?；跀?shù)字高程模型數(shù)據(jù)通過ArcGIS 10.4中的水文分析工具生成的子流域。生物物理表主要用來反映土壤覆蓋度和LULC的屬性,包括LULC編碼、植物蒸散系數(shù)(Kc)和根系深度等,生物物理表中每種LULC類型的系數(shù)來源于文獻(xiàn)[45,48]和《InVEST用戶指南》[49]。

(5)Zhang系數(shù)

Zhang系數(shù)是反映當(dāng)?shù)亟邓窬趾退牡刭|(zhì)特征的氣候季節(jié)因子,其值在1—30之間變化。根據(jù)昌馬鋪水文站多年平均徑流量觀測資料,采用水量平衡法對Zhang系數(shù)進(jìn)行校驗[7],當(dāng)Zhang系數(shù)為3.33時,InVEST產(chǎn)水模塊的模擬效果最好。

(6)NDVI與植被根系深度

NDVI可以準(zhǔn)確反映地表植被覆蓋狀況,是各種植被指數(shù)中最常用的表征地表植被狀況的指標(biāo)[36]。本研究中所涉及的2001—2019年NDVI數(shù)據(jù)源于資源環(huán)境科學(xué)與數(shù)據(jù)中心(http://www.resdc.cn/)。該數(shù)據(jù)集是基于SPOT/VEGETATION NDVI衛(wèi)星遙感數(shù)據(jù),在月數(shù)據(jù)的基礎(chǔ)上采用最大值合成法生成年度NDVI數(shù)據(jù)集。植被根系深度柵格數(shù)據(jù)是InVEST模型基于生物物理參數(shù)表生成的。

2.2.2時空變化趨勢分析

本研究主要通過OriginPro 9.1中的線性擬合方法分析研究目標(biāo)(氣溫、降水量和實際蒸散量、NDVI和產(chǎn)水量)的時間變化趨勢。利用Matlab—R2015B對研究目標(biāo)的空間變化趨勢進(jìn)行線性回歸分析,該方法廣泛應(yīng)用于時序數(shù)據(jù)動態(tài)檢測,線性回歸的斜率是量化動態(tài)變化趨勢的最佳指標(biāo)。具體算法如下[50]：

(8)

式中,n為研究年限(本研究n=19),i為年份序號,Mi為第i年的研究目標(biāo)值,Slope>0表示19年來研究目標(biāo)變化趨勢為正,Slope<0則表示研究目標(biāo)變化趨勢為負(fù)。

利用線性回歸的相關(guān)系數(shù)對研究目標(biāo)動態(tài)趨勢進(jìn)行顯著性檢驗,相關(guān)系數(shù)的計算公式如下：

(9)

式中,i為研究期間的年份序號,Mi為第i年的研究目標(biāo)值,cov為協(xié)方差函數(shù),var為方差函數(shù)。

2.2.3產(chǎn)水功能重要性分級

合理分配、利用和保護(hù)水資源的前提是需對區(qū)域產(chǎn)水功能重要性進(jìn)行分級。本研究基于InVEST模型對研究區(qū)產(chǎn)水量進(jìn)行了定量評估,并依據(jù)《國家生態(tài)保護(hù)紅線—生態(tài)功能紅線劃定技術(shù)指南(試行)》中涵養(yǎng)水源重要生態(tài)功能區(qū)的分級方法[51],利用ArcGIS 10.4中的Quantile(分位數(shù))功能對2001—2019年的年均產(chǎn)水量進(jìn)行分類。按產(chǎn)水量值將其重要性由低到高依次劃分為4個重要性級別,即一般重要[0.00—98.38)、中等重要[98.38—227.33)、高度重要[227.33—344.76)和極重要[344.76—488.72)。

2.2.4統(tǒng)計分析

利用MODIS Reprojection Tools與ArcGIS對所獲取的柵格數(shù)據(jù)進(jìn)行拼接、提取、裁剪等操作,以得到本研究所需的數(shù)據(jù)集。其次,將2001—2019年昌馬堡水文站實測地表年徑流量與模擬年總產(chǎn)水量進(jìn)行線性擬合,并計算調(diào)整R2、P以校驗?zāi)Ｐ徒Y(jié)果；擬合直線斜率與1∶1直線斜率間的差異通過SPSS 22.0中的一般線性模型進(jìn)行分析。另外,利用R v4.0.4對環(huán)境要素與產(chǎn)水量之間相關(guān)性進(jìn)行Pearson相關(guān)矩陣分析。

3 結(jié)果與分析

3.1 模型校驗

將模型預(yù)測結(jié)果與實測數(shù)據(jù)進(jìn)行比較,通過不斷調(diào)整參數(shù),得到最適于研究區(qū)的預(yù)測模型則為模型校驗[14]。本研究利用昌馬堡水文站地表徑流量實測數(shù)據(jù)對InVEST模型產(chǎn)水量的預(yù)測結(jié)果進(jìn)行了校驗。結(jié)果顯示(圖2)：模型模擬的多年平均總產(chǎn)水量為13.66×108m3,水文站實測多年平均總徑流量為13.81×108m3,其相對誤差為1.11%。同時,模型模擬的總產(chǎn)水量與實測地表總徑流量之間存在較強的線性關(guān)系(R2= 0.986P<0.01)。雖然,模擬總產(chǎn)水量相較實測地表總徑流量略低估了0.37×108m3/a,但兩者所擬直線斜率與1沒有顯著差異(P>0.05)?？梢?InVEST模型在本研究區(qū)適用性很好。

圖2 InVEST 模擬總產(chǎn)水量與實測地表總徑流量對比 Fig.2 InVEST modelled total water yield vs measured total surface runoff紅線代表斜率為0,截距為1的直線

3.2 產(chǎn)水量時空分布及變化特征

2001—2019年疏勒河上游產(chǎn)水量顯著增加,增加速率為0.26×108m3/a(圖3)。2003年總產(chǎn)水量最低(9.76×108m3),2017年總產(chǎn)水量最高(17.25×108m3)。從空間分布格局來看,山地高海拔區(qū)產(chǎn)水量較高,河谷低海拔區(qū)產(chǎn)水量較低,尤其是東南部低海拔區(qū)與西北部東側(cè)最低(圖4)。2001—2019年,所有柵格的產(chǎn)水變化量空間趨勢如圖4所示。產(chǎn)水量顯著增加區(qū)域約占整個研究區(qū)的32.01%,主要分布在山地高海拔地區(qū)；僅有8.39%的區(qū)域顯著下降且集中分布在研究區(qū)的西北部。

圖3 2001—2019年產(chǎn)水量時間變化Fig.3 Temporal variation of water yield in 2001—2019

圖4 產(chǎn)水量空間分布及變化Fig.4 Spatial distribution and its variation of water yield

3.3 產(chǎn)水功能重要性分級

研究區(qū)產(chǎn)水功能重要性分級總體呈現(xiàn)出由河谷低海拔區(qū)向山地高海拔區(qū)遞增的階梯分布趨勢(圖5)。其中,一般重要區(qū)域所占面積最大,比例為47.26%,主要分布在東南部低海拔區(qū)與西北部東側(cè)；中等重要區(qū)域所占比例為24.48%,主要分布在西北部；高度重要和極重要區(qū)域所占比例分別為14.51%和13.74%,集中分布于山地高海拔區(qū)。

圖5 產(chǎn)水功能重要性空間分級Fig.5 Importance spatial classification of water yield function

3.4 影響流域產(chǎn)水量的因素分析

疏勒河上游的平均海拔高于3945 m,導(dǎo)致該地區(qū)氣候寒冷干燥,年均氣溫為-6.06℃,年降水量325.96 mm。2001—2019年研究區(qū)年均氣溫、年降水量和年實際蒸散發(fā)的時間變化趨勢如圖6所示。年均氣溫與年降水量分別以0.02 ℃/a (P>0.05)和5.92 mm/a(P<0.05)的速率增加,研究區(qū)氣候呈現(xiàn)暖濕化趨勢。在此影響下,年實際蒸散量和NDVI分別以3.00 mm/a (P>0.05)和0.003/a (P<0.05)的速率遞增。

圖6 年均氣溫、年降水量、年實際蒸散量和年NDVI的時間變化Fig.6 Temporal variations of the mean annual air temperature, annual precipitation, annual actual evapotranspiration and annual NDVI

研究區(qū)氣溫、降水量、實際蒸散量和NDVI的空間分布均具有明顯的空間異質(zhì)性。具體而言,氣溫的空間變化范圍為-16.92—5.33℃,自西北向東南逐漸下降且河谷低海拔區(qū)高于山地高海拔區(qū),2001—2019年研究區(qū)大部分區(qū)域(76.08%)的氣溫呈明顯上升趨勢,但未達(dá)到統(tǒng)計學(xué)顯著水平(P>0.05)。空間上,降水量在127.42—527.13 mm之間變化,總體呈西北低、東南高且河谷低海拔區(qū)低于山地高海拔區(qū)的分布格局(圖7),與氣溫的分布情況恰好相反。研究期間,降水量整體顯著增加(P<0.05)且東南局部地區(qū)(24.41%)極顯著增加(P<0.01)。實際蒸散量的空間變化范圍為0.00—357.54 mm, 其高值區(qū)主要位于東南部低海拔區(qū)與西北部東側(cè)；2001—2019年實際蒸散量顯著增加區(qū)域約為33.39%(P<0.05),主要分布河谷低海拔區(qū),而顯著下降區(qū)域約為4.01%(P<0.05),散落分布在整個研究區(qū)。NDVI空間分布格局與實際蒸散量的分布情況一致,研究期間,NDVI顯著增加區(qū)域約為56.92%(P<0.05),分布在河谷低海拔局部地區(qū),極顯著下降區(qū)域約占11.27%(P<0.05),主要分布在山地高海拔區(qū)。此外,植被根系深度的空間分布與NDVI一致。

圖7 環(huán)境要素空間分布及其變化Fig.7 Spatial distributions and the variations of environmental elements

Pearson相關(guān)分析表明(圖8),產(chǎn)水量與降水量顯著正相關(guān),卻與氣溫、實際蒸散量、NDVI和根系深度顯著負(fù)相關(guān)(P<0.001)。不同環(huán)境因素之間的相關(guān)性表明,實際蒸散量與根系深度、NDVI和氣溫顯著正相關(guān)(P<0.001)。氣溫升高或NDVI、根系深度增加均可能會促進(jìn)實際蒸散量增加,從而對產(chǎn)水量產(chǎn)生負(fù)向影響。

圖8 產(chǎn)水量與影響因子的相關(guān)矩陣Fig.8 Correlation matrix of water yield and driving factors***, P<0.001

4 討論

產(chǎn)水量是一項重要的生態(tài)功能指標(biāo),直接影響著社會經(jīng)濟和生態(tài)系統(tǒng)的可持續(xù)發(fā)展[34]。本研究著重對疏勒河上游產(chǎn)水量的時空變化特征以及影響因素進(jìn)行了探討分析。結(jié)果表明,2001—2019年疏勒河上游產(chǎn)水量顯著上升,這同呂樂婷等[24]對1980—2016年三江源區(qū)的研究結(jié)果類似,卻同王雅等[25]對2000—2014年祁連山黑河中游的研究結(jié)果相反。這可能與各研究區(qū)的氣候條件以及研究的時間尺度不同有關(guān)。產(chǎn)水量是區(qū)域水循環(huán)過程中綜合考慮收(降水量)支(實際蒸散量)平衡的產(chǎn)物[10],因此,區(qū)域降水量和實際蒸散量以及兩者之間平衡關(guān)系是影響產(chǎn)水量和產(chǎn)水能力的主導(dǎo)要素,這是許多產(chǎn)水量研究中的一致結(jié)論[7,52—53]。潘韜等[7]和呂樂婷等[24]對三江源產(chǎn)水量的研究得出,降水和實際蒸散發(fā)是決定生態(tài)系統(tǒng)產(chǎn)水量的兩個關(guān)鍵環(huán)節(jié),其中降水量與產(chǎn)水量顯著正相關(guān),實際蒸散量與其顯著負(fù)相關(guān)。劉軍會等[52]和劉婭等[53]分別對青藏高原和京津冀地區(qū)的產(chǎn)水量進(jìn)行了研究,結(jié)果表明降水量和實際蒸散量綜合影響區(qū)域產(chǎn)水變化量。然而,氣溫、NDVI和植被根系深度會通過影響實際蒸散量來間接影響產(chǎn)水量。其中,NDVI和植被根系深度主要通過下墊面植被變化來影響實際蒸散量,進(jìn)而影響產(chǎn)水量。Lu等[35]對青藏高原拉薩河流域產(chǎn)水量的研究發(fā)現(xiàn),降水量和NDVI變化是影響區(qū)域產(chǎn)水量的關(guān)鍵要素,特別是通過植被變化來推動產(chǎn)水量變化。19年來,疏勒河上游實際蒸散增加趨勢(3.00 mm/a)低于降水增加趨勢(5.92 mm/a),降水量對產(chǎn)水量的補給遠(yuǎn)超過實際蒸散對產(chǎn)水量的消耗,因此,疏勒河上游產(chǎn)水量整體呈顯著增加趨勢。

疏勒河上游發(fā)育大面積多年凍土,其中山地高海拔區(qū)的土層存在凍結(jié)面,這使得降水無法下滲而在凍結(jié)面匯流形成徑流。與此同時,山地高海拔區(qū)的降水量高,氣溫和實際蒸散量低且土地利用類型主要為冰川和裸地。冰川屬于水域,最易融化形成徑流；裸地的植被覆蓋度低且根系淺,截留能力弱,更多降水也直接形成徑流[24,54—55]。相比之下,河谷低海拔區(qū)由于凍土分布面積和深度受限,降水或冰水可以滲入土壤深層,徑流形成的幾率大大降低。另外,河谷低海拔區(qū)的主要土地利用類型為草地,其植被覆蓋度和實際蒸散量都很高,并且植被根系深而密集,可以有效的截留降水,使得產(chǎn)水量偏低。2001—2019年山地高海拔區(qū)產(chǎn)水量顯著增加,主要是因為降水量補給顯著增加而蒸散量無顯著變化,致使生態(tài)產(chǎn)水以增加為主。與其相反,西北地區(qū)的裸地被草地和農(nóng)田大面積取代,相比裸地,草地與農(nóng)田的根系較深且植被生長消耗的水量也比較大,因此產(chǎn)水量顯著下降。

疏勒河上游產(chǎn)水重要性分級與產(chǎn)水量空間分布格局基本一致,這與孫艷偉等[14]和王玉純等[51]研究結(jié)論一致。具體而言,產(chǎn)水極重要和高度重要區(qū)降水量較多、蒸散量較小,并且多分布冰川和裸地,是疏勒河流域水資源集成區(qū),而產(chǎn)水一般重要和中等重要地區(qū)為相對干旱的堆積平原和盆地,為水資源消耗區(qū)。為使研究區(qū)實現(xiàn)健康可持續(xù)發(fā)展,極重要和高度重要區(qū)需作為重點和優(yōu)先保護(hù)的水資源地而納入政策和管理決策之中[56]。

本研究對疏勒河上游產(chǎn)水量的模擬結(jié)果較為理想,基本符合區(qū)域的實際情況。但因模型自身設(shè)定、參數(shù)算法以及數(shù)據(jù)精度等原因,使研究仍存在一些不確定性。首先,InVEST模型簡化了匯流過程,沒有區(qū)分地表徑流、土壤中徑流和基流,并假定所有柵格的產(chǎn)水量都到達(dá)同一出水口[14],并且產(chǎn)水量常與社會經(jīng)濟和人類活動密切相關(guān),而模型的輸入數(shù)據(jù)只有自然數(shù)據(jù)。特別是,疏勒河上游是一個典型的高寒山區(qū),InVEST模型未考慮冰川和凍土融水,這會增加水文模擬的不確定性。因此,InVEST模型在這些方面都有待進(jìn)一步改進(jìn)。其次,降水量和參考蒸散量是精確模擬產(chǎn)水量的重要參數(shù),而西北地區(qū)水文站點較少且氣候變化具有復(fù)雜性,利用ANUSPLIN模型空間插值生成的降水?dāng)?shù)據(jù)以及通過Modified—Hargreaves法計算的參考蒸散發(fā)的準(zhǔn)確性都會受到氣象數(shù)據(jù)精度的影響[57],從而影響模型的模擬精度。盡管存在這些不確定性,但研究結(jié)果仍可以有效反映產(chǎn)水量的時間變化和空間分布,揭示產(chǎn)水量與環(huán)境要素之間的關(guān)系,為疏勒河流域水資源的有效管理、合理利用和保護(hù)提供科學(xué)參考。

5 結(jié)論

基于InVEST模型,結(jié)合環(huán)境要素對疏勒河上游產(chǎn)水量的時空分布格局及影響因素進(jìn)行分析,得出：

(1) 模型中的Zhang系數(shù)為3.3時,研究區(qū)產(chǎn)水量模擬效果最佳；研究區(qū)年均總產(chǎn)水量為13.66×108m3,2001—2019年研究區(qū)產(chǎn)水量以0.26×108m3/a速率顯著增加。從空間分布來看,山地高海拔區(qū)產(chǎn)水量較高,河谷低海拔區(qū)產(chǎn)水量較低,尤其東南部低海拔區(qū)和西北部東側(cè)最低。研究期間,產(chǎn)水量顯著增加區(qū)約占總面積的32.01%,主要分布在山地高海拔區(qū)；僅有8.39%的研究區(qū)顯著下降且集中分布于西北部。

(2) 產(chǎn)水量與降水量顯著正相關(guān),而與氣溫、實際蒸散量、NDVI和根系深度顯著負(fù)相關(guān)；實際蒸散量與氣溫、NDVI和植被根系深度顯著正相關(guān)。

(3) 在對研究區(qū)產(chǎn)水量進(jìn)行定量評估的基礎(chǔ)上,進(jìn)行了產(chǎn)水功能重要性空間分級。得出：一般重要區(qū)域所占面積最大,約占整個研究區(qū)的47.26%,主要位于東南部低海拔區(qū)與西北部東側(cè)；高度重要和極重要區(qū)域共占總面積的28.25%,集中分布于山地高海拔區(qū)。