基于InVEST模型的甘肅白龍江流域生境質量時空分異

鞏 杰, 馬學成, 張玲玲, 柳冬青, 張金茜

(蘭州大學 資源環境學院 西部環境教育部重點實驗室, 蘭州 730000)

生境質量是指生態系統能夠提供給物種個體、種群、群落以及人類持續生存和繁殖所需條件的能力[1]；可用于表征區域地類景觀斑塊生態適宜性好壞，其數值大小能夠反映區域生境斑塊的破碎程度和各景觀斑塊對生境退化的抗干擾性[2-4]。大中尺度上開展生境質量時空分異演變研究，是對自然生態系統健康識別、界定和保育的基礎，也是區域生態環境評價與可持續發展管理的前提條件，對制訂區域保育政策具有重要意義[5]。為此，國內外一些學者開展了不同區域的生境質量分析。Valero等[6]運用河岸帶植被指數(Riparian Forest Quality Index)和河流生境指數(Fluvial Habitat Index)開展了西班牙加利西亞河流生境質量研究，楊建強等[7]采用生態環境綜合指數法分析了萊州灣的生境質量，王瓊等[8]運用層次分析法(AHP)開展了蒲河流域生境質量評估。上述研究多傾向于建立生境評價指標體系來評估生境質量，其結果準確，但時效性相對較差，而在大中尺度區域上進行長時序列生境質量時空分異的報道較少，尤其是易災型山區流域生境質量的研究報道更為鮮見。近年來，隨著3S技術和生態模型的發展，區域生境質量評估逐漸趨向于模型化、空間化、定量化及精細化，如許多學者利用SolVES模型[9]、HSI模型[10]、InVEST模型[11]等，在自然保護區[12]、流域[3]、區域[13]等尺度上對生境質量進行分析研究。其中InVEST模型是目前較為成熟且應用最多的生態系統服務評估模型[14]，其可靠性高[15]。

InVEST模型是由美國斯坦福大學和世界自然基金會等機構研發的一種可以量化多種生態系統服務的評估模型，如產水量、土壤保持、水質凈化等，具備導入數據量少、輸出數據量大的特點，正被用于全球或區域不同生態系統服務評估[11,16]。甘肅白龍江流域位于青藏高原、秦巴山區和黃土高原的脆弱過渡帶，同時也是長江上游重要的水源涵養區和生態屏障，境內地質構造復雜，自然災害頻繁[17]，生態環境脆弱，嚴重威脅流域自然生態與人類生產活動。因此，本文以甘肅白龍江流域為對象，運用RS/GIS技術和InVEST模型對其1977—2014年的生境質量進行分析，擬揭示甘肅白龍江流域生境質量時空演變特征，旨在為流域自然生態系統保育與社會—經濟—自然復合生態系統可持續發展提供參考價值。

1　研究區概況

白龍江發源于四川省與甘肅省交界的碌曲縣郎木寺鎮，為嘉陵江上游的第一大支流，河源海拔4 072 m，河道全長576 km，流域面積為3.18萬km2。甘肅白龍江流域(32°36′—34°24′N，103°00′—105°30′E)位于甘肅省東南部，面積約1.8萬km2(圖1)，地處秦嶺地質構造褶皺帶，新構造活動強烈，溝深坡陡，地勢總體西北高東南低，地形起伏大，主要地貌類型有山地地貌、河谷地貌及黃土地貌，屬于亞熱帶向暖溫帶過渡區，夏季高溫多雨，冬季溫涼少雨，年均氣溫6～15℃，年均降水量400～850 mm。土壤以黃棕壤、棕壤、褐土、黑鈣土為主。植被覆蓋較好，森林廣布，流域源頭屬青藏高原高寒植被區，中上游以暖溫帶落葉闊葉林區為主，下游多為亞熱帶常綠闊葉林[18]，有珙桐(DavidiainvolucrataBaill)、水杉(Metasequoiaglyptostroboides)、銀杉(Cathaya)、香果樹(EmmenopteryshenryiOliv.)、連香樹(Cercidiphyllaceae)、鵝掌樹(AcerpaxiiFranch)、紅椿(ToonaciliataRoem.)、南方紅豆杉(Taxuschinensis(Pilger)Rehd.var.)、水青樹(TetracentronsinenseOliv.)等多種珍稀瀕危植物以及國家一級保護植物7種，二級保護植物19種[19]；當歸(Angelicasinensis)、黨參(Codonopsispilosula(Franch.)Nannf.)、大黃(RheumpalmatumL.)、紅芪(HedysarumpolybotrysHand.Mazz)、柴胡(Radixbupleuri)等藥用植物近700種；擁有大熊貓(Ailuropodamelanoleuca)、金絲猴(Rhinopithecusroxellana)、羚羊(Antidorcasmarsupialis)、短尾猴(Macacaarctoides)、麝(MoschusmoschiferusLinnaeus)、毛冠鹿(Elaphoduscephalophus)、鯢(Megalobatrachusjaponicusdaoidianuas(Blanchard))等20多種保護動物。流域內自然災害種類繁多，其中地質災害主要以滑坡、泥石流和地震為主[17]，其他災害主要表現為旱災、土壤侵蝕、霜凍、冰雹等(甘肅統計年鑒(1984—2014)及地方縣志)。至2014年末，流域內人口為123.77萬人，國內生產總值149.55億元，人均生產總值12 083元。

圖1甘肅白龍江流域位置及高程

2　研究方法

2.1　基礎數據準備

土地覆被數據主要通過Landsat TM/ETM+影像目視解譯獲得，其遙感影像來源于美國地質調查局和國際科學數據服務平臺，分辨率為30 m，時段分別為1977年、1990年、2002年、2014年，月份為7—8月。根據中國土地利用分類的標準和白龍江流域土地資源利用的實際情況，將研究區土地利用類型劃為6大類，即耕地、林地(包括有林地、灌木林地、疏林地、其他林地)、草地(包括高覆蓋草地、中覆蓋草地、低覆蓋草地)、水域(包括河流、湖泊、水庫、坑塘等)、建設用地(包括居民地、工礦用地、城鎮等)、未利用地(包括沙地裸地、荒灘、高山積雪、高寒荒漠等)。對解譯數據進行拓撲檢查和錯誤修改，并利用野外定點數據對其進行精度驗證，符合本研究精度要求。參考流域各縣區土里利用類型變更數據庫、甘肅省土地利用分布圖(1∶10萬，http:∥westdc.westgis.ac.cn/)和白龍江流域植被類型圖，對土地利用類型進一步細分，建立相應數據庫。

2.2　分析方法

(1) InVEST生境質量模塊的原理是利用不同土地覆被類型的脅迫因子敏感度和外界威脅強度，考慮脅迫因子的影響距離、空間權重及土地受法律保護程度等因素，將生境質量視為一個連續變量，由區域內可供生物生存、繁殖和發展所需資源的多少來決定，繼而表征生物多樣性的豐富性，即生境質量好的區域，其生物多樣性水平亦高，反之亦然[2-4,12]。具體計算過程如下：

(1)

(2)

式中：Qxj為土地覆被類型j中柵格單元x的生境質量；Hj為土地覆被類型j的生境適宜度；Dxj為土地覆被類型j中柵格單元x的生境脅迫水平；k為半飽和系數，通常取Dxj最大值的一半；z為常數。R為脅迫因子數；y為脅迫因子r的所有柵格單元；Yr為脅迫因子r所占柵格單元總數；綜合考慮研究區特殊地理環境，結合相關文獻[17,20]與各縣縣志，文章選取對地表景觀影響較大的城鎮(人口大于5 000人)、居民點、人口密度、耕地、道路(國道省道和縣鄉道路)、綜合生態風險源(滑坡、泥石流、地震、干旱、土壤侵蝕)定義為生態脅迫因子，并參考前人研究[2-4,12]對其賦值，見表1—2。ry為土地覆被類型柵格單元y中的脅迫因子r；柵格單元y中的脅迫因子r(ry)對生境柵格單元x的脅迫作用為irxy，公式為：

(3)

(4)

式中：dxy為柵格單元x與柵格單元y之間的直線距離；drmax為脅迫因子r的最大影響距離；wr為脅迫因子r的權重；βx為柵格單元x的合法可達性(或法律準入程度或到達度)，本文因研究區缺失土地保護區數據，暫不考慮白龍江流域特定的土地保護區因素，將流域內各土地覆被類型在國家法律法規上受到保護的程度認為是平等均一的，系統將自動賦值為1，即完全到達；Sjr為土地覆被類型j對脅迫因子r的敏感性，取值0～1。

(2) 土地利用變化動態度(K)是定量描述區域一定時間范圍內某種土地利用類型面積的變化速率，在一定程度上可以反映某種土地資源擴張或萎縮程度[21]。其計算公式為：

(5)

式中：K為研究時段內某土地利用類型的變化動態度；Ua和Ub為研究初期和末期的某土地利用類型面積(km2)；T為研究時段。

表1　脅迫因子的最大影響距離及其權重

3　結果與分析

3.1　土地利用變化

林地和草地一直是甘肅白龍江流域的主要覆被類型(圖2A)，占流域總面積的80%以上。1977—1990年，耕地和建設用地面積的增幅為39.8%和10.68%，其他覆被類型面積減??；1990—2002年，建設用地和耕地呈現增長態勢，年均增長率分別為2.69%和2.28%，林地面積減少量為563.07 km2，草地為171.48 km2；2002—2014年，建設用地和林地面積的增幅為61.81%和30.55%，耕地的減幅為35.59%。耕地在1977—2002年動態度為正值，2002—2014年為負值，其面積的減少速率相當于1977—1990年的增長速率，林地則正好表現出與耕地相反的變化特征(圖2B)，這是因為1999年退耕還林還草政策實施以來，陡坡耕地逐漸退耕導致林草地的增加造成?？傮w上，1977—2014年草地、未利用地和水域3種覆被類型面積減少率約為0.5%，林地和耕地的增率為0.4%和0.36%，而建設用地的增率為3.7%，說明隨著當地人口增加及經濟發展，居民工礦用地、交通用地等建設用地快速擴張，水域面積則因人為填埋及等影響而有所減少。

表2　生境適宜度及其對不同脅迫因子的敏感度

圖2　甘肅白龍江流域土地利用類型面積變化及其動態度

3.2　生境質量變化

甘肅白龍江流域生境質量時間變化見表3，1977年流域整體生境質量較高，高值區(5～7)面積比重較大；1990年以來生境質量整體較1977年有所降低，低值區(1～3)面積比重增加，高值區(5～7)面積比重減少；2014年流域生境質量相對2002年有所改善，在流域西北部及東南部有一定變化(圖3C與圖3D)，相應地高值區(5～7)面積比重有所增加。1977—2014年期間，生境質量呈現出先降低后緩慢增加的趨勢，減幅和增幅為17%和25%，這是因為1977—2002年人們大面積的毀林開荒，使得林草地銳減，自然生態系統長時間的受壓導致景觀破碎化，生境質量下降；隨著甘肅省于1998年啟動實施的天然林資源保護工程以及1999年啟動的退耕還林工程等生態建設政策的實施，加上流域內環境保護與管理程度的不斷加強，到2014年流域生態環境得到一定改善，生境質量相應提高。

表3　甘肅白龍江流域生境質量統計(1977-2014年)

甘肅白龍江流域生境質量空間分布見圖3，其空間格局變化具有一定的規律性，生境質量較高區域(數值大于0.42)面積所占比重較小，主要分布在流域西北部、文縣東南部及舟曲南部，集聚于白水江國家級自然保護區、迭部境內的阿夏和多爾大熊貓自然保護區、達拉溝流域及林業發展區；生境質量較低區域(數值小于0.16)面積所占比重較大，主要分布在人類活動相對頻繁的河谷地帶、城鄉聚落區、耕作區、裸地和低覆蓋草地等區域，集聚于舟曲—武都段白龍江近岸區及文縣和武都交界帶、宕昌西北部等。在縣域尺度上，迭部縣生境質量最好，其次是文縣，武都區生境質量較差。

圖3甘肅白龍江流域1977-2014年生境質量空間變化

4　結論與討論

(1) 綜合分析1977—2014年甘肅白龍江流域不同土地覆被類型變化發現，林地和草地是該流域的主要土地覆被類型，占流域總面積的80%以上；而土地覆被類型的改變能夠影響生境斑塊之間的能量流和物質流過程，繼而影響區域的生境分布格局和功能[22]，說明林草地的增加對流域生境質量的改善和維持具有很重要的作用。

(2) 總體來看，甘肅白龍江流域生境質量在研究期內呈現出先降低后緩慢增加的趨勢。甘肅省政府2013年末發布的《隴南實施退耕還林工程紀實》中指出，二十世紀六七十年代的開荒運動使隴南生態“滿目瘡痍”；1999年啟動實施退耕還林工程后，林地面積累計新增2 166.67 km2及森林覆蓋率提高1.5個百分點的客觀事實間接相符合，與本研究結果相近。甘肅白龍江流域生境質量的空間異質性也呈現出一定的規律性，生境質量較低的區域主要分布在人類活動強度高、自然災害風險較高的白龍江兩岸和宕昌西北部，以城鄉聚落區、耕作區、坡耕地和低覆蓋草地等最為突出。區域景觀格局的高度破碎化，會使生境退化程度加劇，景觀結構穩定性趨于降低，導致生境質量下降；生境質量較高的區域主要分布在流域西北部、文縣東南部及舟曲南部，集聚于自然保護區和林業管護區，如白水江國家自然保護區、阿夏及多爾自然保護區及宕昌縣官鵝溝風景區等，這些區域受人類活動干擾相對較弱，這與其他學者[4,12，23]研究結果相似。

(3) 本文運用InVEST模型時，進行了模型參數本地化處理(包括對模型腳本修改和參數校驗)，并基于RS/GIS開展了流域生境質量評價及其時空演變特征分析，研究結果說明該模型可在山區流域，特別是像甘肅白龍江流域這樣缺乏大量實測數據的易災型流域，能夠替代其他復雜方法，對生境質量變化做出快速而準確的評價。因甘肅白龍江流域面積較大且地勢起伏劇烈，地貌復雜多變，缺乏大型野外觀測臺站，這在一定程度上制約著模型結果的精確驗證。本研究部分參數數據(如不同土地覆被類型對各生態脅迫因子的敏感度)來源于文獻、InVEST模型數據庫以及經驗公式，而生境質量機理相對復雜，加上指標選取的不全面性和模型結構的不確定性，分析結果有一定的誤差，但不會影響流域生境質量的整體變化趨勢。建議后期進一步優化研究方案，理清生境質量的機理，強化評估模型及其參數的適宜性，加強野外實測數據，進行相關參數本地化處理結果的準確驗證，結合宏觀尺度的遙感觀測，更好地補充生境質量在時間和空間尺度上的研究不足[24]。

(4) 綜合考慮流域生境高度脆弱和自然災害頻發的現狀，為進一步保育流域生境，開展生態建設和災害風險防范，實現流域可持續發展，一些對策和建議如下：① 白水江自然保護區南部、武都東南部、拱壩河中上游及博峪河上游區域、宕昌縣南部迭山林區、迭部縣大部分區域主要是流域生境質量優良區域，土地覆被類型以森林為主，加上人類活動強度較弱，自然植被覆蓋良好，對于這些區域應該繼續加強自然保護區的管理和維護，禁止亂砍亂伐，減少或者禁止人類活動對物種生存環境的干擾，積極實施生態工程，恢復自然植被，降低人類活動的擾動，保障流域生態安全。② 甘肅白龍江流域是甘肅省唯一具有河谷亞熱帶氣候及豐富生物資源的森林景觀區，水土流失嚴重，山洪、崩塌、滑坡、泥石流等自然災害頻發，造成強烈侵蝕及堆積物常常堵塞河流和道路，加上多數村莊都建在泥石流堆積扇上，嚴重威脅著居民生命和財產安全。因此，在做好山區災害宣傳教育的同時，應加大災害監測、防治與預警研究，同時實施移民搬遷和環境整治，開展滑坡泥石流等自然災害綜合治理和環境修復[25]，優化人類活動，提升流域生態安全水平和人地和諧關系。

[1]Hall L S, Krausman P R, Morrison M L. The habitat concept and a plea for standard terminology [J]. Wildlife Society Bulletin, 1997,25(1):173-182.

[2]Sharp R, Tallis H T, Ricketts T, et al. InVEST version user′s guide [R]. The Natural Capital Project, Stanford University, University of Minnesota, The Nature Conservancy, World Wildlife Fund, 2016.

[3]肖明. GIS在流域生態環境質量評價中的應用:以昌化江下游為例[D].海口:海南大學,2011.

[4]劉志偉.基于InVEST的濕地景觀格局變化生態響應分析:以杭州灣南岸地區為例[D].杭州:浙江大學,2014.

[5]吳健生,曹祺文,石淑芹,等.基于土地利用變化的京津冀生境質量時空演變[J].應用生態學報,2015,26(11):3457-3466.

[7]楊建強,朱永貴,宋文鵬,等.基于生境質量和生態響應的萊州灣生態環境質量評價[J].生態學報,2014,34(1):105-114.

夏夜，蚊子飛舞在我頭上。蚊子時而落在臉上，時而落于耳邊。我從睡夢中醒來，用手拍打驅趕它。哪知秋后蚊子可不是舉手投足間一招兩式便可消滅的。蚊子飛來飛去，任我怎樣努力，把自己的臉都快打腫了，都不能將其正法。一時間，蚊子侵占了我的世界，整個天地之間似乎都是其嗡嗡聲。隨即睡意全無，心情憤怒。

[8]王瓊,范志平,李法云,等.蒲河流域河流生境質量綜合評價及其與水質響應關系[J].生態學雜志,2015,34(2):516-523.

[9]王玉,傅碧天,呂永鵬,等.基于SolVES模型的生態系統服務社會價值評估:以吳淞炮臺灣濕地森林公園為例[J].應用生態學報,2016,27(6):1767-1774.

[10]王志強,陳志超,郝成元.基于HSI模型的扎龍國家級自然保護區丹頂鶴繁殖生境適宜性評價[J].濕地科學,2009,7(3):197-201.

[11]Goldstein J H, Caldarone G, Duarte T K, et al. Integrating ecosystem-service tradeoffs into land-use decisions [J]. Proceedings of The National Academy of Science of The United States of America, 2012,109(19):7565-7570.

[12]包玉斌,劉康,李婷,等.基于InVEST模型的土地利用變化對生境的影響:以陜西省黃河濕地自然保護區為例[J].干旱區研究,2015,32(3):622-629.

[13]榮月靜,張慧,王巖松.基于Logistic-CA-Markov與InVEST模型對南京市土地利用與生物多樣性功能模擬評價[J].水土保持研究,2016,23(3):82-89.

[14]Wong C P, Jiang B, Kinzig A P, et al. Linking ecosystem characteristics to final ecosystem services for public policy [J]. Ecology Letters, 2015,18(1):108-118.

[15]Terrado M, Sabater S, Chaplin-kramer B, et al. Model development for the assessment of terrestrial and aquatic habitat quality in conservation planning [J]. Science of The Total Environment, 2016,540(1):63-70.

[16]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:1418-1426.

[17]Meng X, Chen G, Guo P, et al. Research of landslides and debris flows in Bailong River Basin:progress and prospect [J]. Marine Geology and Quaternary Geology, 2013,33(4):1-15.

[18]高彥凈.基于CASA模型的植被NPP時空動態研究:以甘肅白龍江流域為例[D].蘭州:蘭州大學,2015.

[19]湯萃文,孫學剛,肖篤寧.甘肅省中國種子植物特有屬物種多樣性保護優先地區分析[J].生態學雜志,2005,24(10):1127-1133.

[20]寧娜,馬金珠,張鵬,等.基于GIS和信息量法的甘肅南部白龍江流域泥石流災害危險性評價[J].資源科學,2013,35(4):892-899.

[21]王秀蘭,包玉海.土地利用動態變化研究方法探討[J].地理科學進展,1999,18(1):81-87.

[22]歐陽志云,鄭華.生態系統服務的生態學機制研究進展[J].生態學報,2009,29(11):6183-6188.

[23]齊姍姍,鞏杰,錢彩云,等.基于SRP模型的甘肅省白龍江流域生態環境脆弱性評價[J].水土保持通報,2017,37(1):224-228.

[24]Lausch A, Bannehr L, Beckmann M, et al. Linking earth observation and taxonomic, structural and functional biodiversity:from local to ecosystem perspectives [J]. Ecological Indicators, 2016,70:317-339.

[25]鞏杰,趙彩霞,王合領,等.基于地質災害的隴南山區生態風險評價:以隴南市武都區為例[J].山地學報,2012,30(5):570-577.