基于土地利用變化的陜北多沙粗沙區(qū)生境質(zhì)量演變及響應(yīng)分析

安文舉， 魚亦凡， 郝姍姍， 王英弟， 李小兵， 麥旭東

（1.黃河水土保持綏德治理監(jiān)督局，陜西 榆林 718000；2.黃土高原水土保持與生態(tài)修復(fù)國家林業(yè)和草原局重點實驗室，陜西 榆林 718000）

土地利用和土地覆被變化（LULC）主要是指土地利用形態(tài)隨時間演變而發(fā)生的動態(tài)轉(zhuǎn)變過程［1］，是人類活動和自然生態(tài)環(huán)境相互作用的最終結(jié)果。近年來，由于人類活動強度加劇，土地利用及覆被變化與生態(tài)效應(yīng)已經(jīng)成為全球環(huán)境變化研究中的重要課題。作為生態(tài)系統(tǒng)的載體之一，土地利用變化會影響土壤侵蝕、能量交換、水循環(huán)和碳循環(huán)等生態(tài)過程，從而對生態(tài)系統(tǒng)服務(wù)產(chǎn)生影響［2］。生態(tài)服務(wù)系統(tǒng)是指通過生態(tài)結(jié)構(gòu)、生態(tài)過程和生態(tài)功能直接或間接得到的、能夠形成和維持人類生存和發(fā)展的環(huán)境條件和效果，關(guān)于該領(lǐng)域的研究，自Ehrlich在1983年系統(tǒng)介紹后，于1997年Costanza等開始根據(jù)效用價值和均衡價值理論等方法核算生態(tài)系統(tǒng)中的經(jīng)濟價值［3］，后逐漸成為生態(tài)領(lǐng)域研究的熱點。生境質(zhì)量作為評價生態(tài)系統(tǒng)服務(wù)的重要指標(biāo)之一受到國內(nèi)外學(xué)者廣泛研究。

生境質(zhì)量是指生態(tài)系統(tǒng)提供適宜個體與種群持續(xù)生存發(fā)展條件的能力，可以在一定程度上反應(yīng)出生態(tài)健康程度以及區(qū)域生物多樣性狀況［4］。目前，評價生境的方法主要有基于評價體系的定性法和利用模型計算的定量法［5］，隨著3S技術(shù)的不斷成熟，采用地理信息系統(tǒng)應(yīng)用程序（SolVES）［6］、棲息地適宜指數(shù)模型（HSI）［7］、生態(tài)系統(tǒng)服務(wù)和權(quán)衡綜合評估（InVEST）模型開展生境研究已經(jīng)成為精細(xì)化評價區(qū)域生態(tài)的重要方法，其中InVEST模型應(yīng)用尺度最為廣泛［8］。土地利用是生物多樣性的基礎(chǔ)，其變化作為人類活動最直接的表現(xiàn)方式，能夠改變生境斑塊間的能量和物質(zhì)轉(zhuǎn)換，最終使得生境質(zhì)量變化，影響生態(tài)環(huán)境和區(qū)域生物多樣性［9］。國內(nèi)外學(xué)者Zheng 等［10］、Chen 等［11］、Yu 等［12］、翟玉鑫等［13］、陳實等［14］、Tang等［15］、de Montaigu 等［16］、程靜等［17］分別從生境質(zhì)量對LULC 響應(yīng)機制及時空演變、生境質(zhì)量與城市化耦合作用以及考慮生物多樣性與生境質(zhì)量相關(guān)性等方面開展研究，為評價目標(biāo)區(qū)域生態(tài)狀況提供了數(shù)據(jù)支撐。

西北地區(qū)生態(tài)脆弱、氣候惡劣，近年來針對該區(qū)域的研究逐漸增多，研究廣度和深度亦有增強，但仍存在問題，一是目前研究集中于城市及周圍，關(guān)于城市發(fā)展中對生態(tài)威脅作用的研究較多［15］；二是對于生態(tài)帶的研究往往尺度較大［8,12-13］，難以準(zhǔn)確反映局部區(qū)域生態(tài)變化情況。陜北多沙粗沙區(qū)作為黃河泥沙主要來源區(qū)之一［18］，黃土深厚，植被稀少，水土流失嚴(yán)重，目前對于該區(qū)域的報道依然較少。近年來，陜北多沙粗沙區(qū)先后實施了退耕還林、飛播造林、高標(biāo)準(zhǔn)農(nóng)田建設(shè)等措施，期間又逢資源開發(fā)熱潮，作為陜西重要的能源化工基地，大量礦區(qū)的建設(shè)和能源開采導(dǎo)致土地利用穩(wěn)定性被打破。因此，本文以陜北多沙粗沙區(qū)為研究區(qū)，基于2000—2020 年土地利用遙感數(shù)據(jù)，分析其土地利用變化、生境質(zhì)量演變和生境質(zhì)量對土地利用轉(zhuǎn)移響應(yīng)規(guī)律，探究土地轉(zhuǎn)移過程中局部生態(tài)變化規(guī)律，以期為陜北多沙粗沙區(qū)生態(tài)恢復(fù)提供數(shù)據(jù)參考。

1 數(shù)據(jù)與方法

1.1 研究區(qū)概況

陜北多沙粗沙區(qū)（107°15′～111°13′E，35°19′～39°35′N）地處晉陜蒙寧甘5 省交界、黃河流域中上游、黃土高原北部，總面積為68266 km2（圖1）。該區(qū)域總體地勢為西高東低，平均海拔約1200～1500 m，多年平均降雨量約400～600 mm，平均蒸發(fā)量1200～1500 mm，屬半干旱氣候。該地區(qū)擁有豐富的礦產(chǎn)資源，對當(dāng)?shù)毓I(yè)和經(jīng)濟發(fā)展起到了重要作用，2022年榆林市GDP達6543.65×108元、延安市達2231.93×108元［19］。

圖1 陜北多沙粗沙區(qū)示意圖Fig.1 Schematic diagram of more sediments and coarse sediments region in northern Shaanxi

1.2 研究方法

1.2.1 土地利用動態(tài)變化2000—2020年土地利用數(shù)據(jù)來自中國科學(xué)院空天信息創(chuàng)新研究院（https://data.casearth.cn/）［20］，分辨率為30 m?；趪量臻g規(guī)劃和國家《土地利用現(xiàn)狀分類》2種分類原則及方法，結(jié)合陜北多沙粗沙區(qū)實際生態(tài)系統(tǒng)服務(wù)狀況，對土地利用數(shù)據(jù)調(diào)整歸類如表1所示。

表1 陜北多沙粗沙區(qū)土地利用分類Tab.1 Land use classification in more sediments and coarse sediments region of northern Shaanxi

利用ArcGIS 10.2 對陜北多沙粗沙區(qū)不同時期土地利用疊加，通過土地利用類型動態(tài)度［1］和土地利用轉(zhuǎn)移矩陣［21］來分析該區(qū)域土地利用變化過程。

土地利用轉(zhuǎn)移矩陣表示階段時間內(nèi)單一土地利用類型變化，其表達式［22］為：

式中：Cxy為第y類土地利用向地類x轉(zhuǎn)移的面積（km2）；x、y為轉(zhuǎn)換前后的土地利用類型，x=1, 2, 3,···,n，y=1,2,3,···,n。

土地利用動態(tài)度可以直觀描述土地利用轉(zhuǎn)型速率，分為單一動態(tài)度和綜合動態(tài)度，單一動態(tài)度是指確定地類轉(zhuǎn)化速率，其表達式為：

式中：Ki為研究區(qū)域內(nèi)土地利用i的動態(tài)度（%）；Ua、Ub分別為研究初、末期土地利用面積（km2）；T為年份初末期時間間隔（a）。

綜合動態(tài)度是描述研究區(qū)域內(nèi)全部地類轉(zhuǎn)移變化速率，用來描述區(qū)域土地變化穩(wěn)定程度，表達式為：

式中：M為綜合動態(tài)度（%）；ΔUi-j為地類i轉(zhuǎn)為j的面積（km2）；Ui為地類i初期面積（km2）；T為年份初末期時間間隔（a）。

積極發(fā)動群眾，對安全施藥技術(shù)進行集中培訓(xùn)學(xué)習(xí)，掌握病蟲害發(fā)生規(guī)律，在關(guān)鍵時期及時施藥，有效減少農(nóng)藥施藥次數(shù)，杜絕重復(fù)施藥及亂施藥現(xiàn)象，并做好對農(nóng)藥包裝袋的統(tǒng)一回收工作，在每個鄉(xiāng)鎮(zhèn)辦事處的作物集中區(qū)通過設(shè)立農(nóng)藥使用宣傳專欄及設(shè)立農(nóng)用物資垃圾回收箱、垃圾池等方式，減少本區(qū)農(nóng)藥施用總量，使本區(qū)農(nóng)業(yè)向生態(tài)、無公害、環(huán)保的方向健康發(fā)展。在新農(nóng)藥的篩選方面，通過田間藥效試驗，針對本區(qū)病蟲害發(fā)生種類篩選出高效、低毒、低殘留農(nóng)藥，并搞好宣傳指導(dǎo)工作，使農(nóng)戶掌握合理的施藥技術(shù)，降低農(nóng)殘及農(nóng)藥施用量。

1.2.2 InVEST 模型利用InVEST 模型中的生境模塊分析陜北多沙粗沙區(qū)生境質(zhì)量及其退化度。In-VEST 模型是由美國斯坦福大學(xué)、大自然保護協(xié)會（TNC）與世界自然基金會（WWF）聯(lián)合開發(fā)的綜合性生態(tài)系統(tǒng)服評價模型，已被中國、美國、瑞士等國家大學(xué)及科研機構(gòu)用于評估城市退化、動植物棲息地等多方面生境質(zhì)量［8］。其表達式為：

式中：Dxj為第j類土地利用中柵格x的生境退化程度；R為脅迫因子數(shù)量；wr為脅迫因子r的權(quán)重；Yr為脅迫因子r中的柵格數(shù)量；ry為圖層范圍脅迫因子強度；βx為柵格x單元的受保護水平；Sjr為生境j對脅迫因子r的敏感度大?。籭rxy為脅迫因子與生境間的距離衰減函數(shù)，其表達式為：

式中：dxy為柵格x、y之間的線性距離；drmax為脅迫因子r的最大線性影響距離。

最終生境質(zhì)量計算公式為：

式中：Qxj為生境質(zhì)量，其值介于[0,1]之間，值越大，生境質(zhì)量越高；Hj為生境適宜度；k為半飽和常數(shù)，一般取值為生境退化度最大值的一半；z為歸一化常數(shù)，一般取0.5。

按照InVEST模型使用手冊，需要設(shè)定生境脅迫因子及權(quán)重、不同生境類型對脅迫因子敏感度表方可計算。脅迫因子是表示該種地類對周圍生境的影響程度；權(quán)重值是指該地類對生境破壞性的相對程度，取值為[0,1]；相對敏感性是指每種地類對每種脅迫因子的相對敏感程度，取值范圍為[0,1]，越接近1即對該脅迫因子敏感度越高。本研究考慮陜北多沙粗沙區(qū)土地利用特征以及人為影響等因素，結(jié)合研究側(cè)重并查閱如陜北地區(qū)、蘭州地區(qū)、黃土高原等類似區(qū)域研究成果，并根據(jù)InVEST 模型手冊、國內(nèi)外相關(guān)研究，將旱地、建設(shè)用地、裸土地設(shè)為脅迫因子，并確定其權(quán)重、最大影響距離、衰退類型及敏感度（表2、表3）。

表2 脅迫因子及權(quán)重Tab.2 Threat factors and weights

表3 不同生境類型對脅迫因子敏感度Tab.3 Sensitivity of different habitat types to threat factors

1.2.3 土地利用轉(zhuǎn)移條件下生境貢獻率為探討生境演變對土地利用變化的響應(yīng)規(guī)律，采用ArcGIS分區(qū)統(tǒng)計將主要土地利用變化柵格與生境質(zhì)量變化柵格相對應(yīng)，提取同一地類變化下生境變化規(guī)律，采用生境貢獻率計算生境質(zhì)量對土地利用變化的響應(yīng)規(guī)律。生境貢獻率是指確定地類轉(zhuǎn)移前后生境質(zhì)量隨之變化比率［23］，其計算公式為：

式中：Rij為地類i轉(zhuǎn)化為j地類時的生境貢獻率（%）；Hi和Hj分別為研究時段內(nèi)地類i、j在變化前后的生境指數(shù)；Si為該變化地類i面積（km2）；St為研究區(qū)總面積（km2）。

2 結(jié)果與分析

2.1 陜北多沙粗沙區(qū)土地利用變化及轉(zhuǎn)移矩陣分析

圖2 2000—2020年土地利用變化Fig.2 Land use change during 2000—2020

2.1.2 土地利用轉(zhuǎn)移情況分析2000—2020 年，耕地、草地和裸地面積轉(zhuǎn)移變化較大（表4）。其中，2005—2010 年草地轉(zhuǎn)出與轉(zhuǎn)入面積分別為1545.5 km2、1484.9 km2，均是該年周期內(nèi)轉(zhuǎn)移面積最大地類，說明陜北多沙粗沙區(qū)草地發(fā)生較大流轉(zhuǎn)；2010—2015年，土地利用轉(zhuǎn)移面積較大的是草地和裸地，轉(zhuǎn)出面積分別為801.1 km2、714.0 km2；2015—2020 年，草地和裸地轉(zhuǎn)出面積最大，耕地轉(zhuǎn)入面積最大，原因可能是2019年大力推行高標(biāo)準(zhǔn)農(nóng)田建設(shè)導(dǎo)致。2000—2020 年，建設(shè)用地主要由草地轉(zhuǎn)入，2000—2005、2005—2010、2010—2015、2015—2020年 轉(zhuǎn) 入 面 積 分 別 為33.2 km2、47.8 km2、50.0 km2、

表4 2000—2020年土地利用轉(zhuǎn)移矩陣Tab.4 Land use transfer matrix during 2000—2020

67.0 km2，其次由裸地和耕地轉(zhuǎn)入。

2.1.3 土地利用變化動態(tài)度分析通過計算轉(zhuǎn)移變化面積較大的旱地、有林地、草地、建設(shè)用地、裸地5種二級地類土地利用變化單一動態(tài)度及所有地類綜合動態(tài)度，繪制變化趨勢95%置信區(qū)間散點圖（圖3），分析陜北多沙粗沙區(qū)土地利用變化速率。2000—2020年除旱地單一動態(tài)度呈增加狀態(tài)，有林地、草地、建設(shè)用地和裸地變化單一動態(tài)度均減小，裸地在近20 a 減少面積超過50%，表明陜北多沙粗沙區(qū)土地利用流轉(zhuǎn)速率逐漸降低。建設(shè)用地單一動態(tài)度在2000—2015年均呈降低狀態(tài)，但在2015—2020 年有所升高，達7.19%，這與2015 年以后大量自然資源的開采和發(fā)掘，城市建設(shè)投入加大的情況相吻合。綜合動態(tài)度從1.32%降至0.61%，說明土地利用變化速率減慢，陜北多沙粗沙區(qū)土地利用逐漸趨于穩(wěn)定。

圖3 土地利用動態(tài)度線性擬合Fig.3 Linear fitting of land use dynamics

2.2 陜北多沙粗沙區(qū)生境質(zhì)量動態(tài)變化情況

采用InVEST 生境質(zhì)量模塊計算陜北多沙粗沙區(qū)2000—2020年的生境質(zhì)量，并采用等間距法對其質(zhì)量從低到高分為低、較低、中、較高和高5個等級，繪制生境質(zhì)量等級分類圖。如圖4 所示，2000 年低等級生境主要分布在榆林市西北各縣，包括橫山區(qū)西北、靖邊北部等區(qū)域，這些區(qū)域是毛烏素沙漠所在地，生境質(zhì)量極差。高等級主要分布在寶塔區(qū)、宜川縣西南等區(qū)域，這些區(qū)域正處于子午嶺國家級自然保護區(qū)東陲，林地、草地面積占總土地利用面積的75.89%，建設(shè)用地僅占0.18%。2000—2010年，高等級生境質(zhì)量帶向志丹縣以北、寶塔以東等地區(qū)擴張。較低生境等級在2000和2020年分布面積最廣，2000年退耕還林前裸土地面積占比6.88%，到2020 年僅占2.95%，但建設(shè)用地面積增加較快，2000—2020 年共增加了447.05 km2，阻礙了生境質(zhì)量的恢復(fù)，使得較低生境質(zhì)量等級仍然分布較廣。

圖4 陜北多沙粗沙區(qū)生境質(zhì)量等級Fig.4 Habitat quality grades in more sediments and coarse sediments region of northern Shaanxi

結(jié)合圖5 可見，陜北多沙粗沙區(qū)生境質(zhì)量指數(shù)在2005 年達到峰值0.43 后逐漸降低至2020 年的0.37，表明該區(qū)域生境質(zhì)量整體較低。低等級面積逐年減少，較低等級面積逐年增加，高等級、較高等級和中等級變化趨勢與生境指數(shù)一致，均在2005年達到峰值后回落，但高于2000年水平。數(shù)據(jù)同樣反映出1999 年開始的生態(tài)恢復(fù)政策顯著提升了陜北多沙粗沙區(qū)的生境質(zhì)量，但后續(xù)城市建設(shè)和農(nóng)村開發(fā)導(dǎo)致這一成效減弱。

圖5 不同等級生境質(zhì)量面積及生境質(zhì)量指數(shù)變化Fig.5 Changes of habitat quality area and habitat quality indices at different levels

2.3 陜北多沙粗沙區(qū)生境退化度動態(tài)變化分析

生境退化度能夠反映生境受脅迫因子的影響程度，其取值在[0, 1]，值越大代表退化程度越高。將生境退化指數(shù)按照等間距法分為弱、較弱、中等、較強、強5個等級，從整體分布來看（圖6），陜北多沙粗沙區(qū)生境退化程度呈現(xiàn)中間強外圍弱的效果，米脂、子洲等中部以旱地為主的區(qū)域退化等級普遍較強，而生境質(zhì)量較高的志丹、宜川等縣區(qū)生境退化度較弱，這些區(qū)域覆蓋面積較大的草地、林地能減弱生境退化程度。

圖6 陜北多沙粗沙區(qū)生境質(zhì)量退化度等級Fig.6 Degradation grades of habitat quality in more sediments and coarse sediments region of northern Shaanxi

從陜北多沙粗沙區(qū)生境退化指數(shù)（表5）來看，2005—2015 年整體退化程度較弱，2005 年最低，為0.35，其后逐漸增加；2000 年和2020 年退化程度較高，說明后期生境下降的概率逐漸增大。結(jié)合生境質(zhì)量退化度空間分布（圖6）和土地利用來看，2015—2020年實行高標(biāo)準(zhǔn)農(nóng)田建設(shè)，草地→耕地面積增加，導(dǎo)致生境退化程度較強。總之生境質(zhì)量等級較高的區(qū)域退化程度一般較低（如寶塔區(qū)西南部等區(qū)域），反之退化程度較高（如綏德縣中部、米脂縣南部等區(qū)域）。

表5 陜北多沙粗沙區(qū)生境退化指數(shù)Tab.5 Habitat degradation indices in more sediments and coarse sediments region of northern Shaanxi

2.4 生境質(zhì)量對土地利用變化的響應(yīng)關(guān)系分析

陜北多沙粗沙區(qū)主要土地利用生境貢獻率（表6）結(jié)合土地利用轉(zhuǎn)移矩陣（表4）反映出生境質(zhì)量變化與土地利用轉(zhuǎn)移關(guān)系密切。按照生境增減來分析，生境增加主要是耕地向林地和草地轉(zhuǎn)移導(dǎo)致植被面積增加，例如2000—2005 年耕地→林地、耕地→草地轉(zhuǎn)移面積分別為156.35 km2、1192.47 km2，生境質(zhì)量分別增加0.0693%、0.3495%。裸土地作為生境質(zhì)量脅迫因子之一，其向耕地、林地、草地轉(zhuǎn)移均造成生境質(zhì)量增加，2000—2005 年，裸地→草地轉(zhuǎn)移面積達1329.65 km2，生境質(zhì)量增加0.90%，其后逐漸減少。生境減少主要是耕地向建設(shè)用地和裸地轉(zhuǎn)移，林地向其他地類轉(zhuǎn)移也會造成生境下降，例如2000—2005、2005—2010、2010—2015、2015—2020年林地→草地面積分別為81.63 km2、84.46 km2、3.10 km2、18.58 km2，其生境質(zhì)量分別下降0.0018%、0.0178%、0.0006%、0.0037%，由此看出林地面積生境質(zhì)量聯(lián)系緊密。草地向耕地、建設(shè)用地和裸地轉(zhuǎn)移引起生境質(zhì)量下降，尤其是向裸地轉(zhuǎn)移面積保持較高水平，削弱了生境質(zhì)量的恢復(fù)水平。綜上，生境質(zhì)量對土地利用類型的轉(zhuǎn)移較為敏感，在研究區(qū)內(nèi)，退耕還林政策有效降低了裸地面積使生境質(zhì)量提升，但同時自然資源的開采也導(dǎo)致了一部分林地和草地轉(zhuǎn)成裸地和建設(shè)用地，為生境恢復(fù)帶來威脅。

表6 陜北多沙粗沙區(qū)主要土地利用生境貢獻率Tab.6 Contribution rates of major land use habitats in more sediments and coarse sediments region of northern Shaanxi

3 討論

3.1 土地利用類型轉(zhuǎn)移分析

陜北多沙粗沙區(qū)土地利用以耕地、有林地、草地、建設(shè)用地和裸地為主，同時也是轉(zhuǎn)移變化過程中的主要地類，旱地極易轉(zhuǎn)為草地，裸地極易轉(zhuǎn)為旱地、草地，草地極易轉(zhuǎn)為旱地和林地，這與該地區(qū)的生產(chǎn)生活習(xí)慣［24］以及長時效的生態(tài)政策［25］密不可分。1999 年實行退耕還林的生態(tài)政策增大了草地和林地面積，也促進了裸地向草地、林地轉(zhuǎn)化，生態(tài)狀況好轉(zhuǎn)。但隨著城市發(fā)展需求不斷增大，林地、草地面積逐年減少，2000—2005 年，耕地→草地、耕地→林地面積分別為1192.5 km2、156.4 km2，林地、草地→耕地面積總計僅755.5 km2，草地→建設(shè)用地面積為33.2 km2，同樣轉(zhuǎn)化關(guān)系的地類在2015—2020 年分別變?yōu)?98.2 km2、25.1 km2、740.7 km2和67.0 km2，這與冀泳帆等［26］通過計算發(fā)現(xiàn)林地和草地向耕地轉(zhuǎn)移占比較大、建設(shè)用地面積激增，導(dǎo)致地區(qū)發(fā)展不均衡的研究一致。陳樂等［27］研究表明耕地擴張雖然能緩解綠洲與荒漠交錯帶的干旱程度，但侵占了草地覆蓋，不利于草地生態(tài)穩(wěn)定。所以，如何解決土地利用協(xié)同發(fā)展，優(yōu)化土地利用布局是陜北多沙粗沙區(qū)目前乃至將來著力解決的問題。

3.2 生境質(zhì)量及退化度演變分析

從生境質(zhì)量來看，毗鄰子午嶺國家生態(tài)保護區(qū)的陜北多沙粗沙區(qū)西南部生境質(zhì)量最高，其次毛烏素沙漠所在地靖邊縣生境質(zhì)量由2000 年的低等級為主、較低等級次之逐漸恢復(fù)到中級，并有零星較高和高等級出現(xiàn)，分析可能是沙漠飛播造林技術(shù)使得毛烏素沙漠區(qū)域的林地、草地面積增加，生境質(zhì)量恢復(fù)。從生境退化度來看，旱地分布廣泛的中部區(qū)域退化度明顯高于外圍地帶，而林地覆蓋面廣的西南部綠地帶生境退化度普遍較低。從生境質(zhì)量對土地利用的響應(yīng)方面來看，生境質(zhì)量對土地利用轉(zhuǎn)移變化較為敏感，轉(zhuǎn)化成林地、草地、濕地和水體能夠有效提升生境質(zhì)量，建設(shè)用地、裸地和耕地則會對生境質(zhì)量的提升造成阻礙，這與胡豐等［28］針對渭河流域生境質(zhì)量的研究結(jié)果一致，后續(xù)建設(shè)用地開發(fā)、耕地擴張應(yīng)該引起關(guān)注。

3.3 后續(xù)生態(tài)保護及研究建議

相較于對大范圍研究區(qū)開展生態(tài)質(zhì)量評價，劃定局部區(qū)域的研究能更精準(zhǔn)地反映其生態(tài)質(zhì)量以及存在的問題，以便于因地制宜開展保護措施。為進一步促進陜北多沙粗沙區(qū)生態(tài)恢復(fù)，應(yīng)從以下2個方面著手：一是堅持當(dāng)前生態(tài)政策，繼續(xù)推行退耕還林還草、飛播造林等措施，擴大有效林地覆蓋持續(xù)增加；二是耕地、建設(shè)用地擴張應(yīng)適度，既要制止耕地向外轉(zhuǎn)移，也要防止耕地侵犯林草、水域等地類，工礦、資源生產(chǎn)用地擴張應(yīng)遵循“邊保護邊開發(fā)”原則，適當(dāng)擴大濕地和水體用地。后期對于具體地域的生態(tài)評價，建議按照生態(tài)邊界、生產(chǎn)生活習(xí)慣等劃定小區(qū)域開展研究；另外，盡管InVEST 模型已經(jīng)進行多次迭代并廣泛運用，但因其中相關(guān)參數(shù)的設(shè)定主要來源于專家評分及參考類似研究區(qū)的文獻設(shè)定，后續(xù)研究應(yīng)當(dāng)進一步優(yōu)化參數(shù)設(shè)定以及全面考慮影響生境質(zhì)量的其他因子作用。

4 結(jié)論

（1）陜北多沙粗沙區(qū)土地利用主要以旱地、有林地、草地、建設(shè)用地和裸土地為主，2000—2020年裸土地、草地、旱地變化面積最大，濕地、其他林地、水體變化面積占比較低。

（2）2000—2020 年裸地轉(zhuǎn)出面積最大，總降幅達57.11%，建設(shè)用地總增幅達272.69%，綜合土地利用動態(tài)度從1.32 降至0.61，區(qū)域土地利用變化趨于穩(wěn)定。

（3）生境質(zhì)量指數(shù)呈先增后減趨勢，峰值達0.43，整體偏低，耕地擴張侵犯綠地面積，威脅生境質(zhì)量恢復(fù)；2000、2020年退化程度較高，退化指數(shù)分別為0.41和0.43，未來生境有惡化趨勢。