石 垚, 袁大鵬, 趙雪杉, 封 乾, 牛志君, 王樹濤

(1.河北農(nóng)業(yè)大學(xué) 資源與環(huán)境科學(xué)學(xué)院, 河北 保定 071000; 2.河北農(nóng)業(yè)大學(xué) 國(guó)土資源學(xué)院, 河北 保定 071000)

生態(tài)系統(tǒng)服務(wù)功能是指人類賴以生存和維持生態(tài)過(guò)程的自然環(huán)境條件和效用[1-2]，它不僅包括各種生態(tài)系統(tǒng)為人類所提供的食物、醫(yī)藥和農(nóng)業(yè)生產(chǎn)的原料，也是支撐和維持地球生命支持系統(tǒng)的重要組成部分[1]。土地利用變化(LUCC)不僅從地表結(jié)構(gòu)上帶來(lái)巨大的變化，而且影響其生態(tài)過(guò)程，如物質(zhì)循環(huán)和能量流動(dòng)，并影響著整個(gè)生態(tài)系統(tǒng)的結(jié)構(gòu)和功能，從而影響生態(tài)系統(tǒng)服務(wù)價(jià)值及其人類福祉[3-4]。伴隨著全球人口增長(zhǎng)和城市化的發(fā)展，環(huán)境問(wèn)題日益突出，生態(tài)系統(tǒng)服務(wù)受到越來(lái)越多的關(guān)注。探討土地利用變化(LUCC)與生態(tài)過(guò)程之間的相互影響和相互作用，對(duì)土地可持續(xù)利用、生態(tài)環(huán)境保護(hù)和區(qū)域生態(tài)安全等方面的研究具有重要意義[5]。

生態(tài)系統(tǒng)服務(wù)評(píng)估有多種方法，而采用經(jīng)濟(jì)學(xué)的手段，從價(jià)值的角度可以更直觀的描述生態(tài)系統(tǒng)服務(wù)價(jià)值[6-7]。Costanza等[8]于1997年提出了一種來(lái)衡量全球生態(tài)系統(tǒng)服務(wù)價(jià)值的方法，謝高地等[9]結(jié)合中國(guó)實(shí)際情況建立了“中國(guó)生態(tài)系統(tǒng)單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量”表，為我國(guó)生態(tài)系統(tǒng)服務(wù)價(jià)值核算以及演化過(guò)程提供重要依據(jù)[10]。國(guó)內(nèi)大多數(shù)學(xué)者結(jié)合謝高地等[9]的研究成果，采用單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值的平均值直接測(cè)算，或者通過(guò)生物多樣性、生物量等修正方法間接去估算區(qū)域生態(tài)系統(tǒng)服務(wù)價(jià)值[11]，通過(guò)建立的區(qū)域生態(tài)系統(tǒng)服務(wù)價(jià)值表，可以更為精確地表達(dá)研究區(qū)域生態(tài)系統(tǒng)服務(wù)價(jià)值的變化。近些年來(lái)國(guó)內(nèi)學(xué)者主要從LUCC對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響等方面來(lái)進(jìn)行研究，然而，很少通過(guò)地形梯度上的土地利用變化對(duì)區(qū)域生態(tài)系統(tǒng)服務(wù)價(jià)值進(jìn)行研究[12]。通過(guò)從不同地形梯度分析生態(tài)系統(tǒng)服務(wù)價(jià)值，能夠更全面的研究LUCC對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響。

懷來(lái)縣地形分異明顯，南部和北部為山地，中部地區(qū)為平原，有較好的地形梯度特征。本文根據(jù)懷來(lái)縣地形梯度上的土地利用變化，通過(guò)對(duì)不同高程梯度和坡度梯度的分析，對(duì)懷來(lái)縣生態(tài)系統(tǒng)服務(wù)價(jià)值進(jìn)行研究，同時(shí)對(duì)懷來(lái)縣敏感性指數(shù)進(jìn)行測(cè)算為懷來(lái)縣生態(tài)環(huán)境建設(shè)及土地資源的合理利用提供理論支持。

1　研究區(qū)概況及數(shù)據(jù)處理

1.1　研究區(qū)概況

懷來(lái)縣位于河北省西北部，屬于張家口市，地處于燕山北側(cè)，永定河上游區(qū)域從中穿過(guò)，介于東經(jīng)115°16′—115°58′、北緯40°4′—40°35′。懷來(lái)縣全境屬燕山山地，南部和北部為山地，中部為河谷平原，燕山支脈從西北和西南兩側(cè)包圍，因中部低，兩側(cè)高而被稱為“懷來(lái)盆地”，懷來(lái)縣主要有山地、河川平原和丘陵等地貌類型。全縣地處于暖溫帶大陸性季風(fēng)氣候區(qū)，受地形及大氣環(huán)流影響，具有光照充足、雨熱同季、四季分明、日溫差大等特點(diǎn)。年平均降水為396 mm，南部和北部?jī)缮狡?，年平均風(fēng)速為2.3 m/s。懷來(lái)縣總面積為1 801 km2，其中河川平原占全縣總面積的35%左右，丘陵占全縣總面積的24.8%左右，山地占全縣總面積的41.7%左右。

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

本研究所用的土地利用現(xiàn)狀數(shù)據(jù)來(lái)源于2003年、2013年懷來(lái)縣的Landsat TM/ETM+遙感影像；地形數(shù)據(jù)來(lái)自中科院國(guó)際數(shù)據(jù)中心下載的DEM和Slop數(shù)據(jù)，分辨率為90 m。懷來(lái)縣2013年Landsat TM/ETM+遙感影像的二次多項(xiàng)式幾何校正通過(guò)ENVI 4.8軟件進(jìn)行，確保校正后誤差小于0.5個(gè)像元。將糾正影像作為參考，采用Image to Image模塊對(duì)2003年遙感影像進(jìn)行自動(dòng)配準(zhǔn)，并以行政區(qū)作為掩膜切割遙感圖像，遙感圖像的分類是通過(guò)監(jiān)督分類和結(jié)合人機(jī)交互方法進(jìn)行的，分別用代碼1，2，3，…，6表示，并通過(guò)地面測(cè)量數(shù)據(jù)(GPS現(xiàn)場(chǎng)測(cè)量)來(lái)測(cè)試解釋精度，隨機(jī)選擇多個(gè)樣本，使用混淆矩陣法進(jìn)行評(píng)價(jià)分析，結(jié)果表明：2期遙感圖像分類的總體精度分別為83.21%和87.69%，解讀精度超過(guò)80%[13]。根據(jù)國(guó)家一級(jí)地類劃分標(biāo)準(zhǔn)，將該區(qū)域劃分為耕地、林地、建設(shè)用地、園地、水域、未利用地6大類?；贏rcGIS 10.2軟件，運(yùn)用ArcGIS中自然間斷法依據(jù)地類特征將其高程分為6類，坡度梯度分為5類(圖1)，高程分別為0～600 m，600～900 m，900～1 200 m，1 200～1 500 m，1 500～1 800 m，1 800～2 100 m；坡度梯度分別為0°～2°，2°～6°，6°～15°，15°～25°，>25°。

從圖中可以看出，高程分布以600～900 m，900～1 200 m，1 500～1 800 m范圍居多，主要分布于懷來(lái)縣中南部和中北部地區(qū)，1 800～2 100 m范圍內(nèi)分布最少；坡度在2°～6°，15°～25°范圍內(nèi)分布較多，同樣分布于中南部和中北部地區(qū)，>25°范圍內(nèi)分布最少。

2　研究方法

2.1　生態(tài)系統(tǒng)服務(wù)價(jià)值動(dòng)態(tài)度

生態(tài)系統(tǒng)服務(wù)價(jià)值動(dòng)態(tài)度指研究區(qū)某一段時(shí)間某種地類生態(tài)系統(tǒng)服務(wù)價(jià)值的變化[14]，它可以在一定程度上反應(yīng)某種地類的土地利用變化，其動(dòng)態(tài)度值越高，表明該地類土地利用變化越頻繁，動(dòng)態(tài)度值越低，表明該地類土地利用變化越少，公式如下：

式中：Va，Vb分別為研究區(qū)某種地類研究期初和期末的生態(tài)系統(tǒng)服務(wù)價(jià)值；T為研究時(shí)間；EV為某種地類生態(tài)系統(tǒng)服務(wù)價(jià)值動(dòng)態(tài)度。

圖1　研究區(qū)高程梯度和坡度梯度

2.2　生態(tài)系統(tǒng)服務(wù)價(jià)值測(cè)算

本文利用Costanza等[8]提出的生態(tài)服務(wù)價(jià)值計(jì)算公式，基于謝高地等[9]人建立的中國(guó)陸地生態(tài)系統(tǒng)單位面積服務(wù)價(jià)值的基礎(chǔ)上，陳影等[15]人結(jié)合研究區(qū)實(shí)際情況對(duì)其修正制定了懷來(lái)縣生態(tài)系統(tǒng)單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值表(表1)，參照懷來(lái)縣生態(tài)系統(tǒng)單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值表計(jì)算懷來(lái)縣生態(tài)系統(tǒng)服務(wù)價(jià)值，生態(tài)系統(tǒng)服務(wù)價(jià)值計(jì)算公式如下：

式中：ESV表示生態(tài)系統(tǒng)價(jià)值；Ai表示研究區(qū)第i類土地類型的面積；VCi為第i種地類生態(tài)系統(tǒng)服務(wù)價(jià)值的系數(shù)。

表1　懷來(lái)縣生態(tài)系統(tǒng)服務(wù)單位面積價(jià)值表　元/hm2

2.3　敏感性分析方法

借鑒相關(guān)研究成果[16-19]，運(yùn)用經(jīng)濟(jì)學(xué)中常用的彈性系數(shù)概念，計(jì)算價(jià)值系數(shù)的敏感性指數(shù)(CS)[20]，確定生態(tài)系統(tǒng)服務(wù)價(jià)值在時(shí)間的變化過(guò)程中對(duì)于價(jià)值系數(shù)的依賴程度。將各種地類的價(jià)值指數(shù)分別調(diào)整50%，來(lái)衡量生態(tài)系統(tǒng)服務(wù)總價(jià)值的變化。敏感指數(shù)計(jì)算公式如下：

式中：ESV表示生態(tài)系統(tǒng)總價(jià)值；VC表示生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)；m，n分別為調(diào)整前后生態(tài)服務(wù)價(jià)值系數(shù)。如果CS>1，則ESV相對(duì)于VC是有彈性的；如果CS<1，則表明ESV是缺乏彈性的，比值越大，表明生態(tài)系統(tǒng)服務(wù)價(jià)值指數(shù)的準(zhǔn)確性越高。

3　結(jié)果與分析

3.1　土地利用變化分析

由圖2分析可知，懷來(lái)縣土地利用主要以林地和未利用地為主，主要以南部和北部山地居多。2003—2013年期間，懷來(lái)縣園地、林地和建設(shè)用地面積呈增加趨勢(shì)，增長(zhǎng)面積分別為7 077.63 hm2，4 331.36 hm2，157.12 hm2；水域、耕地和未利用地呈減少趨勢(shì)減少面積分別為4 190.84 hm2，3 292.74 hm2，207.3 hm2。具體變化區(qū)域?yàn)椋瑘@地增加的主要區(qū)域?yàn)楣購(gòu)d水庫(kù)西部山地，隨著人類社會(huì)經(jīng)濟(jì)活動(dòng)的擴(kuò)大，部分耕地轉(zhuǎn)為經(jīng)濟(jì)效益較好的園地；林地增加的主要區(qū)域?yàn)閼褋?lái)縣南部和北部的山地，由于退耕還林工程的實(shí)施，將耕地坡度較大的區(qū)域退耕為林地，因此林地面積明顯增加，耕地面積減少；建設(shè)用地在中部河谷平原地區(qū)擴(kuò)張明顯，其原因是該地區(qū)交通便利、地勢(shì)平坦和人口集中，隨著懷來(lái)縣經(jīng)濟(jì)迅速發(fā)展，建設(shè)面積不斷增加；水域面積減少區(qū)域主要為懷來(lái)縣中部的官?gòu)d水庫(kù)周邊區(qū)域，主要原因?yàn)楦珊导又祟惿a(chǎn)生活對(duì)水資源的需求增加等因素；近些年人們對(duì)未利用地的開發(fā)及利用使其面積逐漸減少。

3.2　生態(tài)系統(tǒng)服務(wù)價(jià)值的時(shí)空特征

懷來(lái)縣生態(tài)系統(tǒng)服務(wù)價(jià)值分布及變化趨勢(shì)見圖3，其生態(tài)系統(tǒng)服務(wù)價(jià)值總體程南北高中部低的分布趨勢(shì)，生態(tài)系統(tǒng)服務(wù)價(jià)值高的地區(qū)主要分布于南北的山地地區(qū)，中部平原生態(tài)系統(tǒng)服務(wù)價(jià)值較低。2003—2013年懷來(lái)縣生態(tài)系統(tǒng)服務(wù)總價(jià)值增加0.38億元，其中建設(shè)用地和園地的生態(tài)系統(tǒng)服務(wù)價(jià)值動(dòng)態(tài)度較高，增加生態(tài)系統(tǒng)服務(wù)價(jià)值分別為157.12萬(wàn)元、7 077.63萬(wàn)元，年遞增率分別為0.43%，0.19%，增加區(qū)域主要在中部平原區(qū)；耕地和水域的生態(tài)系統(tǒng)服務(wù)價(jià)值分別減少3 292.74萬(wàn)元、4 190.84萬(wàn)元，年遞減率分別為0.17%，0.16%，減少區(qū)域主要分布在北部山地和中部官?gòu)d水域地區(qū)。

圖3　懷來(lái)縣生態(tài)系統(tǒng)服務(wù)價(jià)值分布及變化趨勢(shì)

3.3　生態(tài)系統(tǒng)服務(wù)價(jià)值的地形梯度特征

懷來(lái)縣生態(tài)系統(tǒng)服務(wù)價(jià)值在高程梯度上呈先減少再增加的分布特征，分布不均。在0～600 m的高程梯度范圍內(nèi)，50%的耕地生態(tài)系統(tǒng)服務(wù)價(jià)值集中于該梯度范圍，占該梯度總價(jià)值的13.8%。在600～900 m高程梯度范圍內(nèi)，82.8%的水域、68.6%的建設(shè)用地和57.3%的園地生態(tài)系統(tǒng)服務(wù)價(jià)值主要集中于該高程梯度。高程梯度在900～1 200 m范圍內(nèi)有34.2%林地和33.6%的未利用地生態(tài)系統(tǒng)服務(wù)價(jià)值，其他用地分布較少。在高程梯度為1 200～1 500 m范圍內(nèi)各土地利用類型生態(tài)系統(tǒng)服務(wù)價(jià)值均較少。在1 500～1 800 m高程梯度內(nèi)，主要有42.2%的林地生態(tài)系統(tǒng)服務(wù)價(jià)值集中在該梯度。在1 800～2 100 m高程梯度范圍內(nèi)各土地利用類型生態(tài)服務(wù)價(jià)值最少，主要分布于南部和北部的山地地區(qū)。

根據(jù)高程梯度上生態(tài)系統(tǒng)服務(wù)變化圖(圖4)分析可知。在0～600 m高程梯度上，園地生態(tài)系統(tǒng)服務(wù)價(jià)值增加最多為4 869.87萬(wàn)元，水域生態(tài)系統(tǒng)服務(wù)價(jià)值減少的最多，減少了3 958.10萬(wàn)元，主要由于生產(chǎn)生活對(duì)水資源需求的增加及其干旱等因素，使其水域面積呈減少趨勢(shì)；600～900 m和900～1 200 m高程梯度內(nèi)，耕地生態(tài)系統(tǒng)服務(wù)價(jià)值顯著減少，減少2 682.41萬(wàn)元，主要是由于人類開發(fā)建設(shè)活動(dòng)占用了大量的耕地以及部分耕地轉(zhuǎn)向了經(jīng)濟(jì)效益較好的園地[13]；在1 200～1 500 m，1 500～1 800 m和1 800～2 100 m的高程范圍內(nèi)，林地生態(tài)系統(tǒng)服務(wù)價(jià)值增加顯著為3 247.98萬(wàn)元，生態(tài)保護(hù)政策的實(shí)施和退耕還林工程是其變化的主要原因。

圖4　基于高程梯度的生態(tài)系統(tǒng)服務(wù)價(jià)值變化

懷來(lái)縣生態(tài)系統(tǒng)服務(wù)價(jià)值在坡度上呈波浪型分布，在0°～2°和2°～6°坡度梯度內(nèi)生態(tài)系統(tǒng)服務(wù)價(jià)值隨坡度的增加而增加，6°～15°坡度范圍內(nèi)緩慢減少，在15°～25°坡度梯度內(nèi)達(dá)到最大，>25°坡度范圍又逐漸減少。0°～2°坡度梯度主要分布于中部平原地區(qū)，該梯度內(nèi)生態(tài)系統(tǒng)服務(wù)總價(jià)值為21 907.96萬(wàn)元。15°～25°坡度梯度面積在該區(qū)分布最廣，其生態(tài)系統(tǒng)服務(wù)價(jià)值約為63 011.98萬(wàn)元，達(dá)最大值，主要分布于南部和北部山地地區(qū)。

懷來(lái)縣坡度梯度上生態(tài)系統(tǒng)服務(wù)價(jià)值的變化見圖5。在0°～2°坡度梯度內(nèi)，園地和林地的生態(tài)系統(tǒng)服務(wù)價(jià)值分別增加1 425.18萬(wàn)元和278.71萬(wàn)元，耕地和水域的生態(tài)系統(tǒng)服務(wù)價(jià)值分別減少1 495.27萬(wàn)元和635.21萬(wàn)元；在2°～6°坡度梯度內(nèi)，園地的生態(tài)系統(tǒng)服務(wù)價(jià)值增加最多為4 285.09萬(wàn)元，由于園地的經(jīng)濟(jì)效益較好，懷來(lái)縣的耕地較多的轉(zhuǎn)變?yōu)閳@地，水域的面積在人類活動(dòng)和自然環(huán)境的影響下日益減少，其生態(tài)系統(tǒng)服務(wù)價(jià)值減少最多為2 924.24萬(wàn)元；6°～15°梯度范圍內(nèi)，林地生態(tài)系統(tǒng)服務(wù)價(jià)值增加943.91萬(wàn)元，水域生態(tài)系統(tǒng)服務(wù)價(jià)值減少600.7萬(wàn)元；在15°～25°坡度梯度內(nèi)，林地生態(tài)系統(tǒng)服務(wù)價(jià)值增加最多，為1 902.23萬(wàn)元，該梯度主要分布與南部和北部山地地區(qū)，由于退耕還林工程的實(shí)施，懷來(lái)縣耕地的生態(tài)系統(tǒng)服務(wù)價(jià)值減少，林地生態(tài)系統(tǒng)服務(wù)價(jià)值增加；在大于25°坡度梯度內(nèi)，林地的生態(tài)系統(tǒng)服務(wù)價(jià)值共增加450.66萬(wàn)元，該梯度范圍主要分布于懷來(lái)縣東南部和東北部山地地區(qū)，該地區(qū)是林地生態(tài)服務(wù)價(jià)值增長(zhǎng)的主要原因。

圖5　基于坡度梯度的生態(tài)系統(tǒng)服務(wù)價(jià)值變化

3.4　基于地形梯度的單項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值變化分析

水體具有最高的廢物處理生態(tài)系統(tǒng)服務(wù)功能，而懷來(lái)縣水域面積較少，因此大部分的廢物處理生態(tài)系統(tǒng)服務(wù)價(jià)值由耕地和林地提供。食物生產(chǎn)生態(tài)系統(tǒng)服務(wù)價(jià)值主要由耕地提供，其他土地利用類型的食物生產(chǎn)生態(tài)系統(tǒng)服務(wù)價(jià)值較低。懷來(lái)縣廢物處理和食物生產(chǎn)生態(tài)系統(tǒng)價(jià)值因耕地面積的大量的減少而導(dǎo)致供給服務(wù)大幅減少，其中食物生產(chǎn)生態(tài)系統(tǒng)服務(wù)價(jià)值為369.41萬(wàn)元，廢棄物處理為1 542.36萬(wàn)元。食物生產(chǎn)生態(tài)系統(tǒng)服務(wù)價(jià)值是當(dāng)?shù)鼐用竦纳媾c發(fā)展的物質(zhì)基礎(chǔ)，因此，食物生產(chǎn)價(jià)值的顯著減少應(yīng)當(dāng)引起有關(guān)部門的重視[12]。森林具有較高的土壤形成與保護(hù)、氣體調(diào)節(jié)和生物多樣性保護(hù)等生態(tài)系統(tǒng)服務(wù)價(jià)值，由于林地的增加，氣體調(diào)節(jié)、土壤形成與保護(hù)和生物多樣性保護(hù)等生態(tài)系統(tǒng)服務(wù)價(jià)值增多，分別為1 501.93萬(wàn)元、1 497.85萬(wàn)元、1 099.59萬(wàn)元。

2003—2013年單項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值在高程梯度上變化見圖6。

在0～600 m高程梯度內(nèi)，水源涵養(yǎng)、廢物處理、食物生產(chǎn)和娛樂(lè)文化均有不同程度的減少，分別減少1 096.48萬(wàn)元、1 254.48萬(wàn)元、30.53萬(wàn)元、119.94萬(wàn)元，其他單項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值均有所增加；高程梯度在600～900 m梯度范圍內(nèi)，氣候調(diào)節(jié)、生物多樣性保護(hù)、原材料和娛樂(lè)文化等生態(tài)系統(tǒng)服務(wù)價(jià)值均有所增加，其他5項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值都有不同程度的減少；在900～1 200 m高程梯度內(nèi)，除食物生產(chǎn)項(xiàng)的生態(tài)系統(tǒng)服務(wù)價(jià)值減少外，其他8項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值均增加；在1 200～1 500 m和1 500～1 800 m高程梯度內(nèi)，各項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值均增加；在1 800～2 100 m高程梯度內(nèi)，各單項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值未有明顯變化。

圖6　基于高程梯度的單項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值變化

圖7為各單項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值在坡度梯度上的變化。在坡度梯度0°～2°的范圍內(nèi)，水源涵養(yǎng)、廢物處理和食物生產(chǎn)等的生態(tài)系統(tǒng)服務(wù)價(jià)值均有所減少，其他各單項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值均增加；在2°～6°坡度梯度內(nèi)，水源涵養(yǎng)與廢物處理生態(tài)系統(tǒng)服務(wù)價(jià)值均減少最多，分別為697.42萬(wàn)元、993.84萬(wàn)元，氣體調(diào)節(jié)、氣候條件、土壤形成與保護(hù)、生物多樣性和原材料的生態(tài)系統(tǒng)服務(wù)價(jià)值均增長(zhǎng)最多；在6°～15°，15°～25°和大于25°范圍內(nèi)各單項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值增減程度不一。

圖7　基于坡度梯度的單項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值變化

3.5　敏感性分析

通過(guò)將研究區(qū)各土地利用類型的價(jià)值系數(shù)分別上下調(diào)整50%，并計(jì)算出研究區(qū)2003年、2013年2期各土地利用類型的敏感性指數(shù)(圖8)。研究結(jié)果表明，懷來(lái)縣所有地類生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)(VC)的敏感性指數(shù)(CS)均小于1，并且都接近于0，從高到低分別為耕地、園地、林地、水域、未利用地、建設(shè)用地，說(shuō)明研究區(qū)內(nèi)生態(tài)系統(tǒng)服務(wù)總價(jià)值(ESV)對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)(VC)缺乏彈性，通過(guò)懷來(lái)縣生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)核算的結(jié)果是可信的。

圖8　研究區(qū)敏感性指數(shù)

4　結(jié) 論

懷來(lái)縣在2003—2013年期間園地、林地和建設(shè)用地面積呈增加趨勢(shì)；水域、耕地和未利用地呈減少趨勢(shì)；懷來(lái)縣生態(tài)系統(tǒng)服務(wù)價(jià)值在2003—2013年共增加了0.38億元，增加區(qū)域主要為中部平原區(qū)，以建設(shè)用地和園地為主；生態(tài)系統(tǒng)服務(wù)價(jià)值在高程上呈南北高、中部低的分布趨勢(shì)，主要集中于600～900 m高程梯度內(nèi)，在坡度上隨坡度的增加呈波浪型變化，在15°～25°坡度梯度內(nèi)達(dá)到最大值；單項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值中，由于水域和耕地面積的減少，水源涵養(yǎng)、廢物處理和食物生產(chǎn)等生態(tài)系統(tǒng)服務(wù)價(jià)值損失嚴(yán)重，分別損失539.53萬(wàn)元、1 542.36萬(wàn)元、369.41萬(wàn)元，應(yīng)引起高度重視；懷來(lái)縣各地類生態(tài)系統(tǒng)服務(wù)價(jià)值(ESV)對(duì)于價(jià)值系數(shù)(VC)的敏感性指數(shù)均小于1，表示缺乏彈性，說(shuō)明本文所采用的懷來(lái)縣生態(tài)服務(wù)價(jià)值系數(shù)是可信的。

本文從地形梯度上對(duì)懷來(lái)縣生態(tài)系統(tǒng)服務(wù)價(jià)值進(jìn)行評(píng)估，可為懷來(lái)縣生態(tài)環(huán)境建設(shè)及土地資源的合理利用提供理論支持。在今后的研究中，應(yīng)當(dāng)結(jié)合修正不同梯度范圍內(nèi)的單位面積生態(tài)系統(tǒng)服務(wù)當(dāng)量，從而更加精確的計(jì)算出不同梯度范圍內(nèi)的生態(tài)系統(tǒng)服務(wù)價(jià)值。

參考文獻(xiàn):

[1]歐陽(yáng)志云,王效科.中國(guó)陸地生態(tài)系統(tǒng)服務(wù)功能及其生態(tài)經(jīng)濟(jì)價(jià)值的初步研究[J].生態(tài)學(xué)報(bào),1999,19(5):607-613.

[2]孫慧蘭,李衛(wèi)紅,陳亞鵬,等.新疆伊犁河流域生態(tài)服務(wù)價(jià)值對(duì)土地利用變化的響應(yīng)[J].生態(tài)學(xué)報(bào),2010,30(4):887-894.

[3]Ii B L T, Skole D L, Sanderson S, et al. Land-Use and Land-Cover Change: science/research plan[J]. Global Change Report, 1995,43:669-679.

[4]石龍宇,崔勝輝,尹鍇,等.廈門市土地利用/覆被變化對(duì)生態(tài)系統(tǒng)服務(wù)的影響[J].地理學(xué)報(bào),2010,65(6):708-714.

[5]劉秀麗,張勃,張調(diào)風(fēng),等.黃土高原土石山區(qū)土地利用變化對(duì)生態(tài)系統(tǒng)服務(wù)的影響:以寧武縣為例[J].生態(tài)學(xué)雜志,2013,32(4):1017-1022.

[6]趙景柱,肖寒.生態(tài)系統(tǒng)服務(wù)的物質(zhì)量與價(jià)值量評(píng)價(jià)方法的比較分析[J].應(yīng)用生態(tài)學(xué)報(bào),2000,11(2):290-292.

[7]江波,陳媛媛,饒恩明,等.博斯騰湖生態(tài)系統(tǒng)最終服務(wù)價(jià)值評(píng)估[J].生態(tài)學(xué)雜志,2015,34(4):1113-1120.

[8]Costanza R, D′Arge R, Groot R D, et al. The value of the world′s ecosystem services and natural capital 1[J]. Nature, 1989,25(1):3-15.

[9]謝高地,魯春霞,冷允法,等.青藏高原生態(tài)資產(chǎn)的價(jià)值評(píng)估[J].重慶第二師范學(xué)院學(xué)報(bào),2003,18(3):189-196.

[10]曾杰,李江風(fēng),姚小薇.武漢城市圈生態(tài)系統(tǒng)服務(wù)價(jià)值時(shí)空變化特征[J].應(yīng)用生態(tài)學(xué)報(bào),2014,25(3):883-891.

[11]Feng L, Cheng S, Su H, et al. A theoretical model for assessing the sustainability Of ecosystem services[J].生態(tài)經(jīng)濟(jì)：英文版,2008,4(3):258-265.

[12]王曉峰,薛亞永,張園.基于地形梯度的陜西省生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估[J].冰川凍土,2016,38(5):1432-1439.

[13]哈凱.基于地形梯度的懷來(lái)縣土地利用景觀格局變化及優(yōu)化研究[D].河北保定:河北農(nóng)業(yè)大學(xué),2015.

[14]趙永華,張玲玲,王曉峰.陜西省生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估及時(shí)空差異[J].應(yīng)用生態(tài)學(xué)報(bào),2011,22(10):2662-2672.

[15]陳影,哈凱,賀文龍,等.冀西北間山盆地區(qū)景觀格局變化及優(yōu)化研究:以河北省懷來(lái)縣為例[J].自然資源學(xué)報(bào),2016,31(4):556-569.

[16]周德成,羅格平,許文強(qiáng),等.1960—2008年阿克蘇河流域生態(tài)系統(tǒng)服務(wù)價(jià)值動(dòng)態(tài)[J].應(yīng)用生態(tài)學(xué)報(bào),2010,21(2):399-408.

[17]張偉,張宏業(yè),王秀紅,等.伊犁新墾區(qū)土地利用變化及其對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響[J].資源科學(xué),2009,31(12):2042-2046.

[18]N.格里高利·曼昆,曼昆,梁小民.經(jīng)濟(jì)學(xué)原理上[M].北京:機(jī)械工業(yè)出版社,2003.

[19]趙艷霞,武愛彬,劉欣,等.淺山丘陵區(qū)土地利用地形梯度特征與生態(tài)服務(wù)價(jià)值響應(yīng)[J].水土保持研究,2014,21(3):141-145.

[20]Kreuter U P, Harris H G, Matlock M D, et al. Change in ecosystem service values in the San Antonio area, Texas☆[J]. Ecological Economics, 2001,39(3):333-346.