低丘緩坡建設開發對水源涵養功能的影響

石雅馨 張洪 郭潤展

摘要 運用InVEST模型并結合實地調研結果，對云南10個典型低丘緩坡建設項目區開發前后生態系統水源涵養功能進行定量測算，比較分析不同類型項目區水源涵養功能的差異。結果表明，開發后項目區平均產水量和平均產水深度均低于開發前，城市建設區的產水總量高于工業區，下降程度也高于工業區;大部分項目區水源涵養總量均低于開發前，僅光華項目區水源涵養總量高于開發前，城市建設區的水源涵養總量高于工業區，工業區水源涵養量下降程度高于城市建設區;導致項目區開發前后水源涵養量變化的原因有土地利用結構、植被、地形、徑流、降水、土壤等自然和人為因素，其中受土地利用類型影響較大;低丘緩坡開發導致區域生態環境受到干擾和破壞，但不同的開發模式和功能定位對其影響存在差異。因此在低丘緩坡開發過程中，應采取工程、監管、法律和經濟等綜合措施，實現低丘緩坡開發建設與生態保護的協調。

關鍵詞 低丘緩坡;建設開發;水源涵養;InVEST模型

中圖分類號 S 157.1? 文獻標識碼 A

文章編號 0517-6611（2021）21-0070-09

doi：10.3969/j.issn.0517-6611.2021.21.018

開放科學（資源服務）標識碼（OSID）：

Effect of Construction and Development of Low Hills and Gentle Slopes on Water Conservation Function

SHI Ya-xin， ZHANG Hong， GUO Run-zhan

（School of Urban and Environment， Yunnan University of Finance and Economics， Kunming， Yunnan 650221）

Abstract Using the InVEST model and the results of field investigations， the water conservation functions of the ecosystems before and after the development of 10 typical low-hill and gentle slope construction project areas in Yunnan were quantitatively measured， and the differences in water conservation functions of different types of project areas were compared and analyzed.The results show that the average water production and average water production depth in the project area after development are lower than before development，the total water production of urban construction areas was higher than that of industrial areas， and the decrease degree was also higher than that in industrial areas.The total amount of water conservation in most project areas is lower than before the development. Only the total amount of water conservation in the Guanghua project area was higher than that before the development，the total amount of water conservation in urban construction areas was higher than that in industrial areas， and the decline in water conservation in industrial areas was higher than that in urban construction areas.The reasons for the changes in water conservation before and after the development of the project area are natural and human factors such as land use structure， vegetation， topography， runoff， precipitation and soil， among which the land use

type had a greater impact.The development of low hills and gentle slopes had caused disturbance and destruction of the regional ecological environment， but different development models and functional positioning had different impacts on them.Therefore， in the development process of low hills and gentle slopes， comprehensive measures such as engineering， supervision， law and economy should be adopted to realize the coordination of the development and construction of low hills and gentle slopes with ecological protection.

Key words Low hills and gentle slopes;Construction and development;Water conservation;InVEST model

基金項目 國家自然科學基金項目（71764034）。

作者簡介 石雅馨（1996—），女，湖南新邵人，碩士研究生，研究方向：土地利用規劃。*通信作者，二級教授，博士生導師，從事城市發展與土地利用研究。

收稿日期 2021-06-17

生態系統的水源涵養功能主要體現在改善水文狀況、調節區域水分循環等方面[1]。生態系統中的植被如森林、灌叢、草地對水源涵養功能發揮著巨大的作用，其通過林冠層截留、枯枝落葉層持水和土壤層蓄水來調節地表徑流和補充地下水，在防治旱澇災害、凈化水質等方面具有至關重要的作用[2]。生態系統服務和權衡的綜合評估模型（integrated valuation of ecosystem services and trade-offs，InVEST）從水文角度以中小尺度流域單元為對象，反映不同氣候狀況、土地利用、地形起伏、土壤性質下的產水量大小，功能模型以地圖形式表達了流域的產水能力，間接地對區域的自然價值賦予了適當的衡量標準，而數據與參數的適宜性成為模型結果是否可信的關鍵[3]。近年來國內學者將InVEST模型廣泛運用于北京山區、陜北黃土高原、秦嶺地區、橫斷山區、若爾蓋高原等山地[4-8]，均取得較好的應用效果。國外學者將InVEST模型廣泛應用于美國、英國等國家和地區的產水量計算中[9-11]。

目前，國內關于生態系統服務的研究報道較多，但對整個云南低丘緩坡生態系統服務的研究較少且時序較早。陳武強等[12]研究昆明西華濕地公園水分調節量和價值，趙元藩等[13]估算玉溪市森林生態系統涵養水源量和調節水量價值，楊芳等[14]評估玉龍縣天保工程森林生態系統水源涵養能力，都是從宏觀上估算了區域的水源涵養能力以及為保護環境創造的價值。云南省山地面積占全省土地面積的90%以上，地質環境復雜，生態環境多樣。山地環境具有空間異質性和垂直分異性，海拔、地形起伏和土地覆被在一定的空間尺度上鑲嵌復合。水循環受土地利用、氣候、地形、土壤、植被等自然和人為因素的復雜影響，如何量化這種水源涵養能力，仍是一個懸而未決的問題[15]。鑒于上述原因，筆者針對氣候、地形、土壤、土地利用對山地水循環的影響，采用InVEST產水量模型與地形指數、土壤飽和導水率、流速系數結合的方法，估算研究區水源涵養能力，旨在揭示低丘緩坡綜合開發前后水源涵養功能的時空變化與影響因素，為低丘緩坡合理開發和水資源保護提供科學依據。

1 資料與方法

1.1 研究區概況 該研究在選取研究區時，主要考慮的是功能定位、地理位置、氣候類型和開發情況4個因素，具體理由如下：①功能定位，依據功能定位，項目區大致可分為2種類型，一類是以加工、能源、倉儲物流、高新技術為主的工業區，且具有特色產業;另一類是集旅游、康養、居住等功能的城市建設區。②地理位置，該研究根據地理區位劃分，從滇中選擇6個，滇西選擇2個，滇南選擇2個項目區。③氣候類型，云南省大部分地區屬于亞熱帶季風氣候，滇南部分地區屬于熱帶季風氣候。為了避免氣候類型對結果產生較大影響，該研究所選項目區8個屬于亞熱帶季風氣候，景洪的光華和景大項目區屬于熱帶季風氣候。④開發情況，對項目區進行實地調研，選擇開發程度較高的項目區進行研究。

研究區位于云南省昆明市、曲靖市、玉溪市、大理白族自治州和西雙版納傣族自治州（表1和圖1）。研究區面積在550.98～2 158.40 hm 研究區內土壤類型多樣，有水稻土、暗黃棕壤、黃棕壤、紅壤、山原紅壤、赤紅壤、磚紅壤、褐紅土。植被類型豐富，主要為闊葉林、針葉林、針闊混交林、闊葉灌叢、灌草叢和山地草甸。

1.2 InVEST產水量模型

InVEST模型基于水量平衡原理，通過研究區的降水、地表蒸發、植物蒸騰、土壤深度、根系深度、地形等數據計算每個柵格的產水量。充分考慮不同土地利用類型下土壤滲透的空間差異和地形等因素對徑流的影響，定量估算不同土地利用類型的產水量，結果更為科學。同時，對北京山區、黃土高原和橫斷山區等山區產水或水源涵養評估，也體現了該模型具備參數調整靈活、空間表達性強等優點。

Yjx=1-AETxjPx×Px（1）

式中，Yjx為第j類土地利用/覆被上柵格單元x的年產水量（mm）;Px為柵格單元x的年平均降水量（mm）;AETxj為第j類土地利用/覆被類型上柵格單元x的實際年平均蒸散量（mm），通過公式（2）計算。

AETxjPx=1+ωxRxj1+ωxRxj+1/Rxj（2）

式中，Rxj為土地利用類型j上柵格單元x的干燥指數，無量綱，表示潛在蒸發量與降水量的比值，由公式（3）計算;

k為蒸散系數，由植被葉面積指數LAI計算獲得，見公式（4）;ET0為潛在蒸散量（mm/d），計算公式如式（5）。

ωx為修正植被年可利用水量與降水量的比值，無量綱，由公式（6）計算。

Rxj=k×ET0Px（3）

k=min（ LAI3）（4）

ET0=0.001 3×0.408×RA×（Tavg+17）×（TD-0.012 3P）0.76（5）

ωx=ZAWCxPx（6）

式中，RA為太陽大氣頂層輻射[MJ/（m2·d）];Tavg是日最高溫均值和日最低溫均值的平均值（℃）;TD是日最高溫均值和日最低溫均值的差值（℃）。太陽大氣頂層輻射用氣象站太陽平均總輻射除以50%計算獲得[16]。Z為zhang系數;AWCx為植被有效可利用水，由土壤深度和理化性質決定，計算公式如下：

AWCx=min（maxSDx，RDx）×PAWCx（7）

PAWCx=54.509-0.132Sand-0.003（Sand）2-0.055Silt-0.006（Silt）2-0.738Clay+0.007（Clay）2-2.688OM+0.501（OM）2（8）

式中，maxSDx為最大土壤深度;RDx為根系深度;PAWCx為植被可利用水含量;Sand為土壤砂粒含量（%）;Silt為土壤粉粒含量（%）;Clay為土壤黏粒含量（%）;OM為土壤有機質含量（%）。

1.3 水源涵養模型 該研究選取研究區均為山地，必須考慮地形對水源涵養功能的影響，因此在產水量的基礎上，再利用地形指數、土壤飽和導水率和流速系數對產水量進行修正，得到研究區的水源涵養量。

WR=min（ 249V）×min（ 0.9×TI3）×min（ Ks300）×Yjx（9）

TI=lg（DASD×PS）（10）

lnKs=20.62-0.96×ln（C）-0.66×ln（S）-0.46×ln（OC）-8.43×BD（11）

式中，WR為水源涵養量（mm）;TI為地形指數，無量綱;Ks為土壤飽和導水率（cm/d），該研究選用土壤傳遞函數間接計算得出[17]。V為流速系數，數據無量綱;Yjx為產水量。DA為集水區柵格數量，無量綱;SD為土壤深度（mm）;PS為百分比坡度。C為土壤黏粒含量（%）;S為土壤砂粒含量（%）;OC為土壤有機碳含量（%）;BD為土壤容重（g/cm3）。

1.4 水源涵養功能影響因素 為定量分析影響因素與不同水源涵養功能等級之間的關系，利用地圖代數原理，利用公式（12）對研究區的土地利用類型、坡度與水源涵養功能之間的相關變化信息進行提取。

Cij=Rij×100+Sij（12）

式中，Cij為各個影響因素圖層與水源涵養功能圖層進行空間疊加的綜合圖層，Rij為各個影響因素圖層，Sij為水源涵養功能圖層。通過對2個不同的專題圖層進行空間疊加運算，計算結果既能在空間上反映不同專題圖層的等級變化，同時又能得到不同圖層之間的相互轉換關系。

1.5 數據來源及處理

1.5.1? 產水量數據。

水源涵養功能受到多種因素的影響，其計算過程必然關系到多種因素的信息提取與集成利用。利用InVEST模型進行產水量評估需要的輸入數據包括土地利用/覆被圖、年平均降水量圖、年平均潛在蒸散量圖、土壤根系深度圖、植物可利用水含量（PAWC）、集水區以及生物物理參數表。

（1）年平均降水量。降水數據來源于中國氣象數據網（https：//data.cma.cn/）的中國地面降水月值0.5°×0.5°格點數據集，為對應研究區開發時段，也為了避免單年數據的低代表性，該研究選取2000—2020年的降水數據，運用Python 3.8.3提取項目區周圍格點的月降水量，在ArcGIS中以地形作為協變量進行空間插值獲得研究區多年平均降水量柵格圖層。

（2）年平均潛在蒸散量。太陽輻射數據來源于WorldClim（https：//www.worldclim.org/），氣溫數據來源于中國氣象數據網（https：//data.cma.cn/）的中國地面累年值日值數據集（1981—2010年）的累年平均日最高氣溫和日最低氣溫。估算潛在蒸散量的方法有很多，該研究使用InVEST模型推薦的Modified-Hargreaves。在ArcGIS中進行克里金插值、重采樣和裁剪等處理，并運用柵格計算器對太陽輻射、平均日最高氣溫和最低氣溫計算得到多年平均潛在蒸散量柵格圖。

（3）土地利用類型。土地利用類型由資源三號衛星遙感影像解譯得到。根據數據分析要求將土地利用類型劃分為7個一級地類、14個二級地類，分別是耕地（水田、旱地）、林地（有林地、灌木林地、其他林地）、園地、草地、水域（水庫坑塘、灘涂）、建設用地（農村居民點用地、城鎮用地、其他建設用地）、未利用地（裸地、荒草地）。同時，利用Google Earth 高分辨率影像（2020年），對土地利用解譯結果和植被覆蓋情況進行野外選點驗證并訪談當地居民和工作人員，對土地利用數據進行現場復核。

（4）土壤深度。土壤數據來源于北京大學城市與環境學院地理數據平臺（http：//geodata.pku.edu.cn）的中國土壤數據集。對全國土壤數據進行空間插值后得到云南省土壤根深圖，再運用項目區范圍數據對云南省土壤根深圖進行裁剪得到研究區的土壤根深柵格數據。

（5）植被可利用水含量。結合中國土壤數據集中土壤質地的劃分，將砂粒、粉粒、黏粒、有機質等質量分數導入土壤有效含水量的經驗公式中，在ArcGIS中利用字段算計器求出PAWC，由空間分析和轉換工具處理，得到植被可利用含水量圖。

（6）Z系數。Z系數是表征降水季節性特征的一個常數，其值在1～30，降水主要集中在冬季時，其值越大，降水主要集中在夏季或季節分布較均勻時，其值越小;但在降水量相等的地區，降水事件越頻繁，Z值越大。產水量的大小受到Z系數影響，在熱帶地區Z值取4;在季風氣候區Z值取1;在溫帶海洋性氣候區Z值取9[18];該研究區均屬于季風氣候，故Z值取1。

（7）InVEST模型參數。蒸散系數和流速系數數據是參考相關文獻與結合InVEST用戶指南推薦使用的參考數據以及根據研究區地表植被覆蓋實際情況確定。LULC_veg表示采用何種實際蒸散量公式的信息，含有植被土地利用/覆被賦值為 其他賦值為0（表2）。

1.5.2 水源涵養數據來源及處理。

（1）地形指數。通過項目區資料獲得2012年等高線、高程點、項目區范圍，通過ArcGIS 10.2創建TIN后轉為柵格，得到開發前的DEM;運用2020年資源三號衛星影像立體像對，采用ENVI 5.3提取DEM。通過ArcGIS 10.2的水文分析工具進行填洼、流向、流量等計算步驟，獲得流量柵格，根據項目區情況確定閾值，利用Con條件函數篩選主要徑流，經過盆域分析得到流域圖。再采用面轉柵格工具得到流域單元集水區柵格的數量;通過 ArcGIS中3D分析提取DEM中的百分比坡度，此時會得到坡度為0的柵格，為了不影響后面的計算，將百分比坡度為0的點替換成0.000 0 由于該值遠小于地形指數的計算位數，誤差可忽略不計。最后綜合百分比坡度、流域單元集水區柵格數量以及土壤深度可得到流域的地形指數。

（2）土壤飽和導水率。基于土壤黏粒、砂粒、有機碳、土壤容重含量，運用土壤傳遞函數間接測算出土壤飽和導水率。

（3）流速系數。根據USDA-NRCS提供的國家工程手冊上的流速-坡度-景觀表格，乘以1 000獲得。

2 結果與分析

2.1 產水量時空變化 總體來看（圖2），開發后項目區平均產水量和平均產水深度均低于開發前，城市建設區的產水總量高于工業區，城市建設區的產水量下降程度高于工業區。從平均產水量來看，工業區的平均產水量從高到低依次為中安—后所項目區、北古城項目區、觀音山項目區、高家屯—竇家沖項目區、長坡項目區。城市建設區的平均產水量從高到低依次為光華項目區、景大項目區、沙朗項目區、塔寶山—大湖西項目區、下和—上登項目區，表明平均產水量主要受土地利用類型影響。從產水深度的變化來看，工業區減少量從高到低依次為觀音山項目區、北古城項目區、高家屯—竇家沖項目區、長坡項目區、中安—后所項目區。城市建設區減少量從高到低依次為塔寶山—大湖西項目區、沙朗項目區、下和—上登項目區、光華項目區、景大項目區。從空間上看（圖3～4），產水深度變化受土地利用類型影響較大，耕地、園地、林地、草地、未利用地轉變為建設用地后產水量下降，水域、荒草地轉變為其他地類后產水深度增加。產水量大小與項目區面積具有較強的相關性，產水量變化與降水量變化具有一致性。

2.2 水源涵養功能時空變化

總體來看，大部分項目區水源涵養總量均低于開發前，僅光華項目區水源涵養總量高于開發前，城市建設區的水源涵養總量高于工業區，工業區水源涵養量下降程度高于城市建設區。從水源涵養量變化來看，在工業區中，減少量從高到低依次是北古城項目區、觀音山項目區、高家屯—竇家沖項目區、中安—后所項目區、長坡項目區。在城市建設區中，變化量從高到低依次是下和—上登項目區、景大項目區、沙朗項目區、塔寶山—大湖西項目區、光華項目區。從平均水源涵養能力變化（圖5）來看，工業區減少量從高到低依次是高家屯—竇家沖項目區、中安—后所項目區、北古城項目區、觀音山項目區、長坡項目區。城市建設區減少量從高到低依次是下和—上登項目區、沙朗項目區、塔寶山—大湖西項目區、景大項目區、光華項目區?？臻g趨勢表明（圖6～7），水源涵養能力變化與土地利用類型變化具有一致性。

2.3 影響因素分析

2.3.1 土地利用類型。

從土地利用結構來看（表3、圖8～9），開發后各項目區的耕地減少，大部分項目區園地減少，建設用地面積均有不同程度的擴大，增長量由大到小依次為下和—上登項目區、北古城項目區、高家屯—竇家沖項目區、觀音山項目區、塔寶山—大湖西項目區、中安—后所項目區、光華項目區、沙朗項目區、景大項目區、長坡項目區。

不同土地利用類型的水源涵養能力存在差異。從時間上看，無論開發前還是開發后，林地的水源涵養能力最高，耕地、園地、草地的水源涵養能力在不同的項目區有差異。2012—2020年項目區耕地的水源涵養能力均呈下降趨勢，70%研究區園地和林地的水源涵養能力呈下降趨勢，50%研究區草地的水源涵養能力呈下降趨勢。

2.3.2 植被類型。

在未受人為因素影響的區域，存在植被的自然演替。從表3可以看到，長坡項目區、北古城項目區、中安—后所項目區、塔寶山—大湖西項目區、光華項目區的林地面積均有不同程度的增加。從土地利用類型（圖10～11）可以看出，長坡項目區的部分旱地轉為其他林地;北古城項目區中部分裸地轉變為草地和其他林地，部分草地轉變為有林地;中安—后所項目區中區塊二中部的裸地轉變成草地或其他林地，灌木林地轉變為有林地或其他林地;塔寶山—大湖西項目區中大湖西區塊的草地轉變成其他林地。

2.3.3 地形。時間趨勢表明（圖12～14），2012—2020年各項目區的平均高程均呈現下降趨勢。在工業區中，減少量從高到低依次是中安—后所項目區、長坡項目區、北古城項目區、高家屯—竇家沖項目區、觀音山工業區。城市建設區減少量從高到低依次是塔寶山—大湖西項目區、下和—上登項目區、沙朗項目區、景大項目區、光華項目區。城市建設區高程變化與水源涵養總量變化排序基本一致。

在工業區中，中安—后所項目區的高程變化最大，經實地調研得知該項目區入駐產業為煤礦、金屬加工等重工業，需要大范圍平整土地，因此采用“削峰填谷”的開發方式，對地形的改造程度較高。在城市建設區中，塔寶山—大湖西項目區的高程變化最大，經實地調研得知入駐產業為倉儲物流、汽配廠等，同樣需要大范圍平整土地，且塔寶山—大湖西項目區建設用地面積增加了183.18 hm 為滿足土地利用布局的需要，對地形地貌的改造程度較大。

2.3.4 降水量。

從空間上看，位于滇中的長坡、北古城、觀音山、沙朗項目區降水量較為接近，2020年平均降水量810.80～904.17 mm;位于滇西北下和—上登、塔寶山—大湖西項目區，2020年平均降水量792.48～824.13 mm;滇東北包括高家屯—竇家沖、中安—后所項目區，2020年平均降水量900.77～1 043.31 mm;滇南景大、光華項目區，2020年平均降水量1 307.84～1 447.18 mm。這體現了云南省降水空間分布不均。

降水量是影響水源涵養功能的重要氣候因子，時間趨勢表明（圖15），2012—2020年各項目區的降水量均呈減少狀態。下降量從大到小依次是下和—上登項目區、沙朗項目區、塔寶山—大湖西項目區、觀音山項目區、光華項目區、景大項目區、長坡項目區、中安—后所項目區、高家屯—竇家沖項目區、北古城項目區。降水量與產水量的變化趨勢有一致性。

3 結論與討論

3.1 結論

基于InVEST模型對云南省2012和2020年10個典型項目區水源涵養進行測算，得到以下結論：建設開發后項目區平均產水量均低于開發前;除光華項目區外，建設開發后項目區水源涵養總量均低于開發前。開發后平均水源涵養量從低到高依次是塔寶山—大湖西項目區、沙朗項目區、下和—上登項目區、高家屯—竇家沖工業園、北古城項目區、長坡項目區、中安—后所項目區、光華項目區、觀音山項目區、景大項目區。

該研究發現，造成低丘緩坡項目區建設開發前后水源涵養量變化的原因主要是土地利用結構、植被、地形、徑流、降水、土壤等自然和人為因素。由于土壤9年內變化較小，且受到數據的制約，采用的是一期的數據，故將其剔除。土地利用結構變化顯示，不同地類的水源涵養能力有差異，其中林地的水源涵養能力最高，耕地、園地、草地的水源涵養能力在不同項目區有差異。植被類型表明，在未受人為因素影響的區域，存在植被的自然演替。裸地轉變為草地，草地轉變為其他林地，灌木林地轉變為有林地等，其水源涵養能力也得到一定程度的提升。地形、植被、降水等因素共同影響地表徑流的流速和流量，從而影響水源涵養功能。

3.2 討論

結合實地調研案例，從積極和消極方面分析低丘緩坡建設開發對生態環境產生的影響。①低丘緩坡開發改變原來的土地利用結構，擾動和破壞原有的地表環境，導致生物多樣性的減少。②基礎設施建設過程中需要大量的挖方填方工作，破壞原地表植被，擾動土體，使土壤松動、裸露和堆填，若未及時采取水土保持措施，容易導致地質災害的發生。③人類活動必然會產生污染物排放，如未能恰當處理，會對周邊區域的水源、耕地等生態敏感區產生不良影響，造成生態破壞、環境污染等問題[19]。

根據低丘緩坡建設開發可能造成的生態破壞、基于水源涵養功能提出相應的對策建議。目前低丘緩坡建設開發模式主要有以下3種：削峰填谷式、分臺式、分臺式與緩坡式相結合。該研究依據保護生態優先、因地制宜的原則，提出以下優化措施：首先，科學選址和土地規劃，根據山地立體特征，依據山體不同位置和坡度，布局不同的土地利用類型，宜建則建，宜林則林，減少對原地表的干擾和破壞;其次，調整土地利用結構，保留具有較強生態源的生態用地，加強綠化工程設施建設，提升土地生態功能;第三，落實水源涵養規劃，完善水源涵養生態補償機制，推進低丘緩坡開發利用可持續發展[20-21]。

參考文獻

[1]

周佳雯，高吉喜，高志球，等.森林生態系統水源涵養服務功能解析[J].生態學報，2018，38（5）：1679-1686.

[2] 陳姍姍，劉康，包玉斌，等.商洛市水源涵養服務功能空間格局與影響因素[J].地理科學，2016，36（10）：1546-1554.

[3] 郭麗潔，尹小君，茍貞珍，等.基于InVEST模型的阿克蘇河流域產水量評估及環境因素影響研究[J].石河子大學學報（自然科學版），2020，38（2）：216-224.

[4] 余新曉，周彬，呂錫芝，等.基于InVEST模型的北京山區森林水源涵養功能評估[J].林業科學，201 48（10）：1-5.

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

[6] 寧亞洲，張福平，馮起，等.秦嶺水源涵養功能時空變化及其影響因素[J].生態學雜志，2020，39（9）：3080-3091.

[7] 竇苗，孫建國，陳海鵬.基于InVEST模型的橫斷山區產水量模擬[J].安徽農業科學，2017，45（21）：54-58.

[8] 苑躍，張亮，崔林林.若爾蓋高原生態系統水源涵養功能時空變化特征[J].生態學雜志，2020，39（8）：2713-2723.

[9] REDHEAD J W，STRATFORD C，SHARPS K，et al Empirical validation of the InVEST water yield ecosystem service model at a national scale[J].Science of the total environment， 2016， 569/570：1418-1426.

[10] JULIAN J P，GARDNER R H.Land cover effects on runoff patterns in eastern Piedmont （USA） watersheds[J].Hydrological processes，2014，28（3）：1525-1538.

[11] SCORDO F，LAVENDER T M，SEITZ C，et al.Modeling water yield：Assessing the role of site and region-specific attributes in determining model performance of the InVEST seasonal water yield model[J].Water， 2018，10：1-42.

[12] 陳武強，羅鮮美，季金娥，等.昆明西華濕地公園生態系統服務價值評價研究[J].價值工程，2019，38（10）：167-170.

[13] 趙元藩，宋東華，溫慶忠，等.玉溪市森林生態系統服務功能價值評估[J].林業調查規劃，201 36（1）：12-18，25.

[14] 楊芳，支玲，郭小年，等.云南省玉龍縣天保工程森林生態系統服務功能價值評估[J].西南林業大學學報（社會科學），2017，1（1）：69-76.

[15] 杜軍凱.考慮垂直地帶性的山區分布式水文模擬與應用[D].北京：中國水利水電科學研究院，2019.

[16] 吳丹，邵全琴，劉紀遠，等.中國草地生態系統水源涵養服務時空變化[J].水土保持研究，2016，23（5）：256-260.

[17] 陳麗，郝晉珉，陳愛琪，等.基于二元水循環的黃淮海平原耕地水源涵養功能研究[J].生態學報，2017，37（17）：5871-5881.

[18] ZHANG C Q，LI W H，ZHANG B，et al.Water yield of Xitiaoxi River basin based on InVEST modeling[J].Journal of resources and ecology， 201 3（1）：50-54.

[19] 宋夢意，王盛毅，衛樂樂.浙江省低丘緩坡土地資源開發利用研究：以工業建設用地為視角[J].河南教育學院學報（自然科學版），201 21（3）：52-54.

[20] 王驕.生態視角下城鎮上山實施回顧與規劃對策[D].重慶：重慶大學，2014.

[21] 張健.以提升生態功能為導向的城市綠道系統規劃方法研究：以成都天府綠道為例[J].西部人居環境學刊，2019，34（6）：73-78.