基于InVEST模型農(nóng)牧交錯帶生態(tài)系統(tǒng)價值研究

摘要 ［目的］量化研究農(nóng)牧交錯帶的生態(tài)系統(tǒng)服務(wù)功能及其價值。［方法］以典型農(nóng)牧交錯帶內(nèi)蒙古自治區(qū)興安盟科爾沁右翼中旗為例，基于InVEST模型運行原理、研究區(qū)特點和數(shù)據(jù)可獲得程度，從土壤保持、水源供給和水質(zhì)凈化3個主要方面對科爾沁右旗2000、2010和2019年3個時段的生態(tài)系統(tǒng)服務(wù)功能價值進行評估。［結(jié)果］科爾沁右翼中旗生態(tài)系統(tǒng)服務(wù)功能價值較大，2000、2010和2019年的生態(tài)系統(tǒng)土壤保持服務(wù)功能價值分別為30.54億、31.88億和32.79億元，水源供給服務(wù)功能價值分別為100.69億、105.19億和87.40億元，水質(zhì)凈化服務(wù)功能價值分別為1 685.54萬、1 584.43萬和1 698.02萬元。［結(jié)論］基于InVEST模型對生態(tài)系統(tǒng)服務(wù)功能價值評估，不僅為北方農(nóng)牧交錯帶生態(tài)系統(tǒng)服務(wù)功能重要性區(qū)域和綜合生態(tài)區(qū)劃分奠定了研究基礎(chǔ)，也彌補了InVEST模型在我國北方農(nóng)牧交錯帶生態(tài)系統(tǒng)服務(wù)研究內(nèi)容和方法應(yīng)用的不足。

關(guān)鍵詞 農(nóng)牧交錯帶；生態(tài)系統(tǒng)服務(wù)；價值評估；InVEST模型

中圖分類號 X 171 文獻標識碼 A 文章編號 0517-6611（2024）23-0053-06

doi：10.3969/j.issn.0517-6611.2024.23.012

Research on the Ecosystem Value of the Agriculture-pastoral Ecotone Based on the InVEST Model—Taking Horqin Right Middle Banner as an Example

DING Qing-po1，JIN Wen-juan2，SUN Tian-qi3

（1.Chifeng Natural Resources Reserve and Consolidation Center，Chifeng，Inner Mongolia 024000;2.College of Land and Environment，Shenyang Agricultural University，Shenyang，Liaoning 110000;3.Anshan Natural Resources Bureau No.3 Branch，Anshan，Liaoning 114000）

Abstract ［Objective］To quantitatively study the ecosystem service functions and their values in the agriculture-pastoral ecotone.［Method］Taking the typical agricultural-pastoral ecotone in Horqin Right Middle Banner，Xing′an League，Inner Mongolia Autonomous Region as an example，based on the operation principle of InVEST model，characteristics of the area and data availability，this paper evaluated the value of ecosystem services in 2000，2010 and 2019 from the three main aspects of soil maintenance，water supply and water purification.［Result］The ecosystem service value was large.In 2000，2010 and 2019，the value of soil maintenance service was 3.054 billion，3.188 billion and 3.279 billion，respectively; the value of water supply service was 10.069 billion，10.519 billion and 8.74 billion，respectively; and the value of water purification service was 16.855，4 million，15.844，300 million and 16.980，200 million，respectively.［Conclusion］The value evaluation of ecosystem service function lays a foundation for the division of the important region and comprehensive ecological area of the ecosystem service function in the northern ecosystem，and makes up for the deficiency of the research content and method application of InVEST model in the ecosystem service in northern China.

Key words Agriculture-pasture ecotone;Ecosystem service;Value evaluation;InVEST model

作者簡介 丁青坡（1979—），女，內(nèi)蒙古赤峰人，高級工程師，碩士，從事測繪工程、土地利用相關(guān)研究。*通信作者，博士研究生，研究方向：國土空間生態(tài)修復(fù)。

收稿日期 2024-01-04

我國北方農(nóng)牧交錯帶是一個復(fù)合生態(tài)系統(tǒng)，也是典型的生態(tài)環(huán)境脆弱區(qū)，具有多種農(nóng)牧業(yè)生產(chǎn)方式長期并存的過渡性自然特征。北方農(nóng)牧交錯帶區(qū)域氣候干旱、有效降水少、風(fēng)力侵蝕嚴重，區(qū)域性自然生產(chǎn)力較低，水土流失、水資源短缺及水體污染等一系列生態(tài)問題頻發(fā)，生態(tài)系統(tǒng)服務(wù)功能瀕臨失調(diào)［1］。開展北方農(nóng)牧交錯帶的生態(tài)系統(tǒng)服務(wù)功能價值評估，深入探討多時段、不同土地利用方式下各種生態(tài)系統(tǒng)服務(wù)功能的空間動態(tài)變化特征及其價值，有助于劃分區(qū)域綜合性生態(tài)功能區(qū)，為政府和相關(guān)主管部門科學(xué)合理制定生態(tài)修復(fù)保護政策和決策規(guī)劃提供參考依據(jù)。

當(dāng)前關(guān)于北方農(nóng)牧交錯帶的研究較多集中于土地利用方式演變及其影響［2］、生態(tài)脆弱性［3］、生態(tài)環(huán)境問題［4-5］等方面，缺乏對北方農(nóng)牧交錯帶典型區(qū)域的生態(tài)系統(tǒng)服務(wù)功能及其價值量化研究，尤其在利用InVEST模型評估農(nóng)牧交錯帶的生態(tài)系統(tǒng)服務(wù)功能價值方面。筆者選擇對生態(tài)系統(tǒng)影響較大的土壤保持、水源供給和水質(zhì)凈化3種生態(tài)系統(tǒng)服務(wù)功能，基于InVEST模型，采用經(jīng)驗價值量化法對北方農(nóng)牧交錯帶典型縣域科爾沁右翼中旗的生態(tài)系統(tǒng)服務(wù)功能價值進行評估，以期為農(nóng)牧交錯帶生態(tài)環(huán)境的保護和管理提供理論基礎(chǔ)。

1 資料與方法

1.1 研究區(qū)概況 科爾沁右翼中旗隸屬于內(nèi)蒙古自治區(qū)興安盟市（119°34′～122°18′E，44°15′～46°41′N），地處科爾沁草原中部，是科爾沁沙地向大興安嶺、松嫩平原、松遼平原的過渡地帶。全域總面積15 613 km2，下轄12個蘇木鎮(zhèn)，173個嘎查。西北部發(fā)展畜牧業(yè)，中部至東南進行農(nóng)作物種植，是我國北方農(nóng)牧交錯帶東段一個典型的以農(nóng)牧業(yè)發(fā)展為主的縣域［6］。區(qū)域?qū)僦袦貛О敫珊荡箨懶约撅L(fēng)氣候區(qū)，年均溫3.8～7.2 ℃，無霜期128～289 d，年降水量360～400 mm，年內(nèi)降水分配不均，集中于6—8月，約占全年降水量的47%。境域內(nèi)降水情況南北部存在分異，自北向南逐漸減少。受局部性山洪、干旱、霜凍、沙塵暴和寒潮等的影響，區(qū)域內(nèi)部分地帶生態(tài)環(huán)境尤為敏感脆弱。

1.2 研究方法

1.2.1 InVEST模型。InVEST模型基于GIS平臺開發(fā)，由Natural Capital項目支持，美國斯坦福大學(xué)、世界自然基金和大自然保護協(xié)會等共同研發(fā)，用于評估生態(tài)系統(tǒng)各服務(wù)的空間特征和價值量，衡量維持生態(tài)過程穩(wěn)定性以及人類生存的物質(zhì)需求所提供的各種服務(wù)和效益。InVEST模型的運算過程較靈活，能有針對性地解決主要生態(tài)問題，綜合考慮生態(tài)系統(tǒng)整體循環(huán)過程［7］。基于模擬原理，在模塊的參數(shù)設(shè)置中輸入DEM、降雨、土壤、土地及遙感等數(shù)據(jù)進行測算，從而得出所需數(shù)據(jù)。

1.2.1.1 土壤保持模塊。InVEST模型的土壤保持模塊改進了通用土壤流失方程，綜合計算柵格自身產(chǎn)沙量和對上游柵格產(chǎn)沙的截留量［8］。公式如下：

RKLSx=Rx×Kx×Lx×Sx（1）

USLEx=Rx×Kx×Lx×Sx×Cx×Px（2）

式中：RKLSx是柵格x的土壤潛在侵蝕量（t）；USLEx是實際侵蝕量（t）；Rx為降雨侵蝕力［MJ·mm/（hm2·a）］；Kx為土壤可蝕性［t·a/（MJ·mm）］；Lx是坡長（m）;Sx是坡度（°）；Cx是植被覆蓋度（%）；Px是土壤流失量（t）。

土壤保持量由以下公式得到：

SEDRx=SExx-1y=1USLExx-1s=y+1（1-SEx）（3）

SEDRETxD=RKLSx-USLEx+SEDRx（4）

式中：SEDRx為柵格x截留上游產(chǎn)沙量（t）；SEx為柵格x的截留率；SEDRETxD是柵格x的土壤持留量（t）。

1.2.1.2 水源供給模塊。利用InVEST模型中的產(chǎn)水量模塊，基于水量平衡原理，根據(jù)柵格單元的降雨量與實際蒸散發(fā)量之差計算［9］，公式如下：

Yxk=1-AETxkPx×Px（5）

式中：Yxk為柵格x地類k的平均產(chǎn)水深度（mm）；Px為柵格x的年降雨量（mm）；AETxk為平均實際蒸散發(fā)量（mm）。

1.2.1.3 水質(zhì)凈化模塊。InVEST水質(zhì)凈化模塊主要評估植被、土壤及徑流過程中對N、P物質(zhì)的凈化狀況［10］。N、P保持量越大輸出量越小，水質(zhì)凈化服務(wù)越好［11］。模型只考慮面源污染，其他污染源忽略不計［12］，公式如下：

xexptot=xexpi（6）

xexpi=loadsurf，i×NDRsurf，i+loadsubs，i×NDRsubs，i（7）

式中：xexptot為N、P輸出總量（kg/a）；xexpi為研究區(qū)每個柵格單元內(nèi)的N、P輸出量（kg/a）；loadsurf，i和loadsubs，i分別是地表和地下N、P負荷量［kg/（hm2·a）］；NDRsurf，i和NDRsubs，i分別是地表和地下N、P的傳遞速率（kg/a）。

1.2.2 價值評估法。在確定生態(tài)系統(tǒng)服務(wù)模塊后，采用價值評估法，利用貨幣單位衡量和表征生態(tài)系統(tǒng)服務(wù)功能的價值［13］。該研究主要采用市場價值法、影子價格法和替代成本法評估3種生態(tài)系統(tǒng)服務(wù)功能的價值。其中，土壤保持價值包括保持土壤肥力價值和固土價值，保持土壤肥力價值采用市場價值法計算，固土價值采用替代成本法和影子價格法估算。水源供給服務(wù)價值采用市場價格法根據(jù)居民生活用水價格估算。水質(zhì)凈化服務(wù)價值采用替代成本法估算。

1.2.3 InVEST模型中模塊參數(shù)選擇和價值計算。

1.2.3.1 土壤保持模塊參數(shù)選擇和計算。

（1）降雨侵蝕力（R）。它反映出由于降雨所引起土壤流失的潛在能力，表征由降雨侵蝕條件對土壤潛在流失量的影響。采用月尺度推算法，結(jié)合各氣象站點的月、年降雨量數(shù)據(jù)，利用公式（8）計算降雨侵蝕力（R），并在ArcGIS中利用反距離插值法［14］生成降雨侵蝕力空間分布圖（圖1a）。

R=12i=11.735×101.5lgp2ip-0.818 8×17.02（8）

式中：R為降雨侵蝕力，乘以系數(shù)17.02后轉(zhuǎn)換成國際單位制［MJ·mm/（hm2·a）］；Pi、P分別為月、年降雨量（mm）。

（2）土壤可蝕性（K）。土壤可蝕性表示侵蝕作用下土壤性質(zhì)的敏感性，一般采用EPIC模型［15］進行土壤可蝕性計算（見公式9～11），并根據(jù)研究區(qū)土壤可蝕性特征修正，最后利用反距離插值法形成空間分布圖（圖1b）。

K=（-0.013 83+0.515 75Kepic）×0.131 8 （9）

Kepic=0.2+0.3exp［-0.025 6SAN（1-SIL100）］SILCLA+SIL0.3×1-0.25CC+exp（3.72-2.95C）×1-0.7SN1SN1+exp（22.9SN1-5.51）（10）

SN1=1-SAN100（11）

式中：K為土壤可蝕性［（t·a）/（MJ·mm）］；SAN為土壤砂粒含量（%）；SIL為土壤粉粒含量（%）；CLA為土壤黏粒含量（%）；C為土壤有機碳含量（%）。

（3）植被覆蓋度（C）和土壤流失量（P）。植被覆蓋度與植被指數(shù)、覆蓋率等相關(guān)計算公式如下：

fc=NDVI-NDVIminNDVImax-NDVImin（12）

C=1fc=0C=0.650 8-0.343 61gfc0<fc≤78.3%C=0fc>78.3%（13）

式中：C為植被覆蓋度（%）；fc為植被覆蓋率（%）；NDVI為歸一化植被指數(shù)；NDVImin、NDVImax分別為植被指數(shù)的最小值、最大值。模型的其他參數(shù)包括流域流量累積閾值、K值、IC0值和Max SDR Value均按照模型默認參數(shù)設(shè)置，分別為1 000、2、5和0.8。參考相關(guān)研究結(jié)果［7，16］，該研究的植被覆蓋度和土壤流失量的取值如表1所示。

（4）土壤保持服務(wù)價值。土壤保持價值包括保持土壤肥力價值和固土價值。其中保持土壤肥力價值是指土壤侵蝕導(dǎo)致土壤化學(xué)肥料中施加的N、P、K等元素流失，N、P、K等土壤養(yǎng)分的價值即保持土壤養(yǎng)分價值是由土壤保持量、化肥養(yǎng)分含量以及化肥的市場價格計算得到；保持土壤有機質(zhì)價值按照等量薪材的市場價格計算，具體如下：

Ea=E1+E2 （14）

E1=Q×（mi×ni×pi）（15）

E2=Q×Di×Ps×S（16）

式中：Ea為保持土壤肥力價值（元）；E1為保持土壤養(yǎng)分價值（元）；E2為保持土壤有機質(zhì)價值（元）；Q為土壤保持量（t）；mi為土壤養(yǎng)分中N、P、K的含量，參考相關(guān)研究［17-19］，各元素的含量分別為0.068 0%、0.014 3%、1.826 0%；ni為土壤中各養(yǎng)分含量折算為化肥的系數(shù)，其中y3hAaOOHEcLhPT8ohg4evRktGeR1gSkH6CoyQnzxd3M=堿解氮折算為CH4N2O（尿素）的系數(shù)為2.164，速效磷折算為CaP2H4O8（過磷酸鈣）的系數(shù)為4.056，速效鉀折算為KCl（氯化鉀）的系數(shù)為1.923；pi為3種化肥的平均市場價格，其中CH4N2O（尿素）為1 833 元/t，CaP2H4O8（過磷酸鈣）為566 元/t，KCl（氯化鉀）為2 078 元/t；Di為研究區(qū)土壤有機碳含量，取值1.613%［20］；Ps為薪材平均市場價格，取200元/m3；S為薪材折算為土壤有機質(zhì)的系數(shù)，取0.5。

固土價值是指采用替代成本法和影子價格法估算減少土地廢棄的經(jīng)濟價值。此外，還包括由于泥沙淤積導(dǎo)致水庫蓄水量減少而造成的經(jīng)濟損失［21］。計算公式如下：

Eb=E3+E4（17）

E3=Q×B/（0.6×10 000ρ）（18）

E4=0.24×Q×Vg/ρ（19）

式中：Eb為固土價值（元）；E3為減少土地廢棄的價值（元）；E4為減少泥沙淤積的價值（元）；Q為土壤保持量（t）；B為土地年均生產(chǎn)價值，取科爾沁右翼中旗單位面積農(nóng)林牧漁總產(chǎn)值作為機會成本，根據(jù)《科爾沁右翼中旗統(tǒng)計年鑒》得出其值為2 012.291元/hm2；0.6為研究區(qū)的平均土層厚度（m）；Vg為水庫建設(shè)成本，取0.67元/m3；ρ為土壤容重，取值1.29 t/m3；0.24是土壤侵蝕泥沙淤積率。

1.2.3.2 水源供給模塊參數(shù)選擇和計算。

（1）降雨量（P）。根據(jù)模塊需要，月降雨量根據(jù)各氣象站點的逐年累月平均降雨量取平均值。年降水量在ArcGIS中空間插值生成1 km降雨量空間分布圖（圖2）。

（2）潛在蒸散量（ET0）。潛在蒸散量是指在降水循環(huán)穩(wěn)定狀態(tài)下由于土壤和植物蒸散作用所減少的水量，該研究利用研究區(qū)內(nèi)及周圍的7個氣象站點的溫度、太陽輻射等數(shù)據(jù)，采用Modified-Hargreaves法［22］計算。公式如下：

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

式中：ET0為潛在蒸散量（mm）；RA為太陽大氣頂層輻射［MJ/（m2·d）］，需要根據(jù)所獲取的地面輻射數(shù)據(jù)進行散失量的折算；Tavg為日平均氣溫（℃）；TD為每日最高最低氣溫差值（℃）；P為月降水量（mm）。計算得研究區(qū)各氣象站點的潛在蒸散量，在ArcGIS中Kriging插值得研究區(qū)1 km潛在蒸散量空間分布圖（圖3a）。

（3）土地利用/覆被。將研究時段內(nèi)的LUCC數(shù)據(jù)在ArcGIS中利用空間分析工具按照科爾沁右翼中旗邊界進行提?。▓D3b）。

（4）植物可利用水（PAWC）。參考周文佐等［23］提出的基于土壤質(zhì)地AWC估算模型，數(shù)據(jù)處理后生成植物可利用水空間分布圖（圖3c）。模型如下：

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 （21）

式中：PAWCx為柵格x植被可利用水（%）；sand為土壤砂粒含量（%）；silt是土壤粉粒含量（%）；clay為土壤黏粒含量（%）；OM為土壤有機質(zhì)含量（%）。

（5）生物物理表。參考相關(guān)研究結(jié)果［24-26］、聯(lián)合國糧農(nóng)組織作物系數(shù)在線目錄及《InVEST模型用戶使用手冊》，確定植被蒸散系數(shù)、植物根系深度等，具體見表2。

（6）Zhang系數(shù)。InVEST模型中Zhang系數(shù)表示季節(jié)性降水特征，默認值為9.43。為使水源供給量的評估結(jié)果反映研究區(qū)的實際情況，經(jīng)過多次調(diào)整試驗，該研究的Zhang系數(shù)取值為3.1。

（7）水源供給服務(wù)價值。該研究采用市場價值法，根據(jù)科爾沁右翼中旗居民生活用水價格（2.7元/m3）進行生態(tài)系統(tǒng)水源供給服務(wù)的經(jīng)濟價值估算。計算公式如下：

Ec=G×M&nbsp; （22）

式中：Ec為生態(tài)系統(tǒng)水源供給服務(wù)價值（元）；G為水源供給總量（m3）；M為水價（元/m3）。

1.2.3.3 水質(zhì)凈化模塊參數(shù)選擇和處理。

（1）DEM。DEM可以很好地反映研究區(qū)的海拔與地形特征，該研究采用30 m×30 m精度進行處理并提取研究區(qū)流域空間。

（2）生物物理表。InVEST模型的水質(zhì)凈化模塊的生物物理表詳見表3。其中Load_N、Load_P代表N、P負荷系數(shù)，該研究修正該系數(shù)采用環(huán)境相似性原則，依據(jù)松花江流域［27-29］、松遼流域［30］等相關(guān)研究結(jié)果進行N、P負荷量的修正。

（3）其他參數(shù)。匯水積累量閾值取模型設(shè)定的默認值1 000［12］。地下徑流N、P最大保留率參考相關(guān)文獻取0.4［31］。最大保留率下N、P的輸送距離設(shè)為150 km［32］。Borselli k參數(shù)是水文連接性和養(yǎng)分輸送率之間的校準參數(shù)，默認值2。年降水量、年潛在蒸散發(fā)量、植物可利用水等圖層均與水源供給模塊相同。

（4）水質(zhì)凈化服務(wù)價值。采用替代成本法將生態(tài)系統(tǒng)水質(zhì)凈化服務(wù)價值用處理污染水質(zhì)中的N、P凈化成本替代［33-34］，按照以下公式計算生態(tài)系統(tǒng)水質(zhì)凈化服務(wù)的價值：

Ed=（T×Z）（23）

式中：Ed為生態(tài)系統(tǒng)水質(zhì)凈化價值（元）；T為N、P總持留量（kg）；Z為N、P的凈化單價，分別為1.50和2.50元/kg。

1.3 數(shù)據(jù)來源 該研究中的數(shù)據(jù)來源及獲取方式見表4。其中，土地利用數(shù)據(jù)為2000、2010和2019年的LUCC數(shù)據(jù)。數(shù)字高程模型（DEM）用于提取坡度、坡向、子流域等地理屬性。在模型模擬過程中，為使輸出結(jié)果更為準確，將數(shù)據(jù)重采樣至1 km×1 km。氣象數(shù)據(jù)，根據(jù)科爾沁右翼中旗地理位置，收集研究區(qū)周邊的阿爾山、霍林郭勒、索倫、扎魯特、通遼、白城和通榆共7個氣象站點數(shù)據(jù)。

2 結(jié)果與分析

2.1 生態(tài)系統(tǒng)土壤保持服務(wù)價值 經(jīng)計算，科爾沁右翼中旗2000、2010和2019年的生態(tài)系統(tǒng)土壤保持服務(wù)價值分別為30.54億、31.88億和32.79億元（表5）。其中，各時段的保持土壤肥力價值均占土壤保持服務(wù)價值的90%以上，而固土價值不足土壤保持服務(wù)價值的10%。由此說明，北方農(nóng)牧交錯帶生態(tài)系統(tǒng)的土壤保持服務(wù)主要體現(xiàn)在肥力保持方面。

2.2 生態(tài)系統(tǒng)水源供給價值 經(jīng)計算，科爾沁右翼中旗2000、2010和2019年生態(tài)系統(tǒng)水源供給服務(wù)價值分別為100.69億、105.19億和87.40億元（表6），呈現(xiàn)先上升后下降的趨勢。由此說明，北方農(nóng)牧交錯帶生態(tài)系統(tǒng)的水源供給服務(wù)受氣候影響，尤其是降雨量和蒸散發(fā)量。

2.3 生態(tài)系統(tǒng)水質(zhì)凈化價值 經(jīng)計算，科爾沁右翼中旗2000、2010和2019年生態(tài)系統(tǒng)水質(zhì)凈化服務(wù)價值分別為1 685.54萬、1 584.43萬和1 698.02萬元（表7）。在水質(zhì)凈化服務(wù)價值中，N保持價值約占總價值的95%，而P保持價值僅約占總價值的5%。由此說明，北方農(nóng)牧交錯帶生態(tài)系統(tǒng)的水質(zhì)凈化服務(wù)主要體現(xiàn)在凈化N方面。

3 結(jié)論與討論

該研究基于InVEST模型，采用經(jīng)驗價值量化法對北方農(nóng)牧交錯帶典型縣域科爾沁右翼中旗2000、2010和2019年的生態(tài)系統(tǒng)服務(wù)功能（土壤保持服務(wù)、水源供給服務(wù)、水質(zhì)凈化服務(wù)）價值進行定量計算并評估，結(jié)果表明：①生態(tài)系統(tǒng)中土壤保持服務(wù)、水源供給服務(wù)和水質(zhì)凈化服務(wù)功能時空變化具有一定規(guī)律性。20年來，土壤保持服務(wù)不斷提升；水源供給服務(wù)呈現(xiàn)先升后降趨勢，水源供給服務(wù)能力增強的區(qū)域面積整體大于衰減區(qū)域；水質(zhì)凈化服務(wù)總體穩(wěn)定。②生態(tài)系統(tǒng)中土壤保持服務(wù)、水源供給服務(wù)和水質(zhì)凈化服務(wù)價值較大，2000、2010和2019年的生態(tài)系統(tǒng)土壤保持服務(wù)價值分別為30.54億、31.88億和32.79億元，水源供給服務(wù)功能價值分別為100.69億、105.19億和87.40億元，水質(zhì)凈化服務(wù)功能價值分別為1 685.54萬、1 584.43萬和1 698.02萬元。

該研究基于InVEST模型定量評估了北方農(nóng)牧交錯帶東段典型縣域科爾沁右翼中旗的生態(tài)系統(tǒng)土壤保持、水源供給和水質(zhì)凈服務(wù)功能價值，InVEST模型較其他傳統(tǒng)方法用于生態(tài)系統(tǒng)服務(wù)功能進行評價具有一定優(yōu)勢，尤其是InVEST模型和GIS平臺的有機結(jié)合，對于收集處理降雨侵蝕力、土壤可蝕性、降雨量、蒸發(fā)量、土地利用、DEM等數(shù)據(jù)有明顯的優(yōu)勢。但由于以往研究結(jié)果不足，基礎(chǔ)數(shù)據(jù)相對較少且跨越年度較大，部分環(huán)境特征參數(shù)難以準確獲取，在一定程度上影響評估結(jié)果的精度。該研究僅從科爾沁右翼中旗生態(tài)系統(tǒng)土壤保持、水源供給和水質(zhì)凈化服務(wù)功能的價值評估，科爾沁右翼中旗所處的北方農(nóng)牧交錯帶的地帶特性和農(nóng)牧兼?zhèn)涞漠a(chǎn)業(yè)結(jié)構(gòu)特征致使其他生態(tài)系統(tǒng)服務(wù)能功能也表現(xiàn)出其區(qū)域自身特點，尤其是草原生態(tài)系統(tǒng)的生物多樣性保持、木材等原材料供給、美學(xué)景觀以及游牧民族的文化傳承等服務(wù)能力，應(yīng)繼而進行更為深入細致的研究。今后，基于InVEST模型數(shù)據(jù)收集處理以及模型局限性角度，為更好地服務(wù)于北方農(nóng)牧交錯帶的生態(tài)系統(tǒng)服務(wù)功能價值評估，應(yīng)設(shè)置長時段的實地調(diào)查和觀測活動，完善研究區(qū)自然資源環(huán)境數(shù)據(jù)庫的建立并實時更新，保證數(shù)據(jù)的典型性和準確性。

參考文獻

［1］ 何文清，陳源泉，高旺盛，等.農(nóng)牧交錯帶風(fēng)蝕沙化區(qū)農(nóng)業(yè)生態(tài)系統(tǒng)服務(wù)功能的經(jīng)濟價值評估［J］.生態(tài)學(xué)雜志，2004，23（3）：49-53.

［2］ 楊志鵬，許嘉巍，馮興華，等.基于InVEST模型的東北地區(qū)土地利用變化對生境的影響研究［J］.生態(tài)科學(xué)，2018，37（6）：139-147.

［3］ 黃越，程靜，王鵬.中國北方農(nóng)牧交錯區(qū)生態(tài)脆弱性時空演變格局與驅(qū)動因素：以鹽池縣為例［J］.干旱區(qū)地理，2021，44（4）：1175-1185.

［4］ 張生軍，楊改河，劉和林.北方農(nóng)牧交錯帶水資源與生態(tài)環(huán)境問題初探［J］.安徽農(nóng)業(yè)科學(xué)，2006，34（9）：1945-1947.

［5］ 程序.農(nóng)牧交錯帶研究中的現(xiàn)代生態(tài)學(xué)前沿問題［J］.資源科學(xué)，1999，21 （5）：1-8.

［6］ 陳秀梅.淺談科右中旗中小河流的治理措施［J］.內(nèi)蒙古水利，2015（6）：141-142.

［7］ 王森.延安市土地利用變化及其土壤保持功能效應(yīng)研究：基于InVEST模型的應(yīng)用［D］.楊陵：中國科學(xué)院大學(xué)（中國科學(xué)院教育部水土保持與生態(tài)環(huán)境研究中心），2018.

［8］ 吳哲，陳歆，劉貝貝，等.基于InVEST模型的海南島氮磷營養(yǎng)物負荷的風(fēng)險評估［J］.熱帶作物學(xué)報，2013，34（9）：1791-1797.

［9］ TALLIS H T，RICKETTS T，GUERRY A D，et al.InVEST 2.1 beta user’s guide［R］.Stanford：The Natural Capital Project，2011.

［10］ 潘翔.基于InVEST模型的生態(tài)系統(tǒng)服務(wù)權(quán)衡協(xié)同關(guān)系研究［D］.蘭州：西北師范大學(xué)，2018.

［11］ 王蓓，趙軍，胡秀芳.基于InVEST模型的黑河流域生態(tài)系統(tǒng)服務(wù)空間格局分析［J］.生態(tài)學(xué)雜志，2016，35（10）：2783-2792.

［12］ 吳瑞，劉桂環(huán)，文一惠.基于InVEST模型的官廳水庫流域產(chǎn)水和水質(zhì)凈化服務(wù)時空變化［J］.環(huán)境科學(xué)研究，2017，30（3）：406-414.

［13］ 劉世梁，安南南，王軍.土地整理對生態(tài)系統(tǒng)服務(wù)影響的評價研究進展［J］.中國生態(tài)農(nóng)業(yè)學(xué)報，2014，22（9）：1010-1019.

［14］ 胡勝，曹明明，劉琪，等.不同視角下InVEST模型的土壤保持功能對比［J］.地理研究，2014，33（12）：2393-2406.

［15］ 魏慧，趙文武，王晶.土壤可蝕性研究述評［J］.應(yīng)用生態(tài)學(xué)報，2017，28（8）：2749-2759.

［16］ 張恒瑋.基于InVEST模型的石羊河流域生態(tài)系統(tǒng)服務(wù)評估［D］.蘭州：西北師范大學(xué)，2016.

［17］ 焦燕，趙江紅，徐柱.內(nèi)蒙古農(nóng)牧交錯帶土地利用對土壤性質(zhì)的影響［J］.草地學(xué)報，2009，17（2）：234-238.

［18］ 程冬兵，蔡崇法，左長清.土壤侵蝕退化研究［J］.水土保持研究，2006（5）：252-254，258.

［19］ 楊芙蓉，楊恒智.內(nèi)蒙古巴彥淖爾地區(qū)土壤養(yǎng)分狀況及施肥現(xiàn)狀［J］.內(nèi)蒙古農(nóng)業(yè)科技，2008，36（2）：49.

［20］ 馬玉峰，王文明，賈寶全.半干旱農(nóng)牧交錯區(qū)土壤養(yǎng)分空間變異特征分析：以內(nèi)蒙古伊金霍洛旗為例［J］.干旱區(qū)資源與環(huán)境，2007，21（11）：123-130.

［21］ 許紀泉.武夷山區(qū)森林生態(tài)系統(tǒng)服務(wù)功能價值評估及其空間分布研究［D］.福州：福建師范大學(xué)，2007.

［22］ 任曉東.內(nèi)蒙古不同氣候區(qū)ETo估算方法（FAO-PM溫度法和Hargreaves-Samani法）的適用性評價［D］.呼和浩特：內(nèi)蒙古農(nóng)業(yè)大學(xué)，2014.

［23］ 周文佐，劉高煥，潘劍君.土壤有效含水量的經(jīng)驗估算研究：以東北黑土為例［J］.干旱區(qū)資源與環(huán)境，2003，17（4）：88-95.

［24］ 楊子睿.流域水源涵養(yǎng)價值評估及分區(qū)研究［D］.徐州：中國礦業(yè)大學(xué)，2019.

［25］ 傅斌，徐佩，王玉寬，等.都江堰市水源涵養(yǎng)功能空間格局［J］.生態(tài)學(xué)報，2013，33（3）：789-797.

［26］ 陳姍姍.南水北調(diào)水源區(qū)水源涵養(yǎng)與土壤保持生態(tài)系統(tǒng)服務(wù)功能研究［D］.西安：西北大學(xué)，2016.

［27］ 馬廣文，王業(yè)耀，香寶，等.松花江流域非點源氮磷負荷及其差異特征［J］.農(nóng)業(yè)工程學(xué)報，2011，27（S2）：163-169.

［28］ 張力，徐志金，滕志坤.松花江流域面源污染特征與防治對策［J］.環(huán)境科學(xué)與管理，2008（7）：55-56，61.

［29］ 岳勇，程紅光，楊勝天，等.松花江流域非點源污染負荷估算與評價［J］.地理科學(xué)，2007，27（2）：231-236.

［30］ 王秀娟，劉瑞民，何孟常.松遼流域非點源污染TN時空變化特征研究［J］.水土保持研究，2009，16（4）：192-196，202.

［31］ 潘麗娟.未來土地利用情景下的南京市生態(tài)系統(tǒng)水質(zhì)凈化功能模擬［D］.南京：南京信息工程大學(xué)，2016.

［32］ 范亞寧.秦嶺北麓及周邊生態(tài)系統(tǒng)水源涵養(yǎng)與水質(zhì)凈化功能評估［D］.西安：西北大學(xué)，2017.

［33］ 龔溪，曹銘昌，孫孝平，等.武夷山市生態(tài)系統(tǒng)服務(wù)價值評估［J］.生態(tài)與農(nóng)村環(huán)境學(xué)報，2017，33（12）：1094-1101.

［34］ 張修峰，劉正文，謝貽發(fā)，等.城市湖泊退化過程中水生態(tài)系統(tǒng)服務(wù)功能價值演變評估：以肇慶仙女湖為例［J］.生態(tài)學(xué)報，2007，27（6）：2349-2354.