廈漳泉地區(qū)生態(tài)系統(tǒng)服務(wù)價值時空分布格局演變

胡其玉, 陳松林

(1.福建師范大學(xué) 福建省亞熱帶資源與環(huán)境重點實驗室, 福州350007; 2.福建師范大學(xué) 地理研究所, 福州350007)

當(dāng)今中國的城鎮(zhèn)化進程尚處于加速階段，協(xié)調(diào)生態(tài)和社會經(jīng)濟系統(tǒng)間矛盾是當(dāng)前亟需解決的問題。為實現(xiàn)區(qū)域可持續(xù)發(fā)展，通過探索生態(tài)系統(tǒng)服務(wù)功能變化了解生態(tài)環(huán)境狀況為分區(qū)保育策略的制定提供參考，具有現(xiàn)實意義。生態(tài)系統(tǒng)服務(wù)功能[1]即指生態(tài)系統(tǒng)提供人類生產(chǎn)生活所需的各種自然環(huán)境條件及效用的能力，包括糧食、醫(yī)藥等物質(zhì)生活產(chǎn)品以及水土保持等功能，是人類生存和經(jīng)濟社會發(fā)展的基礎(chǔ)。然而，隨著人類活動對水域、草地、林地等生態(tài)源地的占用，擾動了生態(tài)系統(tǒng)的穩(wěn)定性，造成生態(tài)服務(wù)功能受損。目前，從價值量化角度測算土地利用變化引發(fā)的生態(tài)服務(wù)功能效應(yīng)是建立區(qū)域生態(tài)補償機制的關(guān)鍵所在[2-3]。謝高地等基于中國地區(qū)實際特點修正Costanza等提出的生態(tài)系統(tǒng)服務(wù)價值評估體系，建立適用于中國的陸地生態(tài)服務(wù)價值當(dāng)量表并提出時空動態(tài)修正模型，因該方法適應(yīng)性強、便于測算的特性被普遍用于定量測度不同區(qū)域不同土地利用類型的生態(tài)系統(tǒng)服務(wù)價值。隨著3S技術(shù)的發(fā)展和應(yīng)用，有關(guān)生態(tài)系統(tǒng)服務(wù)價值時空分布的研究愈加受到學(xué)者們的重視[4]，應(yīng)用空間制圖定量刻畫特定時空尺度范圍生態(tài)服務(wù)價值的分布格局成為熱點。而地形作為土地資源構(gòu)成要素之一，也影響著生態(tài)系統(tǒng)服務(wù)的垂直分異特征，通過探索生態(tài)系統(tǒng)服務(wù)多尺度變化和地形梯度效應(yīng)[5]，有利于進一步識別生態(tài)系統(tǒng)服務(wù)價值空間分異特點，推進生態(tài)系統(tǒng)分區(qū)管理。

目前有關(guān)經(jīng)濟社會快速轉(zhuǎn)型的城市群地區(qū)的研究日益突出[6-7]，但對中小城市群的關(guān)注尚有不足。廈漳泉作為典型中小城市群地區(qū)，近10 a來經(jīng)濟社會飛速發(fā)展、人地比例擴大加重了資源承載壓力；另一方面，境內(nèi)平原、丘陵、山地交錯，具有明顯的地形起伏特征。本文采用當(dāng)量因子法評估廈漳泉地區(qū)生態(tài)系統(tǒng)服務(wù)價值，以格網(wǎng)單元的生態(tài)系統(tǒng)服務(wù)價值作為反映地表覆被環(huán)境的指標(biāo)，識別2006—2016年地區(qū)生態(tài)系統(tǒng)服務(wù)價值時空分布格局及地形梯度效應(yīng)，研究成果可為中小城市群地區(qū)維護生態(tài)服務(wù)功能平衡和建立生態(tài)補償機制提供依據(jù)。

1 研究區(qū)概況及數(shù)據(jù)來源

1.1 研究區(qū)概況

廈漳泉地區(qū)位于福建省東南沿海，東臨臺灣海峽，南靠廣東潮汕地區(qū)，處于23°48′—25°56′N和117°—119°05′E，包含廈門、漳州和泉州(不含金門縣)三市，被稱為“閩南金三角”。屬于亞熱帶季風(fēng)氣候，水熱資源充沛，地區(qū)平均森林覆蓋率約為61.22%。地勢總體由西北山地向東南濱海地帶傾斜，地貌構(gòu)成復(fù)雜多樣，山地、丘陵和平原錯落其中。境內(nèi)溪流縱橫，主要有晉江、九龍江、漳江等水系。廈漳泉地區(qū)陸地面積約占福建省陸地面積的20.4%，常住人口占比為45.28%，GDP占全省地區(qū)生產(chǎn)總值的47.87%。

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

2006年、2016年廈漳泉地區(qū)Landsat TM和Landsat8遙感影像資料以及GDEMDEM數(shù)據(jù)(30 m×30 m)源自地理空間數(shù)據(jù)云網(wǎng)站。將中國生態(tài)系統(tǒng)評估與生態(tài)安全數(shù)據(jù)庫2010年福建省土地利用數(shù)據(jù)作為參考,利用ENVI 5.1預(yù)處理后比照Google地圖對遙感影像進行人機交互判讀，根據(jù)中科院土地分類系統(tǒng)結(jié)合區(qū)域特點采用最大似然法解譯出地區(qū)土地利用類型并經(jīng)反復(fù)糾正修改,劃分為林地、草地、耕地、濕地、未利用地和建設(shè)用地(表1)。2006年、2016年全國和廈門、漳州、泉州年均降水量獲取自中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)，各類農(nóng)作物播種面積、單位面積產(chǎn)量分別來源于《廈門經(jīng)濟特區(qū)年鑒2017》、《漳州統(tǒng)計年鑒2017》和《泉州統(tǒng)計年鑒2017》，稻谷、蔬菜、甘薯等作物平均價格依據(jù)福建省糧食局、物價局提供的《價格檢測報告》而確定。

表1 廈漳泉地區(qū)土地利用變化情況 %

2 研究方法

2.1 標(biāo)準(zhǔn)單位生態(tài)系統(tǒng)服務(wù)價值當(dāng)量因子價值量測算

結(jié)合廈漳泉地區(qū)實際，采用1 hm2稻谷、甘薯、花生和蔬菜單位面積的平均價格評估地區(qū)標(biāo)準(zhǔn)生態(tài)系統(tǒng)類型(即1個標(biāo)準(zhǔn)當(dāng)量因子)的價值量。并將2016年地區(qū)4種農(nóng)作物平均價格的1/7作為標(biāo)準(zhǔn)單位生態(tài)服務(wù)價值當(dāng)量因子的經(jīng)濟價值量(公式1)，從而排除人工投入的影響[8-9]，確定2006—2016年標(biāo)準(zhǔn)生態(tài)系統(tǒng)的單位面積生態(tài)服務(wù)功能價值量為5 871.722元/hm2。

(1)

式中：D表示標(biāo)準(zhǔn)生態(tài)系統(tǒng)的單位面積生態(tài)服務(wù)功能價值量(元/hm2)；W1,W2,W3,W4分別為2016年廈漳泉地區(qū)稻谷、甘薯、花生和蔬菜的單位面積平均價格(元/hm2)；S1,S2,S3,S4分別為稻谷、甘薯、花生和蔬菜的播種面積與其播種面積總和之比(%)。

2.2 單位面積生態(tài)系統(tǒng)服務(wù)價值當(dāng)量表修正與價值評估

根據(jù)已有研究成果[10-11]，結(jié)合廈漳泉地區(qū)實際情況，采用地區(qū)平均年均降水量因子修正單位面積生態(tài)服務(wù)價值當(dāng)量表，修正系數(shù)為1.901。對解譯后的土地生態(tài)類型進行合理歸并[12-13]，通過修正后的當(dāng)量因子乘以標(biāo)準(zhǔn)生態(tài)系統(tǒng)的單位面積生態(tài)服務(wù)價值量(D)[14]，得到不同土地生態(tài)類型的價值系數(shù)表(表2)。在此基礎(chǔ)上，根據(jù)公式(2) 計算各類型生態(tài)系統(tǒng)服務(wù)價值及總價值。

(2)

式中：在計算地區(qū)總的ESV(生態(tài)系統(tǒng)服務(wù)價值)時，Ai表示第i種土地生態(tài)類型的面積(hm2)；VCi表示第i種土地生態(tài)類型的單位面積生態(tài)服務(wù)價值(元/hm2)；在計算各種土地生態(tài)類型的ESVi時，Ai表示第i種土地生態(tài)類型的面積；VCk表示第k項的單位面積生態(tài)系服務(wù)價值。

表2 2006-2016年廈漳泉地區(qū)生態(tài)服務(wù)價值系數(shù) 元/hm2

2.3 空間自相關(guān)分析

根據(jù)研究區(qū)土地利用類型圖斑大小，在0.1×0.1～10×10 km2逐步增大間距劃定網(wǎng)格，選取最優(yōu)大小的格網(wǎng)單元展開分析。空間自相關(guān)分析分為全局空間自相關(guān)和局部自相關(guān)[15-16]。前者(公式3)可以測度整個地區(qū)生態(tài)服務(wù)價值的空間集聚性及其特征；后者(公式4)主要是衡量地區(qū)內(nèi)部各個單元間變量的空間集聚性與空間異質(zhì)性，判別“高低”集聚模式。通過公式(5) 可對空間自相關(guān)分析結(jié)果進行顯著性檢驗。

(3)

(4)

Moran′sI的計算結(jié)果>0時表示格網(wǎng)單元具有高—高或者低—低的空間集聚特征；Moran′sI的計算結(jié)果<0時表示格網(wǎng)單元具有低—高或者高—低的空間集聚特征。

(5)

式中：Zi表示空間自相關(guān)的顯著水平；E(I)表示Moran′sI指數(shù)的數(shù)學(xué)期望；V(I)為其方差。

2.4 地形起伏度計算方法

參考前人研究成果[17-18]，采用經(jīng)修訂后的中國人居自然環(huán)境適宜性評價中的算法(式6)測度廈漳泉地區(qū)的地形起伏度，刻畫生態(tài)服務(wù)價值地形梯度效應(yīng)[19]。通過ArcGIS 10.2平臺空間分析模塊中(Zonal Statistics as Table)和(Raster Reclassify)功能實現(xiàn)矢量數(shù)據(jù)和柵格數(shù)據(jù)的疊加分析，基于DEM數(shù)據(jù)提取每個格網(wǎng)的地形起伏度，算得地區(qū)格網(wǎng)單元地形起伏度均值為0.975。

(6)

式中：RDLS表示地形起伏度；ALT為格網(wǎng)單元內(nèi)平均海拔(m)；max(H)、min(H)分別為格網(wǎng)單元內(nèi)最高和最低海拔(m)；P(A)表示格網(wǎng)單元內(nèi)坡度小于等于5°的平地面積(km2)；A表示格網(wǎng)單元的總面積(km2)。

3 結(jié)果與分析

3.1 2006-2016年廈漳泉地區(qū)生態(tài)服務(wù)服務(wù)價值分布總體態(tài)勢

廈漳泉地區(qū)的生態(tài)系統(tǒng)服務(wù)價值測算結(jié)果(表3)顯示，2006—2016年地區(qū)生態(tài)服務(wù)價值總體提升16.51%，提升幅度較大，說明地區(qū)生態(tài)服務(wù)功能有較大程度改善，主要體現(xiàn)在林地、濕地、水域等生態(tài)源地服務(wù)價值的提升。

運用全局趨勢分析廈漳泉地區(qū)生態(tài)服務(wù)價值的總體分布情況(圖1)，圖中y軸代表正北方向，x軸代表正東方向，z軸代表單元生態(tài)服務(wù)價值大小。結(jié)果顯示，x軸(東西方向)呈現(xiàn)下降曲線，y軸(南北方向)呈現(xiàn)先降后升的凹型曲線且2016年顯于2006年，表明地區(qū)生態(tài)服務(wù)價值呈西高東低和南北高中間低的分布格局并且10 a間趨勢效應(yīng)有所增強?？傮w上，廈漳泉地區(qū)生態(tài)服務(wù)價值分布與地形山勢相適應(yīng)，并呈現(xiàn)出自北、西、南向中東部傾斜的趨勢，受制于自然環(huán)境條件和經(jīng)濟發(fā)展引發(fā)的生態(tài)壓力，生態(tài)服務(wù)價值分布于廈漳泉地區(qū)的中東部形成“低谷”，2006—2016年“低谷”效應(yīng)愈加顯著。

表3 2006-2016年各類生態(tài)系統(tǒng)服務(wù)價值變化情況 萬元

3.2 廈漳泉地區(qū)生態(tài)服務(wù)服務(wù)價值空間分布尺度效應(yīng)

以1 km×1 km和2 km×2 km正方形格網(wǎng)為單元的生態(tài)服務(wù)價值空間集聚的特性因粒度太小而并不顯著，進而生成3 km×3 km，4 km×4 km，5 km×5 km 、6 km×6 km(尺度1—4)并將9 km×9 km(尺度5) 格網(wǎng)單元作為參照進行全局空間自相關(guān)分析。劃分等間隔距離對不同尺度下格網(wǎng)單元生態(tài)服務(wù)價值進行空間增量自相關(guān)分析，得到2006年(A)和2016年(B)在不同尺度下的Moran′sI指數(shù)(圖2)，其值大于0。經(jīng)驗Zi>0，p<0.05，表現(xiàn)出正向空間自相關(guān)性，說明兩個時期地區(qū)生態(tài)服務(wù)價值在空間分布上均表現(xiàn)出集聚的特點。并且受間隔距離的制約，生態(tài)服務(wù)價值的集聚程度隨研究尺度增大而呈現(xiàn)波動上升的尺度效應(yīng)。在各個尺度下，2016年廈漳泉地區(qū)生態(tài)服務(wù)價值的Moran′sI指數(shù)基本高于2006年，空間自相關(guān)性更強，集聚程度更高。

注：AB字母后1—5代表不同尺度。

圖2 廈漳泉地區(qū)生態(tài)服務(wù)價值空間自相關(guān)尺度效應(yīng)

3.3 廈漳泉地區(qū)生態(tài)服務(wù)價值時空分異格局特征

全局Moran′sI只是反映出廈漳泉地區(qū)整體的空間自相關(guān)水平和總體趨勢，因而需要通過局部自相關(guān)挖掘地區(qū)內(nèi)部生態(tài)服務(wù)價值的空間分布特點。將3 km×3 km作為主要研究尺度，并以6 km×6 km尺度的研究結(jié)果作為參照。

(1) 將2006年和2016年地區(qū)生態(tài)服務(wù)價值格網(wǎng)數(shù)據(jù)導(dǎo)入GeoDa軟件處理，經(jīng)局部自相關(guān)分析生成廈漳泉地區(qū)生態(tài)服務(wù)價值的Moran指數(shù)分析，通過標(biāo)準(zhǔn)化空間權(quán)重矩陣建立相互聯(lián)系[20]。結(jié)果顯示，兩期的Moran′sI指數(shù)分別為0.107，0.157，且通過顯著性檢驗。觀察圖3并經(jīng)統(tǒng)計可以發(fā)現(xiàn)，2006年和2016年生態(tài)服務(wù)價值散點均集中分布在第一、三象限，其余分布于第二、四象限，表明廈漳泉地區(qū)生態(tài)服務(wù)價值呈現(xiàn)正向空間自相關(guān)性，生態(tài)服務(wù)價值集聚區(qū)域主要集中于高—高和低—低區(qū)域，內(nèi)部差異較小。2006—2016年生態(tài)服務(wù)價值高—高集聚的優(yōu)勢更加明顯，說明地區(qū)生態(tài)服務(wù)功能有所提升和改善。

(2) 通過局部自相關(guān)集聚圖(圖3)，探索2006—2016年廈漳泉地區(qū)生態(tài)服務(wù)價值分布的空間集聚性和空間異質(zhì)性。無論是高—高還是低—低集聚型區(qū)域，p值遠低于0.05逼近于0，說明局部自相關(guān)性顯著。圖3A和圖3B分別為2006年、2016年生態(tài)服務(wù)價值的集聚圖，圖3C表示10 a間單元生態(tài)服務(wù)價值的增減集聚特征。結(jié)果顯示，總體上，生態(tài)服務(wù)價值高、低集聚區(qū)域與上述趨勢面分析結(jié)果一致：生態(tài)服務(wù)價值分布在南北方向呈現(xiàn)出“兩高一低”的格局，北部高生態(tài)服務(wù)價值區(qū)域主要位于德化縣大部和永春縣東部，南部高生態(tài)服務(wù)價值區(qū)域基本位于漳州市西側(cè)，中部低生態(tài)服務(wù)價值區(qū)域主要集中于安溪縣境內(nèi)與以泉州平原和廈門島為中心的沿海一帶；東西方向呈現(xiàn)出“西高東低”分布格局，西部山地、丘陵的生態(tài)服務(wù)價值普遍高于東部平原。2006—2016年廈漳泉地區(qū)生態(tài)服務(wù)價值的增長趨向南北兩極，下降區(qū)域集中于漳州平原和泉州平原，空間格局演變較為顯著，逐漸呈斜三角形分布趨勢。

圖3 廈漳泉地區(qū)生態(tài)服務(wù)價值集聚特征

北部生態(tài)服務(wù)價值高—高集聚區(qū)域山水格局優(yōu)勢明顯，區(qū)域內(nèi)以旅游與農(nóng)業(yè)經(jīng)營為主，年輕勞動力向沿海遷移，土地開發(fā)強度較弱與城鎮(zhèn)化擴張程度較低。南部生態(tài)服務(wù)價值高—高集聚區(qū)域2006年分布于由九龍江北溪、西溪經(jīng)漳州平原延伸至九龍江入?？谝粠?，10 a間由于薌城、龍文、龍?？焖俪擎?zhèn)化的推進，在臨近北溪與西溪的交匯處轉(zhuǎn)向西南方，形成廈漳泉地區(qū)體量最大的生態(tài)服務(wù)功能供給系統(tǒng)。根據(jù)中部生態(tài)服務(wù)價值低—低集聚型區(qū)域特點將其劃分為3個子區(qū)域：一是安溪境內(nèi)水土流失區(qū)，主要分布于安溪縣中西部，山巒起伏，山峰林立，山勢峻峭，河谷狹窄，在此基礎(chǔ)上發(fā)展規(guī)模龐大的茶種植產(chǎn)業(yè)易引發(fā)嚴(yán)重的水土流失，經(jīng)過多年的治理，結(jié)果顯示2006—2016年生態(tài)服務(wù)價值低—低集聚區(qū)域受到限制并縮減，說明水土流失對生態(tài)系統(tǒng)服務(wù)功能的擾動得到改善。二是漳浦縣沿海風(fēng)積地貌區(qū)較為發(fā)育，在海流與風(fēng)力的作用下，逐漸形成沙堤、沙洲，因缺少內(nèi)陸大面積森林等生態(tài)源地的支撐，持續(xù)表現(xiàn)出顯著的低—低集聚狀態(tài)。三是自泉州泉港區(qū)至漳州龍海市沿海一帶的密集城鎮(zhèn)化區(qū)，是廈漳泉地區(qū)面積最大的生態(tài)服務(wù)價值低—低集聚型區(qū)域，也是廈漳泉地區(qū)乃至福建省的經(jīng)濟、人口重心所在，由于經(jīng)濟開發(fā)區(qū)建設(shè)和人口密度增加，10 a間低—低集聚的格局以泉州平原和廈門島為中心向西南—東北瀕海拓展。生態(tài)服務(wù)價值高—低集聚區(qū)域10 a間分布愈加集中，主要在九龍江河口、廈門島、泉州灣以及一些丘陵、平原交錯沿海地帶，應(yīng)注意維護濱海水域和濕地生態(tài)功能的完整性；而低—高集聚區(qū)域與高—低集聚區(qū)域的范圍較小。另外，隨格網(wǎng)尺度的增大使得安溪地區(qū)不再表現(xiàn)出低—低集聚，低—低集聚區(qū)域主要集中于沿海地帶，快速城鎮(zhèn)化區(qū)域生態(tài)服務(wù)價值的損失尤為顯著。

3.4 廈漳泉地區(qū)生態(tài)服務(wù)價值與地形起伏度的相關(guān)分析

應(yīng)用自然斷點法將廈漳泉地區(qū)地形起伏度由低到高劃分為5個等級，并統(tǒng)計各等級生態(tài)服務(wù)價值的均值(表4)[21]，結(jié)果顯示，隨地形起伏度的增加，生態(tài)服務(wù)價值平均水平逐漸提高，但相比于2006年、2016年從第Ⅰ級到第Ⅱ級的增長幅度更為劇烈，并在第Ⅳ，Ⅴ級有所下降，出現(xiàn)波動情況。說明2006—2016年各級別地形起伏度下的生態(tài)服務(wù)功能有所改善，尤其是山地、丘陵等地形起伏較大的區(qū)域，但是地形起伏較小的平原、沿海一帶生態(tài)服務(wù)功能改善趨勢不甚顯著，并且地區(qū)內(nèi)不同地形起伏度下的生態(tài)服務(wù)價值水平表現(xiàn)出分化趨勢。雖然總體上濕地、水域的維護有所成效，但受快速城鎮(zhèn)化擴張的影響，沿海平原的生態(tài)服務(wù)功能仍面臨被擾動的威脅。

通過空間制圖定量刻畫廈漳泉地區(qū)生態(tài)服務(wù)價值和地形起伏度特征，顯示生態(tài)服務(wù)價值與地形起伏度之間存在某種相關(guān)性，參考已有研究[22]并以3 km×3 km格網(wǎng)作為單元，應(yīng)用SPSS軟件生成兩個變量的散點圖并作回歸分析。為完整反映兩者相關(guān)特征，不作剔除異常值處理，而是以生態(tài)服務(wù)價值平均水平為界采用分段回歸分析，進一步分析和驗證其分布特征。當(dāng)?shù)匦纹鸱忍幱?～1時，地區(qū)生態(tài)服務(wù)價值變化劇烈，同時呈現(xiàn)陡增和陡降的趨勢。地形起伏度較小的地區(qū)主要分布在沿海一帶，包含低生態(tài)服務(wù)價值的城鎮(zhèn)建設(shè)區(qū)、沿海風(fēng)沙侵蝕地帶，高生態(tài)服務(wù)價值的園林地、水域與濕地。隨著地形起伏度的增加，生態(tài)服務(wù)價值高于平均水平的濕地、水域逐漸減少而林地向山地、丘陵集中，生態(tài)服務(wù)價值分布漸趨穩(wěn)定。當(dāng)?shù)匦纹鸱忍幱?～2時，單元生態(tài)服務(wù)價值普遍穩(wěn)定于平均水平，單元分布最為集中，但也存在不少生態(tài)服務(wù)價值較低的離群點，通過ArcGIS平臺鎖定這些離群點主要分布于安溪縣中西部和長泰縣東北部。

表4 廈漳泉地區(qū)生態(tài)服務(wù)價值的地形梯度效應(yīng)

4 討論與結(jié)論

高生態(tài)服務(wù)價值區(qū)域集中分布于一些重要山脈(戴云山、大芹山等)、主要溪流、大型水庫等，形成整個地區(qū)的生態(tài)輸出區(qū)[23]，而這些區(qū)域經(jīng)濟發(fā)展較為落后甚至是貧困區(qū)，應(yīng)當(dāng)獲得密集城鎮(zhèn)化區(qū)的生態(tài)補償，通過健全地區(qū)轉(zhuǎn)移支付體系，開展“生態(tài)扶貧”項目給予其經(jīng)濟發(fā)展和生態(tài)維護的支持。針對不同梯度采取相應(yīng)措施，重點保育沿海平原區(qū)的濕地、水域，防治山地、丘陵區(qū)的地表侵蝕。在快速城鎮(zhèn)化地區(qū)探索構(gòu)建“城市立體綠化”格局，以期為高密度城鎮(zhèn)化區(qū)域披上“綠裳”。

生態(tài)服務(wù)價值評估方法的差異可能導(dǎo)致不同的結(jié)果，根據(jù)結(jié)果的實用性選擇適當(dāng)?shù)姆椒ㄟ€有待探討。受限于數(shù)據(jù)獲取，采用兩期數(shù)據(jù)分析2006—2016年這10 a間地區(qū)演變略有不足，繼續(xù)挖掘數(shù)據(jù)以增長研究期和研究樣本，可更客觀科學(xué)地闡釋生態(tài)服務(wù)價值演變規(guī)律。生態(tài)服務(wù)價值時空格局是社會因素和結(jié)構(gòu)性因素作用的復(fù)雜結(jié)果，這些因素對地形梯度效應(yīng)及空間尺度效應(yīng)的驅(qū)動機制需要未來長期深入研究，以彌補本文在這部分的不足。

2006—2016年廈漳泉地區(qū)生態(tài)服務(wù)價值水平總體上有所提升，其分布呈現(xiàn)出三大效應(yīng)。尺度效應(yīng)：隨著格網(wǎng)尺度的增加，地區(qū)生態(tài)服務(wù)價值空間自相關(guān)性具有波動上升的趨勢，而安溪水土流失區(qū)低—低集聚不再明顯，但沿海風(fēng)積地貌區(qū)和密集城鎮(zhèn)化區(qū)低—低集聚態(tài)勢依然存在。時空效應(yīng)：10 a間，生態(tài)服務(wù)價值空間集聚性總體上有所增強；“兩高一低”分布格局愈加明顯，北部高—高集聚區(qū)域有所擴張，南部則發(fā)生大幅度轉(zhuǎn)向漸趨于漳州市西南，而沿海低—低集聚區(qū)域沿東北—西南方向延伸，逐漸接連成片。地形梯度效應(yīng)：隨地形起伏的增加，地區(qū)生態(tài)服務(wù)價值水平總體上升，而各單元價值陡增和陡降的情況并存，并逐漸趨于穩(wěn)定。