陜西關(guān)中地區(qū)土地生態(tài)價(jià)值時(shí)空動(dòng)態(tài)

莫宏偉，全 斌，袁開(kāi)國(guó)，謝金寧，向云波

(湖南科技大學(xué)資源環(huán)境與安全工程學(xué)院， 湖南 湘潭 411201)

陜西關(guān)中地區(qū)土地生態(tài)價(jià)值時(shí)空動(dòng)態(tài)

莫宏偉，全 斌，袁開(kāi)國(guó)，謝金寧，向云波

(湖南科技大學(xué)資源環(huán)境與安全工程學(xué)院， 湖南 湘潭 411201)

為探明土地生態(tài)系統(tǒng)服務(wù)價(jià)值的空間差異性，以關(guān)中地區(qū)遙感影像、氣象氣候、土壤植被及地形地貌數(shù)據(jù)為依據(jù)，借助格網(wǎng)插值技術(shù)和土壤潛力測(cè)算模型獲取土壤潛力柵格圖，運(yùn)用土地利用柵格圖及生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量因子加權(quán)，生成土地生態(tài)系統(tǒng)服務(wù)價(jià)值柵格圖；從而在柵格尺度上揭示土地生態(tài)系統(tǒng)服務(wù)價(jià)值的空間異質(zhì)性。結(jié)果表明：1986—2007年間，關(guān)中地區(qū)土地生態(tài)系統(tǒng)服務(wù)價(jià)值總量由509.61億元增加到561.02億元，增幅約為10.09%，2000—2007年段增速約為1986—2000年段的33.65倍；周邊山地和黃土梁峁區(qū)生態(tài)價(jià)值增加，增量最大的為周至縣；中部平原區(qū)的生態(tài)價(jià)值量減少，減量最大的是西安市區(qū)；研究期內(nèi)的土地利用結(jié)構(gòu)調(diào)整取得了較明顯的生態(tài)效益。

關(guān)中地區(qū)；GIS；土地生態(tài)價(jià)值；空間異質(zhì)性

土地生態(tài)系統(tǒng)服務(wù)是指人類從土地生態(tài)系統(tǒng)獲得的所有惠益，包括供給服務(wù)、調(diào)節(jié)服務(wù)、文化服務(wù)、支持服務(wù)。大多數(shù)土地生態(tài)系統(tǒng)服務(wù)是公共品或準(zhǔn)公共品，具有外部經(jīng)濟(jì)性，無(wú)法進(jìn)入市場(chǎng)，甚至在市場(chǎng)交易中很難發(fā)現(xiàn)對(duì)應(yīng)的補(bǔ)償措施[1]。土地生態(tài)系統(tǒng)服務(wù)價(jià)值研究有利于制定合理的土地生態(tài)資源價(jià)格，促進(jìn)環(huán)境損益進(jìn)入國(guó)民經(jīng)濟(jì)核算體系，為土地生態(tài)建設(shè)與保護(hù)找到合理的資金來(lái)源，達(dá)到利用經(jīng)濟(jì)手段保障土地資源的可持續(xù)利用的目的。自1997年Costanza在《nature》雜志發(fā)表了“全球生態(tài)系統(tǒng)服務(wù)價(jià)值和自然資本”一文以來(lái)[2]，李金昌[3]、歐陽(yáng)志云[4]、陳仲新[5]、肖寒[6]、趙景柱[7]、姜文來(lái)[8]、傅伯杰[9]等眾多學(xué)者從多個(gè)角度在該領(lǐng)域做了許多建設(shè)性的工作；其中，謝高地等在全球生態(tài)系統(tǒng)服務(wù)功能評(píng)價(jià)模型的基礎(chǔ)上，總結(jié)了氣體調(diào)節(jié)、氣候調(diào)節(jié)、水源涵養(yǎng)、土壤形成與保護(hù)、廢物處理、生物多樣性維持、食物生產(chǎn)、原材料生產(chǎn)、休閑娛樂(lè)在內(nèi)的9項(xiàng)生態(tài)系統(tǒng)服務(wù)功能，并對(duì)我國(guó)700位生態(tài)學(xué)者進(jìn)行問(wèn)卷調(diào)查，得到了“中國(guó)生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量因子表”[10-13]；從而把不同生態(tài)系統(tǒng)之間的生態(tài)價(jià)值當(dāng)量區(qū)別開(kāi)來(lái)，但同類型生態(tài)系統(tǒng)價(jià)值的空間異質(zhì)性問(wèn)題仍需更多研究案例來(lái)推進(jìn)。關(guān)中地區(qū)是中國(guó)西部經(jīng)濟(jì)文化中心，涵蓋高原、盆地、山地等多個(gè)地形區(qū)，退耕還林、城市化等全局性的土地利用變化在本區(qū)均有典型體現(xiàn)；隨著西部開(kāi)發(fā)進(jìn)程提速，區(qū)內(nèi)環(huán)境、經(jīng)濟(jì)、人口壓力加大，區(qū)域土地生態(tài)系統(tǒng)的結(jié)構(gòu)和功能變化加??；深入研究土地生態(tài)系統(tǒng)服務(wù)價(jià)值的空間異質(zhì)性，有利于利用經(jīng)濟(jì)手段精準(zhǔn)優(yōu)化土地利用結(jié)構(gòu)，達(dá)到土地資源經(jīng)濟(jì)和生態(tài)雙贏的利用效果，從而為保障區(qū)域土地資源可持續(xù)利用提供決策支持。

1 研究區(qū)概況

關(guān)中地區(qū)位于北緯33°34′～35°52′和東經(jīng)106°18′～110°38′之間，東西長(zhǎng)約360 km，南北寬約170 km，總面積約5.56萬(wàn)km2；轄西安、咸陽(yáng)、寶雞、渭南、銅川五個(gè)地級(jí)市和楊陵區(qū)；東部與河南、山西相鄰，西部與甘肅省接壤，南部是漢中、安康及商洛市，北部為延安市；研究區(qū)中部地勢(shì)低平，南北及西面地勢(shì)較高，有黃土梁峁、黃土臺(tái)塬、黃土塬、平原、山地等五種地貌類型，其面積分別占全區(qū)總面積的6.03%、18.54%、11.91%、20.99%、42.53%；關(guān)中地區(qū)地貌類型及行政區(qū)劃見(jiàn)圖1。本區(qū)屬大陸性季風(fēng)氣候，處于暖溫帶半濕潤(rùn)與半干旱氣候的過(guò)渡地帶，冬冷夏熱、四季分明、降水集中、雨熱同季、易發(fā)生干旱；年平均氣溫為9.9℃～15.8℃，多年平均降雨量500～700 mm，多年平均蒸發(fā)量1 000～1 200 mm。

2 數(shù)據(jù)來(lái)源及研究方法

以遙感影像、氣象、土壤及DEM數(shù)據(jù)為依據(jù)，以ERDAS9.2和ARCGIS9.2作為數(shù)據(jù)處理平臺(tái)，分別獲取關(guān)中地區(qū)土地柵格圖、各地類單位生產(chǎn)力價(jià)值量柵格圖、土壤潛力柵格圖，通過(guò)地圖運(yùn)算得到土地生態(tài)系統(tǒng)服務(wù)價(jià)值柵格圖，從而把土地生態(tài)系統(tǒng)服務(wù)價(jià)值細(xì)化到柵格單元上，初步解決了土地生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)價(jià)的空間異質(zhì)性問(wèn)題，各類柵格圖獲取的具體方法如下：

關(guān)中地區(qū)土地柵格圖是以1986、2000年及2007年的三期TM影像作為基本信息源，采用土地利用二級(jí)分類方法，通過(guò)計(jì)算機(jī)解譯和人工目視解譯相結(jié)合的方式獲取。首先利用1∶5萬(wàn)地形圖對(duì)三期TM影像作幾何精校正，采用回歸分析法并利用DEM對(duì)影像做輻射校正，把校正后的影像轉(zhuǎn)換成Albers投影；進(jìn)行圖像增強(qiáng)處理后，根據(jù)影像特征及野外考察數(shù)據(jù)建立解譯標(biāo)志，運(yùn)用非監(jiān)督和監(jiān)督分類相結(jié)合的辦法對(duì)影像進(jìn)行分類，并對(duì)分類結(jié)果作后處理，最后將土地利用類型歸并為耕地、林地、草地、水域、建設(shè)用地和未利用地；再以野外GPS定位考察點(diǎn)資料對(duì)分類結(jié)果進(jìn)行精度檢驗(yàn)和修改，最終得到關(guān)中地區(qū)三期土地利用現(xiàn)狀圖；再依據(jù)謝高地等[12]2007年編制的新版中國(guó)生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量因子表，并參照陜西省糧食局2009年的糧食價(jià)格(1.60元/kg)，計(jì)算出關(guān)中地區(qū)耕地、林地、草地、水域、建設(shè)用地、未利用地單位生產(chǎn)潛力價(jià)值分別為1 600、5 695、2 364、9 185、0、282 元·t-1，此后依據(jù)土地利用類型圖，生成各地類單位生產(chǎn)力生態(tài)價(jià)值量柵格圖。

關(guān)中地區(qū)土壤潛力分布圖是以氣象數(shù)據(jù)、土壤數(shù)據(jù)、DEM數(shù)據(jù)為依據(jù)，先將研究區(qū)43個(gè)氣象站點(diǎn)各月的輻射、海平面平均氣溫、降水量、蒸發(fā)量近30年的均值進(jìn)行Kriging內(nèi)插，獲取各氣候因子1—12月30 m×30 m的柵格數(shù)據(jù)；然后利用30 m×30 m DEM對(duì)月平均氣溫作高程矯正，獲取實(shí)際地形氣溫分布圖；再依據(jù)溫度、水分及土壤肥力訂正系數(shù)模型獲取關(guān)中地區(qū)土地生產(chǎn)潛力各月的溫度、水分及肥力訂正系數(shù)分布圖；最后依據(jù)光溫階乘模型得出各月土壤潛力圖[14]，再合成研究區(qū)年度土壤潛力柵格圖(見(jiàn)圖2)；具體測(cè)算模型如下：

(1) 光合生產(chǎn)潛力[15-17]

(1)

式中，Y為光合潛力(g·cm-2)；ε為光合有效輻射(0.45～0.55)；α為植物葉面反射率；β為漏射率；γ為光飽和限制率；ω為植物呼吸作用損耗率(約為30%)；ρ為植物非光合器官的無(wú)效吸收，約占光合有效輻射的10%；ψ為量子效率(約為22.4%)；h為單位干物質(zhì)發(fā)熱量(J·g-1)；∑Q為各月太陽(yáng)總輻射；CA為植物體的灰分含量(約為5%)；CM為植物體的含水率(約為14%)。

(2) 光溫生產(chǎn)潛力[19-20]

Yt=Y×f(t)

(2)

(3)

式中，Yt為光溫潛力；f(t)為溫度訂正系數(shù)；t為月平均氣溫(℃)。

(3) 氣候生產(chǎn)潛力[21-23]

Yw=Yt×f(w)

(4)

(5)

E0=0.0018×(t+25)2×(100-f)

(6)

式中，Yw為氣候潛力；f(w)為水分訂正系數(shù)；R為月降水量；E0為月蒸發(fā)量(mm)；t為月平均氣溫(℃)；f為月相對(duì)濕度(%)。

(4) 土壤生產(chǎn)潛力[20,24]

Ys=Yw×f(s)=Yw×U×O

(7)

(8)

(9)

式中，Ys為土壤潛力；f(s)為土壤肥力訂正系數(shù)；O為土壤中有機(jī)質(zhì)訂正系數(shù)；U為土壤中氮磷鉀等養(yǎng)分含量訂正系數(shù)；N、P、K分別為土壤中的堿解氮、速效磷和速效鉀含量(mg·kg-1)；M為土壤中有機(jī)質(zhì)含量(g·kg-1)。

將土壤潛力柵格圖與各地類單位生產(chǎn)力生態(tài)價(jià)值量柵格圖作地圖乘運(yùn)算，可得關(guān)中地區(qū)1986、2000、2007年土地生態(tài)系統(tǒng)服務(wù)價(jià)值分布圖,再做差運(yùn)算可得土地生態(tài)系統(tǒng)服務(wù)價(jià)值變化圖譜(見(jiàn)圖2)，此后依據(jù)關(guān)中地區(qū)土地利用類型圖、地貌分區(qū)圖、行政區(qū)劃圖對(duì)其進(jìn)行分區(qū)統(tǒng)計(jì)可得相應(yīng)的生態(tài)系統(tǒng)服務(wù)價(jià)值分區(qū)統(tǒng)計(jì)表。

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

3.1 土地生態(tài)系統(tǒng)服務(wù)價(jià)值總量構(gòu)成變化分析

由表1可知：1986—2007年，關(guān)中地區(qū)土地生態(tài)系統(tǒng)服務(wù)價(jià)值總量增加了51.41×108元，年度遞增率為0.46%；其中1986—2000年段和2000—2007年段年度遞增率分別為0.04%、1.30%，后時(shí)段增速約為前時(shí)段的33.65倍。耕地、林地、草地、水域、未利用地的生態(tài)系統(tǒng)服務(wù)價(jià)值量占全區(qū)總價(jià)值量百分比在1986年為27.55%、50.57%、17.59%、4.27%、0.02%，2000年為25.44%、53.20%、17.33%、4.01%、0.02%，2007年為22.96%、56.32%、17.30%、3.40%、0.02%；林地生態(tài)系統(tǒng)服務(wù)價(jià)值量占比超過(guò)全區(qū)總量的一半，而且占比呈增加趨勢(shì)，其它生態(tài)系統(tǒng)服務(wù)價(jià)值量占全區(qū)總生態(tài)價(jià)值量百分比下降，其中占比降幅最大的是耕地。綜上可知：1986—2007年間，關(guān)中地區(qū)林地生態(tài)系統(tǒng)服務(wù)價(jià)值量持續(xù)增加，草地生態(tài)系統(tǒng)服務(wù)價(jià)值量是先減后增、總體增加，其它生態(tài)系統(tǒng)服務(wù)價(jià)值量基本是持續(xù)減少的。研究區(qū)土地生態(tài)系統(tǒng)服務(wù)價(jià)值總量增加主要是林地面積增加、質(zhì)量改善的結(jié)果，這一結(jié)果又直接受惠于“退耕還林”政策的實(shí)施；當(dāng)然，經(jīng)濟(jì)情況改善、人口向城鎮(zhèn)流動(dòng)等因素引起的農(nóng)村燃料結(jié)構(gòu)的改變?yōu)榱植葜脖坏淖匀换謴?fù)創(chuàng)造了條件，促進(jìn)了區(qū)域生態(tài)環(huán)境的良性發(fā)展，增強(qiáng)了區(qū)域土地生態(tài)系統(tǒng)的安全性；不過(guò)，作為本區(qū)糧食基本保障的耕地生態(tài)系統(tǒng)服務(wù)價(jià)值不僅總量在快速減少，而且單位面積量也在減少，其質(zhì)與量均多消損于幾乎不具有生態(tài)服務(wù)價(jià)值的建設(shè)用地，這種情況應(yīng)引起我們的警覺(jué)。

3.2 土地生態(tài)系統(tǒng)服務(wù)價(jià)值變化的區(qū)縣差異

研究期內(nèi)，關(guān)中地區(qū)各區(qū)縣中，單位面積土地生態(tài)系統(tǒng)服務(wù)價(jià)值最大的是太白縣，其次是周至縣，最小的是合陽(yáng)縣和澄城縣。西安市區(qū)、興平市、大荔縣、咸陽(yáng)市區(qū)、臨潼區(qū)、高陵縣、武功縣、禮泉縣、閻良區(qū)、涇陽(yáng)縣、乾縣等11個(gè)區(qū)縣的土地生態(tài)系統(tǒng)服務(wù)價(jià)值量減少，其它區(qū)縣均增加，減速最快的是西安市區(qū)，增速最快的是周至縣；其中1986—2000年，岐山縣、扶風(fēng)縣、寶雞市區(qū)、鳳翔縣、隴縣、眉縣、長(zhǎng)安區(qū)、藍(lán)田縣、戶縣、陳倉(cāng)區(qū)、鳳縣、太白縣、周至縣等13個(gè)區(qū)縣的土地生態(tài)系統(tǒng)服務(wù)價(jià)值量增加，其余區(qū)縣減少，增量最大的是周至縣和太白縣，減量最大的是蒲城縣；2000—2007年，西安市區(qū)、興平市、高陵縣、武功縣等4個(gè)區(qū)縣的土地生態(tài)系統(tǒng)服務(wù)價(jià)值量減少，其它區(qū)縣都增加，減量最大的是西安市區(qū)，增量最大的是周至縣和鳳縣(見(jiàn)圖2)。關(guān)中地區(qū)各區(qū)縣土地生態(tài)系統(tǒng)服務(wù)價(jià)值量及其變化的差異性主要受制于其土地利用結(jié)構(gòu)差異，太白縣、周至縣均位于秦嶺，森林覆蓋率很高，人口稀少，高質(zhì)量的自然生態(tài)系統(tǒng)得以保存，土地生態(tài)系統(tǒng)服務(wù)價(jià)值較大；西安市區(qū)等城市生態(tài)系統(tǒng)由于城市擴(kuò)張，周邊的各類自然和半自然生態(tài)系統(tǒng)逐步消失于基本不具有生態(tài)服務(wù)價(jià)值的建設(shè)用地，土地生態(tài)功能價(jià)值量減少甚至消失都是必然的。

圖1 關(guān)中地區(qū)地貌類型及行政區(qū)劃

圖2 關(guān)中地區(qū)1986—2007年土壤潛力分布及土地生態(tài)服務(wù)價(jià)值變化

Fig.2 Soil potential and land ecosystem services value(ESV) changes in Guanzhong area during 1986—2007

表1 關(guān)中地區(qū)1986—2007年土地生態(tài)服務(wù)價(jià)值構(gòu)成及變化/(108元·a-1)

3.3 土地生態(tài)系統(tǒng)服務(wù)價(jià)值變化的地貌差異

從表2可知：1986—2007年，關(guān)中地區(qū)各地貌單元中，單位面積土地生態(tài)系統(tǒng)服務(wù)價(jià)值量最大的是山地區(qū)，黃土塬區(qū)最?。粏挝幻娣e土地生態(tài)系統(tǒng)服務(wù)價(jià)值量持續(xù)增加的只有山地區(qū)，黃土梁峁區(qū)、黃土臺(tái)塬區(qū)、黃土塬區(qū)是先減后增、總體增加，平原區(qū)則是先減后增、總體減少。在總量構(gòu)成方面，黃土梁峁區(qū)、黃土臺(tái)塬區(qū)、黃土塬區(qū)、平原區(qū)、山地區(qū)的生態(tài)系統(tǒng)服務(wù)價(jià)值量占全區(qū)總價(jià)值量百分比在1986年為3.86%、10.25%、6.21%、15.84%、63.83%，2000年為3.64%、9.63%、5.55%、14.34%、66.84%，2007年為3.89%、9.85%、6.12%、13.59%、66.54%；山地區(qū)是研究區(qū)土地生態(tài)系統(tǒng)服務(wù)價(jià)值量占比最大的區(qū)域，而且優(yōu)勢(shì)在繼續(xù)擴(kuò)大，平原區(qū)的生態(tài)服務(wù)價(jià)值量占比則在快速下降。研究期內(nèi)，黃土梁峁區(qū)、黃土臺(tái)塬區(qū)、黃土塬區(qū)、平原區(qū)、山地區(qū)的生態(tài)系統(tǒng)服務(wù)價(jià)值量年度遞增率為0.50%、0.27%、0.39%、-0.27%、0.66%；其中1986—2000年年度遞增率為-0.39%、-0.40%、-0.76%、-0.67%、0.37%；2000—2007年年度遞增率為2.30%、1.63%、2.73%、0.53%、1.24%；由此可見(jiàn)：研究期內(nèi)，關(guān)中地區(qū)的五大地貌單元中，只有平原區(qū)的土地生態(tài)系統(tǒng)服務(wù)價(jià)值量在減少，其余四類地貌單元的土地生態(tài)系統(tǒng)服務(wù)價(jià)值量均增加，且增量最大、增速最快的都是山地區(qū)。各地貌單元土地生態(tài)服務(wù)價(jià)值分布的差異性由其空間位置及地貌形態(tài)所決定，南部的秦嶺山地因降水豐沛而植被茂密，為土地生態(tài)價(jià)值的高值區(qū)，西部和北部山地雖降水較少，但因人類影響較少，故而亦有較好的林草覆被，土地生態(tài)服務(wù)功能也就相對(duì)較強(qiáng)；中部和北部的平原、黃土塬及黃土臺(tái)塬因地勢(shì)平坦而成為人類活動(dòng)的主要場(chǎng)所，土地生態(tài)系統(tǒng)功能破壞嚴(yán)重，為生態(tài)價(jià)值的低值區(qū)。各地貌單元土地生態(tài)價(jià)值變化情況則受社會(huì)觀念、經(jīng)濟(jì)發(fā)展及政策等人文因子的主導(dǎo)，由于改革開(kāi)放造就了內(nèi)地與沿海、平原與山區(qū)之間的經(jīng)濟(jì)發(fā)展差異，使得人類活動(dòng)逐步向山地外轉(zhuǎn)移，山地的生態(tài)服務(wù)功能得到了恢復(fù)和增強(qiáng)；在疊加強(qiáng)制性的生態(tài)措施后，研究區(qū)除平原區(qū)外，各地貌單元的土地生態(tài)價(jià)值都有較大幅度的提升；但因地貌形態(tài)的差異、生態(tài)保護(hù)政策側(cè)重方向及頒布時(shí)間等方面的影響，使得研究區(qū)土地生態(tài)價(jià)值的變化情況在各地貌單元呈現(xiàn)較明顯的時(shí)空差異性。

表2 關(guān)中地區(qū)1986—2007年各地貌單元土地生態(tài)服務(wù)價(jià)值量

4 結(jié)論與討論

1) 利用格網(wǎng)插值方法將輻射、氣溫、降水、蒸發(fā)等氣候因素和有機(jī)質(zhì)、氮、磷、鉀含量等土壤肥力因素以及高程、坡度等地形地貌因素細(xì)化到每一個(gè)柵格單元，然后采用土壤潛力計(jì)算模型進(jìn)行地圖運(yùn)算，得到關(guān)中地區(qū)土壤潛力柵格圖；利用土地利用圖、生態(tài)系統(tǒng)單位面積服務(wù)價(jià)值當(dāng)量因子表及糧食價(jià)格獲取各地類單位生產(chǎn)力服務(wù)價(jià)值分布圖；上述兩圖作地圖乘運(yùn)算得到土地生態(tài)系統(tǒng)服務(wù)價(jià)值柵格圖；由于土壤潛力柵格圖的每一個(gè)柵格均包含了特定空間位置的自然和人文信息，依此而得的土地生態(tài)系統(tǒng)服務(wù)價(jià)值分布圖自然也具有以柵格為單元的空間差異性；從而在柵格尺度上解決了土地生態(tài)系統(tǒng)服務(wù)價(jià)值空間異質(zhì)性問(wèn)題，避免了研究區(qū)內(nèi)相同面積的同一土地利用類型不論質(zhì)量如何、位于何處均取同樣的生態(tài)服務(wù)價(jià)值量的不合理現(xiàn)象，使土地生態(tài)服務(wù)價(jià)值估算更精細(xì)、更科學(xué)。

2) 關(guān)中地區(qū)土地生態(tài)系統(tǒng)服務(wù)價(jià)值總量在1986—2007年間增加了51.41×108元，增幅約為10.09%，2000—2007年段增速約為1986—2000年段的33.65倍；各地類中：林地生態(tài)系統(tǒng)服務(wù)價(jià)值量持續(xù)增加，草地生態(tài)系統(tǒng)服務(wù)價(jià)值量是先減后增、總體增加，其它生態(tài)系統(tǒng)的服務(wù)價(jià)值量基本是持續(xù)減少的；全境有11個(gè)區(qū)縣的土地生態(tài)系統(tǒng)服務(wù)價(jià)值量減少，其它區(qū)縣增加，減速最快的是西安市區(qū)，增速最快的是周至縣；各地貌單元中：平原土地生態(tài)系統(tǒng)服務(wù)價(jià)值量減少，其余四類地貌單元土地生態(tài)系統(tǒng)服務(wù)價(jià)值量均增加，且山地增長(zhǎng)最快。

3) 研究區(qū)土地生態(tài)系統(tǒng)服務(wù)價(jià)值變化的時(shí)空差異是自然和人文因素共同作用的結(jié)果，地理位置、地表形態(tài)等自然因素主導(dǎo)著土地生態(tài)系統(tǒng)服務(wù)價(jià)值的空間分布，而經(jīng)濟(jì)、政策等人文因子則對(duì)其變化起著主控作用。八十年代后期，改革開(kāi)放逐步取得成效，農(nóng)村外出務(wù)工人員增加，經(jīng)濟(jì)條件較差的山區(qū)有較多的青壯年到城市謀生，同時(shí)務(wù)工收入改善了留居人員的生活條件，煤電油等能源替代了薪柴，使山區(qū)的植被得到恢復(fù)、土地生態(tài)價(jià)值量得以提升；但地勢(shì)平坦的平原、黃土塬及黃土臺(tái)塬則因工業(yè)化、城鎮(zhèn)化的推進(jìn)而大量耗損具有較高生態(tài)價(jià)值的農(nóng)用地，土地生態(tài)價(jià)值總量減損嚴(yán)重；2000年以后，大規(guī)模的退耕還林還草工程極大地改善了研究區(qū)生態(tài)環(huán)境、大幅提升了其土地生態(tài)系統(tǒng)服務(wù)價(jià)值總量，特別是地處退耕還林重點(diǎn)區(qū)段的黃土塬、黃土臺(tái)塬及黃土梁峁區(qū)的土地生態(tài)系統(tǒng)服務(wù)價(jià)值急劇增加，但平原區(qū)在經(jīng)濟(jì)建設(shè)和生態(tài)建設(shè)雙方的博弈中，經(jīng)濟(jì)建設(shè)占用生態(tài)用地依然是主流，平原區(qū)的土地生態(tài)價(jià)值總量繼續(xù)下滑。綜上可知：研究區(qū)在研究期內(nèi)，經(jīng)濟(jì)重心向中東部的平原及黃土臺(tái)塬區(qū)轉(zhuǎn)移，生態(tài)重心則向周圍的山區(qū)及黃土梁峁區(qū)轉(zhuǎn)移，并實(shí)現(xiàn)了經(jīng)濟(jì)總量和生態(tài)總量的雙贏；說(shuō)明科學(xué)調(diào)整用地結(jié)構(gòu)、合理利用土地資源能夠一定程度地在有限的土地資源上實(shí)現(xiàn)經(jīng)濟(jì)與環(huán)境的協(xié)調(diào)發(fā)展。

[1] 張志強(qiáng),徐中民,程國(guó)棟.生態(tài)系統(tǒng)服務(wù)與自然資本價(jià)值評(píng)估[J].生態(tài)學(xué)報(bào),2001,21(11):1919-1926.

[2] Costanza R, d'Arge R, deGrootR, et al. The value of the world’s ecosystem services and nature[J]. Nature, 1997,387:253-260.

[3] 李金昌，等.生態(tài)價(jià)值論[M].重慶:重慶大學(xué)出版社,1999.

[4] 歐陽(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.

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

[6] 肖 寒,歐陽(yáng)志云,趙景柱,等.森林生態(tài)系統(tǒng)服務(wù)功能及其生態(tài)經(jīng)濟(jì)價(jià)值評(píng)估初探——以海南島尖峰嶺熱帶森林為例[J].應(yīng)用生態(tài)學(xué)報(bào),2000,11(4):481-484.

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

[8] 姜文來(lái).森林涵養(yǎng)水源的價(jià)值核算研究[J].水土保持學(xué)報(bào),2003,17(2):34-40．

[9] 傅伯杰，呂一河，高光耀.中國(guó)主要陸地生態(tài)系統(tǒng)服務(wù)與生態(tài)安全研究的重要進(jìn)展[J].自然雜志，2012,34(5):261-272.

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

[11] 魯春霞,謝高地,肖 玉,等.青藏高原生態(tài)系統(tǒng)服務(wù)功能的價(jià)值評(píng)估[J].生態(tài)學(xué)報(bào),2004,24(12):2749-2755.

[12] 謝高地,甄 霖,魯春霞,等.一個(gè)基于專家知識(shí)的生態(tài)系統(tǒng)服務(wù)價(jià)值化方法[J].自然資源學(xué)報(bào),2008,23(5):911-919.

[13] 謝高地，張彩霞，張雷明，等.基于單位面積價(jià)值當(dāng)量因子的生態(tài)系統(tǒng)服務(wù)價(jià)值化方法改進(jìn)[J].自然資源學(xué)報(bào)，2015,30(8):1243-1254.

[14] 封志明.資源科學(xué)導(dǎo)論[M].北京:科學(xué)出版社,2004.

[15] 郭秀銳,毛顯強(qiáng).中國(guó)土地承載力計(jì)算方法研究綜述[J].地球科學(xué)進(jìn)展,2000,15(6):705-711.

[16] 劉紀(jì)遠(yuǎn),徐新良,莊大方,等.20世紀(jì)90年代LUCC過(guò)程對(duì)中國(guó)農(nóng)田光溫生產(chǎn)潛力的影響——基于氣候觀測(cè)與遙感土地利用動(dòng)態(tài)觀測(cè)數(shù)據(jù)[J].中國(guó)科學(xué)D輯地球科學(xué),2005,35(6):483-492.

[17] 楊重一,龐士力,孫彥坤.作物生產(chǎn)潛力研究現(xiàn)狀與趨勢(shì)[J].東北農(nóng)業(yè)大學(xué)學(xué)報(bào),2008,39(7):140-144.

[18] 沈思淵,席承藩.淮北渦河流域農(nóng)業(yè)自然生產(chǎn)潛力模型與分析[J].自然資源學(xué)報(bào),1991,6(1):22-33.

[19] 劉 揚(yáng),賈樹(shù)海,那 波.土地生產(chǎn)潛力計(jì)算方法研究[J].中國(guó)農(nóng)學(xué)通報(bào),2005,21(12):376-381.

[20] 朱志輝,張福春.我國(guó)陸地生態(tài)系統(tǒng)的植物太陽(yáng)能利用率[J].生態(tài)學(xué)報(bào),1985,5(4):343-356.

[21] 康西言,馬輝杰.河北省氣候生產(chǎn)潛力的估算與區(qū)劃[J].中國(guó)農(nóng)業(yè)氣象,2008,29(1):37-41.

[22] 張宏利,陳 豫,胡 偉,等.陜西省戶縣氣候生產(chǎn)潛力分析[J].水土保持研究,2009,16(1):111-114

[23] 葛亞寧，劉 洛，徐新良，等.近50a氣候變化背景下我國(guó)玉米生產(chǎn)潛力時(shí)空演變特征[J].自然資源學(xué)報(bào)，2015,30(5):784-795.

The temporal-spatial dynamic of land ecosystem services value in Guanzhong, Shaanxi

MO Hong-wei, QUAN Bin, YUAN Kai-guo, XIE Jin-ning, XIANG Yun-bo

(CollegeofResourceEnvironmentandSafetyEngineering,HunanUniversityofScienceandTechnology,Xiangtan,Hunan411201,China)

On the basis of Costanza's method and expert interviews with structured questionnaire answered by 700 Chinese ecologists, a new method for assessment of ecosystem services was developed, but spatial heterogeneity issue of land ecosystem services value (ESV) assessment is still an open question. Based on remote sensing images, climate data, soil data and DEM data, according to the theory of geographic information images and supported by RS, GIS, and land-use map, soil potential map and equivalent value per unit area of ecosystem services in China to obtain map of land ESV, the spatial heterogeneity issue of land ESV assessment was solved at the grid scale. This method consists of a major flaw of land ESV assessment relying on equivalent factor. After that, the authors analyzed temporal-spatial changes of land ESV in different scales such as grid units, terrain units, and administrative units of Guanzhong area from 1986 to 2007. The results showed that the total ecosystem service value (ESV) of Guanzhong area increased from 50.96 billion yuan in 1986 to 56.10 billion yuan in 2007, increased by 10.09%. The growth rate of the total ESV in 2000—2007 was 33.65 times higher than that in 1986—2000. Among all counties, there were 11 counties whose total ESV increased. Among five geomorphologic units of Guanzhong, the total ESV of plain reduced and the total ESV of the other four geomorphologic units increased, in which both the largest increment and the fastest growing took place in mountains. Adjusting land use structure in the research period has made obvious ecological benefits, providing a good environment platform for economic construction and social development of Guanzhong area.

Guanzhong area; GIS; land ecosystem services value; spatial heterogeneity

1000-7601(2017)04-0167-06

10.7606/j.issn.1000-7601.2017.04.25

2016-04-20

湖南省教育廳科學(xué)研究項(xiàng)目(14C0467)；地理空間信息技術(shù)國(guó)家地方聯(lián)合工程實(shí)驗(yàn)室開(kāi)放基金項(xiàng)目(2014GISNELJJ004)；國(guó)家自然科學(xué)基金項(xiàng)目(41071057)

莫宏偉(1968—)，男，湖南衡陽(yáng)人，博士，副教授，研究方向?yàn)镚IS與城鄉(xiāng)國(guó)土資源評(píng)價(jià)。 E-mail: mhw819@sina.com。

S181

A