華北農(nóng)牧交錯(cuò)區(qū)LUCC對(duì)生態(tài)系統(tǒng)服務(wù)及生態(tài)環(huán)境質(zhì)量的影響
——以內(nèi)蒙古自治區(qū)多倫縣為例

周佳寧, 秦富倉, 鄒 偉, 劉 佳, 楊振奇, 張連根

(1.南京農(nóng)業(yè)大學(xué) 公共管理學(xué)院, 江蘇 南京210095; 2.內(nèi)蒙古農(nóng)業(yè)大學(xué) 沙漠治理學(xué)院, 內(nèi)蒙古 呼和浩特 010019;3.西寧市城市規(guī)劃設(shè)計(jì)研究院, 青海 西寧 810008; 4.內(nèi)蒙古自治區(qū)多倫縣林業(yè)局， 內(nèi)蒙古 多倫 027300)

生態(tài)系統(tǒng)服務(wù)指生態(tài)系統(tǒng)提供人類生存和發(fā)展所需的各種直接或間接惠益的總和[1-2]，生態(tài)環(huán)境質(zhì)量則是描述生態(tài)環(huán)境的優(yōu)劣程度[3],兩者作為評(píng)價(jià)生態(tài)系統(tǒng)服務(wù)功能強(qiáng)弱和優(yōu)劣程度的標(biāo)準(zhǔn)，早已成為國(guó)內(nèi)外學(xué)者關(guān)注的焦點(diǎn)。而土地是自然生態(tài)系統(tǒng)的載體，土地利用/覆被變化(land use/cover change, LUCC)可通過改變生態(tài)系統(tǒng)結(jié)構(gòu)和功能來影響地球生物化學(xué)循環(huán)，進(jìn)而影響生態(tài)系統(tǒng)服務(wù)和生態(tài)環(huán)境質(zhì)量[4]。因此，開展LUCC對(duì)生態(tài)系統(tǒng)服務(wù)及生態(tài)環(huán)境質(zhì)量影響的研究，有利于從國(guó)土資源空間管制和土地利用規(guī)劃的視角正確認(rèn)識(shí)其對(duì)生態(tài)環(huán)境的影響，對(duì)協(xié)調(diào)區(qū)域社會(huì)、經(jīng)濟(jì)、生態(tài)高質(zhì)量發(fā)展具有重要的現(xiàn)實(shí)意義。

黨的十九大報(bào)告作出了“加快生態(tài)文明體制改革，建設(shè)美麗中國(guó)”的重大判斷，縣域空間作為中國(guó)生態(tài)文明建設(shè)第一責(zé)任主體和“一線指揮部”，是國(guó)家和省級(jí)生態(tài)環(huán)境建設(shè)的關(guān)鍵環(huán)節(jié)[5]；目前環(huán)保部亦將縣域作為推進(jìn)全國(guó)生態(tài)文明建設(shè)和優(yōu)化空間開發(fā)格局的基本單元?；谶@一背景，本文擬在以往大多研究主要關(guān)注全球、全國(guó)、省市等[6-11]宏觀或中觀層面環(huán)境變化的基礎(chǔ)上，將研究尺度縮小至某一具體縣域——內(nèi)蒙古自治區(qū)多倫縣。選擇該縣原因其一是中國(guó)幅員遼闊，縣域眾多，而在這些縣域中，地處中國(guó)北方農(nóng)牧交錯(cuò)區(qū)的區(qū)縣往往是人類活動(dòng)強(qiáng)烈影響區(qū)和生態(tài)敏感區(qū)，已成為L(zhǎng)UCC及生態(tài)環(huán)境變化研究熱點(diǎn)區(qū)域[12]；其二是內(nèi)蒙古多倫縣作為華北典型農(nóng)牧交錯(cuò)區(qū)，其地理區(qū)位較為特殊，是京津冀通往塞北的節(jié)點(diǎn)重鎮(zhèn)，也是我國(guó)京津風(fēng)沙源治理的重點(diǎn)區(qū)域，近幾十年來該縣土地利用/覆被變化劇烈，生態(tài)環(huán)境變化也十分明顯。當(dāng)然，國(guó)內(nèi)外針對(duì)農(nóng)牧交錯(cuò)區(qū)LUCC對(duì)生態(tài)環(huán)境影響已有諸多研究[13-16]，但有關(guān)研究仍主要關(guān)注LUCC對(duì)生態(tài)系統(tǒng)服務(wù)總價(jià)值的影響，針對(duì)縣域LUCC如何影響區(qū)域生態(tài)系統(tǒng)服務(wù)價(jià)值及生態(tài)環(huán)境質(zhì)量，并在此基礎(chǔ)上探析生態(tài)環(huán)境變化影響因素的系統(tǒng)性研究仍相對(duì)匱乏。故本文在分析華北典型農(nóng)牧交錯(cuò)區(qū)內(nèi)蒙古多倫縣LUCC的基礎(chǔ)上，結(jié)合Costanza等[17]提出的生態(tài)服務(wù)價(jià)值核算方法與謝高地[18]建立的中國(guó)生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量表，修正并估算了多倫縣生態(tài)服務(wù)價(jià)值；同時(shí)選取生態(tài)系統(tǒng)單項(xiàng)功能服務(wù)價(jià)值和敏感性系數(shù)以及生態(tài)環(huán)境質(zhì)量指數(shù)，建立起LUCC與生態(tài)環(huán)境質(zhì)量之間的關(guān)聯(lián)；最后，對(duì)引起區(qū)域生態(tài)環(huán)境質(zhì)量變化的影響因素進(jìn)行分析。研究既凸顯現(xiàn)階段縣域生態(tài)環(huán)境變化對(duì)區(qū)域乃至全國(guó)生態(tài)文明建設(shè)的重要性，也為類似地區(qū)解決土地利用模式粗放、生態(tài)系統(tǒng)功能退化、資源環(huán)境破壞等問題提供科學(xué)路徑，為因地制宜制定縣域國(guó)土空間用途管制規(guī)劃體系提供理論參考。

1 材料與方法

1.1 研究區(qū)概況

多倫縣(115°51′—116°54′N,41°46′—42°36′E)位于內(nèi)蒙古自治區(qū)中部、錫林郭勒盟東南，地處陰山北麓，小興安嶺余脈及燕山山脈末端；全縣轄4個(gè)鄉(xiāng)鎮(zhèn)(2鄉(xiāng)2鎮(zhèn))，縣域面積3 974.05 km2；地貌類型為平緩起伏的沙地，海拔范圍1 150～1 800 m，平均高度1 350 m；年平均氣溫2.9 ℃，平均降水量385.5 mm，平均大風(fēng)日數(shù)84 d；境內(nèi)有常年性河流47條，年平均流量為1.34×108m3；土壤類型有栗鈣土、風(fēng)沙土、草甸土，成土母質(zhì)多以風(fēng)積、沖積沙為主；植被類型為中溫型草原植被及沙地植被等，植被覆蓋時(shí)間較短。該縣為華北典型的干旱和半干旱農(nóng)牧交錯(cuò)區(qū)，屬生態(tài)脆弱地帶；此外它也是面向京津，通往塞北的節(jié)點(diǎn)重鎮(zhèn)，與首都北京的直線距離僅173 km。

1.2 數(shù)據(jù)來源和方法

本文土地利用數(shù)據(jù)主要來源于多倫縣3期(2004年7月、2009年6月、2014年7月)Landsat-TM(30 m×30 m)遙感影像，并結(jié)合相應(yīng)年份土地利用現(xiàn)狀圖及其變更調(diào)查數(shù)據(jù)修正得到。研究中為保證基礎(chǔ)數(shù)據(jù)的精確性和準(zhǔn)確性，具體提取過程為：首先利用ERDAS 9.0對(duì)3期TM影像進(jìn)行波段合成、幾何精校正及行政區(qū)劃裁剪；再結(jié)合《土地利用現(xiàn)狀分類》標(biāo)準(zhǔn)和國(guó)家《生態(tài)環(huán)境遙感調(diào)查分類規(guī)范》，將多倫縣土地利用/覆被類型劃分為耕地、林地、草地、建設(shè)用地、水域和未利用地6類；最后通過縣域土地利用現(xiàn)狀圖、土地利用年度變更調(diào)查數(shù)據(jù)及野外實(shí)地抽樣對(duì)解譯結(jié)果進(jìn)行糾偏，糾偏結(jié)果顯示各期土地利用數(shù)據(jù)精度均在90%以上，滿足允判精度(≥80%)及研究要求[19](附圖4)。

在此基礎(chǔ)上選取面積變化幅度和單一土地利用動(dòng)態(tài)度指數(shù)，構(gòu)建2004—2014年各土地利用/覆被類型轉(zhuǎn)移率矩陣以分析該區(qū)LUCC的變化數(shù)量及轉(zhuǎn)移情況。其中土地利用/覆被面積變化幅度和單一土地利用動(dòng)態(tài)度可定量描述區(qū)域土地利用/覆被變化的數(shù)量和速率，其公式分別如下：

ΔA=At2-At1

(1)

(2)

式中：ΔA——研究時(shí)段某一土地利用/覆被類型面積變化幅度；k——單一土地利用動(dòng)態(tài)度；At1，At2——研究初期和末期各土地利用/覆被類型面積；t1，t2——變化初期和末期。

土地利用/覆被變化的轉(zhuǎn)移率矩陣可反映區(qū)域內(nèi)不同土地利用/覆被類型間相互轉(zhuǎn)換的過程與規(guī)律。其通用形式為：

(3)

式中：n——轉(zhuǎn)移前后的土地利用/覆被類型數(shù)量；i,j(i，j=1,2，…n)——轉(zhuǎn)移前和轉(zhuǎn)移后的土地利用類型；Pij——i地類轉(zhuǎn)換為j地類的概率。

謝高地等[18]根據(jù)Costanza[17]的研究成果修正并提出了適合中國(guó)國(guó)情的大尺度生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)表，并又進(jìn)一步提出了較為成熟的“基于糧食價(jià)格”和“基于生物量”兩種核算中國(guó)不同地區(qū)小尺度生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量的修訂方法，得到了良好應(yīng)用[4，10]。本研究考慮數(shù)據(jù)可得性，選擇“糧食價(jià)格修訂法”進(jìn)行修正，該方法具體計(jì)算過程為：根據(jù)謝高地等[18]研究指出1個(gè)標(biāo)準(zhǔn)生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量因子相當(dāng)于1 hm2平均產(chǎn)量的農(nóng)田每年自然糧食產(chǎn)量的經(jīng)濟(jì)價(jià)值，其經(jīng)濟(jì)價(jià)值等于當(dāng)年研究區(qū)平均糧食單產(chǎn)市場(chǎng)價(jià)值的1/7，由此可結(jié)合公式(4)—(6)修正得出多倫縣農(nóng)田生態(tài)系統(tǒng)食物生產(chǎn)服務(wù)價(jià)值為943.34元/hm2，最后以此推算出研究區(qū)所有其他土地利用/覆被類型生態(tài)服務(wù)價(jià)值系數(shù)(表2)。

(4)

ESV=∑Ak·VCk

(5)

ESVi=∑Ak·VCik

(6)

式中：En——多倫縣農(nóng)田生態(tài)系統(tǒng)食物生產(chǎn)服務(wù)價(jià)值(元/hm2)；Pi——農(nóng)作物i的價(jià)格(元/kg)；qi——農(nóng)作物i的總產(chǎn)量(kg)；n——區(qū)域農(nóng)作物種類數(shù)量；M——n種農(nóng)作物總面積(hm2)； ESV——生態(tài)系統(tǒng)服務(wù)價(jià)值(元/hm2)；Ak——第k種土地利用/覆被類型面積； VCk——第k種土地利用/覆被類型生態(tài)服務(wù)價(jià)值系數(shù)(元·hm-2/a)； ESVi——生態(tài)系統(tǒng)第i項(xiàng)服務(wù)價(jià)值(元/hm2)； VCik——第k種土地利用/覆被類型第i項(xiàng)服務(wù)價(jià)值系數(shù)(元·hm-2/a)。下同。

在此基礎(chǔ)上，本文采用經(jīng)濟(jì)學(xué)中彈性系數(shù)(CS)作為敏感性評(píng)價(jià)指標(biāo)，來刻畫生態(tài)服務(wù)價(jià)值隨時(shí)間變化對(duì)生態(tài)服務(wù)價(jià)值系數(shù)的依賴程度。若CS>1，表示生態(tài)服務(wù)價(jià)值相對(duì)生態(tài)服務(wù)價(jià)值系數(shù)富有彈性，系數(shù)表需根據(jù)時(shí)序變化加以修正；反之，系數(shù)表無需修正。敏感性系數(shù)的計(jì)算公式為：

(7)

式中：VCi，VCj——變化前后不同土地利用/覆被類型生態(tài)服務(wù)價(jià)值系數(shù)。

相對(duì)生態(tài)服務(wù)價(jià)值則是利用模糊數(shù)學(xué)法將各土地利用/覆被類型生態(tài)服務(wù)價(jià)值在[0，1]區(qū)間內(nèi)賦值。由于單位面積水域生態(tài)服務(wù)價(jià)值最高，賦值為1，其它類型則以此為基準(zhǔn)，按比例確定；最終得到各土地利用/覆被類型相對(duì)生態(tài)服務(wù)價(jià)值(表3)。

生態(tài)環(huán)境質(zhì)量指數(shù)Ei可定量描述區(qū)域生態(tài)環(huán)境總體狀況[20],計(jì)算公式為：

(8)

式中：Ati——T時(shí)段i類土地利用/覆被類型面積；Ci——i類土地利用/覆被類型相對(duì)生態(tài)服務(wù)價(jià)值； TA——區(qū)域總面積。

土地利用/覆被類型轉(zhuǎn)換會(huì)引起區(qū)域生態(tài)環(huán)境質(zhì)量變化；一般來說土地利用/覆被類型相互轉(zhuǎn)換會(huì)使區(qū)域環(huán)境朝著改善和惡化兩種方向發(fā)展；有基于此，可建立LUCC與生態(tài)環(huán)境變化關(guān)聯(lián)體系，通過土地利用/覆被變化生態(tài)貢獻(xiàn)率的計(jì)算，定量研究區(qū)域生態(tài)環(huán)境質(zhì)量變化的影響因素。

LEi=(Ei+1-Ei)×Pi

(9)

式中：LEi——土地利用/覆被變化生態(tài)貢獻(xiàn)率；Ei，Ei+1——某一土地利用/覆被類型變化初期和末期的生態(tài)環(huán)境質(zhì)量指數(shù);Pi——土地利用/覆被類型轉(zhuǎn)移概率。

2 結(jié)果與分析

2.1 土地利用/覆被變化分析

由圖1可以看出，2004—2014年多倫縣林地和建設(shè)用地面積有所增加；而草地、耕地、水域和未利用地面積有所減少。其中，林地和建設(shè)用地在2004—2009年和2009—2014年均為增加，而耕地和未利用地均為減少；草地表現(xiàn)為先增后減，水域與之相反。10 a變化中，草地、耕地和林地面積之和始終占區(qū)域面積比重的86%以上，且以草地占據(jù)優(yōu)勢(shì)地位。

圖1 多倫縣2004-2014年各土地利用/覆被類型面積

從面積變化幅度和單一土地利用動(dòng)態(tài)度的分析可知(圖2)： ①10 a間林地面積增加幅度最大，且2004—2009年和2009—2014年單一土地利用動(dòng)態(tài)度均超過10%，說明其變化速率較高； ②建設(shè)用地在2004—2009年大幅度增加，而2009—2014年增加幅度明顯降低；其前一時(shí)段單一土地利用動(dòng)態(tài)度為后一時(shí)段27倍，且明顯高于同時(shí)期其他土地利用/覆被類型。 ③耕地和未利用地在2004—2009年和2009—2014年均呈減少趨勢(shì)，且都為前一時(shí)段減少幅度最大；但未利用地在兩時(shí)段單一土地利用動(dòng)態(tài)度均超過-15%，分別為同時(shí)期耕地面積變化率的5倍和17倍，可見未利用地減少之快。 ④10 a間草地和水域面積變化最??；研究時(shí)段各土地利用/覆被類型變化率從高到低依次為林地、未利用地、建設(shè)用地、耕地、草地和水域。

圖2 多倫縣2004-2014年土地利用/ 覆被類型面積變化幅度和變化速率

根據(jù)多倫縣近10 a來LUCC轉(zhuǎn)移率矩陣分析可知(表1)： ①2004—2014年，林地和建設(shè)用地以轉(zhuǎn)入為主；草地、耕地、水域和未利用地以轉(zhuǎn)出為主。 ②從各土地利用/覆被類型相互轉(zhuǎn)移概率看：草地和林地相互轉(zhuǎn)換較為頻繁，二者之間相互轉(zhuǎn)換概率均較高；建設(shè)用地以轉(zhuǎn)向草地和耕地為主，轉(zhuǎn)移率分別為8.00%和7.61%，相差不大；水域轉(zhuǎn)向草地概率達(dá)24.81%，其次為轉(zhuǎn)向林地；未利用地以轉(zhuǎn)向草地和林地為主，且轉(zhuǎn)向草地概率最大，為51.89%；耕地轉(zhuǎn)向林地和草地概率均超過10%，且由于研究區(qū)退耕還林補(bǔ)助高于還草，使得耕地轉(zhuǎn)向前者概率略高于后者，分別是12.77%和11.89%。

表1 多倫縣2004-2014年不同土地利用/覆被類型轉(zhuǎn)移率矩陣 %

2.2 LUCC與生態(tài)環(huán)境變化分析

2.2.1 生態(tài)系統(tǒng)服務(wù)價(jià)值變化分析 根據(jù)上述各項(xiàng)農(nóng)田生態(tài)系統(tǒng)服務(wù)價(jià)值指標(biāo)計(jì)算公式(4)—(6)修正計(jì)算得到多倫縣不同地類生態(tài)服務(wù)價(jià)值系數(shù)(表2)及區(qū)域生態(tài)系統(tǒng)服務(wù)總價(jià)值和各地類生態(tài)服務(wù)價(jià)值及其變化(表3)。

表2 多倫縣不同土地利用/覆被類型生態(tài)服務(wù)價(jià)值系數(shù) 元/(hm2·a)

分析可知，草地、耕地和林地是研究區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值組成中貢獻(xiàn)率最大的3種地類；各期比重之和均在95%以上。2004—2014年，多倫縣生態(tài)系統(tǒng)服務(wù)價(jià)值總量不斷提高；由2004年的30.90億元增加至2009年的38.28億元，再增加至2014年的46.08億元；兩時(shí)段增加幅度分別為7.37和7.81億元，增加速率分別是23.86%和20.39%，年均增加速率均超過4%。從各時(shí)段變化看，2004—2009年，草地和林地生態(tài)系統(tǒng)服務(wù)價(jià)值呈增加趨勢(shì)，且林地增加幅度最大，為7.36億元；耕地、水域和未利用地生態(tài)系統(tǒng)服務(wù)價(jià)值有不同程度減小，減小幅度最大的是耕地，為0.77億元，這與該時(shí)段耕地被大量占用有關(guān)。由于水域面積變化不明顯，其生態(tài)系統(tǒng)服務(wù)價(jià)值減少幅度僅為0.18億元；而未利用地面積雖然在該時(shí)段減少最多，但由于其生態(tài)服務(wù)價(jià)值系數(shù)相對(duì)較小，故其生態(tài)系統(tǒng)服務(wù)價(jià)值減少幅度也很小，但減少率卻高達(dá)94.45%，為該時(shí)段減少最快。2009—2014年，由于草地面積大量減少，其生態(tài)服務(wù)價(jià)值也呈減少趨勢(shì)，減少幅度為該時(shí)段最大(3.50億元)；而林地和水域生態(tài)服務(wù)價(jià)值有所增加，且林地增加最大，為11.34億元，增加率達(dá)66.56%；未利用地生態(tài)服務(wù)價(jià)值基數(shù)小，小幅度變化即會(huì)導(dǎo)致變化率很大，因此其生態(tài)系統(tǒng)服務(wù)價(jià)值減少率為該時(shí)段最大，為74.23%，年均減少速度較快。

表3 多倫縣2004-2014年生態(tài)系統(tǒng)服務(wù)價(jià)值變化

2.2.2 敏感性分析 根據(jù)公式(7)，分別將研究區(qū)不同土地利用/覆被類型的生態(tài)服務(wù)價(jià)值系數(shù)上下調(diào)整50%，可以計(jì)算出該區(qū)生態(tài)服務(wù)價(jià)值敏感性系數(shù)，結(jié)果詳見表4。

分析可知，各土地利用/覆被類型生態(tài)服務(wù)價(jià)值敏感性系數(shù)均小于1，說明研究區(qū)生態(tài)服務(wù)價(jià)值與生態(tài)服務(wù)價(jià)值系數(shù)之間缺乏彈性；因此根據(jù)該系數(shù)計(jì)算出的多倫縣生態(tài)系統(tǒng)服務(wù)價(jià)值結(jié)果可信、模型穩(wěn)健。不同土地利用/覆被類型生態(tài)服務(wù)價(jià)值敏感性系數(shù)變化規(guī)律與各自面積變化趨勢(shì)相一致；其中由于草地面積基數(shù)大，生態(tài)服務(wù)價(jià)值高，敏感性指數(shù)值也較高；其次為耕地和林地，二者生態(tài)服務(wù)價(jià)值敏感性系數(shù)均在0.04以上，但都低于0.2；雖然水域生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)最高，但由于其面積變化幅度小，導(dǎo)致其敏感性指數(shù)值變化并不明顯；而未利用地面積減少幅度大使得其敏感性指數(shù)值也顯著減小。

表4 多倫縣生態(tài)服務(wù)價(jià)值敏感性指數(shù)

2.2.3 生態(tài)環(huán)境質(zhì)量指數(shù)變化分析 分析表5可知，整體而言，2004—2014年多倫縣生態(tài)環(huán)境質(zhì)量指數(shù)不斷提高。在2004—2009年，生態(tài)環(huán)境指數(shù)由0.179 32增加到0.222 11；增加了23.86%；2009—2014年增加幅度略高于前一階段，為0.045 29，生態(tài)環(huán)境指數(shù)上升到0.267 40；但從變化細(xì)節(jié)上來分析，研究區(qū)生態(tài)環(huán)境存在改善和惡化雙向驅(qū)動(dòng)機(jī)制。

表5 多倫縣2004-2014年生態(tài)環(huán)境質(zhì)量指數(shù)

2.2.4 LUCC與生態(tài)環(huán)境效應(yīng)的影響因素分析 結(jié)合公式(9)和表6，分析多倫縣生態(tài)環(huán)境改善和惡化的兩種土地利用/覆被變化的影響因素體系，同時(shí)計(jì)算出相應(yīng)貢獻(xiàn)率；結(jié)果表明： ①2004—2014年，不同土地利用/覆被類型相互轉(zhuǎn)換對(duì)研究區(qū)生態(tài)環(huán)境改善的貢獻(xiàn)率高于對(duì)生態(tài)環(huán)境惡化的貢獻(xiàn)率；說明研究時(shí)段各土地利用/覆被變化使多倫縣生態(tài)環(huán)境質(zhì)量有所提高。 ②引起生態(tài)環(huán)境質(zhì)量改善的影響因素中，未利用地向草地和林地轉(zhuǎn)換是引起區(qū)域生態(tài)環(huán)境質(zhì)量提高的主導(dǎo)因素，二者貢獻(xiàn)率之和占總貢獻(xiàn)率的72.89%；其次為耕地向林地轉(zhuǎn)換的生態(tài)貢獻(xiàn)率，為1.774%；未利用地轉(zhuǎn)向耕地生態(tài)貢獻(xiàn)率最小，為0.109%；由此可以看出，草地和林地變化顯著影響區(qū)域生態(tài)環(huán)境質(zhì)量，退耕還林還草和提高未利用地向林地和草地遷移率是改善區(qū)域生態(tài)環(huán)境質(zhì)量的有效手段。 ③草地轉(zhuǎn)向耕地和建設(shè)用地是造成區(qū)域生態(tài)環(huán)境質(zhì)量惡化的主要影響因素，二者對(duì)區(qū)域生態(tài)環(huán)境惡化的貢獻(xiàn)率之和達(dá)89.45%，而其他土地利用/覆被類型轉(zhuǎn)換造成的生態(tài)環(huán)境惡化貢獻(xiàn)率均低于0.01%；可見今后加強(qiáng)生態(tài)脆弱區(qū)草原保護(hù)任重而道遠(yuǎn)，且意義深遠(yuǎn)。

表6 多倫縣2004-2014年生態(tài)環(huán)境 變化影響因素類型及其貢獻(xiàn)率

3 結(jié) 論

(1) 2004—2014年，多倫縣林地和建設(shè)用地不斷增加，草地、耕地、水域和未利用地不斷減少；其中林地增加幅度最大，而未利用地減少最多，二者的變化率均較高。

(2) 土地利用轉(zhuǎn)移率矩陣顯示，區(qū)域耕地轉(zhuǎn)向林地的概率高于轉(zhuǎn)向草地的概率，這是由于退耕還林補(bǔ)助高于還草所致；而建設(shè)用地以轉(zhuǎn)向草地和耕地為主；未利用地轉(zhuǎn)向草地和林地的概率相對(duì)較大。

(3) 多倫縣生態(tài)系統(tǒng)服務(wù)價(jià)值總量不斷提高。草地、耕地和林地是研究區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值組成中貢獻(xiàn)率最大的3種地類；各土地利用/覆被類型生態(tài)服務(wù)價(jià)值變化與各自面積變化規(guī)律相一致。

(4) 就生態(tài)環(huán)境質(zhì)量變化看，一方面，農(nóng)牧交錯(cuò)區(qū)生態(tài)環(huán)境質(zhì)量不斷提高；另一方面，生態(tài)環(huán)境質(zhì)量變化亦存在改善和惡化兩種趨勢(shì)。其中未利用地向草地和林地轉(zhuǎn)換是生態(tài)環(huán)境質(zhì)量改善的主導(dǎo)因素；而草地轉(zhuǎn)向耕地和建設(shè)用地則對(duì)生態(tài)環(huán)境質(zhì)量惡化的貢獻(xiàn)率較高。

(5) 本文研究表明，在全國(guó)深入推進(jìn)生態(tài)文明建設(shè)的體制改革中，要?jiǎng)?wù)必重視生態(tài)脆弱區(qū)LUCC對(duì)生態(tài)環(huán)境的影響；針對(duì)華北農(nóng)牧交錯(cuò)區(qū)而言，應(yīng)切實(shí)提高其土地利用效率，加強(qiáng)地區(qū)草原和森林保護(hù)；統(tǒng)籌規(guī)劃，建立合理的區(qū)域國(guó)土空間管控體系，促進(jìn)區(qū)域人口、經(jīng)濟(jì)、資源和環(huán)境均衡持續(xù)發(fā)展。當(dāng)然，由于生態(tài)系統(tǒng)是一個(gè)具有物質(zhì)循環(huán)、能量流動(dòng)和信息傳遞等復(fù)雜動(dòng)態(tài)過程的系統(tǒng)，對(duì)該系統(tǒng)生態(tài)環(huán)境效應(yīng)的分析，在未來研究中應(yīng)進(jìn)一步結(jié)合具體生態(tài)過程加以深入探討。