快速城鎮(zhèn)化背景下山地城市生態(tài)系統(tǒng)服務(wù)價(jià)值變化梯度特征
——以貴陽市為例

韓會(huì)慶,劉 悅,蔡廣鵬,白玉梅,馬淑亮,陳思盈，羅瑞堯

(1.貴州理工學(xué)院 建筑與城市規(guī)劃學(xué)院,貴陽 550003;2.貴州師范大學(xué) 地理與環(huán)境科學(xué)學(xué)院,貴陽 550001)

城市作為人類活動(dòng)最為強(qiáng)烈地區(qū)，其生態(tài)系統(tǒng)服務(wù)功能深受人為活動(dòng)干擾[1]。城市生態(tài)系統(tǒng)服務(wù)功能的高低對(duì)城市可持續(xù)發(fā)展至關(guān)重要[2]。然而，在快速城鎮(zhèn)化背景下，通過改變地表景觀格局、生態(tài)過程、生物棲息地、生物地球化學(xué)循環(huán)等，影響城市生態(tài)系統(tǒng)結(jié)構(gòu)和組分，從而深刻影響城市生態(tài)系統(tǒng)服務(wù)供給水平[3]。因此，城鎮(zhèn)化影響下生態(tài)系統(tǒng)服務(wù)研究成為國(guó)內(nèi)外研究的焦點(diǎn)。

目前，城鎮(zhèn)化與生態(tài)系統(tǒng)服務(wù)的關(guān)系研究主要集中于城鎮(zhèn)化對(duì)生態(tài)系統(tǒng)服務(wù)的影響以及城鎮(zhèn)化與生態(tài)系統(tǒng)服務(wù)之間相關(guān)性研究等方面，其中城鎮(zhèn)化對(duì)生態(tài)系統(tǒng)服務(wù)的影響研究涉及城市生態(tài)系統(tǒng)[4]、水生態(tài)系統(tǒng)[5]、土壤生態(tài)系統(tǒng)[6]、農(nóng)林生態(tài)系統(tǒng)[7]、文化生態(tài)系統(tǒng)[8]、生物多樣性[9]等。研究區(qū)域包括城市圈[10]、海岸帶[11]、特定空間尺度(如市域、地區(qū))[12-13]等。生態(tài)系統(tǒng)服務(wù)之間權(quán)衡與協(xié)同關(guān)系對(duì)城鎮(zhèn)化的響應(yīng)亦受到關(guān)注[14]。研究時(shí)間既有過去時(shí)點(diǎn)，也涉及未來預(yù)估時(shí)點(diǎn)，如孟悅[15]分析了1991—2014年寧波市城鎮(zhèn)化對(duì)生態(tài)系統(tǒng)服務(wù)的影響，夏敏等[16]通過預(yù)估宜興市官林鎮(zhèn)2021年和2024年生態(tài)用地，分析了未來生態(tài)用地變化趨勢(shì)對(duì)生態(tài)系統(tǒng)服務(wù)的影響。城鎮(zhèn)化與生態(tài)系統(tǒng)服務(wù)之間相關(guān)性研究主要關(guān)注城鎮(zhèn)化與生態(tài)系統(tǒng)服務(wù)的空間相關(guān)性[17]、城鎮(zhèn)化對(duì)生態(tài)系統(tǒng)服務(wù)的作用機(jī)制[18]等。研究方法有地理加權(quán)回歸[18]、空間統(tǒng)計(jì)法(如莫蘭指數(shù)、空間自相關(guān))[17,19]、耦合協(xié)調(diào)度模型[20]等。然而，這些研究較多關(guān)注我國(guó)東部和中部社會(huì)經(jīng)濟(jì)發(fā)達(dá)，且地形平坦的城市或地區(qū)，而對(duì)西部經(jīng)濟(jì)相對(duì)落后的山地城市關(guān)注較少。

貴陽市是貴州省省會(huì)，素有“林城”美譽(yù)，是我國(guó)西南部重要的中心城市之一，同時(shí)貴陽市位于云貴高原東部，屬喀斯特丘原山地地貌，是典型的山地城市[21]。2000年以來，在高速經(jīng)濟(jì)發(fā)展推動(dòng)下，快速城鎮(zhèn)化使得貴陽市土地利用發(fā)生劇烈變化[22]，進(jìn)而影響該市生態(tài)系統(tǒng)服務(wù)供給能力。然而，在快速城鎮(zhèn)化影響下山地城市生態(tài)系統(tǒng)服務(wù)變化規(guī)律尚不清楚，這成為當(dāng)前亟需解決的科學(xué)問題。鑒于此，以貴陽市為例，分析快速城鎮(zhèn)化背景下山地城市生態(tài)系統(tǒng)服務(wù)價(jià)值變化的梯度特征，以期為山地城市土地資源合理利用，生態(tài)環(huán)境保護(hù)提供科學(xué)依據(jù)。

1 數(shù)據(jù)來源與研究方法

1.1 數(shù)據(jù)來源

本研究數(shù)據(jù)包括土地利用、地形、社會(huì)經(jīng)濟(jì)3種數(shù)據(jù)。土地利用數(shù)據(jù)是利用1995年、2005年、2015年Landsat遙感影像，采用人機(jī)交互方式解譯獲得。為驗(yàn)證解譯精度，以貴陽市翁朵村為典型案例地，基于0.5 m分辨率的Pleiades高分衛(wèi)星解譯得到2015年翁朵村土地利用數(shù)據(jù)(采用人工目視解譯，經(jīng)實(shí)地驗(yàn)證，精度大于93%)，與基于Landsat遙感影像解譯的土地利用數(shù)據(jù)進(jìn)行比較，發(fā)現(xiàn)基于Landsat的土地利用數(shù)據(jù)精度較高(大于83%)，這說明解譯的土地利用數(shù)據(jù)可以滿足研究需要。依據(jù)研究區(qū)特點(diǎn)，將土地利用類型劃分為水田、旱地、有林地、灌木林、草地、水域、建設(shè)用地、未利用地共8種地類;地形數(shù)據(jù)有30 m DEM和坡度數(shù)據(jù)，均下載于國(guó)家地球系統(tǒng)科學(xué)數(shù)據(jù)中心;社會(huì)經(jīng)濟(jì)數(shù)據(jù)包括來源于歷年貴陽市統(tǒng)計(jì)年鑒的人口密度數(shù)據(jù)和糧食產(chǎn)量數(shù)據(jù)、來源于矢量化的貴陽交通圖和行政區(qū)劃圖。

1.2 研究方法

(1) 生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估方法。生態(tài)系統(tǒng)服務(wù)價(jià)值計(jì)算是依據(jù)謝高地等[23]提出的中國(guó)生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)法。由于該方法是適用于全國(guó)尺度，在貴陽市生態(tài)系統(tǒng)服務(wù)價(jià)值計(jì)算過程中需進(jìn)行修正。具體修正步驟:貴陽市生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)確定是基于2010—2015年貴陽市單位面積糧食產(chǎn)值。貴陽市主要糧食作物為玉米和水稻，種植結(jié)構(gòu)約為6∶4,2010—2015年玉米和水稻平均單產(chǎn)分別是4 200 kg/hm2,5 550 kg/hm2，收購(gòu)價(jià)格1.9元/kg和2.2元/kg，根據(jù)種植比例，計(jì)算得到單位面積糧食社會(huì)經(jīng)濟(jì)產(chǎn)值9 672元/hm2，單位面積農(nóng)田生態(tài)系統(tǒng)服務(wù)價(jià)值等于其社會(huì)經(jīng)濟(jì)產(chǎn)值1/7，計(jì)算得到1個(gè)當(dāng)量生態(tài)系統(tǒng)服務(wù)價(jià)值為1 382元。在系數(shù)修正基礎(chǔ)上，利用各生態(tài)系統(tǒng)服務(wù)的當(dāng)量表[23]以及1個(gè)當(dāng)量生態(tài)系統(tǒng)服務(wù)價(jià)值計(jì)算得到貴陽市生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)(表1)。

表1 貴陽市單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù) 元/hm2

(2) 生態(tài)系統(tǒng)服務(wù)價(jià)值空間異質(zhì)性和空間變化分析方法。利用ArcGIS軟件的Getis-Ord Gi工具分析生態(tài)系統(tǒng)服務(wù)價(jià)值的熱點(diǎn)和冷點(diǎn)。依據(jù)其計(jì)算的Z值和p值將熱點(diǎn)和冷點(diǎn)分為高熱點(diǎn)(Z≥2,p>90%)、低熱點(diǎn)(2>Z≥1,p>90%)，高冷點(diǎn)(Z<-2,p>90%)，低冷點(diǎn)(-1>Z≥-2,p>90%)和不顯著(1>Z≥-1,p<90%)5種類型。其計(jì)算過程:基于三期單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值空間圖，分別轉(zhuǎn)為生態(tài)系統(tǒng)服務(wù)價(jià)值Point文件，進(jìn)而利用Kriging插值得到生態(tài)系統(tǒng)服務(wù)價(jià)值的熱點(diǎn)和冷點(diǎn)空間圖。生態(tài)系統(tǒng)服務(wù)價(jià)值空間變化分析方法是利用ArcGIS軟件的柵格處理器，進(jìn)行1995—2005年和2005—2015年生態(tài)系統(tǒng)服務(wù)價(jià)值疊加處理，進(jìn)而得到10 a間生態(tài)系統(tǒng)服務(wù)價(jià)值變化的空間格局，負(fù)值表示減少區(qū)，正值表示增加區(qū)。

(3) 地形及社會(huì)經(jīng)濟(jì)梯度分級(jí)。地形起伏度梯度劃分是利用柵格鄰域計(jì)算工具，在一定窗口下通過計(jì)算最大高程與最低高程差值得到地形起伏度，并利用ArcGIS軟件的自然斷點(diǎn)法將其分為6個(gè)梯度帶;地形位指數(shù)梯度劃分是在計(jì)算任意坡度與平均坡度以及任意高程與平均高程之間關(guān)系，進(jìn)而得到地形位指數(shù)，并利用自然斷點(diǎn)法其分為6個(gè)梯度帶;交通梯度劃分是以貴陽市交通網(wǎng)為基礎(chǔ)，從道路中心線往兩側(cè)生成不同距離的緩沖區(qū)，并依據(jù)距離道路遠(yuǎn)近分為6個(gè)交通梯度帶;人口密度梯度劃分是以2015年貴陽市各區(qū)縣人口密度圖，將其分為不同人口密度等級(jí)，進(jìn)而得到6個(gè)人口密度梯度帶;城鄉(xiāng)梯度劃分是以各縣區(qū)行政中心為基礎(chǔ)，往周邊生成不同距離的緩沖區(qū)，并依據(jù)距離行政中心遠(yuǎn)近分為6個(gè)城鄉(xiāng)梯度帶(表2)。緩沖區(qū)的生成是利用ArcGIS軟件多環(huán)緩沖工具。

表2 貴陽市地形與社會(huì)經(jīng)濟(jì)梯度分級(jí)

2 結(jié)果與分析

2.1 貴陽市生態(tài)系統(tǒng)服務(wù)價(jià)值變化整體特征

1995—2005年貴陽市供給服務(wù)及總服務(wù)價(jià)值呈上升趨勢(shì)，調(diào)節(jié)服務(wù)、支持服務(wù)及文化服務(wù)價(jià)值均呈下降趨勢(shì)，其中文化服務(wù)價(jià)值下降幅度較大。2005—2015年貴陽市各類型服務(wù)和總服務(wù)價(jià)值均呈上升趨勢(shì)，其中供給服務(wù)價(jià)值上升幅度較大(表3)。

表3 貴陽市生態(tài)系統(tǒng)服務(wù)價(jià)值的整體變化

1995年貴陽市生態(tài)系統(tǒng)服務(wù)價(jià)值冷點(diǎn)集中分布于南部地區(qū)和西北部，其他地區(qū)零星分布，且以低冷點(diǎn)為主，中南部為高冷點(diǎn)主要分布區(qū)。熱點(diǎn)主要分布于西北部邊緣地區(qū)、西南部以及東北部的零星地區(qū);2005年和2015年貴陽市生態(tài)系統(tǒng)服務(wù)價(jià)值冷熱點(diǎn)與1995年較為接近，其中2005年的差異之處主要體現(xiàn)在西部的熱點(diǎn)區(qū)、2015年的差異區(qū)域主要有中南部的高冷點(diǎn)、東部的低熱點(diǎn)(圖1)。

圖1 1995年、2005年、2015年貴陽市生態(tài)系統(tǒng)服務(wù)價(jià)值熱點(diǎn)分析

1995—2005年貴陽市大部分地區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值無變化。中西部、北部的零星地區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值增加，中部和北部的小部分地區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值減少。2005—2015年貴陽市中南部地區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值變化較為明顯，主要表現(xiàn)為生態(tài)系統(tǒng)服務(wù)價(jià)值減少，西南部和中部的小部分地區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值增加較為明顯(圖2)。

圖2 1995-2015年貴陽市生態(tài)系統(tǒng)服務(wù)價(jià)值變化空間格局

2.2 貴陽市生態(tài)系統(tǒng)服務(wù)價(jià)值變化的地形梯度特征

隨著地形起伏度的增加，1995年、2005年、2015年調(diào)節(jié)服務(wù)價(jià)值呈現(xiàn)先下降后增加趨勢(shì)，供給服務(wù)、支持服務(wù)、文化服務(wù)和總服務(wù)價(jià)值呈增加趨勢(shì)，其中供給服務(wù)、支持服務(wù)、文化服務(wù)的變化特點(diǎn)更加接近，總服務(wù)價(jià)值呈線性增加趨勢(shì)(圖3)。

圖3 貴陽市生態(tài)系統(tǒng)服務(wù)價(jià)值在地形起伏度梯度上特征

1995—2005年貴陽市供給服務(wù)價(jià)值在地形起伏度Ⅰ至Ⅵ梯度均呈上升趨勢(shì)，其中在Ⅵ梯度的上升幅度較大，其他服務(wù)及總服務(wù)價(jià)值在地形起伏度Ⅰ至Ⅴ梯度均呈下降趨勢(shì)，在Ⅵ梯度均呈上升趨勢(shì)，其中均在Ⅰ梯度的下降幅度較大。2005—2015年貴陽市供給服務(wù)價(jià)值在地形起伏度Ⅰ至Ⅵ梯度均呈上升趨勢(shì)，其中在Ⅰ梯度的上升幅度較大;其他服務(wù)及總服務(wù)價(jià)值在地形起伏度Ⅰ梯度均呈下降趨勢(shì)，在Ⅱ至Ⅵ梯度均呈上升趨勢(shì)，其中均在Ⅲ梯度的上升幅度較大(表4)。

表4 貴陽市生態(tài)系統(tǒng)服務(wù)價(jià)值變化在地形起伏度梯度上的變化率 %

隨著地形位指數(shù)梯度的增加，1995年、2005年、2015年貴陽市供給服務(wù)、支持服務(wù)及文化服務(wù)價(jià)值均持續(xù)上升，而調(diào)節(jié)服務(wù)及總服務(wù)價(jià)值先下降后上升，Ⅱ梯度是最低點(diǎn)(圖4)。

圖4 貴陽市生態(tài)系統(tǒng)服務(wù)價(jià)值在地形位指數(shù)梯度上特征

1995—2005年貴陽市供給服務(wù)價(jià)值在地形位指數(shù)Ⅰ至Ⅵ梯度均呈上升趨勢(shì)，其中在Ⅰ梯度的上升幅度較大;其他服務(wù)及總服務(wù)價(jià)值在地形位指數(shù)Ⅰ至Ⅴ梯度均呈下降趨勢(shì)，在Ⅵ梯度均呈上升趨勢(shì)，其中調(diào)節(jié)服務(wù)及總服務(wù)價(jià)值在Ⅱ梯度的下降幅度較大，支持服務(wù)價(jià)值在Ⅰ梯度的下降幅度較大，文化服務(wù)價(jià)值在Ⅰ梯度的下降幅度較大。2005—2015年貴陽市供給服務(wù)價(jià)值在地形位指數(shù)Ⅰ至Ⅵ梯度均呈上升趨勢(shì)，其中在Ⅰ梯度的上升幅度較大;其他服務(wù)及總服務(wù)價(jià)值在地形位指數(shù)Ⅰ梯度均呈下降趨勢(shì)，在Ⅱ至Ⅵ梯度均呈上升趨勢(shì)，其中均在Ⅵ梯度的上升幅度較大(表5)。

表5 貴陽市生態(tài)系統(tǒng)服務(wù)價(jià)值變化在地形位指數(shù)梯度上的變化率 %

2.3 貴陽市生態(tài)系統(tǒng)服務(wù)價(jià)值變化的社會(huì)經(jīng)濟(jì)梯度特征

隨著交通梯度的增加，1995年、2005年、2015年貴陽市各類型服務(wù)和總服務(wù)價(jià)值均持續(xù)上升，且各類型的服務(wù)價(jià)值變化特征也較為接近(圖5)。

圖5 貴陽市生態(tài)系統(tǒng)服務(wù)價(jià)值在交通、人口密度、城鄉(xiāng)梯度上特征

1995—2005年貴陽市供給服務(wù)價(jià)值在交通Ⅰ至Ⅵ梯度均呈上升趨勢(shì)，其中在Ⅵ梯度的上升幅度較大;其他服務(wù)及總服務(wù)價(jià)值在交通Ⅰ至Ⅴ梯度均呈下降趨勢(shì)，在Ⅵ梯度均呈上升趨勢(shì)，其中均在Ⅴ梯度的下降幅度較大。2005—2015年貴陽市供給服務(wù)價(jià)值在交通Ⅰ至Ⅵ梯度均呈上升趨勢(shì)，其中在Ⅵ梯度的上升幅度較大;其他服務(wù)及總服務(wù)價(jià)值在交通Ⅰ梯度均呈下降趨勢(shì)，在Ⅱ至Ⅵ梯度均呈上升趨勢(shì)，其中均在Ⅵ梯度的上升幅度較大(表6)。

表6 貴陽市生態(tài)系統(tǒng)服務(wù)價(jià)值在交通梯度上的變化率 %

隨著人口密度梯度增加，1995年、2005年、2015年貴陽市各類型服務(wù)和總服務(wù)價(jià)值均呈現(xiàn)波動(dòng)起伏特點(diǎn)，供給服務(wù)價(jià)值的高值位于Ⅴ和Ⅵ梯度，低值位于Ⅱ和Ⅳ梯度，調(diào)節(jié)服務(wù)價(jià)值的高值位于Ⅲ和Ⅴ梯度，低值位于Ⅰ和Ⅱ梯度，支持服務(wù)價(jià)值的高值位于Ⅴ和Ⅵ梯度，低值位于Ⅱ梯度，文化服務(wù)價(jià)值的高值位于Ⅴ和Ⅵ梯度，低值位于Ⅱ梯度，總服務(wù)價(jià)值的高值位于Ⅲ和Ⅴ梯度，低值位于Ⅰ和Ⅱ梯度(圖6)。

1995—2005年貴陽市供給服務(wù)價(jià)值在人口密度Ⅰ,Ⅲ至Ⅵ梯度呈上升趨勢(shì)，在Ⅱ梯度呈下降趨勢(shì)，其中在Ⅵ梯度的上升幅度較大;其他服務(wù)及總服務(wù)價(jià)值在人口密度Ⅰ至Ⅲ,Ⅴ,Ⅵ梯度均呈下降趨勢(shì)，在Ⅳ梯度均呈上升趨勢(shì)，其中調(diào)節(jié)服務(wù)價(jià)值在Ⅱ梯度的下降幅度較大，而支持服務(wù)、文化服務(wù)及總服務(wù)價(jià)值在Ⅰ梯度的下降幅度較大;2005—2015年貴陽市供給服務(wù)價(jià)值在人口密度Ⅱ至Ⅵ梯度呈上升趨勢(shì)，在Ⅰ梯度呈下降趨勢(shì)，其中在Ⅲ梯度的上升幅度較大;調(diào)節(jié)服務(wù)價(jià)值在人口密度Ⅰ,Ⅴ,Ⅵ梯度呈上升趨勢(shì)，在Ⅱ至Ⅳ梯度呈下降趨勢(shì);支持服務(wù)、文化服務(wù)及總服務(wù)價(jià)值在人口密度Ⅰ,Ⅱ,Ⅳ梯度均呈下降趨勢(shì)，在Ⅲ,Ⅴ,Ⅵ梯度均呈上升趨勢(shì)，其中均在Ⅱ梯度的下降幅度較大，在Ⅴ梯度的上升幅度較大(表7)。

表7 貴陽市生態(tài)系統(tǒng)服務(wù)價(jià)值在人口密度梯度上的變化率 %

隨著城鄉(xiāng)梯度的增加，1995年、2005年、2015年貴陽市供給服務(wù)價(jià)值持續(xù)上升，調(diào)節(jié)服務(wù)及總服務(wù)價(jià)值呈先上升后下降再上升最后下降趨勢(shì)，其中梯度Ⅱ和Ⅴ的服務(wù)價(jià)值較高，支持服務(wù)及文化服務(wù)價(jià)值呈先上升后下降趨勢(shì)，其中梯度Ⅳ和Ⅴ的服務(wù)價(jià)值較高(圖7)。

1995—2005年貴陽市供給服務(wù)價(jià)值在城鄉(xiāng)Ⅰ至Ⅵ梯度均呈上升趨勢(shì)，其中在Ⅰ梯度的上升幅度較大;其他服務(wù)及總服務(wù)價(jià)值在城鄉(xiāng)Ⅰ至Ⅴ梯度均呈下降趨勢(shì)，在Ⅵ梯度均呈上升趨勢(shì)，其中調(diào)節(jié)服務(wù)及總服務(wù)價(jià)值均在Ⅴ梯度的下降幅度較大，而支持服務(wù)價(jià)值在Ⅰ梯度的下降幅度較大，文化服務(wù)價(jià)值在Ⅳ梯度的下降幅度較大。2005—2015年貴陽市供給服務(wù)價(jià)值在城鄉(xiāng)Ⅰ至Ⅵ梯度均呈上升趨勢(shì)，其中在Ⅰ梯度的上升幅度較大;其他服務(wù)價(jià)值及總服務(wù)價(jià)值在城鄉(xiāng)Ⅰ梯度均呈下降趨勢(shì)，在Ⅱ至Ⅵ梯度均呈上升趨勢(shì)，其中支持服務(wù)、文化服務(wù)及總服務(wù)價(jià)值在均在Ⅴ梯度的上升幅度較大，而調(diào)節(jié)服務(wù)價(jià)值在Ⅵ梯度的上升幅度較大(表8)。

表8 貴陽市生態(tài)系統(tǒng)服務(wù)價(jià)值在城鄉(xiāng)梯度上的變化率 %

3 討論與結(jié)論

3.1 討 論

在山地地形約束下，貴陽城鎮(zhèn)化過程主要沿地形平坦的平壩地?cái)U(kuò)張[24]，這使得城鎮(zhèn)化集中于地形起伏度和地形位指數(shù)較小的梯度帶，進(jìn)而導(dǎo)致這些區(qū)域生態(tài)系統(tǒng)服務(wù)價(jià)值下降。同時(shí)由于山地地形的約束，貴陽市交通、行政中心往往集中于低緩平壩地，隨著城鎮(zhèn)化擴(kuò)張，這些距離交通線和行政中心近的梯度帶生態(tài)系統(tǒng)服務(wù)價(jià)值下降。此外，受地形限制，地形起伏和地形位指數(shù)較大的梯度、遠(yuǎn)離交通和行政中心的梯度帶主要以山地為主，這些區(qū)域?qū)俚湫痛嗳醯目λ固厣降厣鷳B(tài)系統(tǒng)，近年來貴陽市實(shí)施的退耕還林還草、天然林保護(hù)工程等[25]使得這些山地為主的梯度帶生態(tài)系統(tǒng)服務(wù)價(jià)值呈增加趨勢(shì)。從空間格局看，貴陽市中南部為盆壩地，周邊為山地。在快速城鎮(zhèn)化影響下，貴陽市中南部地形平坦的盆壩地(亦是中心城區(qū))人類活動(dòng)強(qiáng)烈，因此中南部地區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值有所降低，這在高冷點(diǎn)的分布變化有所體現(xiàn)，而周邊山地在一系列生態(tài)保護(hù)措施下，其生態(tài)系統(tǒng)服務(wù)價(jià)值有所提高。值得注意的是，在山地地形影響下，為了疏散中心城區(qū)人口壓力，貴陽市城鎮(zhèn)化并不集中于人口稠密區(qū)，往往選取有一定人口基礎(chǔ)，地形相對(duì)平坦的地區(qū)[26]，這使得人口密度對(duì)城鎮(zhèn)化的影響較小，故而在生態(tài)系統(tǒng)服務(wù)價(jià)值變化在人口密度梯度上呈現(xiàn)波動(dòng)起伏特點(diǎn)。

貴陽市生態(tài)系統(tǒng)服務(wù)價(jià)值變化的地形梯度特征與陳奕竹等[27]對(duì)湘西地區(qū)的研究結(jié)果較為相似，這說明相似兩個(gè)地區(qū)的人類活動(dòng)在山地上表現(xiàn)具有一致性。然而，由于人類對(duì)不同地域山地的利用存在差異，進(jìn)而對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響不同，這在本研究與李理等[28]對(duì)華北地區(qū)淇河研究結(jié)果的差異性有所體現(xiàn)。

3.2 結(jié) 論

貴陽市1995—2005年和2005—2015年供給服務(wù)和總服務(wù)價(jià)值增加，調(diào)節(jié)服務(wù)、支持服務(wù)及文化服務(wù)價(jià)值在1995—2005年呈下降趨勢(shì)，在2005—2015年呈增加趨勢(shì);在地形影響下山地城市生態(tài)系統(tǒng)服務(wù)價(jià)值變化空間格局破碎化突出，地形相對(duì)平坦區(qū)域是生態(tài)系統(tǒng)服務(wù)價(jià)值變化的主要區(qū)域。隨著地形起伏度、地形位指數(shù)、交通梯度、城鄉(xiāng)梯度的增加，供給服務(wù)呈增加趨勢(shì)。除高梯度外，調(diào)節(jié)服務(wù)、支持服務(wù)、文化服務(wù)和總服務(wù)價(jià)值在1995—2005年地形起伏度、地形位指數(shù)、交通、城鄉(xiāng)的中低梯度上呈下降趨勢(shì)。除低梯度外，在中高梯度上調(diào)節(jié)服務(wù)、支持服務(wù)、文化服務(wù)和總服務(wù)價(jià)值2005—2015年呈上升趨勢(shì)。各種服務(wù)及總服務(wù)價(jià)值在人口密度梯度上呈現(xiàn)下降—上升—下降的波動(dòng)特點(diǎn)。