封建民, 文 琦, 郭玲霞

(1.咸陽(yáng)師范學(xué)院 資源環(huán)境與歷史文化學(xué)院, 陜西 咸陽(yáng) 712000; 2.寧夏大學(xué) 資源環(huán)境學(xué)院, 銀川 750021)

土地是人類賴以生產(chǎn)和生活的物質(zhì)基礎(chǔ)，是一個(gè)自然經(jīng)濟(jì)的綜合體。土地利用狀態(tài)及變化趨勢(shì)是人類基于社會(huì)經(jīng)濟(jì)目的，對(duì)土地進(jìn)行長(zhǎng)期經(jīng)營(yíng)的各類活動(dòng)與自然環(huán)境相互作用的集中表現(xiàn)[1-2]，是區(qū)域生態(tài)環(huán)境質(zhì)量的綜合表現(xiàn)形式，并且能在一定程度上反映生態(tài)環(huán)境優(yōu)劣程度，同時(shí)也是全球變化的重要組成部分和驅(qū)動(dòng)因子之一[3-5]。土地利用類型的變化將直接影響區(qū)域生態(tài)系統(tǒng)向社會(huì)提供物質(zhì)和能量，進(jìn)而影響生態(tài)系統(tǒng)直接向人類社會(huì)提供各項(xiàng)服務(wù)的能力，因而近年來區(qū)域土地利用變化及其生態(tài)系統(tǒng)服務(wù)價(jià)值的響應(yīng)研究已成為生態(tài)學(xué)和生態(tài)經(jīng)濟(jì)學(xué)的一個(gè)熱點(diǎn)[6-10]。研究土地利用背景下的生態(tài)系統(tǒng)服務(wù)價(jià)值變化對(duì)促進(jìn)區(qū)域生態(tài)建設(shè)、研究區(qū)域可持續(xù)發(fā)展具有重要意義。

榆林市地處毛烏素沙漠和黃土高原過渡地帶，風(fēng)蝕沙化和水土流失嚴(yán)重，東南部丘陵溝壑區(qū)，是黃河中游水土流失最嚴(yán)重的地區(qū)。進(jìn)入“九五”以來，國(guó)家把榆林列為防沙治沙重點(diǎn)地區(qū)。近年來隨著土地沙漠化治理力度的加大、退耕還林政策的實(shí)施以及社會(huì)經(jīng)濟(jì)的發(fā)展，榆林市土地利用結(jié)構(gòu)發(fā)生了深刻的變化，這必然會(huì)引起生態(tài)服務(wù)價(jià)值的變化。本研究利用RS和GIS技術(shù)，分析近25 a來榆林市土地利用變化特征，并根據(jù)當(dāng)?shù)貙?shí)際情況修正各生態(tài)系統(tǒng)的生態(tài)服務(wù)價(jià)值當(dāng)量，計(jì)算區(qū)域生態(tài)服務(wù)價(jià)值及變化特征，旨在為區(qū)域土地利用規(guī)劃、生態(tài)環(huán)境保護(hù)等決策提供科學(xué)依據(jù)。

1 研究區(qū)概況

榆林市位于陜西省北部，地處鄂爾多斯盆地東南部的黃土高原和毛烏素沙地過渡帶，區(qū)內(nèi)平均海拔1 000～1 500 m，呈西北高、東南低的特點(diǎn)，范圍位于東經(jīng)107°28′—111°15′，北緯36°57′—39°34′，東臨黃河與山西相望，西連寧夏、甘肅，北鄰內(nèi)蒙，南接陜西省延安市。地域東西長(zhǎng)385 km，南北寬約263 km，總面積43 578.68 km2。地貌大體以長(zhǎng)城為界，北部為風(fēng)沙草灘區(qū)，占全市總面積的42%，該區(qū)地勢(shì)較為平坦，海撥較高，地表多為流動(dòng)和半固定沙地，土地利用以沙地、灌木林地、草地和耕地為主；南部為黃土丘陵溝壑區(qū)，占全市總面積的52%，該區(qū)相對(duì)高差較大，山大溝深，石多土薄，水土流失嚴(yán)重，土地利用以草地、耕地、林地為主。榆林氣候?qū)倥瘻貛Ш蜏貛О敫珊荡箨懶约撅L(fēng)氣候，寒暑劇烈，氣候干燥，災(zāi)害頻發(fā)，四級(jí)分明，冬季漫長(zhǎng)寒冷，夏季短促，溫差大，多西北風(fēng)，風(fēng)沙頻繁，無霜期短，年平均氣溫10℃，平均降水400 mm左右，無霜期150 d左右。全市轄榆陽(yáng)、神木、府谷、定邊、靖邊、橫山、佳縣、米脂、吳堡、綏德、清澗、子洲共1區(qū)11縣，全市總?cè)丝?37萬人。

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

2.1 數(shù)據(jù)來源及土地利用信息提取

以研究區(qū)1990年、2000年、2010年、2015年4期Landsat TM/ETM遙感影像為數(shù)據(jù)源，運(yùn)用ERDAS 9.3軟件對(duì)4期遙感影像進(jìn)行波段組合、幾何校正、投影轉(zhuǎn)換、圖像鑲嵌、裁剪等預(yù)處理的基礎(chǔ)上，按照2007年8月頒布的《土地利用現(xiàn)狀分類》國(guó)家標(biāo)準(zhǔn)，結(jié)合研究區(qū)土地資源利用特征、遙感影像特點(diǎn)及研究需要，將研究區(qū)土地利用類型歸納劃分為耕地(水澆地和旱地)、林地(有林地和灌木林地)、草地、水域、建設(shè)用地(居住地、工業(yè)用地和交通用地)、未利用土地(沙地和裸地)6個(gè)類型。利用監(jiān)督分類和人工目視判讀相結(jié)合的方法進(jìn)行遙感影像的解譯，生成4期土地利用數(shù)據(jù)，利用榆林市1∶5萬地形圖以及Google Earth高分辨率影像進(jìn)行精度評(píng)價(jià)，評(píng)價(jià)結(jié)果均超過85%。最后應(yīng)用GIS技術(shù)對(duì)4期數(shù)據(jù)進(jìn)行空間疊置，最終得出3個(gè)時(shí)段各個(gè)地類變化數(shù)據(jù)。

2.2 研究方法

2.2.1 土地利用變化動(dòng)態(tài)分析 本研究中，通過單一土地利用動(dòng)態(tài)度和土地利用變化轉(zhuǎn)移矩陣來分析研究區(qū)土地利用結(jié)構(gòu)變化。其中，單一土地利用動(dòng)態(tài)度表征研究區(qū)某段時(shí)間內(nèi)某種土地利用的數(shù)量變化速度[10]，其表達(dá)式為：

(1)

式中：K為研究區(qū)某一土地利用類型動(dòng)態(tài)度；Ua，Ub分別為研究初期和末期某一種土地類型的面積；T為研究時(shí)段，當(dāng)T為年份時(shí)，K值表示研究區(qū)某種土地利用類型的年變化率。

K代表某一種土地利用類型面積的年際變化，反映了該類土地變化的劇烈程度，然而當(dāng)該類土地與其他類型土地轉(zhuǎn)入、轉(zhuǎn)出雙向變化較為頻繁時(shí)，則不能很好地反映其轉(zhuǎn)化的劇烈程度。為了反映某種土地利用類型轉(zhuǎn)化的劇烈程度，采用公式(2)表示某一土地利用類型變化的動(dòng)態(tài)度[11]：

(2)

式中：ΔUout為研究時(shí)段T內(nèi)某一類型土地轉(zhuǎn)變?yōu)槠渌愋偷拿娣e之和；ΔUin為其他類型轉(zhuǎn)變?yōu)樵擃愋偷拿娣e之和；Ks為某一土地利用類型變化的動(dòng)態(tài)度，其值趨于0時(shí)，說明該類型土地在研究時(shí)段內(nèi)趨于穩(wěn)定，其值越大，表明該類型土地在研究時(shí)段內(nèi)變化愈劇烈。

2.2.2 生態(tài)服務(wù)價(jià)值評(píng)估 1997年Costanza等[12]提出了生態(tài)系統(tǒng)服務(wù)價(jià)值估算的原理和方法，為區(qū)域土地利用變化與生態(tài)服務(wù)價(jià)值估算搭建了定量評(píng)價(jià)的橋梁。計(jì)算公式如下：

(3)

式中：ESV為區(qū)域生態(tài)系統(tǒng)服務(wù)總價(jià)值；Ai為第i類生態(tài)系統(tǒng)的面積；Eij為第i類生態(tài)系統(tǒng)第j項(xiàng)服務(wù)功能的價(jià)值；m為生態(tài)系統(tǒng)類型數(shù)；n為生態(tài)系統(tǒng)服務(wù)類別數(shù)。

2003年謝高地等[13]，根據(jù)中國(guó)實(shí)情提出了中國(guó)不同陸地生態(tài)系統(tǒng)單位面積服務(wù)價(jià)值當(dāng)量表。以榆林市多年平均糧食產(chǎn)量1 916.04 kg/hm2為基準(zhǔn)單產(chǎn)，按陜西省多年糧食平均單價(jià)1.47元/kg，得出榆林市農(nóng)田自然糧食產(chǎn)量的經(jīng)濟(jì)價(jià)值為402.37元/(hm2·a)。根據(jù)以上數(shù)據(jù)測(cè)算出榆林市不同生態(tài)系統(tǒng)單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值表[14]。

2.2.3 敏感性分析 為了揭示生態(tài)系統(tǒng)服務(wù)價(jià)值隨時(shí)間的變化對(duì)于價(jià)值指數(shù)的依賴程度，選取彈性系數(shù)來計(jì)算研究區(qū)價(jià)值系數(shù)對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的敏感性指數(shù)(coefficient of sensitive,CS)[15-16]，敏感性指數(shù)是指生態(tài)服務(wù)價(jià)值系數(shù)變動(dòng)1%引起生態(tài)系統(tǒng)服務(wù)價(jià)值的變化情況。本文將研究區(qū)各類土地利用類型的價(jià)值指數(shù)分別上下調(diào)整50%，進(jìn)行檢驗(yàn)。所用公式如下：

(4)

式中：VCik和VCjk代表調(diào)整前后的第k類生態(tài)系統(tǒng)單位面積生態(tài)服務(wù)價(jià)值系數(shù)；ESVi和ESVj分別代表調(diào)整前后的生態(tài)服務(wù)總價(jià)值。CS為研究區(qū)各生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)的敏感度，如果CS<1，表明ESV相對(duì)于VC是缺乏彈性的，即生態(tài)服務(wù)價(jià)值對(duì)價(jià)值系數(shù)不敏感；反之，如果CS>l，表明ESV相對(duì)于VC是富有彈性的，即生態(tài)服務(wù)價(jià)值對(duì)價(jià)值系數(shù)敏感。CS越大，說明VC對(duì)ESV的影響越大，其取值的準(zhǔn)確性也越重要，表明生態(tài)服務(wù)價(jià)值系數(shù)的準(zhǔn)確性越關(guān)鍵。

3 結(jié)果與分析

3.1 土地利用變化分析

3.1.1 土地利用數(shù)量變化分析 從榆林市土地利用總體情況來看(圖1)，草地所占面積最大，占到全市面積的60%以上，是研究區(qū)最主要的土地利用類型；其次是耕地，占到全市面積20%以上；水域所占面積比例最小，4期都在0.7%以下。從生態(tài)系統(tǒng)類型來看，榆林市以草地、耕地、沙地、林地生態(tài)系統(tǒng)為主要特征，這也是風(fēng)沙過渡區(qū)特有的景觀現(xiàn)象。

1990—2015年榆林市土地利用變化很大，總體趨勢(shì)表現(xiàn)為耕地、未利用地減少，而林地、草地、水域、建設(shè)用地面積增加，其中耕地減少最多，達(dá)181 807.54 hm2，由1990年的26.36%下降到2015年的22.19%，未利用地減少33 151.88 hm2，減幅為2.30%；面積增加最多的為草地，25 a間增加了159 315.89 hm2，增幅0.94%，建設(shè)用地面積增加58 134.18 hm2，增加的比例最大，達(dá)84.09%。1990—2000年，林地、草地、建設(shè)用地面積增加，耕地、水域、未利用地面積減少，增幅最大的為建設(shè)用地87.03%，減幅最大的為水域31.22%。2000—2010年、2010—2015年變化趨勢(shì)基本相同，表現(xiàn)為耕地、未利用地減少，而建設(shè)用地面積顯著增加，林地、草地、水域略有增加。在退耕還林政策實(shí)施和城市化進(jìn)程中，榆林市耕地減少、建設(shè)用地顯著增加成為必然趨勢(shì)，耕地轉(zhuǎn)化為草地和建設(shè)用地占用耕地是耕地減少的主要因素。

3.1.2 土地利用動(dòng)態(tài)度分析 根據(jù)公式(2)計(jì)算得出榆林市土地利用變化動(dòng)態(tài)度。榆林市2000年以后土地利用變化的速度有所降低(圖2)。在3個(gè)研究時(shí)段內(nèi)，6類土地利用類型中，耕地的動(dòng)態(tài)度始終最大，1990—2000年、2000—2010年、210—2015年分別達(dá)到8.98，5.46，6.76，水域次之，分別為5.77，2.20，1.89，林地和草地最小，均在0.7以下。

圖1 1990-2015年榆林市土地利用變化

圖2 1990-2015年榆林市土地利用變化動(dòng)態(tài)度

3.2 生態(tài)服務(wù)價(jià)值分析

3.2.1 土地類型的服務(wù)功能價(jià)值變化分析 榆林市1990年、2000年、2010年和2015年的生態(tài)服務(wù)價(jià)值分別為133.84，132.88，132.58，132.35億元，總體上呈現(xiàn)逐年減少的趨勢(shì)，25年間減少了1.49億元(表1)。耕地生態(tài)價(jià)值減少了5.05億元；未利用土地面積持續(xù)減少，生態(tài)價(jià)值減少了0.06億元；林地和草地的生態(tài)價(jià)值分別增加了0.47，4.64億元，不足以彌補(bǔ)耕地和未利用土地生態(tài)價(jià)值的損失，故榆林市生態(tài)價(jià)值呈現(xiàn)了逐年減少的態(tài)勢(shì)。在各種土地利用類型中草地產(chǎn)生的生態(tài)價(jià)值最大，約為總價(jià)值的60%左右，其次是耕地和林地，分別占總價(jià)值的20.32%～23.87%和12.77%～13.27，最少的為未利用地，其產(chǎn)生的生態(tài)價(jià)值僅約占總價(jià)值的0.3%。林地和水域面積約占總土地面積的4.5%和0.5%，但其產(chǎn)生的生態(tài)價(jià)值卻約占總生態(tài)價(jià)值的13%和3%。可見，在干旱風(fēng)沙過渡區(qū)的榆林市林地和水域的生態(tài)效應(yīng)與景觀效應(yīng)是十分明顯的。

表1 1990-2015年榆林市生態(tài)服務(wù)價(jià)值

3.2.2 生態(tài)系統(tǒng)服務(wù)價(jià)值的組成分析 由表2可以看出，生態(tài)系統(tǒng)各項(xiàng)服務(wù)功能產(chǎn)生的價(jià)值大小排序?yàn)椋和寥佬纬膳c保護(hù)>廢物處理>生物多樣性保護(hù)>氣候調(diào)節(jié)>水源涵養(yǎng)>氣體調(diào)節(jié)>食物生產(chǎn)>原材料>娛樂文化。在9項(xiàng)服務(wù)功能中，土壤形成與保護(hù)服務(wù)功能占總價(jià)值比例為23.25%～23.71%，廢物處理服務(wù)價(jià)值所占比例為18.17%～18.65%，生物多樣性保護(hù)服務(wù)價(jià)值所占比例為13.72%～13.97%，氣候調(diào)節(jié)服務(wù)價(jià)值所占比例為12.03%～12.17%，水源涵養(yǎng)服務(wù)價(jià)值所占比例為11.86%～12.20%，氣體調(diào)節(jié)服務(wù)價(jià)值所占比例為10.30%～10.59%，其他各服務(wù)功能的價(jià)值均在6%以下。這也正好反映了處于干旱風(fēng)沙過渡區(qū)的榆林市生態(tài)系統(tǒng)類型的分布格局和服務(wù)功能的構(gòu)成情況，榆林市以草地和耕地生態(tài)系統(tǒng)為主，這兩者的面積占到總面積的88%以上，而草地和耕地的土壤形成與保護(hù)、廢物處理、生物多樣性保護(hù)、氣體調(diào)節(jié)的價(jià)值較大，因此榆林市的生態(tài)服務(wù)價(jià)值主要體現(xiàn)在土壤形成與保護(hù)、廢物處理、生物多樣性保護(hù)上，三者占到總價(jià)值的55%以上。從時(shí)間變化上來看，1990—2015年，氣體調(diào)節(jié)和土壤形成與保護(hù)的服務(wù)價(jià)值明顯增加，而食物生產(chǎn)、廢物處理和水源涵養(yǎng)的服務(wù)價(jià)值顯著減少，其他服務(wù)功能價(jià)值變化不大。造成這一結(jié)果的主要原因是產(chǎn)生較大氣體調(diào)節(jié)和土壤形成與保護(hù)服務(wù)功能的林地和草地面積不斷擴(kuò)大，而產(chǎn)生最大食物生產(chǎn)、廢物處理和水源涵養(yǎng)服務(wù)價(jià)值的耕地和水域面積持續(xù)減少。

表2 1990-2015年榆林市生態(tài)系統(tǒng)各類服務(wù)價(jià)值構(gòu)成

3.2.3 生態(tài)系統(tǒng)服務(wù)價(jià)值的敏感性分析 運(yùn)用公式(4)，將各類生態(tài)類型用地的價(jià)值指數(shù)分別上下調(diào)整50%，計(jì)算了研究區(qū)1990年、2000年、2010年和2015年的生態(tài)系統(tǒng)服務(wù)價(jià)值敏感性指數(shù)(表3)。結(jié)果表明，各類生態(tài)用地的敏感性指數(shù)都小于1，ESV對(duì)VC缺乏彈性，說明研究結(jié)果是可信的，所采用的生態(tài)服務(wù)價(jià)值系數(shù)適合當(dāng)?shù)厍闆r的。1990—2015年林地和草地敏感性指數(shù)呈增長(zhǎng)趨勢(shì)，耕地敏感性指數(shù)呈減少趨勢(shì)，這就表明林地和草地的生態(tài)系統(tǒng)價(jià)值系數(shù)變化會(huì)對(duì)研究區(qū)生態(tài)系統(tǒng)服務(wù)總價(jià)值產(chǎn)生放大作用，耕地生態(tài)系統(tǒng)價(jià)值系數(shù)變化對(duì)生態(tài)系統(tǒng)服務(wù)總價(jià)值產(chǎn)生了縮小的作用。未利用地的敏感性指數(shù)接近于0，位于0.002 8～0.003 2，表明未利用地生態(tài)系統(tǒng)價(jià)值系數(shù)變化對(duì)榆林市生態(tài)系統(tǒng)服務(wù)總價(jià)值的變化影響不大。水域敏感性指數(shù)在1990—2000年呈減少趨勢(shì)，而在2000—2015年呈增加趨勢(shì)，表明水域的生態(tài)系統(tǒng)價(jià)值系數(shù)變化在1990—2000年對(duì)榆林市生態(tài)系統(tǒng)服務(wù)總價(jià)值產(chǎn)生縮小作用，而在2000—2015年對(duì)榆林市生態(tài)系統(tǒng)服務(wù)總價(jià)值產(chǎn)生放大作用。敏感性指數(shù)最高值為草地，為0.590 9～0.632 6，說明草地對(duì)服務(wù)價(jià)值的貢獻(xiàn)率最大，當(dāng)草地生態(tài)系統(tǒng)價(jià)值系數(shù)增加1%時(shí)，生態(tài)系統(tǒng)服務(wù)價(jià)值增加0.5909%～0.6326%，其次為耕地和林地，分別為0.203 2～0.238 7，0.127 7～0.132 7，水域和未利用地均在0.04以下。

表3 榆林市生態(tài)系統(tǒng)服務(wù)價(jià)值的敏感性

4 討論與結(jié)論

近25 a來，受土地沙漠化治理和退耕還林政策直接影響，榆林市土地利用結(jié)構(gòu)發(fā)生了明顯的改變，表現(xiàn)為草地和建設(shè)用地面積顯著增長(zhǎng)，而水域、耕地和未利用地面積大幅減少。這種變化，特別是水域面積減少和建設(shè)用地面積的增加，造成了生態(tài)服務(wù)價(jià)值的減小，尤其是廢物處理、水源涵養(yǎng)、食物生產(chǎn)3項(xiàng)服務(wù)功能顯著減少。水體生態(tài)系統(tǒng)的恢復(fù)和改善，對(duì)生態(tài)系統(tǒng)服務(wù)功能的改善有著非常重要的影響，因而在未來的土地利用中，既要盡量減小氣候干旱化對(duì)生態(tài)系統(tǒng)的影響，更要采取嚴(yán)厲措施防止礦產(chǎn)資源開發(fā)對(duì)水體的破壞和污染。本文在Costanza等[12]生態(tài)系統(tǒng)服務(wù)價(jià)值研究的基礎(chǔ)上，根據(jù)榆林市的實(shí)際情況，對(duì)研究區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量因子進(jìn)行修正，計(jì)算了1990—2015年榆林市的生態(tài)系統(tǒng)服務(wù)價(jià)值。敏感度分析表明生態(tài)系統(tǒng)服務(wù)價(jià)值對(duì)生態(tài)服務(wù)功能價(jià)值指數(shù)是缺乏彈性的，研究結(jié)果可信。本文未計(jì)算建設(shè)用地產(chǎn)生的生態(tài)服務(wù)價(jià)值，根據(jù)余嘉琦等[16]的研究，建設(shè)用地單位面積的生態(tài)服務(wù)價(jià)值當(dāng)量小于0，研究區(qū)近25 a來建設(shè)用地面積增長(zhǎng)幅度很大，因而實(shí)際榆林市生態(tài)系統(tǒng)服務(wù)價(jià)值的損失可能比研究結(jié)果要大。另外，本研究未對(duì)榆林市各縣區(qū)的去年情況進(jìn)行具體分析，這些問題和因素將會(huì)在接下來的工作中繼續(xù)進(jìn)行研究。

參考文獻(xiàn):

[1] Foley J A, DeFries R, Asner G P, et al. Global consequences of land use[J]. Science, 2005,309(5734):570-574.

[2] 李慧燕,常慶瑞,馬廷剛,等.銅川市土地利用/土地覆被動(dòng)態(tài)變化研究[J].水土保持通報(bào),2011,31(2):211-214.

[3] 周秋文,楊勝天,羅亞,等.2001—2009年中國(guó)及周邊國(guó)家土地利用/覆被變化分析[J].世界地理研究,2012,21(2):64-71.

[4] 封建民,王黎,文琦.榆神府礦區(qū)土地利用變化及合理性分析[J].水土保持研究,2015,22(1):188-193.

[5] Denis Ruelland, Antoine Tribotte, Christian Puech, et al. Comparison of methods for LUCC monitoring over 50 years from aerial photograph sand satellite images in a Sahelian catchment[J]. International Journal of Remote Sensing, 2011,32(6):1747-1777.

[6] 張楊,劉艷芳,顧漸萍,等.武漢市土地利用覆被變化與生態(tài)環(huán)境效應(yīng)研究[J].地理科學(xué),2011,31(10):1280-1285.

[7] 封建民,郭玲霞.陜西省神木縣土地利用格局和生態(tài)服務(wù)價(jià)值和生態(tài)服務(wù)價(jià)值變化[J].水土保持通報(bào),2014,34(6):293-298.

[8] 高練,周勇.武漢市土地利用/土地覆蓋變化的生態(tài)環(huán)境效應(yīng)分析[J].農(nóng)業(yè)工程學(xué)報(bào),2008,24(1):73-77.

[9] 貢璐,呂光輝,丁建麗,等.塔里木河上游土地利用變化中的生態(tài)價(jià)值損益分析[J].生態(tài)學(xué)雜志,2006,25(5):526-530.

[10] 吳海珍,阿如旱,郭田保,等.基于RS和GIS的內(nèi)蒙古多倫縣土地利用變化對(duì)生態(tài)服務(wù)價(jià)值的影響[J].地理科學(xué),2011,31(1):110-116.

[11] 周厚俠,望勇.基于RS和GIS的黑河中游生態(tài)服務(wù)價(jià)值對(duì)土地利用變化的影響[J].生態(tài)經(jīng)濟(jì),2016,32(4):183-188.

[12] Costanza R, d′Arge R, de Groot R, et al. The value of the world′s ecosystem service and natural capital[J]. Nature, 1997,386:253-260.

[13] 謝高地,魯春霞.青藏高原生態(tài)資產(chǎn)的價(jià)值評(píng)估[J].自然資源學(xué)報(bào),2003,18(2):189-196.

[14] 文琦,劉彥隨,王建興.生態(tài)脆弱區(qū)土地利用格局演變及其生態(tài)響應(yīng):以榆林市為例[J].地域研究與開發(fā),2010,29(2):104-109.

[15] 祖皮艷木·買買提,玉米提·哈力克,肉孜·阿基,等.基于生態(tài)系統(tǒng)服務(wù)價(jià)值變化的焉耆盆地環(huán)境與經(jīng)濟(jì)協(xié)調(diào)發(fā)展[J].應(yīng)用生態(tài)學(xué)報(bào),2015,26(3):875-883.

[16] 余嘉琦,李鋼,趙華,等.江蘇省沛縣土地利用變化及其生態(tài)服務(wù)價(jià)值研究[J].江蘇農(nóng)業(yè)科學(xué),2015,43(6):371-376.