申建秀，王秀紅，劉 羽，張 偉，姚麗娜

（1.中國(guó)科學(xué)院 地理科學(xué)與資源研究所，北京100101；2.中國(guó)科學(xué)院研究生院，北京100049；3.杭州師范大學(xué)理學(xué)院／遙感與地球科學(xué)研究院，杭州311121；4.浙江省城市濕地與區(qū)域變化研究重點(diǎn)實(shí)驗(yàn)室，杭州311121）

土地利用／土地覆被變化是人與自然關(guān)系最為密切的環(huán)節(jié)，可以引起許多自然現(xiàn)象和生態(tài)過(guò)程的變化，在全球環(huán)境變化中具有重要的作用，對(duì)生態(tài)系統(tǒng)服務(wù)功能起著決定性的作用[1]。分析和評(píng)價(jià)土地利用變化所誘發(fā)的生態(tài)與環(huán)境問(wèn)題既是國(guó)際可持續(xù)發(fā)展研究的熱點(diǎn)之一，也是當(dāng)前土地管理科學(xué)、環(huán)境科學(xué)及生態(tài)科學(xué)的交叉前沿領(lǐng)域[2]。生態(tài)系統(tǒng)服務(wù)是指通過(guò)生態(tài)系統(tǒng)的結(jié)構(gòu)、過(guò)程和功能，直接或間接得到的生命支持產(chǎn)品和服務(wù)[3]，其價(jià)值化已被廣泛運(yùn)用在理解和探討人類的可持續(xù)發(fā)展研究中[4]。我國(guó)的西北退耕還林區(qū)多為生態(tài)環(huán)境脆弱的區(qū)域，自1999年開(kāi)始實(shí)施退耕還林工程以來(lái)，土地利用變化劇烈。深入研究該區(qū)域土地利用變化引起的生態(tài)系統(tǒng)服務(wù)功能價(jià)值的變化，有助于提高我們對(duì)退耕還林政策實(shí)施效果的認(rèn)識(shí)，進(jìn)而為宏觀政策的完善和微觀土地利用行為的選擇提供科學(xué)依據(jù)?，F(xiàn)有對(duì)退耕還林區(qū)生態(tài)系統(tǒng)的研究多從農(nóng)作物生產(chǎn)穩(wěn)定性[5]、生態(tài)格局安全[6]、人口密度變化[7－8]等方面入手進(jìn)行研究，而對(duì)退耕還林工程實(shí)施前后區(qū)域生態(tài)服務(wù)價(jià)值時(shí)空變化研究不足。本文以甘肅省正寧縣為例，對(duì)其退耕還林工程實(shí)施前后（1995—2010年）土地利用變化引起的生態(tài)服務(wù)價(jià)值的時(shí)間變化和空間變化進(jìn)行分析，以期為我國(guó)西北地區(qū)退耕還林政策的完善和區(qū)域土地資源的可持續(xù)利用提供有益的參考。

1 研究區(qū)概況

甘肅省正寧縣位于甘肅省東南端，地處我國(guó)黃土高原溝壑區(qū)，總面積1 319.5km2，介于東經(jīng)107°56′20″—108°38′08″、北緯35°14′40″—35°36′18″。地勢(shì)東北高，西南低，平均海拔1 460m，地形東寬西窄，略呈三角形，東西長(zhǎng)63.5km，支黨河和四郎河流經(jīng)縣域中西部。正寧縣屬溫帶大陸性半濕潤(rùn)氣候，全縣年平均降水量598.9mm，年平均日照時(shí)數(shù)2 682.3 h。全縣有大小不等溝壑1 200多條，約占總面積的82%，主要集中在東北部及子午嶺一帶，這些溝壑侵蝕十分嚴(yán)重。為了進(jìn)一步分析正寧縣ESV空間分異特點(diǎn)，本文結(jié)合正寧縣地形特點(diǎn)和土地利用狀況，以行政鄉(xiāng)界為依據(jù)，將東部子午嶺一帶溝壑嚴(yán)重的月明鄉(xiāng)、西坡鄉(xiāng)、五頃塬回族鄉(xiāng)和三嘉鄉(xiāng)劃分為東部地區(qū)；將農(nóng)村居民點(diǎn)分布廣泛的永正鄉(xiāng)、山河鎮(zhèn)、永和鎮(zhèn)、湫頭鄉(xiāng)以及羅川鄉(xiāng)東部劃分為中部地區(qū)；將地處平坦區(qū)的周家鄉(xiāng)、宮河鎮(zhèn)、榆林子鎮(zhèn)以及羅川鄉(xiāng)西部劃分為西部地區(qū)。

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

研究區(qū)1995年、2000年兩期土地利用圖來(lái)自中國(guó)科學(xué)院資源與環(huán)境數(shù)據(jù)庫(kù)中1∶10萬(wàn)土地利用數(shù)據(jù)；2010年Landsat TM影像來(lái)自美國(guó)地質(zhì)調(diào)查局（USGS）網(wǎng)站。在ENVI和ArcGIS軟件平臺(tái)的支持下，將2010年遙感影像與已有的1995年和2000年土地利用類型圖進(jìn)行坐標(biāo)匹配，實(shí)現(xiàn)縣界與圖像套合。并參考1995年、2000年土地利用分類方案，結(jié)合正寧縣的實(shí)際情況，將研究區(qū)土地利用類型劃分為耕地、林地、草地、水域、工礦及居民點(diǎn)用地（建設(shè)用地）、未利用地6種一級(jí)地類以及相應(yīng)的16種二級(jí)地類。根據(jù)土地利用調(diào)查規(guī)程的遙感解譯最小圖上面積規(guī)定，通過(guò)專家目視人工解譯成圖、建庫(kù)，得到研究區(qū)2010年土地利用現(xiàn)狀；研究區(qū)1995年、2000年、2010年糧食單產(chǎn)量的統(tǒng)計(jì)資料來(lái)源于甘肅農(nóng)村年鑒（1996—2009）；全國(guó)2000年凈初級(jí)凈生產(chǎn)力（Net Primary Productivity，NPP）數(shù)據(jù)來(lái)自于中國(guó)科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心。

3 研究方法

3.1 研究區(qū)生態(tài)系統(tǒng)服務(wù)功能價(jià)值的核算

區(qū)域生態(tài)服務(wù)價(jià)值是定量表征區(qū)域生態(tài)系統(tǒng)服務(wù)功能的重要指標(biāo)[9]。本文以Costanza等的研究成果為基礎(chǔ)[3]，根據(jù)相關(guān)專家對(duì)我國(guó)陸地生態(tài)系統(tǒng)服務(wù)價(jià)值的實(shí)際計(jì)算結(jié)果[10]，計(jì)算研究區(qū)各土地覆被類型單位面積的生態(tài)價(jià)值。由于中國(guó)不同陸地生態(tài)系統(tǒng)單位面積生態(tài)服務(wù)價(jià)值表是以全國(guó)平均水平狀況下的農(nóng)田每年自然糧食產(chǎn)量的經(jīng)濟(jì)價(jià)值為基準(zhǔn)的，直接用于研究區(qū)3期生態(tài)服務(wù)價(jià)值的計(jì)算并不合理，因此本研究在謝高地等[11]的研究成果基礎(chǔ)上進(jìn)行了進(jìn)一步的空間、時(shí)間和地類的校正，使之更適合研究區(qū)域不同年份相應(yīng)地類的實(shí)際情況。

（1）生態(tài)系統(tǒng)服務(wù)價(jià)值核算方法。計(jì)算區(qū)域生態(tài)系統(tǒng)服務(wù)價(jià)值的公式如下[12－13]：

式中：ESV——生態(tài)系統(tǒng)服務(wù)價(jià)值；Ak——研究區(qū)第k種土地利用類型的分布面積；VCk——生態(tài)服務(wù)功能價(jià)值指數(shù)，即經(jīng)過(guò)時(shí)空和地類校正后的單位面積的生態(tài)系統(tǒng)服務(wù)價(jià)值。

（2）生態(tài)系統(tǒng)服務(wù)價(jià)值參數(shù)。第一是空間校正。根據(jù)師慶三等的研究成果，NPP與生態(tài)系統(tǒng)服務(wù)價(jià)值之間具有良好的線性關(guān)系[14]，因此，本文以空間校正是以全國(guó)2000年NPP為基礎(chǔ)，在ArcGIS支持下，在配準(zhǔn)后的全國(guó)數(shù)據(jù)基礎(chǔ)上進(jìn)行裁切獲得研究區(qū)數(shù)據(jù)，進(jìn)而進(jìn)行全國(guó)和研究區(qū)生態(tài)系統(tǒng)平均生物量的計(jì)算；第二是時(shí)間校正。因糧食產(chǎn)出量與生態(tài)系統(tǒng)各要素關(guān)系密切，具有生態(tài)服務(wù)價(jià)值代表性，所以本文通過(guò)1995—2000年、1995—2005年、2000—2010年3個(gè)時(shí)間段糧食單產(chǎn)均值的比值來(lái)確定各期生態(tài)系統(tǒng)單位面積生態(tài)服務(wù)價(jià)值的校正系數(shù)；第三是地類校正。由于研究區(qū)是我國(guó)西北黃土溝壑區(qū)退耕還林的代表區(qū)域，林地、草地的二級(jí)地類對(duì)生態(tài)系統(tǒng)生態(tài)服務(wù)價(jià)值的影響較大，本文適當(dāng)調(diào)整了相關(guān)文獻(xiàn)中的計(jì)算參數(shù)。其中，有林地、灌木林地、其他林地的調(diào)整系數(shù)依據(jù)唐秀美等[15]經(jīng)過(guò)專家問(wèn)卷調(diào)查法得到的林地生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)調(diào)整表進(jìn)行校正；高覆蓋率草地、中覆蓋率草地、低覆蓋率草地的調(diào)整是在把高覆蓋率草地系數(shù)設(shè)為1的基礎(chǔ)上，根據(jù)相應(yīng)草地的地類覆蓋率系數(shù)[16]比值得到。

3.2 生態(tài)服務(wù)價(jià)值質(zhì)心計(jì)算

由于地理事物分布的分散性，無(wú)法以它們的實(shí)際位置計(jì)算其重心，通常將研究區(qū)劃分為若干小區(qū)，以各小區(qū)代表點(diǎn)坐標(biāo)作為該小區(qū)的位置。本文中根據(jù)研究區(qū)土地利用劃分為耕地、林地、草地、水域、建設(shè)用地5個(gè)區(qū)，利用ArcGIS中Spatial Analyst Tool模塊的Zonal Geometry as Table命令分別計(jì)算出3期各地類的景觀幾何中心坐標(biāo)和研究區(qū)的幾何中心，根據(jù)地理事物重心的計(jì)算公式[17－18]計(jì)算研究區(qū)ESV重心：

3.3 敏感性分析

參照相關(guān)研究成果[2]，應(yīng)用經(jīng)濟(jì)學(xué)中常用的彈性系數(shù)概念來(lái)計(jì)算價(jià)值系數(shù)的敏感性指數(shù)（coefficient of sensitivity，CS），以確定生態(tài)系統(tǒng)服務(wù)價(jià)值隨時(shí)間的變化情況對(duì)價(jià)值系數(shù)變化的依賴程度。如果CS＞1，表明ESV相對(duì)于VC是富有彈性的；如果CS＜1，ESV則被認(rèn)為是缺乏彈性的，比值越大，表明生態(tài)服務(wù)功能價(jià)值指數(shù)的準(zhǔn)確性越關(guān)鍵。計(jì)算公式為：

式中：ESV——估算的總生態(tài)系統(tǒng)服務(wù)價(jià)值；VC——生態(tài)價(jià)值系數(shù)；i，j——初始總價(jià)值和生態(tài)價(jià)值系數(shù)調(diào)整以后的總價(jià)值；k——各土地利用類型。

4 結(jié)果與分析

4.1 土地利用變化分析

由正寧縣退耕還林前后15a的土地利用結(jié)構(gòu)變化情況來(lái)看（表1），耕地、林地、草地為研究區(qū)的主要地類，其中耕地占總面積比重最大，說(shuō)明正寧縣是以種植業(yè)為主導(dǎo)發(fā)展的。耕地和水域面積一直在減小，且2000—2010年減小的幅度明顯大于1995—2000年的減小幅度。其中，耕地在退耕前表現(xiàn)為輕微減少趨勢(shì)，由46.5%減少到45.5%，在退耕后耕地面積又大幅度降低到2010年的41.3%，前后共減少了7 072.9hm2。在退耕還林工程實(shí)施前的1995年到2000年，面積減小的地類為有林地、灌木林地、水域和耕地，灌木林地減小速度最大，為－3.46%。同時(shí)，面積增加的地類為其他林地、草地、建設(shè)用地。因黃土溝壑區(qū)草地主要來(lái)源于未成林的林地和退化后的耕地、水域，建設(shè)用地的增加占用的是平原耕地，說(shuō)明這期間研究區(qū)的生態(tài)環(huán)境受到了一定程度的干擾。退耕還林工程實(shí)施后的2000年到2010年，有林地以9.16%的年均速度大幅度增加，表明退耕還林效果明顯。建設(shè)用地比例較小，但一直保持一定的增長(zhǎng)速度。

表1 正寧縣土地利用面積變化和年變化率 %

綜上，從整個(gè)研究期來(lái)看，退耕還林工程的實(shí)施和人口增長(zhǎng)、人們生活生產(chǎn)條件的不斷改善促使當(dāng)?shù)馗睾退蛎娣e不斷減少，建設(shè)用地不斷增加，同時(shí)林地面積的變化表現(xiàn)為先減后增。

4.2 生態(tài)系統(tǒng)服務(wù)價(jià)值的時(shí)序變化分析

（1）各地類生態(tài)系統(tǒng)服務(wù)價(jià)值的變化。在退耕還林工程實(shí)施前的1995—2000年間，林地、耕地、水域、的生態(tài)服務(wù)價(jià)值明顯減少，其中灌木林地變化最為顯著，5a間減少了19.97%，且在研究區(qū)內(nèi)的ESV價(jià)值比例下降了4.07%。耕地面積的減少導(dǎo)致了耕地生態(tài)服務(wù)價(jià)值的減小，減小幅度為5.39%，但其在比重上有所增加，歸因于毀林造田過(guò)程中林地生態(tài)服務(wù)價(jià)值比重的減小。正寧縣ESV總價(jià)值由1995年的15.523億元／a減小到2000年的14.512億元／a，共減少了6.512%，說(shuō)明在退耕還林工程實(shí)施前，正寧縣的土地利用變化不利于整體生態(tài)服務(wù)價(jià)值的提高。退耕還林工程實(shí)施后的2000—2010年間，有林地ESV增加明顯，價(jià)值幅度上增加了117.98%，比重上增加了17.23%。生態(tài)條件的改善保證了耕地和水域ESV的增加。由于耕地和水域面積的不斷減少，以及林地ESV的顯著增加，二者在正寧縣總ESV中的比重也有所下降。從整個(gè)研究期來(lái)看，正寧縣總生態(tài)服務(wù)價(jià)值先減少（－6.512%）后增加（29.369%），15a間的增加幅度為20.94%，說(shuō)明正寧縣在實(shí)施退耕還林工程以后，土地利用變化向著可持續(xù)性增強(qiáng)的方向發(fā)展。

我國(guó)競(jìng)爭(zhēng)中立制度的實(shí)施策略應(yīng)采取一種逐步推進(jìn)的模式:首先在自由貿(mào)易試驗(yàn)區(qū)進(jìn)行試點(diǎn)實(shí)踐，然后在全國(guó)范圍內(nèi)逐步推廣，最后在國(guó)際經(jīng)貿(mào)協(xié)定中廣泛適用。

表2 1995年、2000年、2010年正寧縣各地類生態(tài)服務(wù)價(jià)值變化

圖1 1995年、2000年、2010年正寧縣各單項(xiàng)生態(tài)服務(wù)價(jià)值構(gòu)成情況

（2）單項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值（ESVj）的變化。在退耕還林工程實(shí)施前的1995—2000年間，所有單項(xiàng)生態(tài)服務(wù)價(jià)值均不同程度的減少，主要由林地面積的減少和耕地退化所致。退耕還林工程實(shí)施后的2000—2010年間，除食物生產(chǎn)功能外，所有單項(xiàng)功能價(jià)值均明顯增加，比重最大的是土壤形成與保護(hù)功能價(jià)值，為18.67%；幅度最大的是娛樂(lè)文化功能價(jià)值，為48.35%（圖1）。在正寧縣3期單項(xiàng)生態(tài)服務(wù)價(jià)值中，土壤形成與保護(hù)功能和生物多樣性功能比重最大，變化幅度劇烈；娛樂(lè)文化功能比重最小但變化幅度很大；食物生產(chǎn)功能價(jià)值比重小且隨土地利用變化而變化的幅度也小。說(shuō)明地處黃土溝壑區(qū)的正寧縣在進(jìn)行土地利用時(shí)已經(jīng)把生態(tài)價(jià)值放在重要位置。

4.3 生態(tài)系統(tǒng)服務(wù)價(jià)值的空間分異

由2010年正寧縣單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值空間分異圖（圖2）可以看出：正寧縣東部溝壑嚴(yán)重地區(qū)生態(tài)服務(wù)主要由林地和耕地提供。相對(duì)于中部地區(qū)和西部地區(qū)，其有林地和灌木林地的比較優(yōu)勢(shì)最為突出，這說(shuō)明該區(qū)今后應(yīng)以生態(tài)保護(hù)為主，發(fā)展成為正寧縣的生態(tài)涵養(yǎng)區(qū)。中部地區(qū)與西南部地區(qū)的情況比較相似，主要表現(xiàn)為耕地的生態(tài)服務(wù)優(yōu)勢(shì)較大。因中部地區(qū)人口密度相對(duì)大，又有支黨河、四郎河過(guò)境，應(yīng)加強(qiáng)中部地區(qū)的農(nóng)田水利基礎(chǔ)設(shè)施建設(shè)，保證農(nóng)業(yè)縣的糧食安全。西部地區(qū)地勢(shì)相對(duì)平坦，水源充足，但因正寧縣西南地區(qū)耕地被開(kāi)墾因而無(wú)法充分利用，林地又無(wú)法被有效保護(hù)的現(xiàn)象很多，致使各項(xiàng)生態(tài)服務(wù)價(jià)值均偏低，所以，應(yīng)注重對(duì)該區(qū)域內(nèi)土地利用結(jié)構(gòu)的優(yōu)化配置，加大退耕還林工程的實(shí)施力度，從而發(fā)揮西南部最大綜合效益。

正寧縣西、東、中部地區(qū)不同的土地利用空間組合造成了3個(gè)地區(qū)不同的功能分異，也使得不同地區(qū)單位面積的生態(tài)服務(wù)價(jià)值差異明顯。退耕還林前后的3個(gè)年期，正寧縣單位面積的ESV按區(qū)排序均為東部地區(qū)＞中部地區(qū)＞西部地區(qū)（表3），因?yàn)檠芯繀^(qū)東北山地溝壑區(qū)的地類受人為干擾程度要小于西南相對(duì)平坦地區(qū)，東部地區(qū)有林地、灌木林地若能得到有效保護(hù)，就能夠充分發(fā)揮其生態(tài)服務(wù)功能。中部地區(qū)人口相對(duì)集中，對(duì)地類的干擾程度較大，單位面積ESV低于受人類干擾程度小的東部地區(qū)，通過(guò)確保該地區(qū)耕地的充分利用和合理開(kāi)發(fā)，才能充分發(fā)揮該地區(qū)的生態(tài)服務(wù)功能。西部地區(qū)土地利用結(jié)構(gòu)沒(méi)有得到科學(xué)配置，削減了其單位面積的生態(tài)服務(wù)價(jià)值。

圖2 2010年正寧縣生態(tài)系統(tǒng)服務(wù)價(jià)值的空間分異

表3 正寧縣西、中、東部地區(qū)單位面積ESV變化

4.4 生態(tài)系統(tǒng)服務(wù)價(jià)值的質(zhì)心變化

由正寧縣退耕還林前后生態(tài)服務(wù)價(jià)值質(zhì)心轉(zhuǎn)移跟蹤方位示意圖（圖3）可知，正寧縣生態(tài)服務(wù)價(jià)值質(zhì)心在整個(gè)研究期內(nèi)先向西南方向遷移，后又向東南方向遷移。1995—2000年向西南方向遷移了323.78 m，2000年至2010年向東南方向偏移了486.50m，1995—2000年向東南偏移了272.45m，整體表現(xiàn)出“東進(jìn)南下”的趨勢(shì)。

圖3 退耕還林前后正寧縣生態(tài)服務(wù)價(jià)值質(zhì)心轉(zhuǎn)移跟蹤方位示意圖

研究區(qū)的生態(tài)服務(wù)價(jià)值質(zhì)心與其幾何中心并不匹配，這表明研究區(qū)的生態(tài)服務(wù)價(jià)值在空間分布上不均衡，可以通過(guò)正寧縣西南部、中、東部地區(qū)的3a期單位面積ESV變化來(lái)體現(xiàn)（表3）。

追蹤正寧縣ESV質(zhì)心的遷移，在退耕還林工程實(shí)施前，因不合理的土地利用，使得正寧縣各區(qū)單位面積生態(tài)服務(wù)價(jià)值均下降，東部地區(qū)減小8.24%，明顯大于西部地區(qū)和中部地區(qū)的凈減小幅度，使得正寧縣的ESV質(zhì)心向西南方向移動(dòng)，是一種不可持續(xù)的遷移。退耕還林工程實(shí)施后，各區(qū)單位面積生態(tài)服務(wù)價(jià)值增加明顯，東部響應(yīng)最大，凈增加32.74%，西部?jī)粼黾?6.04%，因各區(qū)單位面積ESV增加幅度上的不平衡，使得質(zhì)心在退耕還林工程實(shí)施后向增加幅度大的東南方向遷移，土地利用配置并沒(méi)有達(dá)到最優(yōu)。

4.5 敏感性指數(shù)分析

根據(jù)CS計(jì)算公式，本文將各地類土地利用類型的價(jià)值指數(shù)分別調(diào)整50%，計(jì)算出研究區(qū)1995，2000，2010年3期各地類的敏感性指數(shù)，來(lái)衡量總生態(tài)系統(tǒng)服務(wù)價(jià)值的變化。結(jié)果表明，ESV對(duì)VC的敏感性指數(shù)均小于1，最高值是2010年期的水域0.088。這表明在2010年，當(dāng)水域的VC增加1%時(shí)，相應(yīng)的ESV增加0.008%。敏感性分析的結(jié)果表明，研究區(qū)內(nèi)ESV對(duì)VC缺乏彈性，基于謝高地[11]價(jià)值系數(shù)校正的生態(tài)系統(tǒng)服務(wù)功能價(jià)值核算結(jié)果在本區(qū)域內(nèi)是可信的。

5 結(jié)論

（1）在退耕還林工程實(shí)施前的1995—2000年，面積減少的地類為有林地、灌木林地、水域和耕地，面積增加的地類為其他林地、草地、建設(shè)用地。退耕還林工程實(shí)施后，有林地面積大幅度增加，是退耕還林的明顯成效。

（2）正寧縣在實(shí)施退耕還林工程前后，總生態(tài)服務(wù)價(jià)值表現(xiàn)為先減少（－6.512%）后增加（29.369%）。表現(xiàn)在地類上為：退耕前各地類ESV比重大小排序?yàn)楦兀竟嗄玖值兀居辛值兀局懈采w率草地＞其他林地＞高覆蓋率草地＞水域＞低覆蓋率草地＞建設(shè)用地，退耕后各地類ESV比重大小排序?yàn)橛辛值兀靖兀竟嗄玖值兀酒渌值兀局懈采w率草地＞水域＞低覆蓋率草地＞建設(shè)用地，有林地在比重和幅度上的顯著提高是退耕還林工程的體現(xiàn)。由于地處黃土溝壑區(qū)，正寧縣所提供的生態(tài)服務(wù)功能以生物多樣性保護(hù)、土壤形成與保護(hù)、氣體調(diào)節(jié)、水源涵養(yǎng)為主，食物生產(chǎn)和娛樂(lè)文化功能比重最小。

（3）追蹤正寧縣生態(tài)服務(wù)價(jià)值質(zhì)心的遷移，在退耕還林工程實(shí)施前的1995—2000年向西南方向偏移，退耕工程實(shí)施后的2000—2010年向東南方向偏移，縣域單位面積的ESV明顯增加。正寧縣應(yīng)在保持已有的退耕還林成果基礎(chǔ)上，進(jìn)一步加強(qiáng)中部地區(qū)農(nóng)業(yè)基礎(chǔ)設(shè)施的建設(shè)，控制建設(shè)用地?cái)U(kuò)張，優(yōu)化種植結(jié)構(gòu)，擴(kuò)大中部地區(qū)耕地的生態(tài)服務(wù)價(jià)值。此外，應(yīng)重視西南部邊際耕地的退耕和防護(hù)林工程的建設(shè)，優(yōu)化土地利用結(jié)構(gòu)，發(fā)揮西南部最大的綜合效益。

（4）敏感性分析的結(jié)果表明，研究區(qū)內(nèi)的生態(tài)系統(tǒng)服務(wù)價(jià)值（ESV）對(duì)于生態(tài)價(jià)值系數(shù)（VC）的敏感性指數(shù)均小于1，說(shuō)明基于謝高地價(jià)值系數(shù)校正后的生態(tài)系統(tǒng)服務(wù)功能價(jià)值核算結(jié)果在本區(qū)域內(nèi)具有較高的可信度。

[1]Semwal R L，Nautiyal S，Sen K K，et al.Patterns and ecological implications of agricultural land－use changes：a case study from central Himalaya，India[J].Agriculture，Ecosystems＆Environment，2004，102（1）：81－92.

[2]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.

[3]Costanza R，d′Arge R，De Groot R，et al.The value of the world＇s ecosystem services and natural capital[J].Nature，1997，387（6630）：253－260.

[4]Bingham G，Bishop R，Brody M，et al.Issues in ecosystem valuation：improving information for decision making[J].Ecological Economics，1995，14（2）：73－90.

[5]張應(yīng)龍，謝永生，李曉，等.黃土溝壑區(qū)主要糧食作物生產(chǎn)穩(wěn)定性及成本效益分析[J].水土保持通報(bào)，2010，30（4）：201－204.

[6]付艷玲，張曉萍，陳鳳娟，等.黃土區(qū)退耕背景下土地利用／覆被格局現(xiàn)狀調(diào)查研究[J].水土保持研究，2010，17（6）：81－86.

[7]王海鴻，馬瓊，付士波，等.西北干旱半干旱區(qū)耕地面積變化與人口—經(jīng)濟(jì)發(fā)展的相關(guān)關(guān)系研究：以甘肅省為例[J].干旱區(qū)資源與環(huán)境，2011，25（1）：74－79.

[8]Shoshany M，Goldshleger N.Land－use and population density changes in Israel—1950to 1990：analysis of regional and local trends[J].Land Use Policy，2002，19（2）：123－133.

[9]冉圣宏，李秀彬，呂昌河.土地覆被及生態(tài)服務(wù)價(jià)值變化的多時(shí)間尺度模擬：以四川省漁子溪流域?yàn)槔齕J].地理學(xué)報(bào)，2006，61（10）：1113－1120.

[10]陳仲新，張新時(shí).中國(guó)生態(tài)系統(tǒng)效益的價(jià)值[J].科學(xué)通報(bào)，2000，45（1）：17－22.

[11]謝高地，魯春霞，冷允法，等.青藏高原生態(tài)資產(chǎn)的價(jià)值評(píng)估[J].自然資源學(xué)報(bào)，2003，18（2）：189－196.

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

[13]宋宏利，張曉楠，倫更永.冀南土地利用變化對(duì)區(qū)域生態(tài)服務(wù)價(jià)值的影響分析[J].水土保持研究，2011，18（1）：236－238.

[14]師慶三，王智，吳友均，等.新疆生態(tài)系統(tǒng)服務(wù)價(jià)值測(cè)算與 NPP的相關(guān)性分析[J].干旱區(qū)地理，2010，33（3）：427－432.

[15]唐秀美.北京市土地利用生態(tài)分類與格局變化研究[D].北京：中國(guó)科學(xué)院研究生院，2011.

[16]中華人民共和國(guó)土資源部.第二次全國(guó)土地調(diào)查技術(shù)規(guī)程（中華人民共和國(guó)土地管理行業(yè)標(biāo)準(zhǔn)[S].TD／T1014—2007）.北京：2007.

[17]徐映雪，任志遠(yuǎn)，張文梅.基于GIS的陜西省耕地與人口空間變化分析[J].干旱地區(qū)農(nóng)業(yè)研究，2008，26（3）：216－221.

[18]許月卿，李雙成.我國(guó)人口與社會(huì)經(jīng)濟(jì)重心的動(dòng)態(tài)演變[J].人文地理，2005（1）：117－120.

[19]李秀彬.地區(qū)發(fā)展均衡性的可視化測(cè)度[J].地理科學(xué)，1999，19（3）：254－257.

[20]He Y，Chen Y，Tang H，et al.Exploring spatial change and gravity center movement for ecosystem services value using a spatially explicit ecosystem services value index and gravity model[J].Environmental Monitoring and Assessment，2011，175（1）：563－571.