屈雯，謝保鵬，陳英，賀瑩，馬麗婭

（甘肅農(nóng)業(yè)大學(xué)管理學(xué)院，甘肅 蘭州 730070）

生境質(zhì)量是生態(tài)系統(tǒng)為個(gè)體和種群提供適宜、持續(xù)性生存條件的能力［1］，是區(qū)域生物多樣性及生態(tài)服務(wù)水平的重要反映［2-3］。近年來(lái)，隨著人為活動(dòng)的不斷加強(qiáng)，使得土地利用方式、利用強(qiáng)度以及景觀格局發(fā)生了巨大的變化［4］，導(dǎo)致區(qū)域生態(tài)系統(tǒng)結(jié)構(gòu)和功能的變化，進(jìn)而引發(fā)了生物多樣性銳減、生境質(zhì)量下降、生境退化等一系列生態(tài)環(huán)境問(wèn)題［3］。地形作為自然環(huán)境中的重要影響因子，對(duì)地表物質(zhì)的遷移與能量的轉(zhuǎn)換有著顯著影響［2］，研究表明地形因子是影響土地利用景觀格局變化的重要因素［3，5］，同時(shí)也是影響生境系統(tǒng)物質(zhì)交換和能量循環(huán)的重要因素［4］。因此，從地形梯度效應(yīng)角度出發(fā)，對(duì)生境質(zhì)量及生境退化進(jìn)行研究，將有助于準(zhǔn)確的識(shí)別生境質(zhì)量的變化效應(yīng)，對(duì)于挖掘區(qū)域生境質(zhì)量的地形信息圖譜具有重要的意義。目前，國(guó)內(nèi)外學(xué)者進(jìn)行了大量關(guān)于生境質(zhì)量的評(píng)估研究，國(guó)外學(xué)者側(cè)重于研究生境斑塊對(duì)生物多樣性的影響［6-9］，國(guó)內(nèi)學(xué)者側(cè)重于研究生境質(zhì)量的時(shí)空變化規(guī)律［10］?？傮w而言，生境質(zhì)量的評(píng)估方法分為兩大類(lèi)［1］，一是實(shí)地調(diào)查法，通過(guò)對(duì)研究區(qū)進(jìn)行實(shí)地調(diào)研，構(gòu)建評(píng)價(jià)指標(biāo)對(duì)生境質(zhì)量進(jìn)行綜合評(píng)價(jià)；二是基于土地利用數(shù)據(jù)，采用模型對(duì)研究區(qū)生境質(zhì)量進(jìn)行評(píng)估，使用最普遍的是InVEST模型，此類(lèi)研究尺度涉及全國(guó)［10］、?。?1-13］、市［14-15］縣［16-17］、流域［18-19］、城市群［20］、自然保護(hù)區(qū)［21-22］等，研究?jī)?nèi)容多從生境質(zhì)量時(shí)空變化［4，23］、影響因素［24-25］、生境質(zhì)量與土地整治［26］、生境質(zhì)量與生態(tài)系統(tǒng)服務(wù)價(jià)值［27］等方面展開(kāi)。

位于甘肅省中部地區(qū)的榆中縣是東部季風(fēng)區(qū)、西北干旱半干旱區(qū)、青藏高原區(qū)三大自然區(qū)的交匯區(qū)，其地形地勢(shì)和土地利用方式具有明顯的過(guò)渡性特征［28］，特殊的地理位置使其生境具有顯著的地域性特征。近年來(lái)，隨著蘭州市的不斷發(fā)展，其作為蘭州市“東擴(kuò)”的主要承載區(qū)，經(jīng)濟(jì)發(fā)展勢(shì)頭迅猛，人口不斷增加，資源與環(huán)境約束日益凸顯，因此，以榆中縣作為研究區(qū)，基于土地利用數(shù)據(jù)、DEM數(shù)據(jù)，借助InVEST模型評(píng)估榆中縣生境退化和生境質(zhì)量，以期為榆中縣生境保護(hù)、生態(tài)保護(hù)及國(guó)土空間規(guī)劃編制提供理論參考。

1 材料與方法

1.1 研究區(qū)概況

榆中縣位于甘肅省中部地區(qū)，是省會(huì)蘭州市的東部郊區(qū)（E 103°49′15″~104°34′40″，N 35°34′20″~36°26′30″），縣域海拔1 369~3 670 m，總面積為32 946.69 km2。地處黃土高原西部，同時(shí)也是西北干旱區(qū)、青藏高寒區(qū)和東部季風(fēng)區(qū)的交匯地帶［28］?？h域地形復(fù)雜，南部為高寒陰濕區(qū)，馬銜山和興隆山從西北向東南橫列整個(gè)縣域的南部，該區(qū)降水充沛，植被茂盛，森林資源豐富，是榆中縣主要的生態(tài)保障區(qū)；中部為興隆山山前沖積、洪積平原和宛川河谷地組成的川塬河谷平原區(qū)，地勢(shì)平坦，土壤肥沃，交通便利，是榆中縣主要的工業(yè)和農(nóng)業(yè)生產(chǎn)區(qū)；北部為干旱山區(qū)，降水稀少，植被稀疏，是榆中縣水土流失最為嚴(yán)重的區(qū)域。

1.2 數(shù)據(jù)來(lái)源

（1）榆中縣DEM數(shù)據(jù)來(lái)源于中國(guó)科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心，其分辨率為30 m。

（2）榆中縣2008年和2018年土地利用柵格數(shù)據(jù)來(lái)源于榆中縣土地利用變更數(shù)據(jù)庫(kù)，分辨率為30 m。根據(jù)研究需要，按照《土地利用現(xiàn)狀分類(lèi)GB/T 21010—2017》的分類(lèi)依據(jù)［29］和土地利用/土地覆被遙感監(jiān)測(cè)數(shù)據(jù)分類(lèi)標(biāo)準(zhǔn)［10］，利用ArcGIS 10.2軟件平臺(tái)對(duì)變更數(shù)據(jù)中的地類(lèi)進(jìn)行進(jìn)一步劃分，分為耕地、園地、林地、草地、水域、建設(shè)用地和未利用地7大類(lèi)。

1.3 研究方法

1.3.1 基于In VEST模型的生境質(zhì)量和生境退化評(píng)估 In VEST模型中開(kāi)展生境質(zhì)量評(píng)價(jià)首先要計(jì)算區(qū)域的生境退化度，其計(jì)算公式如下：

式中：Dxj是生境退化指數(shù)，R是脅迫因子的數(shù)量；Y r是脅迫因子r在土地利用柵格圖中的總個(gè)數(shù)；ry是土地利用柵格單元上的脅迫因子個(gè)數(shù)；ωr是脅迫因子的權(quán)重；βx是柵格x的可達(dá)性水平；Sjr是土地利用類(lèi)型j對(duì)脅迫因子r的敏感性程度，取值范圍為0~1；irxy是脅迫因子的影響距離；d xy為柵格x與y之間的距離；d rmax是脅迫因子r的最大影響范圍。

生境質(zhì)量的計(jì)算公式為：

式中：Qxj為第j種土地利用類(lèi)型中柵格單元x的生境質(zhì)量；H j為第j種土地利用類(lèi)型的生境適宜；z為尺度常數(shù)是模型的默認(rèn)參數(shù)；k為半飽和常數(shù)。

In VEST模型運(yùn)行時(shí)需要輸入?yún)?shù)，本研究在對(duì)InVEST模型手冊(cè)及相關(guān)研究［11-12，27-29］分析的基礎(chǔ)之上，參數(shù)設(shè)置具體參考劉春芳等［30］、賀瑩等［31］的研究，參數(shù)如表1~2所示。

表1 研究區(qū)內(nèi)各威脅因子的參數(shù)Table 1 Parameters of threat factors in the study area

1.3.2 圖譜分析 文章借助土地利用圖譜分析的相關(guān)原理［12，32-34］，構(gòu)建生境圖譜分析模型，以記錄和分析不同時(shí)空條件下生境質(zhì)量及生境退化的復(fù)合變化信息。計(jì)算公式如下：

式中：A為整個(gè)研究期內(nèi)生境質(zhì)量或生境退化圖譜單元變化類(lèi)型；A 1、A 2為研究初期和研究末期的生境質(zhì)量或生境退化圖譜柵格狀態(tài)編碼。

1.3.3 地形梯度分析 地形條件在很大程度上會(huì)影響土地利用類(lèi)型的分布［35，2］，本研究引進(jìn)復(fù)合地形因子，即地形位梯度［36］，以綜合反映地形條件與生境質(zhì)量、生境退化之間的關(guān)系。計(jì)算公式如下：

式中：T表示地形位梯度指數(shù)，E表示研究區(qū)內(nèi)某一點(diǎn)的高程，S表示研究區(qū)某一點(diǎn)的坡度，Eˉ表示研究區(qū)的平均高程，Sˉ表示研究區(qū)的平均坡度。

由于不同地形位梯度下生境質(zhì)量、生境退化面積存在較大的差異，本文為消除這種面積大小差異帶來(lái)的量綱影響［36］，引進(jìn)地形分布指數(shù)，對(duì)榆中縣生境質(zhì)量及生境退化的地形位梯度效應(yīng)進(jìn)行研究。參照周潤(rùn)芳等［37］對(duì)榆中縣地形地貌的研究，將地形位梯度分為5級(jí)，其計(jì)算公式如下：

式中：P為分布指數(shù)，Sie為第i級(jí)生境質(zhì)量或生境退化度在e級(jí)地形因子等級(jí)下的面積，Si為第i級(jí)生境質(zhì)量或生境退化度的總面積，Se為第e級(jí)地形因子所占的面積，S為研究區(qū)的總面積。

表2 各地類(lèi)對(duì)脅迫因子的敏感度參數(shù)Table 2 Sensitivity parameters of each category to stress factors

2 結(jié)果與分析

2.1 生境質(zhì)量和生境退化度的時(shí)空變化特征

2.1.1 生境質(zhì)量時(shí)空變化特征 為了準(zhǔn)確識(shí)別榆中縣生境質(zhì)量的變化規(guī)律，研究利用ArcGIS軟件中的自然斷點(diǎn)法將生境質(zhì)量分為低值區(qū)（1級(jí)）、較低值區(qū)（2級(jí)）、中值區(qū)（3級(jí)）、較高值區(qū)（4級(jí)）、高值區(qū)（5級(jí)）5個(gè)等級(jí)。

由圖1可以看出，2008~2018年，榆中縣生境質(zhì)量的空間格局整體上較為穩(wěn)定，從南到北總體呈現(xiàn)出由西南向東北遞減的分布特點(diǎn)。生境質(zhì)量高值區(qū)域主要分布在縣域南部的興隆山、馬銜山沿線，區(qū)內(nèi)人為活動(dòng)稀少，常年降水充足，植被茂密，土地利用多以草地、林地為主，生境質(zhì)量整體上優(yōu)于其他區(qū)域。生境質(zhì)量較高值區(qū)域主要分布在生境質(zhì)量高值區(qū)域的外圍，除此之外在貢井林場(chǎng)有小部分分布。生境質(zhì)量中值區(qū)域整體上分布零散，在園子岔鄉(xiāng)、貢井鄉(xiāng)、中連川鄉(xiāng)、清水驛鄉(xiāng)、韋營(yíng)鄉(xiāng)等個(gè)別鄉(xiāng)鎮(zhèn)分布較為集中，除此之外在宛川河沿岸河谷兩側(cè)也有零星分布。生境質(zhì)量較低值區(qū)域位于北部山區(qū)和中部平原地區(qū)，該區(qū)土地利用主要是耕地和其他草地、農(nóng)村居民點(diǎn)。低值區(qū)域位于縣域中部和北部，該區(qū)域分布大量的居民點(diǎn)用地、城鎮(zhèn)用地以及工礦用地，人類(lèi)活動(dòng)頻繁，因此生境質(zhì)量普遍偏低。由表3可以看出，2008年和2018年研究區(qū)平均生境質(zhì)量分別為0.228 9和0.226 6，生境質(zhì)量總體上略有下降。研究期間1等生境占比由46.29%變?yōu)?7.01%，結(jié)合圖2可以發(fā)現(xiàn)，變化區(qū)域主要是由于受到建設(shè)用地需求增加的影響，人類(lèi)對(duì)區(qū)域生態(tài)環(huán)境有所干預(yù)，破壞了原有的自然生境，使得生境質(zhì)量略有退化；2等生境占比由38.74%變?yōu)?8.06%，主要是由于研究期間耕地面積減少的影響；3等、4等、5等生境總體上保持相對(duì)穩(wěn)定的狀態(tài)。

圖1 2008~2018年榆中縣生境質(zhì)量空間分布圖Figure 1 The spatial distribution of habitat quality in Yuzhong County from 2008 to 2018

表3 2008~2018年榆中縣生境質(zhì)量變化情況Table 3 Habitat quality changes in Yuzhong County from 2008 to 2018

圖2 榆中縣土地利用類(lèi)型Figure 2 Land use types in Yuzhong County

2.1.2 生境退化度時(shí)空變化特征 由圖3可以看出，2008~2018年，榆中縣生境退化度的空間格局整體上較為穩(wěn)定，并且表現(xiàn)出明顯的空間差異，具體表現(xiàn)為以宛川河為界，呈現(xiàn)“河北低而河南高”的分布特點(diǎn)。生境退化度最嚴(yán)重的區(qū)域?yàn)橹胁看▍^(qū)的和平鎮(zhèn)、定遠(yuǎn)鎮(zhèn)等重點(diǎn)鄉(xiāng)鎮(zhèn)，這些區(qū)域是榆中縣主要的城鎮(zhèn)及工農(nóng)業(yè)發(fā)展區(qū)，如圖2所示，隨著縣域人為活動(dòng)日益增強(qiáng)，對(duì)居住地、工礦用地以及交通運(yùn)輸用地等建設(shè)用地的需求快速增加，使原有的自然生境遭到破壞，導(dǎo)致生境破碎度加劇、連通性變差，最終使得這些區(qū)域的生境質(zhì)量下降；整個(gè)研究期間縣域生境退化程度有所增加，平均生境退化度由2008年的0.007 370變?yōu)?018年的0.007 389。具體而言，研究期間生境退化度高值區(qū)域進(jìn)一步擴(kuò)散，2008年和2018年的占比分別為0.15%、0.49%，說(shuō)明榆中縣中部川區(qū)存在生態(tài)退化的風(fēng)險(xiǎn)；生境退化度較低值區(qū)有所增加，占比由2008年的67.35%變?yōu)?018年的74.26%，主要分布在縣域北部石質(zhì)低山丘陵區(qū)，主要是由于榆中縣近年來(lái)實(shí)施了一系列“退耕還林”工程和“生態(tài)北山”建設(shè)工程，加強(qiáng)了對(duì)北部山區(qū)生態(tài)環(huán)境的保護(hù)；生境退化較高值區(qū)域有所縮減，由2008年的9.20%變?yōu)?018年的7.54%，變化主要發(fā)生在縣域中部，轉(zhuǎn)化成了生境退化中值區(qū)；生境退化中值區(qū)域也有所縮減，由2008年的23.30%變?yōu)?018年的7.54%，變化主要集中在縣域北部，主要轉(zhuǎn)化成了生境退化低值區(qū)域。

2.2 生境質(zhì)量及生境退化的圖譜分析

按照前文分析，將生境質(zhì)量5級(jí)區(qū)間的區(qū)間編碼設(shè)置為1、2、3、4、5，分別對(duì)應(yīng)生境質(zhì)量低值、較低值、中值、較高值和高值5種狀態(tài)。由圖4可以看出，榆中縣生境質(zhì)量轉(zhuǎn)移圖譜共有25種類(lèi)型，其中，有5類(lèi)生境質(zhì)量圖譜單元未發(fā)生變化，其他20類(lèi)生境質(zhì)量圖譜單元或多或少都發(fā)生了變化。整體而言，未發(fā)生變化的圖譜編碼為11（表示研究期初和期末的生境質(zhì)量狀態(tài)都為低值）、22、33、44、55總占比為71.45%，其中，低-低（編碼為11）的圖譜單元最為穩(wěn)定，占未變化單元的49.75%，主要分布在縣域北部和中部偏西的區(qū)域，和草地、未利用地的分布基本一致，中-中（編碼為33）的圖譜單元最不穩(wěn)定，占未變化單元的2%，分布比較零散，主要分布在縣域中東部區(qū)域。發(fā)生變化的圖譜占總圖譜的比例為28.55%，其中，“低-較低”和“較低-低”（編碼為12和21）圖譜變化最為明顯，由圖3可以看出，這兩類(lèi)圖譜主要分布在榆中縣北部、東部，以及宛川河以南興隆山以北的區(qū)域，縣域北部和東部分布著大量的耕地，隨著人類(lèi)耕地利用程度的不斷加深，使得原有的耕地景觀遭到破壞，生境質(zhì)量有所下降，宛川河以南興隆山以北的區(qū)域地勢(shì)平坦，土壤肥沃，是主要的生產(chǎn)生活場(chǎng)所，尤其是近幾年，縣域人口增加導(dǎo)致建設(shè)用地需求增長(zhǎng)，從而影響了生境變化。由圖4可以看出，榆中縣生境退化轉(zhuǎn)移圖譜共有19種類(lèi)型，其中，有4類(lèi)生境退化圖譜單元未發(fā)生變化，其他15類(lèi)生境退化圖譜單元或多或少都發(fā)生了變化。整體而言，未發(fā)生變化的圖譜編碼為22、33、44、55，總占比為81.10%，其中，較低-較低（編碼為22）的圖譜單元最為穩(wěn)定，占未變化單元的76.93%。發(fā)生變化的圖譜占總圖譜的比例為18.90%，其中，中-較低（編碼為32）圖譜變化最為明顯，集中分布在縣域北部，主要是“北山生態(tài)工程”的建設(shè)使得該區(qū)中等退化轉(zhuǎn)化成了較低程度的退化。

圖4 2008~2018年榆中縣生境質(zhì)量及生境退化的圖譜空間分布Figure 4 The spatial distribution map of habitat quality and habitat degradation in Yuzhong County from 2008 to 2018

表4 2008~2018年榆中縣生境退化變化情況Table 4 Habitat degradation changes in Yuzhong County from 2008 to 2018

圖3 2008~2018年榆中縣生境退化度空間分布Figure 3 The spatial distribution of habitat degradation in Yuzhong County from 2008 to 2018

2.3 生境質(zhì)量的地形梯度效應(yīng)

2.3.1榆中縣的地形梯度特征 研究上述公式（5）計(jì)算榆中縣地形位梯度，結(jié)果區(qū)間為0.22~1.09，采用自然斷點(diǎn)法將其分成5個(gè)等級(jí)，等級(jí)1表示地形位指數(shù)在0.22~0.42，等級(jí)2表示地形位指數(shù)在0.42~0.54，等級(jí)3表示地形位指數(shù)在0.54~0.64，等級(jí)4表示地形位指數(shù)在0.64~0.75，等級(jí)5表示地形位指數(shù)在0.75~1.09，等級(jí)越低表示地形越平坦。

由圖5可以看出，榆中縣地形位指數(shù)呈現(xiàn)出“南部高、中部低、北部高”的分布格局，與縣域南高北低，中間低凹的地勢(shì)基本一致。榆中縣地形位以3級(jí)為主，占比為30.80%左右，其次是4級(jí)、2級(jí)、1級(jí)，占比分別為24.23%、16.94%、15.55%左右，5級(jí)地形位占比最小，為12.48%左右。由表5可知，榆中縣1等地形位和2等地形位的土地利用類(lèi)型以耕地為主，占比達(dá)60%左右，3等、4等及5等地形位以未利用地為主，具體而言耕地、園地、水域、建設(shè)用地在1等地形位上占比最大，林地、草地在5等地形位上占比最大，未利用地在4等地形位上占比最大。隨著地形位指數(shù)的增加，耕地、園地、水域、建設(shè)用地的占比出現(xiàn)下降趨勢(shì)，相反，林地、草地、未利用地的占比出現(xiàn)了增加的趨勢(shì)。主要是由于人類(lèi)活動(dòng)強(qiáng)度最大的區(qū)域一般都是地形位低的位置，地形位指數(shù)值低表明區(qū)域自然環(huán)境越有利于人類(lèi)生存，人類(lèi)對(duì)土地的利用程度越強(qiáng)，相反，隨著地形位指數(shù)的增加，人類(lèi)對(duì)自然環(huán)境的改造程度越小，因此土地利用多以林地、草地以及未利用地為主。

表5 不同地形位梯度的土地利用特征Table 5 Land use characteristics of different terrain position gradients

圖5 榆中縣地形位梯度分布情況Figure 5 Topographical position gradient distribution in Yuzhong County

2.3.2 生境質(zhì)量及生境退化的地形位梯度分布特征地形要素是影響生物種群分布格局以及維持生物多樣性的一個(gè)重要因素［4］，地形分異性會(huì)影響土地利用類(lèi)型進(jìn)而也會(huì)在很大程度上影響生境質(zhì)量及生境退化的空間分布格局。本文在計(jì)算了研究區(qū)地形位指數(shù)的基礎(chǔ)上，采用分布指數(shù)來(lái)描述生境質(zhì)量、生境退化在不同地形位梯度上的分布特征，以揭示綜合地形變化對(duì)生境質(zhì)量和生境退化的影響。如圖6、7所示，生境質(zhì)量和生境退化在地形位梯度上呈現(xiàn)出較大的差異。具體而言，低質(zhì)量生境在第3、第4級(jí)地形位梯度上的分布指數(shù)大于1，占絕對(duì)優(yōu)勢(shì)，較低質(zhì)量生境在第1、第2、第3級(jí)地形位梯度上的分布指數(shù)大于1，占絕對(duì)優(yōu)勢(shì)，中等質(zhì)量生境在第3、第4級(jí)地形位梯度上的分布指數(shù)大于1，占絕對(duì)優(yōu)勢(shì)。較高質(zhì)量生境及高質(zhì)量生境在第4、第5級(jí)地形位梯度上的分布指數(shù)大于1，占絕對(duì)優(yōu)勢(shì)。生境退化低值在第1級(jí)地形位梯度上的分布指數(shù)大于1，占絕對(duì)優(yōu)勢(shì)，生境退化較低值在第3、第4、第5級(jí)地形位梯度上的分布指數(shù)大于1，占絕對(duì)優(yōu)勢(shì)，生境退化中等區(qū)域在第1、第2、第3級(jí)地形位梯度上的分布指數(shù)大于1，占絕對(duì)優(yōu)勢(shì)，較高生境退化和高生境退化第1級(jí)地形位梯度上的分布指數(shù)大于1，占絕對(duì)優(yōu)勢(shì)，因此，第3級(jí)地形位梯度是生境質(zhì)量及生境退化的混合梯度。

圖6 生境質(zhì)量在地形梯度上的分布指數(shù)Figure 6 Distribution index of habitat quality on topographic gradient

3 討論

本研究應(yīng)用InVEST模型和地形位指數(shù)對(duì)榆中縣生境質(zhì)量進(jìn)行了比較充分的分析，其中生境質(zhì)量的時(shí)空分布結(jié)果與劉春芳等［11，28，30］的研究結(jié)果相近，說(shuō)明本文的研究結(jié)果比較符合榆中縣的實(shí)際情況，應(yīng)該能反映近年來(lái)榆中縣生境質(zhì)量及生境退化的變化，而且本研究首次從地形視角對(duì)榆中縣的生境質(zhì)量及生境退化的分布特點(diǎn)進(jìn)行了分析，較好地揭示了研究區(qū)生境質(zhì)量的地形梯度效應(yīng)，研究成果可為榆中縣生態(tài)環(huán)境保護(hù)、生物多樣性維護(hù)提供理論依據(jù)。但研究也存在以下的不足：（1）研究在對(duì)生境質(zhì)量進(jìn)行評(píng)估時(shí)，應(yīng)用的是傳統(tǒng)的評(píng)估模型，缺乏方法的創(chuàng)新性，而且模型計(jì)算中參數(shù)設(shè)置可能存在一定的主觀性；（2）生境質(zhì)量和生境退化會(huì)受到多種因素的影響，本文只是選取了較為重要的地形因素，對(duì)人為因素的研究還需進(jìn)一步探討。

圖7 生境退化度在地形梯度上的分布指數(shù)Figure 7 Distribution index of habitat degradation degree on topographic gradient

榆中縣是西北地區(qū)重要的生態(tài)屏障區(qū)，作為三大自然帶的交匯地帶具有明顯的地形分異特征，區(qū)域生境質(zhì)量及生境退化受到地形變化特征的影響?？h域南部地勢(shì)高聳，分布大量的天然植被，林地、草地是該區(qū)的土地利用類(lèi)型，是縣域生境質(zhì)量的高值區(qū)；中部宛川河流域地勢(shì)低平，是近年來(lái)蘭州市東擴(kuò)的主要承載區(qū)，同時(shí)該區(qū)域是主要的城鄉(xiāng)建設(shè)用地增加區(qū)，隨著人類(lèi)活動(dòng)日益加強(qiáng)，導(dǎo)致原有的土地利用結(jié)構(gòu)、生態(tài)結(jié)構(gòu)和功能遭到破壞，從而使得此處生境退化程度不斷加深，生境質(zhì)量不斷下降；北部地區(qū)系主要的黃土丘陵區(qū)，區(qū)域自然環(huán)境較差，分布大量的草地和未利用地以及部分農(nóng)村居民點(diǎn)。因此，在南部山區(qū)，應(yīng)該加強(qiáng)對(duì)森林的保護(hù)，適當(dāng)?shù)脑黾恿植莞采w度，同時(shí)有序控制人類(lèi)活動(dòng)，以保持其生態(tài)功能；在中部川區(qū)，應(yīng)該限制建設(shè)用地增長(zhǎng)，同時(shí)要遏制生境退化，對(duì)用地結(jié)構(gòu)進(jìn)行適當(dāng)調(diào)整，尤其是園地要盡可能的降低其利用強(qiáng)度；北部山區(qū)，一方面要適當(dāng)?shù)脑黾恿植莞采w率，另一方面要控制現(xiàn)有的居民點(diǎn)用地，盡可能使該區(qū)生態(tài)保持自然的狀態(tài)，以使其進(jìn)行必要的自我調(diào)節(jié)和恢復(fù)。

4 結(jié)論

本研究利用DEM數(shù)據(jù)、2008年和2018年土地利用柵格數(shù)據(jù)，基于InVEST模型和地形位指數(shù)，評(píng)估了榆中縣生境質(zhì)量時(shí)空的變化特征及其地形梯度效應(yīng)。結(jié)果表明：1）近10 a來(lái)，榆中縣生境質(zhì)量整體上呈下降趨勢(shì)，生境退化度有小幅度的增加趨勢(shì)，但二者整體上變化幅度較小，同時(shí)生境質(zhì)量及生境退化未發(fā)生變化的圖譜單元占比遠(yuǎn)遠(yuǎn)大于發(fā)生變化的圖譜單元占比，說(shuō)明榆中縣生境整體上存在惡化的風(fēng)險(xiǎn)，但這種風(fēng)險(xiǎn)比較小。

2）研究期內(nèi)生境質(zhì)量總體上呈現(xiàn)出由西南向東北遞減的分布特點(diǎn)，生境退化度呈現(xiàn)出“河北低而河南高”的分布格局。

3）榆中縣生境質(zhì)量及生境退化分布指數(shù)具有明顯的地形位梯度變化效應(yīng)，二者都隨著地形位起伏發(fā)生變化，但生境質(zhì)量和生境退化的優(yōu)勢(shì)地形位梯度差異較大，高質(zhì)量生境在第4、5等地形位梯度上分布最多，占絕對(duì)優(yōu)勢(shì)，生境退化高值在1等地形位梯度上占絕對(duì)優(yōu)勢(shì)，二者的混合梯度都出現(xiàn)在第3級(jí)地形位梯度。