李存秀劉峰貴周源濤陳 瓊周 強何文鑫田萬霞

(1.青海師范大學(xué) 地理科學(xué)學(xué)院,青海 西寧810008;2.高原科學(xué)與可持續(xù)發(fā)展研究院,青海 西寧810008)

土地覆被/土地利用作為人類活動響應(yīng)氣候變化的重要途徑之一[1],在全球變化中扮演著重要的角色,是全球變化的重要驅(qū)動力。工業(yè)革命以來,隨著經(jīng)濟和科學(xué)技術(shù)的迅速發(fā)展,以及在全球變化的高速推動下,人地矛盾日益尖銳[2]。近年來中國特別強調(diào),保障生態(tài)文明有效建設(shè)和可持續(xù)發(fā)展的關(guān)鍵是維持人地關(guān)系協(xié)調(diào)發(fā)展[3]。在生態(tài)環(huán)境脆弱、自然環(huán)境復(fù)雜的西藏地區(qū),人地關(guān)系協(xié)調(diào)演進無疑是一項重大挑戰(zhàn)。隨著人地關(guān)系研究的不斷推進,西藏地區(qū)人地關(guān)系優(yōu)質(zhì)協(xié)調(diào)發(fā)展的勢頭良好[4]。珠峰北麓朋曲流域是中國珠穆朗瑪峰自然保護區(qū)的核心腹地,高山草甸作為環(huán)珠峰地區(qū)的主要植被類型,20世紀(jì)末至21世紀(jì)初,隨著牧業(yè)比重的增大,環(huán)珠峰地區(qū)過牧現(xiàn)象嚴(yán)重,草地退化加劇,森林資源也處于明顯的退化狀態(tài),退化面積大于恢復(fù)面積,區(qū)域生態(tài)環(huán)境遭到破壞,人地關(guān)系沖突日益明顯[5],因此,加強環(huán)珠峰地區(qū)生態(tài)環(huán)境保護問題被提到重要地位。

從20世紀(jì)60年代第一次展開珠峰科學(xué)考察開始,珠峰自然保護區(qū)就成為高原人地關(guān)系的研究熱點,研究聚焦人與自然相互作用的評估[6]、氣候變化對環(huán)珠峰地區(qū)植被動態(tài)變化規(guī)律及影響等研究[7-8]。研究結(jié)果表明,自2002年以來,隨著西藏退耕還林、還草工程的實施以及環(huán)珠峰地區(qū)生態(tài)保護工程的實施,珠峰南坡植被類型恢復(fù)良好,珠峰北坡植被穩(wěn)定發(fā)展[9],但耕地面積明顯減少。環(huán)珠峰地區(qū)土地覆被研究主要集中于植被現(xiàn)狀和恢復(fù)研究,但對未來情景下的土地覆被動態(tài)的研究不足。

當(dāng)前用于土地覆被變化模擬預(yù)測的模型主要有:CLUE-S模型、Logistic回歸模型、Markov模型和FLUS模型等[10]。FLUS模型是對傳統(tǒng)預(yù)測模型的改進,能有效解決傳統(tǒng)預(yù)測模型中出現(xiàn)的參數(shù)設(shè)置和類型之間轉(zhuǎn)換困難等問題,有效處理人為和自然因素雙重影響下的土地覆被/土地利用變化不確定性;Markov模型對預(yù)測未來情景下的土地覆被/土地利用像元數(shù)量具有較高的實用性,但難以對變化進行空間表達(dá)[11]。為此,基于Markov-FLUS模型中和了Markov模型的缺點,以達(dá)到像元數(shù)量預(yù)測和空間表達(dá)的雙重優(yōu)勢,適用于中小流域未來土地覆被情況變化的研究。

在青藏高原生態(tài)保護和可持續(xù)發(fā)展的大背景下,對耕地和生態(tài)安全保護將對未來情境下的流域水土響應(yīng)與人地關(guān)系的和諧發(fā)展產(chǎn)生重要影響。因此,針對流域內(nèi)土地覆被變化未來發(fā)展的研究非常之必要,既可以對保障糧食安全和保護生態(tài)安全提供建設(shè)性意見,又為下一步生態(tài)文明發(fā)展提供基礎(chǔ)理論支撐。朋曲流域作為環(huán)珠峰北麓地區(qū)的人口中心,是珠穆朗瑪峰國家級自然保護區(qū)的核心區(qū)與人地關(guān)系相對復(fù)雜的地帶,同時作為珠峰地區(qū)重要的水文基本單元,統(tǒng)籌流域的相關(guān)研究,是構(gòu)建國家生態(tài)安全屏障的重要研究對象。對朋曲流域未來土地覆被情況的研究,可以反映高海拔高脆弱性的特殊自然環(huán)境下未來土地利用情況的動態(tài)變化規(guī)律,以及不同土地利用限制情景下,土地覆被在未來近期和中期可能發(fā)生的演變。本文基于Markov-FLUS模型,運用不同情景約束因子,對朋曲流域短期(2030年)和中遠(yuǎn)期(2050年)的土地覆被類型進行模擬預(yù)測分析,對比兩種類型下的耕地發(fā)展和生態(tài)保護發(fā)展,為促進朋曲流域耕地保護工作和實現(xiàn)生態(tài)空間的科學(xué)布局提供思路,為實現(xiàn)高原地區(qū)和諧的人地關(guān)系服務(wù)。

1 研究區(qū)概況

朋曲發(fā)源于喜馬拉雅中段的希夏邦馬峰,源頭以冰川融水為主,大氣降水補給為輔,河流自西向東在定結(jié)縣境內(nèi)折向南,切穿喜馬拉雅山脈入尼泊爾,最后注入印度恒河,上游中國境內(nèi)流域范圍涉及西藏日喀則市南部的崗巴縣、定日縣、定結(jié)縣、聶拉木縣和薩迦縣,全長376 km,流域總面積2.50×104km2,平均海拔5 259 m,支流主要有扎嘎曲、葉如藏布江、西寧藏布和卡得藏布等[12]。喜馬拉雅中段的珠峰地區(qū)是青藏高原人類活動歷史悠久且較為頻繁的區(qū)域,該河流屬于跨境河流,是喜馬拉雅中段人類南北遷徙和貿(mào)易的主要通道,也是珠穆朗瑪峰國家自然保護區(qū)的核心區(qū)域,因此,該流域受到廣泛關(guān)注。

朋曲流域焚風(fēng)效應(yīng)明顯,印度洋季風(fēng)暖濕氣流在喜馬拉雅山脈南坡形成大量降水,越過山脊到達(dá)北坡后下沉,導(dǎo)致該地區(qū)常年氣溫較低,降水偏少,年均降雨量僅僅維持在600 mm左右,屬于典型的高原大陸性氣候,寒冷而干旱,由此北坡植被多為草甸和灌叢。位于珠峰南坡的少部分區(qū)域,受印度洋季風(fēng)影響,降雨充沛[13],氣候溫暖濕潤,形成大片原始森林。朋曲水系包括熱曲藏布、扎嘎曲、羅別藏布、彭作浦曲、西糾曲、魯魯曲等小河流和主河道朋曲河[14]。據(jù)統(tǒng)計,2020年朋曲流域總?cè)丝跒?6.5萬人,人口密度約為4.7人/km2,半農(nóng)半牧,耕地主要分布在河谷地帶,農(nóng)作為以種植青稞、小麥、油菜為主[15]。

2 數(shù)據(jù)來源與方法

2.1 數(shù)據(jù)來源與處理

本文數(shù)據(jù)包括朋曲流域縣界數(shù)據(jù),空間驅(qū)動因子數(shù)據(jù)集和土地覆被數(shù)據(jù)集,研究區(qū)縣界數(shù)據(jù)來源于2019年西藏自治區(qū)縣級行政邊界提取的數(shù)據(jù);空間驅(qū)動因子數(shù)據(jù)集包含自然因子、社會經(jīng)濟因子和距離因子等數(shù)據(jù);2010,2015和2020年3期土地覆被數(shù)據(jù)來源“地球大數(shù)據(jù)工程數(shù)據(jù)共享服務(wù)系統(tǒng)”[16](表1)。本文以2020年最新研發(fā)的高精度、全球性土地分類數(shù)據(jù)(GLC_FCS30-2020)為基準(zhǔn)數(shù)據(jù),該數(shù)據(jù)集共包含29個地表覆蓋類型,更新周期為5 a。本文通過篩選所需要的土地覆被類型,將29個地表類型進行重分類操作,最后提取需要的5個土地覆被類型。數(shù)據(jù)來源詳見表1。

表1 研究所用各類數(shù)據(jù)來源

2.2 FLUS模型

FLUS模型是主要用于當(dāng)前和未來情景下的土地覆被模擬。該模型的原理源自元胞自動(selfadaptive inertia and competition mechanism CA),并且在其基礎(chǔ)上做了進一步修改,增加了人工神經(jīng)網(wǎng)絡(luò)模型算法(ANN-based probability-of-occurrence estimation)和馬爾科夫鏈計算(Markov chain)兩大模塊[18],模型的改善和發(fā)展進一步耦合了人類活動和自然要素之間相互競爭和轉(zhuǎn)化的模擬優(yōu)勢[19]。

2.2.1 適宜性概率統(tǒng)計 人工神經(jīng)網(wǎng)絡(luò)模型可用以表征未來情景土地覆被類型和多項空間驅(qū)動因子之間的對應(yīng)關(guān)系,將初始年份的土地覆被數(shù)據(jù)和各項驅(qū)動因子作為輸入數(shù)據(jù)源,采用神經(jīng)網(wǎng)絡(luò)算法[20](ANN)整合并計算朋曲流域內(nèi)5種土地覆被類型出現(xiàn)在每個像元上的概率[21],計算公式為:

式中:Nj(p,q)為隱藏層神經(jīng)元j接收的輸入層信號;wij為訓(xùn)練中輸入層和隱藏層之間設(shè)定的權(quán)重[22]。

2.2.2 馬爾科夫鏈預(yù)測未來像元數(shù) 利用馬爾科夫鏈模塊預(yù)測未來情景下的土地覆被類型像元數(shù)量[23],在2010,2015和2020年的數(shù)據(jù)基礎(chǔ)上,預(yù)測2020年和2025年的土地覆被類型像元數(shù)量(表2),計算公式為:

表2 研究區(qū)2010—2050年多期土地覆被類型像元總數(shù)

式中:Tn+1為n+1時刻土地覆被狀態(tài);Tn為n時刻土地覆被狀態(tài);pij為各地類的轉(zhuǎn)變適宜性概率。

2.2.3 領(lǐng)域因子權(quán)重計算 領(lǐng)域因子用以表征不同土地覆被類型的擴張能力[18],擴張能力范圍為[0,1],范圍數(shù)值大小與擴張強度成正比,當(dāng)擴張數(shù)值越大接近于1時,表示該土地覆被類型在其領(lǐng)域范圍內(nèi)的擴張能力越強。當(dāng)前用于確定權(quán)重參數(shù)的方法主要有:因子分析法、Moore鄰域模型法[24]、無量綱化。無量綱化方法將原始數(shù)據(jù)進行無量綱化處理,根據(jù)每個指標(biāo)值的差異,客觀確定指標(biāo)權(quán)重[25]。本文考慮土地類型擴張規(guī)律的客觀性,使用朋曲流域各土地類型歷史總面積變化數(shù)量進行無量綱化處理,來確定各地類的領(lǐng)域因子權(quán)重值[26](表3)。

表3 朋曲流域無量綱化得到的領(lǐng)域因子權(quán)重值

2.2.4 轉(zhuǎn)換成本計算 轉(zhuǎn)換成本反映各土地類型之間的變化規(guī)則,通常用矩陣表示,0表示兩種地類之間轉(zhuǎn)換成本高,不容易轉(zhuǎn)換,1表示兩地類之間轉(zhuǎn)換成本最低,最容易轉(zhuǎn)換。當(dāng)前關(guān)于計算轉(zhuǎn)換成本矩陣的方法主要有專家意見咨詢、經(jīng)驗轉(zhuǎn)換和政策導(dǎo)向轉(zhuǎn)換規(guī)則,本文利用朋曲流域2010—2020年土地類型轉(zhuǎn)移矩陣為基礎(chǔ),以專家咨詢意見為輔,確定各地類的轉(zhuǎn)換成本系數(shù)。

2.3 模型精度驗證

本文將驗證模型模擬精度高低的評價指標(biāo)分為兩類:分類別精度評價指標(biāo)〔生產(chǎn)者精度PA(producer’s accuracy)〕和用戶精度UA(user’s accuracy)檢驗[27])和綜合精度評價指標(biāo)〔總體精度OA(overall accuracy)和kappa系數(shù)(kappa coefficient)檢驗〕。其中,總體精度OA和Kappa系數(shù)主要檢驗?zāi)M整體結(jié)果與實際結(jié)果之間的一致性,用以評價模擬圖像的精度問題,兩個結(jié)果的檢測值越接近于1,表明兩個圖像之間具有較高的相似性[28],本文主要采用該方法進行模擬精度的驗證。以朋曲流域2015年為初始年份數(shù)據(jù),對影響因子歸一化處理,采用FLUS模型的馬爾科夫鏈模塊預(yù)測2020年的土地類型數(shù)量,計算得到朋曲地區(qū)2015—2020年的土地覆被轉(zhuǎn)移適宜性概率數(shù)據(jù),然后再對2020年的朋曲流域土地覆被進行模擬,將模擬結(jié)果與2020年實測的土地覆被類型進行精度驗證,利用FLUS模型模擬出的2020年數(shù)據(jù)結(jié)果總體精度(OA)為90%,kappa系數(shù)為0.82。其中草地的分類精度最優(yōu),PA和UA均在90%以上;其他用地類型的分類精度次之,PA和UA均超過88%;林地的UA高達(dá)93%,但是PA僅有66%;水域的分類精度相較于前3種類型較低,PA和UA僅有65%和77%;耕地的分類精度是最低的,PA和UA均不超過60%,分類精度一般(圖1)?？傮w而言,FLUS模型對朋曲流域土地覆被類型的模擬結(jié)果精度較高,達(dá)到預(yù)期目標(biāo),模型可以用于本文關(guān)于未來情景土地覆被模擬預(yù)測研究。

圖1 朋曲流域2020年土地覆被現(xiàn)狀與模擬結(jié)果比較

3 結(jié)果與分析

3.1 歷史時期朋曲流域土地覆被變化

根據(jù)2010—2020年朋曲流域土地覆被面積變化表和空間分布圖(圖2)可知,從2010—2020年期間,林地的面積變化最高,其次是草地,再次是耕地,其他用地和水域的面積變化是最小的。林地主要分布在聶拉木縣、定結(jié)縣和定日縣南部地區(qū),其面積變化最為明顯,10 a間增加了778.16 km2,主要增加的區(qū)域為聶拉木西北部和定日縣中部地區(qū);其次是草地,草地是研究區(qū)內(nèi)面積最大的地物類型,廣泛分布于研究區(qū)內(nèi),面積呈顯著較大的減少變化趨勢,10 a以來共減少726.89 km2,草地面積減少區(qū)域主要為林地的擴張;再次為耕地面積變化,河谷地區(qū)零星呈帶狀分布的耕地,面積呈減少趨勢,但是減少面積較小,為67.78 km2;其他用地面積稍有增加,為8.46 km2;朋曲流域的河流受到冰川融雪和大氣降水的多重補給,常年徑流量穩(wěn)定,水域面積基本無變化;造成上述這種變化的主要原因是朋曲流域是以牧業(yè)為主農(nóng)業(yè)為輔的地區(qū),畜牧業(yè)是當(dāng)?shù)剞r(nóng)牧民最重要的經(jīng)濟來源,草地面積在該流域占很重要的份額(表4)。隨著氣候的變化和社會經(jīng)濟的發(fā)展,越來越多的農(nóng)牧民開始脫離草場,從事非農(nóng)牧活動[14],人類活動區(qū)域逐漸向城鎮(zhèn)集中,耕地面積逐漸減少,同時隨國家“退耕還林”與珠峰國家級自然保護區(qū)的建立,撂荒地逐漸被林地取代。

圖2 朋曲流域2010—2020年土地覆被空間分布

表4 朋曲流域2010—2020年土地覆被面積變化

3.2 未來情景土地覆被變化分析

利用FLUS模型模擬預(yù)測未來短期(2030年)和未來中遠(yuǎn)期(2050年)朋曲流域土地覆被狀況,但由于馬爾科夫模塊預(yù)測的土地覆被類型像元數(shù)量結(jié)果,可分為3種限制情景的模擬預(yù)測,未來情景的土地覆被類型預(yù)測需要同時考慮不同限制因子對土地擴張的約束。因此,設(shè)置多個約束用地變化的限制因子(表1),對未來短期與未來中遠(yuǎn)期不同限制情景下的流域土地覆被變化進行模擬預(yù)測,其結(jié)果如下:

3.2.1 自由發(fā)展情景 自由發(fā)展情景旨在不考慮相關(guān)約束用地發(fā)展的政策因子,僅考慮在自然和人為影響下土地覆被的變化情況,圖3和圖4反映了在此情境下未來短期(2030年)和未來中遠(yuǎn)期(2050年)朋曲流域土地覆被空間分布狀況與面積變化,較2020年:①未來短期土地覆被變化中,林地和其他用地面積呈增加趨勢,而耕地、草地和水域較歷史時期呈減少趨勢,其中,其他用地的面積變化是所有土地覆被類型中變化最為明顯的類型;②未來中遠(yuǎn)期土地覆被變化相較于2020年有些地類有明顯的向外擴張趨勢,林地面積的變化最為明顯,較2020年增加3 498.35 km2,從空間分布圖上看出2050年林地面積擴張明顯,擴張范圍主要集中在崗巴縣和定結(jié)縣范圍內(nèi),同時陳穎峰等[29]人的研究結(jié)果也表示,在2000—2009年間,珠峰東南區(qū)域(崗巴縣和定結(jié)縣)林地等植

圖3 朋曲流域自由發(fā)展情景下土地覆被類型面積變化比較

圖4 朋曲流域自由發(fā)展情景下土地覆被類型面積變化

被類型增加明顯,研究結(jié)果與以上結(jié)論吻合。在不受政策等限制因子去約束土地類型擴張的條件下,其他用地類型將快速擴張,擴張的面積主要占據(jù)了草地、耕地和水域等用地類型,若不加以相關(guān)政策因子限制,將直接威脅研究區(qū)內(nèi)的生態(tài)環(huán)境和糧食安全。

3.2.2 耕地保護情景 基于耕地保護情景的模擬預(yù)測,以基本農(nóng)田保護區(qū)為限制性因素,在該情景下的耕地發(fā)展穩(wěn)定性較強,耕地類型向其他類型轉(zhuǎn)化的條件受到限制。未來短期和未來中遠(yuǎn)期的土地覆被變化在此限制情景下發(fā)生了較大的變化(圖5—6)。①在此情景下,未來近期耕地面積126.8 km2,較2020年耕地面積無明顯變化,耕地減少的速率變緩,耕地保護的措施得到有效發(fā)揮。面積變化最明顯的是其他用地類型,較2020年面積增加393.19 km2,這部分增加的面積主要來自于對草地面積的轉(zhuǎn)換,致使草地面積減少幅度較大。朋曲流域的土地覆被類型復(fù)雜,其他用地類型里面包含的二級類型又充滿了復(fù)雜性,其他用地類型的擴張,一方面有朋曲流域少部分的建設(shè)用地和冰川積雪地向外擴張的原因,而最重要的一方面主要為草地轉(zhuǎn)化為裸土地,成為其他用地外擴的主要原因;②未來中遠(yuǎn)期林地的面積變化最大,較2020年面積增加3 485.55 km2,平均每年增加116.2 km2,這部分增加的面積主要來源于對草地和其他用地類型的轉(zhuǎn)化。耕地變化是其余地類中面積變化最小的一類,耕地保護情景約束了未來中遠(yuǎn)期耕地向其他地類轉(zhuǎn)出的趨勢。從面積變化趨勢分析,加入農(nóng)田保護區(qū)作為限制因子的朋曲流域土地覆被變化,減緩了耕地轉(zhuǎn)化的速率,可以有效保護耕地。

圖5 朋曲流域耕地保護情景下土地覆被類型面積變化比較

圖6 朋曲流域耕地保護情景下土地覆被類型面積變化

3.2.3 生態(tài)保護情景 朋曲流域大部分區(qū)域處于珠穆朗瑪峰自然保護區(qū)范圍內(nèi),生態(tài)保護是該區(qū)域可持續(xù)發(fā)展的基本方向。在生態(tài)保護情景下,從各用地類型數(shù)量變化見圖7—8。①此情景下未來短期土地覆被類型變化較大的仍集中在草地和其他用地類型,相較于2020年的草地發(fā)展,未來近期草地面積減少362.67 km2,其他用地面積增加389.29 km2。此情景下的土地覆被類型面積變化與耕地保護情景下的面積變化類似,但是在空間分布上卻出現(xiàn)了較大的不同。受自然保護區(qū)的影響,林地范圍在研究區(qū)東南方向擴張迅速,并且呈大塊連片分布;②未來中遠(yuǎn)期覆被類型面積變化也是主要草地、林地和其他用地,林地仍然是變化最大最為明顯的地類,而草地和其他用地仍然呈減少趨勢,面積受到進一步壓縮,草地和其他用地的轉(zhuǎn)出的主要方向為林地。但值得一提的是,此時期水域面積較2020年增加了117.538 km2,由于對氣候變化的響應(yīng),珠峰自然保護區(qū)變暖趨勢要比中國其他區(qū)域要更加顯著[30],致使珠峰地區(qū)湖泊冰雪融水補給量增加,面積擴大,這與王毅等[31]人的研究結(jié)果相吻合。

圖7 朋曲流域生態(tài)保護情景下土地覆被類型面積變化比較

圖8 朋曲流域生態(tài)保護情景下土地覆被類型面積變化

生態(tài)保護情景下,土地覆被變化特征主要為草地的減少以及其他用地類型的增加,耕地面積變化在整體變化中所占比例不高,但依舊呈減少趨勢。此情景下土地覆被類型變化較大主要集中在草地和其他用地類型,但是相較于自由發(fā)展情景和耕地保護情景下的土地擴張,生態(tài)保護下的土地覆被總體而言變化最小,土地轉(zhuǎn)化速率趨于放緩,這對保護生物多樣性,保障生態(tài)安全區(qū)域生態(tài)安全具有重要價值,同時我們也關(guān)注到,其他利用地與草地在此情景下變化仍最明顯,究其原因可能與模型中使用的馬爾可夫鏈有關(guān)。執(zhí)行農(nóng)田保護區(qū)和自然保護區(qū)的同時限制的相關(guān)政策,將對本區(qū)域的生態(tài)恢復(fù)、發(fā)展和經(jīng)濟建設(shè)提供重要的參考,有助于實現(xiàn)區(qū)域的可持續(xù)發(fā)展。

4 討論與結(jié)論

4.1 結(jié)論

(1)利用Markov-FLUS模型,模擬2020年朋曲流域的土地覆被狀況,和2020年實際的土地覆被現(xiàn)狀進行對比。經(jīng)過驗證,模型分類的總體精度(OA)為90%,kappa系數(shù)為0.82,精度較高,達(dá)到預(yù)期模擬精度目標(biāo),所以可以采用此模型進行未來情景的土地覆被模擬預(yù)測。

(2)通過統(tǒng)計分析2010,2015和2020年土地覆被狀況,林地和草地的面積變化最為明顯,草地面積在10 a間減少726.89 km2,林地面積增加778.16 km2,主要是由草地轉(zhuǎn)換為林地,主要是因為西藏地區(qū)近些年的植樹造林政策的有效實施。

(3)在3種情境下,未來近期(2030年)和未來中期(2050年)的土地覆被均發(fā)生明顯變化。未來近期和未來中期的土地覆被變化分別在連續(xù)3種情境下的表現(xiàn)一致,未來短期土地覆被變化最明顯均集中于草地和其他用地兩個類型,未來中遠(yuǎn)期覆被面積變化最明顯地集中在草地、林地和其他用地。也就是說,在隨著氣候變化、經(jīng)濟發(fā)展和政策的逐步推行之下,朋曲流域的林地面積將迅速擴大,而這部分的轉(zhuǎn)入是草地和其他用地的轉(zhuǎn)出獲得。

4.2 討論

不同發(fā)展情景下的土地覆被變化研究對區(qū)域的經(jīng)濟和可持續(xù)發(fā)展具有特別重要的影響,也對當(dāng)前朋曲流域的土地覆被空間格局變化具有良好的指導(dǎo)意義。朋曲流域作為珠峰自然保護區(qū)的主體部分,既要保障流域內(nèi)糧食生產(chǎn)安全,又要保證流域內(nèi)生物多樣性的有序發(fā)展。盡管流域內(nèi)的人口比重不大,當(dāng)前的糧食生產(chǎn)能滿足當(dāng)前的人口規(guī)模,但是未來人口的不斷增長會在一定程度上加重對糧食需求的壓力,所以耕地保護情景下的土地覆被變化研究具有非常之必要性。朋曲流域生態(tài)保護作為珠峰自然保護區(qū)的核心問題,生態(tài)意義重大,生態(tài)保護情景下的朋曲流域土地覆被研究可以進一步彌補生態(tài)建設(shè)基于現(xiàn)狀分析上的短板問題。雖然關(guān)于朋曲流域未來短期(2030年)和未來中遠(yuǎn)期(2050年)的模擬研究成果較好地擬合了未來土地發(fā)展?fàn)顩r,但是仍存在一些不足:①在使用模型中,未利用地與草地始終是變化最明顯的兩種覆被類型,這可能與模型中的馬爾可夫鏈的特性有關(guān),馬爾可夫鏈的引用雖解決了隨時間變化的土地覆被類型的預(yù)估問題,但作為一種基于概率論和數(shù)理統(tǒng)計的模型,在受全球變化與人類社會發(fā)展等一系列復(fù)雜機制影響的土地利用變化研究中,模型的預(yù)測效力還需進一步討論;②模型驗證中,實測結(jié)果與模擬結(jié)果kappa值為82%,雖滿足了模型要求的75%以上的精度,但對于小流域的預(yù)測來說,仍可能與實際產(chǎn)生較大的差距。這種偏差的產(chǎn)生可能與青藏高原特殊的自然環(huán)境條件和區(qū)域政策有關(guān),模型針對青藏高原進行區(qū)域性的改進和更多驅(qū)動因子選擇與其權(quán)重的研究將是下一步的工作方向;③FLUS模型在進行未來情景模擬預(yù)測時,所選取的驅(qū)動因子指標(biāo)具有一定的主觀性,驅(qū)動因子指標(biāo)的選取多參照前人研究成果中所選取的指標(biāo),因為所選取的研究區(qū)不同,受地勢、海拔、氣溫、降水、人口、經(jīng)濟、宗教等因素的影響,不同研究區(qū)內(nèi)對于土地利用的方式不同,不同地類發(fā)展的趨向也就各不相同,這就會在之后的模擬研究中增加模擬誤差,模擬結(jié)果與實際發(fā)展方向背道而馳。另外在選擇驅(qū)動因子時,無法準(zhǔn)確分析驅(qū)動因子選取的多少是否對模擬結(jié)果產(chǎn)生重要影響,一系列問題還需進一步模擬驗證。