韓 蕊，孫思琦，郭 濼①，陳奕竹

(1.中央民族大學生命與環(huán)境科學學院，北京 100081；2.北京大學城市與環(huán)境學院，北京 100871)

生態(tài)系統(tǒng)服務價值(ecosystem service value， ESV)是生態(tài)系統(tǒng)與過程所形成并維持人類賴以生存的自然環(huán)境條件與功能[1]，包括自然資本的物流、能流和信息流構成的生態(tài)系統(tǒng)服務和非自然資本生產(chǎn)的人類福利[2]。生態(tài)系統(tǒng)不僅是人類賴以生存之本，更是區(qū)域可持續(xù)發(fā)展核心[3]，掌握生態(tài)系統(tǒng)服務功能時空演變規(guī)律是開展科學生態(tài)系統(tǒng)管理的重要前提[4]。1997年COSTANZA等[2]基于土地利用類型開展了生態(tài)系統(tǒng)服務價值核算研究，推動了生態(tài)系統(tǒng)結構和功能研究的發(fā)展。如KREUTER等[5]對美國德州貝克薩爾縣3個主要流域進行評估，KONARSKA等[6]也對全美生態(tài)系統(tǒng)功能進行量化研究。謝高地等[7-8]根據(jù)我國實際情況建立并修訂了中國陸地二級生態(tài)系統(tǒng)單位面積服務價值當量表，已被用于區(qū)域生態(tài)系統(tǒng)服務價值評估[9-10]。土地利用/覆被變化(land-use and land-cover change， LUCC)是全球生態(tài)系統(tǒng)服務研究中關鍵切入點，環(huán)境問題很可能引起大尺度、轉換性土地利用類型變化[11]，造成土地利用格局改變，進而對生態(tài)系統(tǒng)產(chǎn)生重大影響。因此將生態(tài)系統(tǒng)服務價值與土地利用類型進行相關關系分析，有助于了解土地利用類型與生態(tài)系統(tǒng)服務價值的關聯(lián)度，揭示人類活動對地表生態(tài)過程及效應的影響機制。

川東地區(qū)是我國中西部地區(qū)重要經(jīng)濟建設區(qū)和糧食生產(chǎn)基地，經(jīng)濟和城市化快速發(fā)展對研究區(qū)生態(tài)系統(tǒng)構成嚴重威脅。因此，探究研究區(qū)生態(tài)系統(tǒng)服務價值時空演變驅動因素，有助于深入了解社會在滿足發(fā)展需求的同時對生態(tài)系統(tǒng)服務功能的改變程度，對川東地區(qū)生態(tài)環(huán)境保護和可持續(xù)發(fā)展具有重要意義。

1 材料與方法

1.1 研究區(qū)概況

川東地區(qū)位于四川盆地東部，處于27°48′～32°58′ N，105°12′～110°22′ E，毗鄰陜西、湖北、湖南和貴州省，地處青藏高原與長江中下游平原過渡帶，全區(qū)包含重慶市、南充市、達州市、廣安市和巴中市的5市54縣(區(qū)、市)，總面積為126 714 km2(圖1)。研究區(qū)屬亞熱帶季風性濕潤氣候區(qū)，年平均氣溫為16～18 ℃，年平均降水量為1 000 mm，年平均相對濕度為70%～80%，日照時數(shù)為1 200 h。研究區(qū)海拔范圍為73.1～2 796.8 m之間，地形差異大，地勢表現(xiàn)為由南北部向中西部地區(qū)逐級降低。研究區(qū)水系眾多，縱橫交錯，分布著眾多國家重點生態(tài)功能區(qū)，擁有豐富的水、礦產(chǎn)和動植物資源。研究區(qū)是長江上游及西南地區(qū)最大的經(jīng)濟中心，2015年戶籍總人口為5 222.52萬，全年生產(chǎn)總值約為2.01萬億元。

圖1 川東地區(qū)行政區(qū)劃

1.2 研究方法

1.2.1土地利用動態(tài)度

土地利用類型動態(tài)度是定量評價土地利用變化速度的指標，包括單一土地利用類型動態(tài)度和綜合土地利用動態(tài)度[12]，計算公式如下：

單一土地利用類型動態(tài)度：

(1)

式(1)中，K為研究期某種土地利用類型動態(tài)度，%；Ua和Ub分別為研究期初和期末某種土地利用類型面積，hm2；T為研究時長，a。

綜合土地利用動態(tài)度(C，%)：

(2)

式(2)中，Ui為研究期初第i類土地利用類型面積，hm2；Ui，j為研究期第i類土地利用類型轉為非i類土地利用類型面積的絕對值，hm2；T為研究時長，a。

1.2.2生態(tài)系統(tǒng)服務價值評估

參照生態(tài)系統(tǒng)服務價值研究[13-15]，采用謝高地等[7-8]的價值當量換算法確定川東地區(qū)生態(tài)系統(tǒng)服務價值。將生態(tài)系統(tǒng)生產(chǎn)凈利潤作為該生態(tài)系統(tǒng)所能提供的生產(chǎn)價值，將單位面積農(nóng)田生態(tài)系統(tǒng)糧食生產(chǎn)凈利潤作為1個標準當量因子的生態(tài)系統(tǒng)服務價值量。中國2010年生態(tài)系統(tǒng)價值當量的經(jīng)濟價值為3 406.50元·hm-2，考慮到數(shù)據(jù)可獲得性，以重慶市單位面積糧食產(chǎn)量代替川東地區(qū)單位面積糧食產(chǎn)量，根據(jù)《重慶統(tǒng)計年鑒2010》中重慶市單位面積糧食產(chǎn)量為5 152.35 kg·hm-2，而同期全國單位面積糧食產(chǎn)量為4 973.58 kg·hm-2，確定修正系數(shù)為1.036，得到研究區(qū)生態(tài)系統(tǒng)服務價值當量為3 529.13元·hm-2。根據(jù)單位面積生態(tài)系統(tǒng)服務價值當量因子表[7]和修正后的價值當量得到川東地區(qū)各土地利用類型生態(tài)系統(tǒng)服務價值當量。川東地區(qū)生態(tài)系統(tǒng)服務價值計算公式為

(3)

(4)

式(3)～(4)中，V為研究區(qū)總生態(tài)系統(tǒng)服務價值，元；Ai為研究區(qū)第i種土地利用類型面積，hm2；VC，i為研究區(qū)第i種土地利用類型生態(tài)系統(tǒng)服務價值當量，元·hm-2；Vf為第f項生態(tài)系統(tǒng)服務價值，元；VCf，i為研究區(qū)第i種土地利用類型第f項生態(tài)系統(tǒng)服務價值當量，元·hm-2。

采用SPSS 19.0軟件計算不同土地利用類型面積變化與總生態(tài)系統(tǒng)服務價值變化的相關性。當相關系數(shù)為>0.8～1.0時，表示兩者相關性極強，當為>0.6～0.8時，表示相關性強，當為>0.4～0.6時，表示相關性中等，當為>0.2～0.4時，表示相關性弱，當為0～0.2時，表示相關性極弱或無相關。

1.2.3驅動因素分析

生態(tài)系統(tǒng)服務價值時空差異的驅動因素包括自然和人為因素兩大類，在短時間尺度上主要受人為因素的干擾。參照文獻[16]，并考慮到數(shù)據(jù)可獲取性，選擇總人口、城鎮(zhèn)人口、農(nóng)村人口、人口密度和城鎮(zhèn)化率5項人口因素指標，GDP總量、人均國內(nèi)生產(chǎn)總值、第一產(chǎn)業(yè)總產(chǎn)值、第二產(chǎn)業(yè)總產(chǎn)值、第三產(chǎn)業(yè)總產(chǎn)值、農(nóng)業(yè)總產(chǎn)值、林業(yè)總產(chǎn)值、牧業(yè)總產(chǎn)值、漁業(yè)總產(chǎn)值和工業(yè)總產(chǎn)值10項經(jīng)濟因素指標以及國內(nèi)旅游收入、入境旅游外匯收入和境外游客人數(shù)3項旅游因素指標，采用SPSS 19.0軟件進行主成分及相關性分析，選出其中對研究區(qū)生態(tài)系統(tǒng)服務價值時空變化影響較大的指標作為生態(tài)系統(tǒng)服務價值時空變化的主要驅動因素，并進行多元線性回歸分析。

1.2.4空間自相關分析

參照文獻[17-18]，綜合驅動因素分析結果，最終選取表征川東地區(qū)各縣市區(qū)人口經(jīng)濟要素的人均國內(nèi)生產(chǎn)總值、人口密度和城市化率3個指標，分別與生態(tài)系統(tǒng)服務價值變化率進行空間相關性分析。采用GeoDa模型中雙變量莫蘭指數(shù)(bivariate Moran′sI)分析研究區(qū)生態(tài)系統(tǒng)服務價值與社會經(jīng)濟指標的空間相關性，并繪制LISA集聚圖，將空間自相關達到顯著性水平α=0.05的區(qū)域分為4類：高-高型為社會經(jīng)濟指標變化率和生態(tài)系統(tǒng)服務價值變化率都高；低-低型為社會經(jīng)濟指標變化率和生態(tài)系統(tǒng)服務價值變化率都低；低-高型為社會經(jīng)濟指標變化率低，而生態(tài)系統(tǒng)服務價值變化率高；高-低型為社會經(jīng)濟指標變化率高，而生態(tài)系統(tǒng)服務價值變化率低。

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

川東地區(qū)1990、1995、2000、2005、2010和2015年共6期土地利用類型數(shù)據(jù)來源于中國土地利用現(xiàn)狀遙感監(jiān)測數(shù)據(jù)庫(http：∥www.resdc.cn)，由Landsat TM/ETM遙感影像人工目視解譯生成，空間分辨率為1 km。經(jīng)室外調(diào)查及隨機抽樣檢查，土地利用一級類型平均解譯精度達94.3%[19-20]。川東地區(qū)土地利用一級類型為林地、草地、水域、耕地、建筑用地和未利用地(圖2)，其中二級分類如下：耕地分為水田和旱地；林地分為針葉林、闊葉林、針闊混交林和灌木林；草地分為草甸、草原和灌草叢。

18項社會經(jīng)濟數(shù)據(jù)來源于四川統(tǒng)計年鑒和重慶統(tǒng)計年鑒。中國單位面積糧食產(chǎn)量和重慶市單位面積糧食產(chǎn)量源于同時期全國統(tǒng)計數(shù)據(jù)年鑒和重慶統(tǒng)計年鑒。

2 結果與分析

2.1 川東地區(qū)土地利用變化特征

1990—2015年川東地區(qū)最主要景觀類型為耕地，占總面積的70%～74%，其次為林地，占總面積的22%(圖2)。由表1可知，25 a來綜合土地利用動態(tài)度<0.1%，1990—2010年綜合土地利用動態(tài)度逐漸降低，而2010—2015年增加，但總體來說，研究區(qū)土地利用綜合動態(tài)度變化較小。

1990—2015年林地、水體和建筑用地單一土地利用類型動態(tài)度呈正向變化，而耕地、草地和未利用地單一土地利用類型動態(tài)度呈負向變化。1990—2015年建筑用地面積變化極為明顯，增加15.92×104hm2，而未利用地和草地面積各減少100×104和44.49×104hm2，這可能與研究區(qū)城鎮(zhèn)化率增加有關。研究區(qū)作為長江上游地區(qū)重要生態(tài)屏障，分布著眾多生態(tài)功能區(qū)，水源涵養(yǎng)功能保持較好，25 a來水域面積增加2.47×104hm2，林地面積增加3.35×104hm2，而耕地面積減少17.25×104hm2。

圖2 川東地區(qū)土地利用分類Fig.2 Land use type in Chuandong region

表1 1990—2015年川東地區(qū)土地利用景觀類型動態(tài)度

Table 1 Dynamic degree and the information entropy of land use in Chuandong region from 1990 to 2015

時期單一土地利用類型動態(tài)度/%耕地林地草地水域建設用地未利用地綜合土地利用動態(tài)度/%1990—1995年-0.260.310.792.112.512.860.161995—2000年0.290.05-0.86-0.574.11-2.000.142000—2005年-0.120.28-0.28-1.183.92-0.170.112005—2010年-0.150.050.09-1.377.2600.112010—2015年-0.23-0.02-0.080.4312.810.610.181990—2015年-0.100.14-0.070.9710.51-0.110.09

2.2 川東地區(qū)生態(tài)系統(tǒng)服務價值時空變化

2.2.1生態(tài)系統(tǒng)服務價值時間維度差異性分析

1990—2015年研究區(qū)總生態(tài)系統(tǒng)服務價值增長143.14億元(表2)。各單項土地利用類型中，林地和水域生態(tài)系統(tǒng)服務價值呈增長趨勢，耕地和草地生態(tài)系統(tǒng)服務價值呈下降趨勢?？傮w來看，研究區(qū)生態(tài)系統(tǒng)服務價值變化趨勢與土地利用動態(tài)度變化趨勢相似。

表2 1990—2015年川東地區(qū)各土地利用類型生態(tài)系統(tǒng)服務價值變化

Table 2 Changes of ecosystem service value in Chuandong region from 1990 to 2015

年份耕地林地草地水域未利用地ESV/億元占比/%ESV/億元占比/%ESV/億元占比/%ESV/億元占比/%ESV/億元占比/%總計/億元1990973.1432.281 279.5742.45151.215.02610.6120.250.020.003 014.551995960.5431.091 299.7542.07156.695.07672.5221.770.030.003 089.532000974.4231.641 301.7542.27150.124.87653.6121.220.030.003 079.932005968.6131.091 319.3942.34147.804.74680.3021.830.020.003 116.122010961.3530.471 322.5141.91148.664.71722.7622.910.020.003 155.302015950.3030.101 321.3041.84148.044.69738.0123.370.020.003 157.67

ESV為生態(tài)系統(tǒng)服務價值。

從不同土地利用類型(表2)來看，林地對研究區(qū)生態(tài)系統(tǒng)服務價值貢獻率最大，耕地次之，1990—2015年兩者貢獻率分別超過41%和30%，分別約占總面積的22%和72%。林地屬于自然生態(tài)系統(tǒng)，耕地屬于人工生態(tài)系統(tǒng)，林地對生態(tài)系統(tǒng)服務功能發(fā)揮著關鍵作用。水域、草地和未利用地面積相對較小，生態(tài)系統(tǒng)服務價值也明顯小于耕地和林地?？傮w來看，研究區(qū)不同土地利用類型對總生態(tài)系統(tǒng)服務價值貢獻率由大到小依次為林地、耕地、水域、草地和未利用地。

1990—2015年川東地區(qū)各項生態(tài)系統(tǒng)服務價值見表3。川東生態(tài)系統(tǒng)各服務類型中價值最突出的為調(diào)節(jié)服務，占總生態(tài)系統(tǒng)服務價值的70%，支持服務和供給服務對總價值的貢獻率分別為20%和7%，文化服務能力價值最低，僅占總價值的3%左右。同時，25 a來調(diào)節(jié)服務和供給服務生態(tài)系統(tǒng)服務價值呈增長趨勢，支持服務和文化服務變化不大，這表明研究區(qū)生態(tài)系統(tǒng)服務價值格局比較穩(wěn)定。

表3 1990—2015年川東地區(qū)各項生態(tài)系統(tǒng)服務功能價值估算

Table 3 Ecosystem service value of each kind of the service function in Chuandong region from 1990 to 2015

生態(tài)系統(tǒng)服務類型生態(tài)系統(tǒng)服務價值/億元一級二級1990年1995年2000年2005年2010年2015年增長率/%供給服務食物生產(chǎn)276.83273.59276.22275.06273.62270.54-2.27原料生產(chǎn)113.23113.59114.85114.82114.38113.710.42水資源供給-169.39-158.43-159.24-155.92-152.00-147.0813.17調(diào)節(jié)服務氣體調(diào)節(jié)344.50344.16345.79345.92345.06342.44-0.60氣候調(diào)節(jié)522.34529.12528.75532.52533.47531.921.83凈化環(huán)境179.79184.25183.32185.37187.29187.434.25水文調(diào)節(jié)1 062.941 111.041 095.671 122.651 157.501 165.039.60支持服務土壤保持333.62336.36339.33339.88339.06337.611.20維持養(yǎng)分48.7548.4648.8048.7248.5248.07-1.39生物多樣性206.71210.22209.64210.06210.87210.531.85文化服務美學景觀95.2397.1796.8097.0497.5397.472.35總計3 014.553 089.533 079.933 116.123 155.303 157.67

從生態(tài)系統(tǒng)服務功能二級分類來看，水文調(diào)節(jié)和氣候調(diào)節(jié)功能生態(tài)系統(tǒng)服務價值占比最大，分別為35%和17%，氣體調(diào)節(jié)、土壤保持、食物供給和生物多樣性功能占比也均大于6%。其中，25 a來食物供給價值降低與耕地面積減少密切相關，土壤保持和生物多樣性功能價值增加與林地面積增加密切相關。維持養(yǎng)分循環(huán)功能價值占比最小，約為2%，且呈減少趨勢，這與草地和耕地生態(tài)系統(tǒng)服務價值變化趨勢一致，這表明養(yǎng)分循環(huán)功能可能受耕地和草地面積減少的影響。

2.2.2土地利用類型面積變化與生態(tài)系統(tǒng)服務價值相關性分析

生態(tài)系統(tǒng)服務價值與各土地利用類型面積有相關關系，不同土地利用類型面積變化與生態(tài)系統(tǒng)服務價值間相關性有所不同。由表4可知，川東地區(qū)不同土地利用類型對生態(tài)系統(tǒng)服務價值影響程度不同。林地和水域與生態(tài)系統(tǒng)服務價值相關性較顯著；其次為耕地，與生態(tài)系統(tǒng)服務價值呈明顯負相關關系。

二級生態(tài)系統(tǒng)服務類型中食物生產(chǎn)、水資源供給、氣體調(diào)節(jié)、水文調(diào)節(jié)和維持養(yǎng)分與耕地呈顯著相關；水資源供給、氣候調(diào)節(jié)、凈化環(huán)境、水文調(diào)節(jié)、土壤保持、生物多樣性和美學景觀與林地呈顯著相關；食物生產(chǎn)、水資源供給、氣候調(diào)節(jié)、凈化環(huán)境、水文調(diào)節(jié)、生物多樣性和美學景觀與水域呈顯著相關(P<0.05或P<0.01)。

表4 1990—2015年川東地區(qū)土地利用類型面積變化與生態(tài)系統(tǒng)服務價值相關性系數(shù)矩陣

Table 4 Correlation coefficient matrix between land use type area change and ecosystem service value in Chuandong region from 1990 to 2015

生態(tài)系統(tǒng)服務類型一級二級 耕地林地草地水域未利用地供給服務食物生產(chǎn)0.99??-0.67 0.11 -0.91?? 0.65原料生產(chǎn)0.300.55-0.460.18-0.38水資源供給-0.81?0.93??-0.380.98??-0.32調(diào)節(jié)服務氣體調(diào)節(jié)0.80?-0.07-0.10-0.45-0.15氣候調(diào)節(jié)-0.560.98??-0.370.87?-0.36凈化環(huán)境-0.760.96??-0.360.97??-0.34水文調(diào)節(jié)-0.83?0.93??-0.400.98??-0.38支持服務土壤保持-0.090.81?-0.470.54-0.39維持養(yǎng)分0.98??-0.480.04-0.79-0.00生物多樣性-0.610.89?-0.130.88?-0.09文化服務美學景觀-0.680.88?-0.120.92??-0.07

*表示在0.05水平(雙側)上顯著相關，**表示在0.01水平(雙側)上顯著相關。

2.2.3生態(tài)系統(tǒng)服務價值空間維度差異性分析

川東地區(qū)單位面積生態(tài)系統(tǒng)服務價值總體分布趨勢為由東北和東南兩翼地區(qū)向中部及西部地區(qū)逐級遞減，這主要是因為東北和東南地區(qū)是研究區(qū)重要生態(tài)功能區(qū)及自然保護區(qū)的集中分布地，主要土地利用類型為林地，生態(tài)系統(tǒng)服務功能好，而中部及西部低海拔區(qū)主要土地利用類型為耕地，生態(tài)系統(tǒng)服務功能較差。1990—2015年研究區(qū)19個縣(區(qū)、市)單位面積生態(tài)系統(tǒng)服務價值呈下降趨勢，而36個縣(區(qū)、市)呈增加趨勢，主要分布在東北和中部以及南部部分地區(qū)(圖3)。

圖3 1990—2015年川東地區(qū)單位面積生態(tài)系統(tǒng)服務價值空間變化

2.3 川東地區(qū)生態(tài)系統(tǒng)服務價值時空差異驅動因素分析

人類活動對土地利用產(chǎn)生影響，進而影響區(qū)域生態(tài)系統(tǒng)服務價值。為分析生態(tài)系統(tǒng)服務價值時空差異的驅動因素，選取對川東地區(qū)生態(tài)系統(tǒng)服務價值可能產(chǎn)生影響的社會、經(jīng)濟和旅游3個方面指標進行相關性和主成分分析，剔除對其影響較小且不顯著的指標，構建的回歸方程為

Y=-0.246X1-1.082X2-0.022X3+

0.000 006X4+4.696 2。

(5)

式(5)中，Y為生態(tài)系統(tǒng)服務價值，億元；X1為農(nóng)業(yè)人口數(shù)，人；X2為人口密度，人·hm-2；X3為工業(yè)總產(chǎn)值，億元；X4為國內(nèi)旅游收入，萬元。

回歸方程擬合度達到1，擬合效果極好。分析可知，農(nóng)業(yè)人口數(shù)越小，人口密度越小，工業(yè)總產(chǎn)值越小，國內(nèi)旅游收入越高，生態(tài)系統(tǒng)服務價值就越高。

2.4 川東地區(qū)社會經(jīng)濟因素與生態(tài)系統(tǒng)服務價值空間相關性分析

根據(jù)文獻[17-18]，綜合多元回歸分析統(tǒng)計結果，對人口社會經(jīng)濟指標與生態(tài)系統(tǒng)服務價值變化率進行空間相關性分析。25 a來，研究區(qū)生態(tài)系統(tǒng)服務價值變化與城鎮(zhèn)率變化、人口密度變化和人均國內(nèi)生產(chǎn)總值變化的全局空間相關指數(shù)分別為-0.23、-0.14和-0.04，均呈顯著負相關(P<0.05)。由圖4可知，人均國內(nèi)生產(chǎn)總值、人口密度和城市化率與生態(tài)系統(tǒng)服務價值的局部空間相關性具有高度一致性，且空間相關性分析均達到顯著性水平，這說明這3個因素與生態(tài)系統(tǒng)服務價值變化率的局部空間相關顯著性良好。城鎮(zhèn)化率、人口密度和人均國內(nèi)生產(chǎn)總值變化率與生態(tài)系統(tǒng)服務價值變化率中低-低型和低-高型均聚集在研究區(qū)西部的西充、南充、岳池和武勝縣一帶；高-低型聚集在南部的江津市和南川市一帶；高-高型則分布在西南部的榮昌縣和雙橋區(qū)。

3 討論

生態(tài)系統(tǒng)服務功能復雜性決定其評價方法是一個不斷完善的過程，當量因子法具有數(shù)據(jù)穩(wěn)定、表達直觀等特點，適用于區(qū)域及全球尺度生態(tài)系統(tǒng)服務價值評估[21]。采用當量因子法研究川東地區(qū)生態(tài)系統(tǒng)服務價值時空演變特征及其驅動因素對于研究區(qū)可持續(xù)發(fā)展具有指導意義。

川東地區(qū)土地利用類型以林地和耕地為主，是西南地區(qū)重要經(jīng)濟中心和糧食基地，同時也是長江上游重要生態(tài)屏障。由于對生態(tài)環(huán)境治理重視度加強，經(jīng)濟發(fā)展迅速，1990—2015年研究區(qū)土地利用變化表現(xiàn)出兩個趨勢，一方面受“退耕還林(草)”“天然林資源保護工程”等林業(yè)政策影響而呈耕地縮減和林地擴張的趨勢，另一方面受快速城鎮(zhèn)化影響呈現(xiàn)建筑用地擴張的趨勢。人類活動對土地利用產(chǎn)生影響，進而影響區(qū)域生態(tài)系統(tǒng)服務價值，研究發(fā)現(xiàn)25 a來川東地區(qū)土地開發(fā)力度的加強對生態(tài)系統(tǒng)服務價值具有一定負面效應，該結論也與其他區(qū)域研究結果[22]一致。同時，筆者研究發(fā)現(xiàn)川東地區(qū)生態(tài)系統(tǒng)服務分配格局的空間差異明顯，土地開發(fā)強度低的東北、東南兩翼林地聚集區(qū)生態(tài)系統(tǒng)服務單位面積價值較高，而土地開發(fā)強度高的中西部耕地及經(jīng)濟開發(fā)聚集區(qū)價值較低。這表明在繼續(xù)退耕還林(草)和重視生態(tài)用地建設與保護的同時，還要盡量注重合理規(guī)劃城市區(qū)域發(fā)展，重視耕作技術改良和農(nóng)業(yè)種植結構調(diào)整，變革低產(chǎn)田，再造優(yōu)質(zhì)田，以期提高研究區(qū)生態(tài)系統(tǒng)服務能力。

圖4 1990—2015年川東地區(qū)社會經(jīng)濟與生態(tài)系統(tǒng)服務價值局部空間相關性

從社會經(jīng)濟方面研究影響川東生態(tài)系統(tǒng)服務價值的驅動因素及空間相關性分析，結果表明生態(tài)系統(tǒng)服務價值變化受人為因素影響顯著，空間差異明顯。因此，當環(huán)境壓力增大進而導致生態(tài)系統(tǒng)服務功能降低時，政府應制定有效政策，提升生態(tài)系統(tǒng)服務功能，使自然資本能得以保值增值。隨著研究區(qū)人口增長，工業(yè)水平提升，生態(tài)系統(tǒng)供給和調(diào)節(jié)服務需求會持續(xù)增加，提高森林覆蓋和保護東南、東北兩翼的生態(tài)功能區(qū)對川東地區(qū)優(yōu)化生態(tài)環(huán)境、促進區(qū)域可持續(xù)發(fā)展具有重要意義。

4 結論

基于1990—2015年川東地區(qū)土地利用數(shù)據(jù)，分析土地利用格局和生態(tài)系統(tǒng)服務價值時空演變特征，并研究兩者相關性，進一步探究生態(tài)系統(tǒng)服務價值變化的驅動因素，為研究區(qū)生態(tài)環(huán)境保護提供決策支持。研究表明：(1)1990—2015年川東地區(qū)土地利用類型主要為耕地和林地，生態(tài)系統(tǒng)服務價值呈增長趨勢，其中生態(tài)系統(tǒng)服務中水文調(diào)節(jié)功能增長最突出。(2)單位面積生態(tài)系統(tǒng)服務價值呈由東北、東南兩翼向中西部地區(qū)逐級遞減的空間格局。(3)不同土地利用類型對生態(tài)系統(tǒng)服務價值影響程度不同，其中林地、水域與生態(tài)系統(tǒng)服務價值呈明顯正相關。(4)生態(tài)系統(tǒng)服務價值時空差異主要驅動因素為農(nóng)業(yè)人口數(shù)、人口密度、工業(yè)總產(chǎn)值和國內(nèi)旅游收入。(5)空間相關性分析結果表明，社會經(jīng)濟因素與生態(tài)系統(tǒng)服務價值變化率呈負相關，并且低-低型和高-低型聚集在西部地區(qū)，低-高型和高-高型聚集在南部地區(qū)。