賈 琦

(鄭州輕工業(yè)大學 環(huán)境設(shè)計系, 河南 鄭州 450002)

生態(tài)系統(tǒng)服務(wù)是指人類在生產(chǎn)生活過程中，通過生態(tài)系統(tǒng)代謝過程和功能服務(wù)進而直接或者間接的方式得到所需的產(chǎn)品和相關(guān)服務(wù)等，以便全方位滿足人類生存和可持續(xù)發(fā)展的需求[1]。早在2005年聯(lián)合國發(fā)布的《千年生態(tài)系統(tǒng)評估綜合報告》中便已明確指出[2]，半個世紀以來，由于人口急劇增長，人類過度開發(fā)和使用地球資源，從而導致生態(tài)功能中有60%在退化和被人類以不可持續(xù)的方式利用，并且退化程度在未來50 a內(nèi)將進一步加劇。

土地利用演化作為人類各種活動影響區(qū)域生態(tài)系統(tǒng)的主要方式之一，在很大程度上作用于生態(tài)系統(tǒng)的格局分布與變化強度，將直接影響生態(tài)系統(tǒng)服務(wù)價值(value of ecosystem service， ESV)的空間分布[3]。以人類活動作用于土地利用變化過程為視角評估生態(tài)系統(tǒng)服務(wù)價值響應(yīng)特征儼然已成為全球景觀生態(tài)學領(lǐng)域研究的熱點方向之一[1]。當前關(guān)于土地利用變化及其生態(tài)系統(tǒng)服務(wù)價值響應(yīng)研究相關(guān)成果已較為豐富，諸如寧佳等[4]基于降水貯存量法、通用土壤流失方程和修正土壤風蝕方程等方法探討黃土高原土地利用和生態(tài)系統(tǒng)服務(wù)時空變化特征；Geange等[5]基于生態(tài)系統(tǒng)服務(wù)混合矩陣方法探討新西蘭海洋生態(tài)功能演化過程及其影響范圍；王大菊等[6]采用替代花費法、市場價值法和水量平衡法對三峽庫區(qū)蓄水前后多時期的土地利用變化和生態(tài)系統(tǒng)服務(wù)價值進行量化評估。但是同類研究多探討兩者的演化過程和分布規(guī)律，尚未深入揭示各土地利用類型對生態(tài)系統(tǒng)的影響機制。又如王磊等[7]、趙麗等[8]、Kalfas等[9]、鄧偉等[10]通過探究土地利用與生態(tài)系統(tǒng)服務(wù)價值演化過程，揭示其驅(qū)動機理和影響強度，但是從相關(guān)研究結(jié)論來看，多數(shù)學者側(cè)重從土地利用變化視角入手，涉及空間相互作用視角的驅(qū)動因子考慮依舊不足。土地利用變化受區(qū)域自然條件、社會經(jīng)濟發(fā)展、規(guī)劃政策等多種因素影響，在空間上具有很強的差異性和分布規(guī)律[11-12]。因此，不同空間分布的土地利用類型及其變化對區(qū)域生態(tài)系統(tǒng)影響均不相同。目前已有相關(guān)學者開始關(guān)注這一研究方向，如Tanner等[13]、李子君等[14]、朱增云等[15]采用生態(tài)系統(tǒng)服務(wù)價值評估法、空間自相關(guān)分析和Getis-Ord熱點區(qū)分析等方法探討了土地利用演化及其生態(tài)系統(tǒng)響應(yīng)機制，從而能夠較好地從空間層面有效揭示兩者的影響機制。但在格網(wǎng)尺度下，空間演化過程、分異特征及其作用機理等方面的研究仍有待進一步挖掘。

鞏義市作為河南省鄭州都市區(qū)副中心城市，地形條件較為復(fù)雜，區(qū)域內(nèi)平原、山地、丘陵占比分別為1∶4∶5，在鄭州國家中心城市建設(shè)中承載著極重要的生態(tài)功能。自20世紀90年代以來的快速城鎮(zhèn)化進程擴張影響下，鞏義市土地利用變化程度較高，進而影響全域的生態(tài)系統(tǒng)價值變化，生態(tài)功能不斷退化[16]。因此，研究區(qū)域土地利用變化及其生態(tài)系統(tǒng)價值的時空響應(yīng)機制對于以鞏義市為代表的快速城鎮(zhèn)化地區(qū)國土資源合理化配置和生態(tài)環(huán)境保護具有極其重要的指導意義。本研究在基于鞏義市1990—2019年共4期的土地利用演化及其生態(tài)系統(tǒng)服務(wù)價值時空分布特征分析的基礎(chǔ)之上，采用轉(zhuǎn)移矩陣法、空間自相關(guān)分析等方法，從格網(wǎng)尺度下分析區(qū)域土地利用演化過程對生態(tài)系統(tǒng)價值的影響機制，以期為以鞏義市為代表的快速城鎮(zhèn)化地區(qū)生態(tài)系統(tǒng)保護和國土空間規(guī)劃提出切實有效的決策依據(jù)。

1 研究區(qū)概況

鞏義市位于河南省西部，在鄭州市和洛陽市交界地帶，市域總面積為1 043 km2。地勢呈現(xiàn)出東南高西北低的特點。東南部主要以高山區(qū)為主，中部則為淺山丘陵，北部為鄺嶺。氣候?qū)偌撅L暖溫帶氣候，年均氣溫11.5～22.3 ℃，年均降水量394.5 mm左右，常年降水量偏少，且時空分布差異明顯，富水區(qū)分布在沿黃河、伊洛河附近，而東南部山區(qū)為貧水區(qū)。全市土壤以潮土、褐土、棕壤土為主，在老廟山及周邊封山育林等低山丘陵區(qū)植被率高達50%～60%，占全市植物種類分布的80%。土地利用方面，農(nóng)業(yè)用地占比較高(62.72%)，其次為建設(shè)用地，占比20.52%，其他土地中，水域和自然保留地面積也占有一定比例，分別為3.94%和12.82%。截至2019年，鞏義市轄5個街道，15個鎮(zhèn)，全市人口83.83萬人，實現(xiàn)地區(qū)生產(chǎn)總值801.2億元。鞏義市憑借其重要的生態(tài)地資源優(yōu)勢和經(jīng)濟發(fā)展實力，經(jīng)過多年的城鎮(zhèn)化建設(shè)，2017年以來長期入選全國工業(yè)百強縣(市)、綜合競爭力全國百強縣(市)和全國制造業(yè)百強縣(市)等稱號，但是在取得上述發(fā)展成績的同時，鞏義市生態(tài)環(huán)境問題日益加劇。

2 研究數(shù)據(jù)與方法

2.1 數(shù)據(jù)來源及預(yù)處理

選取1990，2000，2010和2019年4期的鞏義市域6—9月的Landsat TM，ETM+，OLI遙感衛(wèi)星影像數(shù)據(jù)，數(shù)據(jù)均來源于美國地質(zhì)勘探局網(wǎng)站(https:∥glovis.usgs.gov/)。隨后基于ENVI5.0軟件平臺下，分別對各時段的影像數(shù)據(jù)進行融合、校正、裁切等數(shù)據(jù)的預(yù)處理工作。在進行土地利用分類時采用最大似然法監(jiān)督分類與google earth歷史影像方式相結(jié)合的方式，對4個年份研究區(qū)進行土地利用分類。分類方法上結(jié)合《土地利用現(xiàn)狀分類(GB/T21010-2017)》和謝高地等[17]構(gòu)建的單位面積生態(tài)系統(tǒng)服務(wù)價值分類體系，同時結(jié)合鞏義市土地資源分布現(xiàn)狀，將研究區(qū)土地利用分為耕地、喬木林、灌木林、高覆蓋度草地、低覆蓋度草地、水域和建設(shè)用地等7種類型(圖1)。隨后采用ENVI 5.0軟件中的精度評價模塊，對隨機采樣獲取的250個樣本點，結(jié)合4個時段土地利用解譯數(shù)據(jù)進行評價，進而得出各期分類結(jié)果的精度分別為87.13%，92.75%，90.37%和89.02%，可較好地滿足研究需要。此外，研究涉及數(shù)據(jù)中的DEM高程數(shù)據(jù)來源于地理空間數(shù)據(jù)云平臺(http:∥www.gscloud.cn/)。

2.2 研究方法

2.2.1 土地利用轉(zhuǎn)移分析 學界關(guān)于土地利用類型及面積轉(zhuǎn)移分析方法主要以土地利用轉(zhuǎn)移矩陣測度法為主，它多用來表征特定時段內(nèi)研究區(qū)域在研究期初至研究期末各土地利用類型之間的面積及其比例轉(zhuǎn)移情況，可有效揭示區(qū)域土地利用演化規(guī)律，計算方法如公式(1)所示[18]：

(1)

式中：S為各土地利用類型的總面積；n為區(qū)域土地利用類型的數(shù)量；i，j分別表示研究時段內(nèi)期初和期末的各土地利用類型。

2.2.2 基于格網(wǎng)尺度的數(shù)據(jù)處理方法 GIS格網(wǎng)分析法是指基于格網(wǎng)面狀單元類型，將其作為區(qū)域單元指標因子的數(shù)據(jù)評價分析單元，進而揭示區(qū)域空間內(nèi)部分異特征的方法[19]。在本研究中，基于格網(wǎng)單元對鞏義縣域的ESV空間分布及演化特征分析可在空間層面進一步細化研究尺度，從而使得其生態(tài)系統(tǒng)分析和優(yōu)化變得更易操作和精細化管理。由于縣域面積較大，研究基于ArcGIS 10.5軟件，運用Create Fishnet工具創(chuàng)建漁網(wǎng)樣本，為提高ESV時空分異精度，經(jīng)試驗得到最優(yōu)格網(wǎng)單元尺度為300 m×300 m，并最終得到覆蓋縣域范圍單元數(shù)量共計12 019個(邊界處格網(wǎng)面積稍有差異)。分別統(tǒng)計各格網(wǎng)內(nèi)生態(tài)系統(tǒng)服務(wù)價值，并將ESV值賦值至格網(wǎng)面矢量文件中，從而完成空間數(shù)據(jù)的重構(gòu)，進而得到基于格網(wǎng)尺度下鞏義市ESV分布及其自相關(guān)格局。格網(wǎng)化計算方法如公式(2)所示[19]：

H=(1/2)×min(Ai)1/2

(2)

式中：H為單元格網(wǎng)的邊長;Ai為格網(wǎng)單元內(nèi)最小圖斑的面積。

2.2.3 生態(tài)系統(tǒng)服務(wù)價值評估 研究基于Costanza等[17]對不同生物群區(qū)的貨幣化評估方法[20]以及謝高地等學者所建立的基于功能單元的單位面積生態(tài)系統(tǒng)服務(wù)價值評價表，結(jié)合鞏義市土地利用現(xiàn)狀實際(因分類后未利用地面積較少且調(diào)研后發(fā)現(xiàn)多以農(nóng)戶季節(jié)性種植菜地為主，故將其納入耕地類型)，以旱地—耕地、闊葉林—喬木林、灌草叢—灌木林、草地(平均值)—高覆蓋度草地、草地(平均值的0.7倍)—低覆蓋度草地、水系—水域及城市—建設(shè)用地等將各生物群區(qū)與研究區(qū)生態(tài)類型相對應(yīng)，由于本文不考慮人工生態(tài)系統(tǒng)的調(diào)節(jié)功能(相比縣域整體的生態(tài)系統(tǒng)作用可忽略不計)，因此將建設(shè)用地的生態(tài)系統(tǒng)服務(wù)價值默認為0。由于全國各省市自有生態(tài)系統(tǒng)差異較大，為避免直接使用謝高地等學者建立的單位面積生態(tài)系統(tǒng)服務(wù)價值當量因子法[20]產(chǎn)生的誤差，研究進一步參考了謝高地等[21]因地制宜修正后的各省份農(nóng)田生態(tài)系統(tǒng)生物量權(quán)重因子系數(shù)。筆者依據(jù)區(qū)域修正系數(shù)(河南省為1.39)和相關(guān)文獻中基于伏牛山等[22]同類型文獻研究成果選定的當量因子對鞏義市生態(tài)服務(wù)價值當量進行修正(表1)。測算方法如公式(3)所示[21]：

(3)

式中：ESV為生態(tài)系統(tǒng)服務(wù)價值;Pi為當量因子系數(shù);Si為第i種土地利用類型的面積;n為縣域所有土地利用類型的數(shù)量。

圖1 鞏義市1990-2019年土地利用類型分布

表1 鞏義市單位面積生態(tài)服務(wù)價值當量權(quán)重因子

2.2.4 空間探索性分析 空間探索性分析主要以空間全局自相關(guān)和局部自相關(guān)分析為主，主要通過某種現(xiàn)象的空間分布特征或規(guī)律進行可視化分析，進而探索空間內(nèi)部相關(guān)指標屬性的相似性或者差異性[23]。上述分析主要結(jié)合ArcGIS和GeoDa空間探索性分析工具來完成。

(1) 全局空間自相關(guān)。全局空間自相關(guān)以Global Moran’sI指數(shù)表示，可有效揭示鞏義市ESV空間相關(guān)性的總體趨勢。指數(shù)變化范圍為[-1,1]，若Moran’sI>0，則表示空間呈正相關(guān)，越接近1表示相關(guān)性越顯著，相似特征明顯；相反，若Moran’sI<0，則表示空間呈負相關(guān)，且越接近-1，則表示差異越大，分布越分散；若Moran’sI=0，則表示全域空間分布特征沒有相關(guān)性，呈隨機性分布。Global Moran’sI計算方法如公式(4)所示[23]：

(4)

(2) 局部空間自相關(guān)。局部空間自相關(guān)以Moran’sIi指數(shù)表示，可有效揭示鞏義市ESV空間相關(guān)性的局部分布特征。如果Moran’sIi指數(shù)>0，則表示局部空間單元內(nèi)部的ESV值存在高—高或者低—低局部聚集的特征；如果Moran’sIi指數(shù)<0，則表示局部空間單元內(nèi)部的ESV值存在高—低或者低—高局部聚集的特征。Moran’sIi計算方法如公式(5)所示[23]：

(5)

3 結(jié)果與分析

3.1 土地利用變化分析

基于GIS軟件對鞏義市1990—2019年4個時段內(nèi)各土地利用類型面積分別統(tǒng)計，可得到各時期內(nèi)土地類型面積統(tǒng)計(表2)。從土地利用結(jié)構(gòu)來看，鞏義市域最主要的用地類型為耕地，大約占鞏義市域面積的40%，其次為建設(shè)用地，占比約20.39%。從土地利用變化程度來看，以建設(shè)用地增幅最高(15.09%)，面積增加了2 965.23 hm2，其次為喬木林，共計增加了828.45 hm2；與此相反，研究區(qū)水域和灌木林則出現(xiàn)大幅減少，兩者降幅分別為-45.17%和-10.09%。究其主要原因，主要歸因于自20世紀90年代以來我國城鎮(zhèn)化進程持續(xù)加快，同時鞏義市工業(yè)經(jīng)濟增長迅速，進而導致區(qū)域城鎮(zhèn)擴張明顯，區(qū)域水域、灌木林地、耕地等不斷遭到侵占。此外，受市域北部黃河灘區(qū)不斷人為開發(fā)需要，大量水域空間轉(zhuǎn)變?yōu)楦?，這也是導致區(qū)域土地利用發(fā)生變化的主要原因。將1990，2019年土地利用矢量數(shù)據(jù)導入GIS軟件中，通過交集分析工具獲得1990—2019年鞏義市土地利用轉(zhuǎn)移矩陣(表3)。

表2 鞏義市1990-2019年土地利用及生態(tài)系統(tǒng)服務(wù)價值變化

表3 鞏義市1990-2019年土地利用轉(zhuǎn)移矩陣 hm2

由表3可知，鞏義市各土地利用類型轉(zhuǎn)移特征較為明顯，其中水域和耕地轉(zhuǎn)出面積最多，分別高達3 010.23和2 928.63 hm2，主要轉(zhuǎn)出類型分別為耕地和建設(shè)用地；灌木林轉(zhuǎn)向喬木林面積相對較高，轉(zhuǎn)移面積為523.16 hm2；其他用地類型轉(zhuǎn)入轉(zhuǎn)出程度相對較低。表明鞏義市擴張主要來自于周邊耕地類型，而耕地面積增長區(qū)域則來自于沿黃周邊的水域空間。同時受鞏義市植樹造林等生態(tài)保護政策影響，區(qū)域喬木林地呈現(xiàn)出逐年增長的趨勢。

3.2 生態(tài)系統(tǒng)服務(wù)價值分析

3.2.1 生態(tài)系統(tǒng)服務(wù)價值測度 結(jié)合單位面積生態(tài)服務(wù)價值當量權(quán)重因子和各時段土地利用面積數(shù)據(jù)，進而得出鞏義市各土地類型和不同生態(tài)功能的生態(tài)系統(tǒng)服務(wù)價值統(tǒng)計及變化情況(表2，4)。由表2和表4可知，鞏義市生態(tài)系統(tǒng)服務(wù)價值整體呈現(xiàn)出波動式下跌的趨勢，近30 a間共計減少了34.85億元(表2)。結(jié)合表2中不同土地利用類型的ESV來看，水域價值為各土地類型中最高，占比高達42.22%；其次為喬木林和耕地，分別占總ESV的32.01%和23.34%；其他用地類型的ESV則均低于9億元。從研究期間的變化特征來看，不同土地利用類型的ESV表現(xiàn)各異。其中，喬木林和高覆蓋度草地的ESV出現(xiàn)小幅增長的趨勢，ESV分別增加了1.901和0.116億元；而水域、灌木林、耕地和低覆蓋度草地則出現(xiàn)了不同程度的下降趨勢，尤以水域ESV下降最為顯著，共計減少了35.549億元，其減幅為各用地類型中最高，其他用地類型ESV減幅均在1億元以內(nèi)。

從表4所示不同服務(wù)功能下的生態(tài)價值分布來看，水文調(diào)節(jié)服務(wù)價值最高(年均58.51億元)，ESV占總價值量的43.29%；其次為氣候調(diào)節(jié)和水土保持服務(wù)，分別占13.69%和9.63%。其他類型ESV則占比較小，且ESV均低于10億元。從各時段變化特征來看，研究期內(nèi)所有功能ESV均處于持續(xù)下降趨勢，尤以水資源供給和水文調(diào)節(jié)降幅最為明顯，分別為-38.2%和-37.572%，兩者分別減少了-2.333億元和-28.764億元；其次為凈化環(huán)境和美學景觀功能，兩者降幅分別為-17.454%和-12.267%?？傮w而言，以水文調(diào)節(jié)和氣候調(diào)節(jié)為主的調(diào)節(jié)服務(wù)在鞏義市生態(tài)系統(tǒng)服務(wù)價值構(gòu)成中占據(jù)主導地位，同時其近30 a來的損失率也處于各功能中較高位置。因此，在未來相當長一段時期內(nèi)應(yīng)當嚴格維持區(qū)域生態(tài)調(diào)節(jié)服務(wù)，以避免快速城鎮(zhèn)擴張和人為活動對區(qū)域生態(tài)系統(tǒng)造成負面影響。

表4 1990-2019年鞏義市不同生態(tài)功能的生態(tài)服務(wù)價值變化 億元

3.2.2 生態(tài)系統(tǒng)服務(wù)價值時空演化特征 基于ArcGIS 10.5軟件中的創(chuàng)建漁網(wǎng)功能，通過對各網(wǎng)格下土地利用類型分布進行生態(tài)系統(tǒng)權(quán)重因子賦值，進而得出格網(wǎng)尺度下鞏義市生態(tài)系統(tǒng)服務(wù)價值，采用Jenks最佳自然斷裂法對區(qū)域ESV空間變化進行分類，共計劃分為7個等級，分別為：0～213萬元，213～470萬元，470～1001萬元，1 001～1868萬元，1 868～4 440萬元，4 440～9 542萬元，9 542～12 600萬元(圖2)。從空間分布及變化特征來看，1990年鞏義市ESV較高區(qū)域主要位于市域北部黃河、伊洛河沿線以及南部伏羲山等山脈，而ESV較低區(qū)域則主要以鞏義市市區(qū)及東西地形較為平坦區(qū)域一帶；進入2000年以后，鞏義市ESV高值片區(qū)開始出現(xiàn)快速減少趨勢，且多集中在黃河沿線，多歸因于黃河濕地持續(xù)被改造為耕地導致ESV持續(xù)減少，此時低值ESV區(qū)域則整體變化不大；而在2000—2010年，位于鞏義市主城區(qū)的低值ESV區(qū)域受快速城鎮(zhèn)擴張影響，開始出現(xiàn)蔓延趨勢，且該趨勢在2010年以后不斷加劇，除圈層式蔓延外，在鞏義市區(qū)至西向的回郭鎮(zhèn)之間已形成帶狀ESV低值集聚區(qū)，該片區(qū)同時也是鄭州—洛陽工業(yè)走廊中心地帶，城鎮(zhèn)化進程相對較高。從1990—2019年ESV整體分布特征來看，鞏義市ESV近30 a出現(xiàn)了較大幅度的下降，區(qū)域城鎮(zhèn)化進程對生態(tài)系統(tǒng)影響的負面效應(yīng)較為顯著，且ESV變化的分布基本上與水域、林地和建設(shè)用地的變化相吻合，鞏義市北部以水域、耕地和建設(shè)用地為主的ESV高值區(qū)和低值區(qū)變化較明顯，南部以林地和高覆蓋度草地為主的生態(tài)資源集中區(qū)則變化相對穩(wěn)定。

圖2 鞏義市1990-2019年生態(tài)服務(wù)價值(ESV)空間演化格局

為深入探索格網(wǎng)單元下鞏義市ESV的內(nèi)在演化機制和分布特征，研究基于圖2中的ESV柵格數(shù)據(jù)進行矢量化處理，隨后導入GeoDa1.14空間探索性分析軟件中，對其全局Moran’sI指數(shù)、散點圖以及LISA聚集特征進行分析。全局Moran’sI指數(shù)基于Queen鄰接原則(秩設(shè)定為1)將其空間自相關(guān)特征進行定量化分析，并最終得到1990—2019年4期鞏義市全局Moran’sI指數(shù)分別為0.881，0.723，0.851和0.741，表明鞏義市ESV存在較為顯著的空間自相關(guān)特征。結(jié)合圖3中1990—2019年鞏義市ESV的LISA空間集聚圖可明顯發(fā)現(xiàn)，ESV在空間集聚分布上近30 a整體處于相對穩(wěn)定的格局，呈現(xiàn)出北部ESV高值聚集而中部及西南側(cè)低值聚集的特征，上述兩種類型的空間聚集效應(yīng)較為明顯，而高—低和低—高類型的聚集效應(yīng)則并不突出。

(1) 高—高ESV等級聚集區(qū)。鞏義市ESV處于高—高聚集區(qū)范圍主要以北部黃河、伊洛河流域為主，同時隨時間推移，東南部青龍山、九龍峽、浮戲山等山脈開始出現(xiàn)高—高聚集區(qū)。結(jié)合用地類型來看，上述區(qū)域多集中分布在生態(tài)系統(tǒng)保存較好的片區(qū)，主要以水域、林地等類型為主，同時也是鞏義市ESV的主要來源地。

(2) 低—低ESV等級聚集區(qū)。鞏義市低—低ESV等級集聚區(qū)域多分布在鞏義市主城區(qū)以及鞏義—芝田鎮(zhèn)—回郭鎮(zhèn)、鞏義—竹林鎮(zhèn)—米河鎮(zhèn)等西向和東向兩條帶狀片區(qū)，且隨著時間變化呈現(xiàn)出不斷蔓延的趨勢，尤以鞏義市區(qū)南向蔓延特征相對明顯。其他分散分布的低—低ESV集聚區(qū)則多以縣域各鄉(xiāng)鎮(zhèn)為主。從地形特征來看，鞏義市區(qū)及其東西向沿線工業(yè)走廊涉及鄉(xiāng)鎮(zhèn)多以平原為主，且用地類型主要為建設(shè)用地。鞏義市在近年來生態(tài)保護政策取得一定成效的同時，在低—低ESV聚集區(qū)域受城鎮(zhèn)持續(xù)擴張影響，未來如何優(yōu)化該片區(qū)生態(tài)系統(tǒng)仍值得重點關(guān)注。

(3) 高—低和低—高ESV等級聚集區(qū)。鞏義市ESV處于高—低和低—高等級的聚集區(qū)域整體分布相對較少。其中，低—高等級聚集區(qū)多集中分布在北部的黃河流域附近，結(jié)合Moran散點圖可知，處于第一象限的格網(wǎng)分散化趨勢顯著，而處于第4象限的高—低等級聚集區(qū)格網(wǎng)數(shù)則年均不超過4個，且零星分布在南部山區(qū)。

圖3 鞏義市1990-2019年LISA空間集聚格局

4 討論與結(jié)論

4.1 討 論

自20世紀90年代以來，涉及土地利用演化及其生態(tài)系統(tǒng)服務(wù)價值的影響研究已有大量實證研究成果，但從空間視角探討其影響特征及其驅(qū)動機制仍未出現(xiàn)較為完善的研究方法[24]。本研究以工業(yè)型快速城鎮(zhèn)化為代表的河南省鞏義市為研究對象，在分析區(qū)域土地利用及其ESV損益的基礎(chǔ)上，通過格網(wǎng)尺度從微觀視角下探討土地利用對生態(tài)系統(tǒng)服務(wù)價值的空間影響特征。研究結(jié)果表明，以耕地為主的鞏義市中部地區(qū)應(yīng)在基本農(nóng)田得到保障的基礎(chǔ)之上，合理控制區(qū)域工業(yè)化進程導致的建設(shè)用地擴張態(tài)勢；在市域南部尤其是東南山地區(qū)域，應(yīng)制定嚴格的生態(tài)保護措施，有效保障該片區(qū)生態(tài)系統(tǒng)價值穩(wěn)定發(fā)展。區(qū)域ESV的演變不僅受土地利用和自然地形影響，同時還受到氣候條件、產(chǎn)業(yè)分布、經(jīng)濟水平、人口密度等諸多自然和社會經(jīng)濟因素的影響。因此，在未來研究中可重點關(guān)注： ①基于地理探測器方法下研究大數(shù)據(jù)驅(qū)動下的區(qū)域生態(tài)系統(tǒng)服務(wù)價值空間分布及預(yù)測模擬研究； ②行政區(qū)劃范圍內(nèi)的ESV演化機制通常情況下還會受到周邊城市的影響，因此在研究范圍的確定上應(yīng)更關(guān)注特定經(jīng)濟影響區(qū)域范圍下的ESV空間特征更具有說服力。

受持續(xù)城鎮(zhèn)擴張和工業(yè)化進程影響，鞏義市土地利用變化直接作用于區(qū)域生態(tài)系統(tǒng)。故可通過相應(yīng)的規(guī)劃調(diào)控等相關(guān)措施不斷優(yōu)化區(qū)域土地資源配置進而維護生態(tài)平衡。相關(guān)政策建議包括： ①鑒于林地與水域生態(tài)用地自身價值較高，因而在未來生態(tài)環(huán)境修復(fù)與保護過程中應(yīng)重點加強縣域內(nèi)黃河流域與伏羲山等生態(tài)重點保護區(qū)的生態(tài)紅線控制與保護； ②針對現(xiàn)有的工礦產(chǎn)業(yè)布局與開發(fā)現(xiàn)狀，應(yīng)加快開展綠色礦山生態(tài)修復(fù)工程，優(yōu)先治理因工礦產(chǎn)業(yè)開發(fā)引起的生態(tài)環(huán)境污染問題； ③縣域內(nèi)各鄉(xiāng)鎮(zhèn)應(yīng)統(tǒng)一建立起建設(shè)用地動態(tài)監(jiān)測與預(yù)警機制，同時在保護耕地紅線的基礎(chǔ)上開展退耕還林政策，以期進一步提升區(qū)域生態(tài)系統(tǒng)價值。

4.2 結(jié) 論

(1) 研究期內(nèi)工業(yè)城鎮(zhèn)化進程持續(xù)推進，不斷促使建設(shè)用地快速擴張(增加2 965.23 hm2)，與此相反，水域和灌木林等被建設(shè)用地不斷吞噬而出現(xiàn)大幅減少，兩者降幅分別為-45.17%和-10.09%。用地類型轉(zhuǎn)移方面，水域和耕地轉(zhuǎn)出面積最多，主要轉(zhuǎn)出類型以耕地和建設(shè)用地為主。

(2) 鞏義市生態(tài)系統(tǒng)服務(wù)價值整體呈波動式下跌，近30 a共減少了34.85億元。從變化特征來看，喬木林和高覆蓋度草地小幅增長，而水域、灌木林、耕地和低覆蓋度草地則出現(xiàn)不同程度下降，尤以水域ESV跌幅明顯。從不同服務(wù)功能的生態(tài)價值來看，水文調(diào)節(jié)服務(wù)價值最高，其次為氣候調(diào)節(jié)和水土保持服務(wù)。以水文調(diào)節(jié)和氣候調(diào)節(jié)為主的調(diào)節(jié)服務(wù)在鞏義市生態(tài)系統(tǒng)服務(wù)價值構(gòu)成中占據(jù)主導地位，同時其近30 a來的損失率也處于高位。

(3) 鞏義市城鎮(zhèn)化進程對生態(tài)系統(tǒng)影響的負面效應(yīng)顯著，北部以水域、耕地和建設(shè)用地為主的ESV高值區(qū)和低值區(qū)變化較大，南部以林地和高覆蓋度草地為主的生態(tài)資源集中區(qū)由于受社會經(jīng)濟活動影響較少則相對穩(wěn)定。1990年ESV較高區(qū)域集中在北部黃河、伊洛河沿線及南部伏羲山等山脈，而ESV較低區(qū)域則分布在鞏義市區(qū)及東西兩側(cè)；2000年以后，主城區(qū)低值ESV區(qū)域受快速城鎮(zhèn)擴張影響逐步蔓延，在城區(qū)及鞏義—芝田鎮(zhèn)—回郭鎮(zhèn)、鞏義—竹林鎮(zhèn)—米河鎮(zhèn)等西向和東向兩條帶狀片區(qū)已形成ESV低值集聚區(qū)。