中國(guó)土地利用與生態(tài)服務(wù)價(jià)值空間動(dòng)態(tài)研究
——以地級(jí)及以上城市為例

危小建, 趙 莉, 程朋根, 謝亞娟, 王慧敏

(1.東華理工大學(xué) 測(cè)繪工程學(xué)院, 南昌 330013; 2.自然資源部 環(huán)鄱陽(yáng)湖區(qū)域礦山環(huán)境監(jiān)測(cè)與治理重點(diǎn)實(shí)驗(yàn)室, 南昌 330013; 3.仲愷農(nóng)業(yè)工程學(xué)院 城鄉(xiāng)建設(shè)學(xué)院, 廣州 510225)

土地利用是人與自然進(jìn)行緊密交流的重要紐帶，土地利用/土地覆被變化影響著生態(tài)系統(tǒng)的結(jié)構(gòu)和功能[1]，帶來(lái)了生態(tài)系統(tǒng)服務(wù)價(jià)值的變化。生態(tài)系統(tǒng)服務(wù)(Ecosystem Services，ES)是指生態(tài)系統(tǒng)與生態(tài)過(guò)程直接或間接產(chǎn)生的，保證人類(lèi)賴(lài)以生存的不可缺少的生命支撐產(chǎn)品和服務(wù)[2]，它是人類(lèi)生存與現(xiàn)代文明發(fā)展的重要保障，同時(shí)也是科學(xué)技術(shù)無(wú)法取代的產(chǎn)物[3]。土地利用與生態(tài)系統(tǒng)相互作用、密不可分，分析土地利用與生態(tài)服務(wù)價(jià)值兩者的相互關(guān)系，對(duì)合理規(guī)劃土地資源，保證生態(tài)系統(tǒng)的協(xié)調(diào)發(fā)展具有重要意義。

近年來(lái)，國(guó)內(nèi)外不少學(xué)者從多方面對(duì)土地利用與生態(tài)服務(wù)價(jià)值的關(guān)系開(kāi)展了一系列研究[4-9]，已然成為生態(tài)學(xué)、地理學(xué)等領(lǐng)域的研究熱點(diǎn)之一[10-11]。從研究驅(qū)動(dòng)因素來(lái)看，茆長(zhǎng)寶[12]、韓蕊[13]等基于土地利用數(shù)據(jù)和相關(guān)社會(huì)經(jīng)濟(jì)統(tǒng)計(jì)數(shù)據(jù)分別深入探究了促使江蘇省和川東地區(qū)土地利用及生態(tài)服務(wù)價(jià)值產(chǎn)生變化的驅(qū)動(dòng)因素；從研究方法和手段來(lái)看，張鵬巖等[14]采用熱點(diǎn)分析、重心模型以及地理探測(cè)器模型研究了黃河流域生態(tài)服務(wù)價(jià)值的空間分異特征；姚小薇等[15]應(yīng)用空間自相關(guān)理論分析了人口、空間、經(jīng)濟(jì)、生活4個(gè)不同層面城鎮(zhèn)化水平對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響。從研究尺度來(lái)看，目前開(kāi)展的研究中主要是以單個(gè)省市縣[16-18]、流域[19]或生態(tài)經(jīng)濟(jì)區(qū)[20]為研究對(duì)象，而針對(duì)中觀尺度研究鮮少，以地級(jí)城市為研究單元也更為少見(jiàn)。同時(shí)，研究不同地級(jí)城市土地利用、社會(huì)經(jīng)濟(jì)發(fā)展與ESV的空間相關(guān)性，能夠進(jìn)一步揭示人類(lèi)活動(dòng)與社會(huì)經(jīng)濟(jì)發(fā)展對(duì)不同地區(qū)生態(tài)環(huán)境的影響，而當(dāng)前從土地利用與社會(huì)經(jīng)濟(jì)兩方面共同出發(fā)，對(duì)二者與生態(tài)服務(wù)價(jià)值的空間相關(guān)性方面也少有研究。因此，本文基于1990—2018年中國(guó)地級(jí)及以上城市土地利用數(shù)據(jù)和社會(huì)經(jīng)濟(jì)數(shù)據(jù)，分析土地利用與生態(tài)系統(tǒng)服務(wù)價(jià)值的動(dòng)態(tài)變化情況，同時(shí)對(duì)生態(tài)服務(wù)價(jià)值進(jìn)行敏感性分析，最后再通過(guò)研究生態(tài)服務(wù)價(jià)值與土地利用動(dòng)態(tài)度、社會(huì)經(jīng)濟(jì)指標(biāo)綜合值的空間相關(guān)性，為優(yōu)化土地利用結(jié)構(gòu)、保證中國(guó)土地資源可持續(xù)利用以及促進(jìn)生態(tài)文明社會(huì)協(xié)調(diào)發(fā)展提供一定的決策支持。

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

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

本研究采用的中國(guó)地級(jí)及以上城市1990年、2000年、2010年、2018年的4期空間分辨率為1 km×1 km的土地利用數(shù)據(jù)由資源環(huán)境科學(xué)數(shù)據(jù)中心(http:∥www.resdc.cn)整理而得，土地利用類(lèi)型主要通過(guò)耕地、林地、草地、水域、建設(shè)用地和未利用土地6種一級(jí)地類(lèi)的形式呈現(xiàn)。相應(yīng)年份的社會(huì)經(jīng)濟(jì)數(shù)據(jù)主要來(lái)源于中國(guó)統(tǒng)計(jì)年鑒和中國(guó)城市統(tǒng)計(jì)年鑒。本文研究對(duì)象是中國(guó)地級(jí)及以上城市單元，包括除香港、澳門(mén)和臺(tái)灣地區(qū)以外的地級(jí)及以上城市和省直轄縣級(jí)行政單元。西藏自治區(qū)無(wú)各個(gè)地級(jí)市GDP、城鎮(zhèn)及農(nóng)村人口指標(biāo)等詳細(xì)數(shù)據(jù)，本文以整個(gè)西藏地區(qū)的各項(xiàng)社會(huì)經(jīng)濟(jì)指標(biāo)的平均值替代。另外，部分研究單元如鎮(zhèn)江市、淮陰市及宿遷市等社會(huì)經(jīng)濟(jì)數(shù)據(jù)缺失，新疆維吾爾自治區(qū)無(wú)1990年城市統(tǒng)計(jì)年鑒等情況，使得部分?jǐn)?shù)據(jù)難以獲取，本文以相鄰年份統(tǒng)計(jì)年鑒或同省周邊相鄰市縣的數(shù)據(jù)插值替代。

1.2 研究方法

1.2.1 土地利用動(dòng)態(tài)研究方法 土地利用動(dòng)態(tài)度是分析土地利用結(jié)構(gòu)變化動(dòng)態(tài)的重要指標(biāo)，主要包括單一土地利用動(dòng)態(tài)度和綜合土地利用動(dòng)態(tài)度[14,21]。

(1) 單一土地利用動(dòng)態(tài)度。單一土地利用動(dòng)態(tài)度(K)是指在一定時(shí)間范圍內(nèi)某種土地利用類(lèi)型的數(shù)量變化情況。其公式為：

(1)

式中:Ua，Ub分別為研究初期、末期某種土地利用類(lèi)型的面積(km2)；T為研究時(shí)段。

(2) 綜合土地利用動(dòng)態(tài)度。綜合土地利用動(dòng)態(tài)度(LC)是指在某一時(shí)段所在研究區(qū)內(nèi)，結(jié)合所有土地利用類(lèi)型的變化情況而計(jì)算的綜合土地利用變化的年變化率。其公式為：

(2)

式中：Ui為研究初期第i類(lèi)土地利用類(lèi)型的面積；ΔUi-j為研究時(shí)段第i類(lèi)土地利用類(lèi)型向其他非i類(lèi)土地利用類(lèi)型轉(zhuǎn)化的面積的絕對(duì)值；n為土地利用類(lèi)型的種類(lèi)。

1.2.2 生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估模型構(gòu)建 Costanza等[22]在效用價(jià)值理論和均衡價(jià)值理論的基礎(chǔ)上，首先對(duì)全球生態(tài)系統(tǒng)服務(wù)價(jià)值進(jìn)行了評(píng)估，隨后謝高地等[23-25]通過(guò)分析Costanza等[22]得出的全球生態(tài)系統(tǒng)服務(wù)價(jià)值，結(jié)合中國(guó)生態(tài)系統(tǒng)結(jié)構(gòu)特點(diǎn)，進(jìn)一步制定了中國(guó)生態(tài)系統(tǒng)單位面積的生態(tài)服務(wù)價(jià)值當(dāng)量。本文借助謝高地等[24]的“中國(guó)陸地生態(tài)系統(tǒng)單位面積生態(tài)服務(wù)價(jià)值當(dāng)量表”測(cè)算中國(guó)地級(jí)及以上城市的生態(tài)服務(wù)價(jià)值。其中，生態(tài)服務(wù)價(jià)值當(dāng)量因子的經(jīng)濟(jì)價(jià)值約為全國(guó)平均糧食單產(chǎn)價(jià)值的1/7。參照前人研究[6,15,22]，建設(shè)用地的生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估過(guò)程中需要考慮因素眾多，測(cè)算較為復(fù)雜，本研究對(duì)其將不作估算。未利用地的生態(tài)服務(wù)價(jià)值估算時(shí)，其當(dāng)量因子參數(shù)以草地與荒漠當(dāng)量因子的平均值為基礎(chǔ)進(jìn)行估算[26]。因此，根據(jù)中國(guó)各省、自治區(qū)、直轄市不同糧食單產(chǎn)的實(shí)際情況計(jì)算當(dāng)?shù)貑挝幻娣e糧食產(chǎn)量的經(jīng)濟(jì)價(jià)值，計(jì)算公式如下：

(3)

式中:Ea為研究區(qū)單位面積糧食產(chǎn)量的經(jīng)濟(jì)價(jià)值(元/hm2);P為2018年全國(guó)糧食平均價(jià)值(元/kg),本研究以國(guó)家發(fā)改委價(jià)格監(jiān)測(cè)中心2018年的全國(guó)早中晚秈稻最低收購(gòu)價(jià)格平均值與小麥最低收購(gòu)價(jià)格的平均值作為2018年的糧食平均價(jià)格,取值為2.403元/kg;Q為根據(jù)1990年、2000年、2010年、2018年中國(guó)統(tǒng)計(jì)年鑒中31個(gè)省、自治區(qū)、直轄市糧食產(chǎn)量數(shù)據(jù)計(jì)算得到的1990—2018年各研究區(qū)域平均糧食單產(chǎn)(kg/hm2)(表1)。

表1 1990-2018年各地區(qū)糧食平均單產(chǎn) kg/hm2

然而，隨著社會(huì)不斷發(fā)展，不同時(shí)間各研究區(qū)域社會(huì)發(fā)展程度有所不同，生態(tài)系統(tǒng)服務(wù)價(jià)值也會(huì)隨之產(chǎn)生變化。為提高評(píng)價(jià)模型的精確性，本研究將以2018年為基準(zhǔn)年，引入與社會(huì)發(fā)展相關(guān)的支付意愿系數(shù)(Dt)，構(gòu)建適用于中國(guó)地級(jí)及以上城市的生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估模型[27]，其公式如下：

(4)

式中:ESV為研究區(qū)生態(tài)系統(tǒng)服務(wù)總價(jià)值(億元);Ai為第i類(lèi)生態(tài)系統(tǒng)的面積(km2);Mij為第i類(lèi)土地第j種生態(tài)系統(tǒng)服務(wù)功能所對(duì)應(yīng)當(dāng)量因子表數(shù)值;m為生態(tài)系統(tǒng)類(lèi)型數(shù);n為生態(tài)系統(tǒng)服務(wù)類(lèi)功能項(xiàng)數(shù);Ea表示單位面積糧食產(chǎn)量的經(jīng)濟(jì)價(jià)值(元/hm2)；Dt為引入恩格爾系數(shù)通過(guò)Logistic回歸模型修正后的社會(huì)發(fā)展階段系數(shù),詳見(jiàn)參考文獻(xiàn)[27]。

1.2.3 生態(tài)系統(tǒng)敏感性指數(shù) 敏感性指數(shù)(Coeficient of Sensitive,CS)是指當(dāng)生態(tài)服務(wù)價(jià)值系數(shù)變動(dòng)1%時(shí)所引起的生態(tài)系統(tǒng)服務(wù)價(jià)值的變化情況[28]。其公式如下：

(5)

式中:ESVi和ESVj分別為調(diào)整1%前、調(diào)整1%后的生態(tài)服務(wù)價(jià)值;VCik和VCjk則分別為調(diào)整前后的生態(tài)服務(wù)價(jià)值系數(shù)。

1.2.4 空間自相關(guān)分析 空間自相關(guān)主要是用來(lái)研究某一變量在一定研究范圍內(nèi)是否成相關(guān)特征的空間分析方法，主要包括全局空間自相關(guān)(Moran′sI指數(shù))和局部空間自相關(guān)(Local Moran′sI指數(shù))[15]。

全局Moran′sI指數(shù)計(jì)算公式如下：

(6)

對(duì)于局部空間單元i，其Local Moran′sI指數(shù)計(jì)算公式為：

(7)

為了揭示多個(gè)要素的空間相關(guān)性，相關(guān)學(xué)者在Moran′sI指數(shù)的基礎(chǔ)上進(jìn)一步拓展了雙變量全局自相關(guān)和局部自相關(guān)，其公式為：

(8)

2 結(jié)果與分析

2.1 土地利用變化分析

根據(jù)各年份土地利用數(shù)據(jù)，利用公式(1)和公式(2)分別計(jì)算1990—2000年、2000—2010年、2010—2018年以及1990—2018年的單一土地利用動(dòng)態(tài)度和綜合土地利用動(dòng)態(tài)度。如圖1—2和表2所示，中國(guó)地級(jí)及以上城市不同土地利用類(lèi)型在空間上存在著明顯的區(qū)域分異特征，林地、草地分布面積較大。從整個(gè)時(shí)間跨度來(lái)看，除草地呈減少的趨勢(shì)外，其他土地利用類(lèi)型均呈增加的趨勢(shì)。其中，建設(shè)用地在研究區(qū)內(nèi)變化速率最大，達(dá)2.50%，整個(gè)研究期共增加10.72萬(wàn)km2；水域和未利用地次之，整個(gè)研究期內(nèi)土地利用動(dòng)態(tài)度為0.29%，0.40%；耕地和林地變化速率則相對(duì)較小，整體上呈增加趨勢(shì)，分別增加了1.50萬(wàn)km2，2.33萬(wàn)km2；同時(shí)，整個(gè)研究時(shí)段內(nèi)綜合土地利用動(dòng)態(tài)度大于0.15%的區(qū)域主要集中于東北、西北和西部方向。

圖1 1990-2018年中國(guó)地級(jí)及以上城市土地利用類(lèi)型變化情況

注：基于標(biāo)準(zhǔn)地圖服務(wù)系統(tǒng)下載的審圖號(hào)GS(2016)1585號(hào)的標(biāo)準(zhǔn)地圖制作，底圖未做修改，下圖同。

從單個(gè)時(shí)間段跨度來(lái)看，草地的土地利用動(dòng)態(tài)度一直為負(fù)值，其減少幅度在2010—2018年最明顯，為-1.17%；1990—2000年，耕地的土地利用動(dòng)態(tài)度主要表現(xiàn)為增加的趨勢(shì)，在2000—2010年、2010—2018年兩個(gè)時(shí)段內(nèi)均以較小的幅度減少，減少幅度不斷增加；林地的土地利用動(dòng)態(tài)度在1990—2000年為負(fù)值，在2000—2018年為正值，且在各個(gè)時(shí)間段始終小于0.1%，變化幅度較不明顯；水域一直保持增加的趨勢(shì)，尤其是2000—2010年期間變化最為強(qiáng)烈；未利用地在1990—2010年時(shí)段內(nèi)以減少的趨勢(shì)變化，變化幅度很小，在2010—2018年土地利用的變化幅度略微增大，呈增加趨勢(shì)；建設(shè)用地面積則在研究時(shí)段內(nèi)一直保持增加的趨勢(shì)，變化幅度逐漸強(qiáng)烈，這主要是由于人口增加、城市擴(kuò)張等原因，造成建設(shè)用地面積也隨之逐漸增加；同時(shí)，綜合土地利用動(dòng)態(tài)度在2010—2018年變化最為明顯，相對(duì)于其他類(lèi)型土地而言，草地變化面積占總變化面積比重較大，且主要分布于西部地區(qū)，使得綜合土地利用動(dòng)態(tài)度變化明顯區(qū)域主要集中于西部地區(qū)。

表2 研究區(qū)1990-2018年單一土地利用變化年動(dòng)態(tài)度 %

2.2 生態(tài)系統(tǒng)服務(wù)價(jià)值變化分析

從時(shí)間尺度來(lái)看，研究區(qū)在1990年、2000年、2010年、2018年的生態(tài)服務(wù)價(jià)值ESV分別為：2.354×1013，2.253×1013，2.264×1013，2.197×1013元，整體呈減少的趨勢(shì)。由表3可知，整個(gè)研究期內(nèi)，生態(tài)服務(wù)總價(jià)值減少了1.57×1012元。1990—2010年期間，耕地的生態(tài)服務(wù)價(jià)值增加幅度逐漸減小，在2010—2018年則呈減少的趨勢(shì)。林地和草地是我國(guó)最主要的兩種土地利用類(lèi)型，是構(gòu)成生態(tài)系統(tǒng)服務(wù)總價(jià)值的重要組成部分，林地的生態(tài)服務(wù)價(jià)值在生態(tài)服務(wù)總價(jià)值中占比最大，在1990—2018年整個(gè)研究時(shí)段內(nèi)，林地和草地兩者的生態(tài)服務(wù)價(jià)值均呈減少特征，其中，林地減少2.8×1011元，草地減少了1.47×1012元。與生態(tài)服務(wù)總價(jià)值相比，由于水域占總面積比例小、未利用地產(chǎn)生的生態(tài)效用較小，使得兩者產(chǎn)生的生態(tài)服務(wù)價(jià)值較小，變化幅度較為穩(wěn)定。

從空間尺度來(lái)看(圖3)，生態(tài)服務(wù)總價(jià)值超過(guò)1.5×1011元的區(qū)域主要分布于西藏、新疆和內(nèi)蒙古等西北地區(qū)一帶，這些地區(qū)海拔較高，草地分布較為廣泛。生態(tài)服務(wù)價(jià)值在5×1010～1.5×1011元的區(qū)域主要集中于東北、東南和南方一帶，這是由于所處地區(qū)氣候和地理?xiàng)l件原因，使得林地面積分布范圍廣泛，產(chǎn)生生態(tài)服務(wù)價(jià)值較高。中東部則主要以耕地為主，該地帶生態(tài)服務(wù)價(jià)值普遍低于5×1010元。

表3 1990-2018年研究區(qū)土地利用類(lèi)型生態(tài)服務(wù)價(jià)值的變化 1010元

圖3 1990-2018年研究區(qū)土地利用類(lèi)型生態(tài)服務(wù)價(jià)值在空間上的變化

2.3 敏感性分析

通過(guò)公式(5)，將整個(gè)研究區(qū)域的各項(xiàng)生態(tài)服務(wù)價(jià)值系數(shù)分別變動(dòng)1%，計(jì)算不同研究時(shí)期各土地利用類(lèi)型的生態(tài)服務(wù)價(jià)值敏感性指數(shù)(表4)。整個(gè)研究期內(nèi)，不同土地利用類(lèi)型的敏感性指數(shù)一直處于穩(wěn)定狀態(tài)，且均小于1，說(shuō)明生態(tài)服務(wù)價(jià)值系數(shù)在研究區(qū)內(nèi)存在合理。另外，研究區(qū)內(nèi)各土地利用類(lèi)型敏感性指數(shù)從高到低依次為林地>草地>耕地>未利用地>水域。以上結(jié)果表明，林地對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的貢獻(xiàn)最大，未利用地和水域貢獻(xiàn)最小。因此，為了保證生態(tài)系統(tǒng)服務(wù)價(jià)值的穩(wěn)定協(xié)調(diào)，應(yīng)加大對(duì)林地的保護(hù)。

表4 生態(tài)系統(tǒng)服務(wù)價(jià)值的敏感性指數(shù)

2.4 空間相關(guān)分析

2.4.1 生態(tài)服務(wù)價(jià)值變動(dòng)與土地利用變化動(dòng)態(tài)度的相關(guān)性分析 采用Geoda軟件，構(gòu)建空間權(quán)重矩陣，得到1990—2000年、2000—2010年、2010—2018年以及1990—2018年4個(gè)時(shí)間段內(nèi)生態(tài)服務(wù)價(jià)值變化量與綜合土地利用動(dòng)態(tài)度的全局Moran′sI值，其結(jié)果為0.11，-0.08，-0.27，-0.18。結(jié)果表明，除1990—2000年，其余幾個(gè)時(shí)間段內(nèi)生態(tài)服務(wù)價(jià)值變化量與綜合土地利用動(dòng)態(tài)度均呈負(fù)相關(guān)。同時(shí)根據(jù)z檢驗(yàn)結(jié)果，繪制1990—2018年整個(gè)研究期間內(nèi)生態(tài)服務(wù)價(jià)值與綜合土地利用變化動(dòng)態(tài)度的雙變量局部空間自相關(guān)LISA聚集圖。由圖4可知(高—高型表示生態(tài)服務(wù)價(jià)值與綜合土地利用動(dòng)態(tài)度均表現(xiàn)為較高值；高—低型表示生態(tài)服務(wù)價(jià)值較高，綜合土地利用動(dòng)態(tài)度表現(xiàn)為較低值；低—高型表示生態(tài)服務(wù)價(jià)值較低，綜合土地利用動(dòng)態(tài)度表現(xiàn)為較高值；低—低型表示生態(tài)服務(wù)價(jià)值與綜合土地利用動(dòng)態(tài)度均表現(xiàn)為較低值)，高—低型主要分布在中部地區(qū)，低—高型主要分布于西藏地區(qū)，低—低型主要分布于南部，高—高型則零星分布。

2.4.2 生態(tài)服務(wù)價(jià)值與社會(huì)經(jīng)濟(jì)的相關(guān)性分析 通過(guò)收集1990年、2000年、2010年、2018年的中國(guó)地級(jí)及以上城市總?cè)丝跀?shù)、城鎮(zhèn)人口數(shù)、城鎮(zhèn)化率以及GDP這4個(gè)指標(biāo)，利用歸一化和熵值法求取社會(huì)經(jīng)濟(jì)指標(biāo)綜合值，再結(jié)合生態(tài)服務(wù)價(jià)值進(jìn)行空間自相關(guān)分析，得到4個(gè)時(shí)期研究區(qū)生態(tài)服務(wù)價(jià)值與社會(huì)經(jīng)濟(jì)指標(biāo)綜合值的全局Moran′sI值，其結(jié)果分別為-0.14，-0.09，-0.10，-0.13。結(jié)果可知，生態(tài)服務(wù)價(jià)值與社會(huì)經(jīng)濟(jì)之間存在顯著的空間負(fù)相關(guān)性，且兩者在研究期內(nèi)z檢驗(yàn)值均<-1.96。根據(jù)生態(tài)服務(wù)價(jià)值與社會(huì)經(jīng)濟(jì)指標(biāo)綜合值的雙變量空間自相關(guān)LISA聚集圖(圖5)，1990—2018年研究區(qū)生態(tài)服務(wù)價(jià)值與社會(huì)經(jīng)濟(jì)的局部空間特點(diǎn)總體上較為統(tǒng)一。其中，高—高型主要分布于中部和東部偏北地區(qū)，該區(qū)域經(jīng)濟(jì)較為發(fā)達(dá)，且林地所占面積較大，其生態(tài)價(jià)值較高；低—高型主要分布于中、東部地區(qū)，這部分地區(qū)經(jīng)濟(jì)發(fā)展較快，對(duì)生態(tài)資源占用現(xiàn)象嚴(yán)重；低—低型主要分布于西南地區(qū)，該區(qū)域山多地少、交通相對(duì)閉塞、地理?xiàng)l件較差，其經(jīng)濟(jì)和生態(tài)條件均處于落后階段；高—低型則主要分布于昌都、噶爾縣、拉薩和喀什市等西部地區(qū)，這是因?yàn)槲鞑康貐^(qū)草地面積分布范圍廣泛，生態(tài)服務(wù)價(jià)值較高，但其經(jīng)濟(jì)發(fā)展較為落后。

圖4 1990-2018年生態(tài)服務(wù)價(jià)值與土地利用變化動(dòng)態(tài)度的

圖5 1990-2018年生態(tài)服務(wù)價(jià)值與社會(huì)經(jīng)濟(jì)的雙變量空間自相關(guān)LISA集聚圖

3 討論與結(jié)論

(1) 1990—2018年，草地在整個(gè)研究期內(nèi)一直呈減少的趨勢(shì)，共減少了38.74萬(wàn)km2；耕地在1990—2000年表現(xiàn)為增加的趨勢(shì)，在2000—2018年呈減少的趨勢(shì)，減少幅度略微增大，但整個(gè)研究期內(nèi)仍保持增加趨勢(shì)，凈面積共增加了1.50萬(wàn)km2；林地在各個(gè)研究階段變化幅度不明顯；水域的土地利用動(dòng)態(tài)度一直為正值，整個(gè)研究期內(nèi)增加了2.06萬(wàn)km2；未利用地在1990—2010年以減少的趨勢(shì)變化，2010—2018年有所增加，整體仍呈增加的趨勢(shì)；建設(shè)用地在整個(gè)研究期內(nèi)變化最為強(qiáng)烈，土地利用動(dòng)態(tài)度始終為正值，增加了10.72萬(wàn)km2；綜合土地利用動(dòng)態(tài)度隨著時(shí)間變化也有所增加，在2010—2018年變化尤為明顯。

(2) 1990—2018年，中國(guó)地級(jí)及以上城市生態(tài)服務(wù)價(jià)值較高的區(qū)域主要集中于西北地區(qū)，總體上呈減少的趨勢(shì)，28 a間共減少了1.57×1012元，而林地和草地產(chǎn)生的生態(tài)服務(wù)價(jià)值是構(gòu)成研究區(qū)生態(tài)系統(tǒng)服務(wù)總價(jià)值的主要部分，其在研究期內(nèi)分別減少了2.8×1011元、1.47×1012元。

(3) 1990—2018年，各土地利用類(lèi)型的生態(tài)服務(wù)價(jià)值敏感性指數(shù)在研究區(qū)內(nèi)均表現(xiàn)為<1，研究期內(nèi)不同土地利用類(lèi)型敏感性指數(shù)由高到低依次為林地、草地、耕地、未利用地以及水域。其中，林地敏感性指數(shù)最高，產(chǎn)生的生態(tài)服務(wù)價(jià)值最大。

(4) 1990—2018年，生態(tài)服務(wù)價(jià)值與綜合土地利用動(dòng)態(tài)度整體上呈負(fù)相關(guān)，各年份生態(tài)服務(wù)價(jià)值與社會(huì)經(jīng)濟(jì)指標(biāo)綜合值呈負(fù)相關(guān)，且雙變量LISA圖空間特點(diǎn)較為一致。低—高型主要分布于經(jīng)濟(jì)發(fā)展較快的中、東部地區(qū)；高—低型則主要分布于草地面積廣泛經(jīng)濟(jì)發(fā)展較為落后的西部地區(qū)。這表明社會(huì)經(jīng)濟(jì)的發(fā)展、土地利用格局的變化一定程度上影響著生態(tài)系統(tǒng)結(jié)構(gòu)和功能的轉(zhuǎn)變。當(dāng)?shù)卣块T(mén)應(yīng)結(jié)合實(shí)際情況，通過(guò)加大對(duì)生態(tài)系統(tǒng)的保護(hù)力度，采取有效措施保證土地合理規(guī)劃，促進(jìn)社會(huì)—經(jīng)濟(jì)—生態(tài)和諧發(fā)展。

與此同時(shí)，文中仍存在一些缺陷。本文僅根據(jù)謝高地等[23-25]制作的中國(guó)陸地生態(tài)系統(tǒng)單位面積生態(tài)服務(wù)價(jià)值表測(cè)算中國(guó)各地級(jí)城市生態(tài)服務(wù)價(jià)值，未根據(jù)各地區(qū)實(shí)際情況制定出適用于各地區(qū)的陸地生態(tài)系統(tǒng)單位面積生態(tài)服務(wù)價(jià)值表，導(dǎo)致測(cè)算有所偏差，在今后需要加強(qiáng)這一方面的研究。