?？谑泻０稁鷳B(tài)系統(tǒng)服務(wù)及其時(shí)空變異

喻露露, 張曉祥, 李楊帆, 趙杏杏, 馬生全, Huhua Cao,曹均闊

1 河海大學(xué)地理信息科學(xué)與工程研究所，南京　210098 2 廈門大學(xué)海洋與海岸帶發(fā)展研究院，廈門　361102 3 海南師范大學(xué)信息學(xué)院，?？凇?71158 4 Department of Geography, School of International Development and Global Studies, University of Ottawa, Ottawa K1N 6N5, Canada

?

海口市海岸帶生態(tài)系統(tǒng)服務(wù)及其時(shí)空變異

喻露露1, 張曉祥1, 李楊帆2,*, 趙杏杏1, 馬生全3, Huhua Cao4,曹均闊3

1 河海大學(xué)地理信息科學(xué)與工程研究所，南京210098 2 廈門大學(xué)海洋與海岸帶發(fā)展研究院，廈門361102 3 海南師范大學(xué)信息學(xué)院，海口571158 4 Department of Geography, School of International Development and Global Studies, University of Ottawa, Ottawa K1N 6N5, Canada

摘要:海岸帶位于陸地、海洋交匯地帶，有著重要的生態(tài)價(jià)值，但也是生態(tài)環(huán)境較為脆弱的地區(qū)。以?？谑袨檠芯繀^(qū)，研究?？谑泻０稁鷳B(tài)系統(tǒng)服務(wù)及其時(shí)空變異。基于1990、1995、2004和2012年4期Landsat TM影像得到土地利用覆蓋數(shù)據(jù)；采用基于當(dāng)量因子的生態(tài)系統(tǒng)服務(wù)價(jià)值評價(jià)方法，結(jié)合?？谑袑?shí)際情況，構(gòu)建了?？谑泻０稁鷳B(tài)系統(tǒng)服務(wù)價(jià)值動(dòng)態(tài)評估模型；借助GIS空間分析技術(shù)，應(yīng)用局部空間自相關(guān)模型，定量研究了生態(tài)系統(tǒng)服務(wù)價(jià)值(ESV)的時(shí)空變異特征。研究結(jié)果表明：(1)在生態(tài)系統(tǒng)服務(wù)的總值方面，過去20多年的ESV整體變化不大，但近年來ESV有所降低，生態(tài)環(huán)境出現(xiàn)退化趨勢。(2)在生態(tài)系統(tǒng)服務(wù)的時(shí)空變異方面，?？谑泻０稁SV的空間正相關(guān)性顯著，空間集聚程度較高，高值區(qū)主要集中于中心城區(qū)西部和東寨港周圍，低值區(qū)主要集中于中心城區(qū)；近來，在發(fā)展城市經(jīng)濟(jì)和保護(hù)生態(tài)環(huán)境的雙重作用下，位于中心城區(qū)的低值區(qū)和位于東寨港紅樹林自然保護(hù)區(qū)的高值區(qū)都有所輻射擴(kuò)張。

關(guān)鍵詞：生態(tài)系統(tǒng)服務(wù)；GIS；空間相關(guān)性；?？谑泻０稁?/p>

生態(tài)系統(tǒng)服務(wù)(Ecosystem Services)是指生態(tài)系統(tǒng)與生態(tài)過程所形成及維持的人類賴以生存的自然環(huán)境條件與效用[1]，近年來，全球土地利用的變化導(dǎo)致全球生態(tài)系統(tǒng)服務(wù)功能的大量損失[2]。生態(tài)系統(tǒng)服務(wù)體現(xiàn)了生態(tài)系統(tǒng)的功能狀況，對于像海岸帶這種人類活動(dòng)劇烈、脆弱性程度高的敏感地帶，進(jìn)行生態(tài)系統(tǒng)服務(wù)的價(jià)值評價(jià)與空間分析具有重要意義[3- 4]。1997年Costanza等對全球生態(tài)系統(tǒng)服務(wù)功能進(jìn)行劃分與評估[5]，成為定量評價(jià)生態(tài)系統(tǒng)服務(wù)的里程碑。近年來，國外學(xué)者從生態(tài)系統(tǒng)服務(wù)的定義、度量指標(biāo)、評估方法等方面進(jìn)行了更深入的研究[6- 8]，在評估過程中越來越重視模型方法的使用[9- 10]，并開始注重研究單項(xiàng)生態(tài)系統(tǒng)服務(wù)之間的相互作用[11- 12]。而國內(nèi)學(xué)者主要在以下幾方面取得較大發(fā)展：對生態(tài)系統(tǒng)服務(wù)價(jià)值評估方法的比較[13- 14]，對不同區(qū)域和不同自然地理區(qū)的生態(tài)系統(tǒng)服務(wù)價(jià)值的評估[15]，對森林[16]、海洋[17]、草地[18]、湖泊[19]等不同類型生態(tài)系統(tǒng)服務(wù)價(jià)值的評估，對生態(tài)多樣性保護(hù)[20]、土壤保持[21]、環(huán)境凈化[22]等單項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值的評估，對生態(tài)系統(tǒng)服務(wù)價(jià)值響應(yīng)區(qū)域土地利用[23]或景觀格局變化[24]的研究，對生態(tài)系統(tǒng)服務(wù)評估新方法的探索[25- 26]。但總體而言，國內(nèi)學(xué)者在評估生態(tài)服務(wù)價(jià)值時(shí)，大多依據(jù)Costanza等[5]或謝高地等[27]的研究成果，直接使用單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值的平均值進(jìn)行計(jì)算，或者通過生物量、生物多樣性等修訂方法間接估算區(qū)域生態(tài)系統(tǒng)服務(wù)的價(jià)值。這種方法主觀性較強(qiáng)，不能充分地反映特定區(qū)域生態(tài)系統(tǒng)的結(jié)構(gòu)和特性，因此，構(gòu)建符合研究區(qū)域客觀情況的評價(jià)模型很有必要。

空間相關(guān)性和空間異質(zhì)性普遍存在于地理對象的空間格局問題中，對地理數(shù)據(jù)進(jìn)行空間分析是地理學(xué)領(lǐng)域的重要方向?？臻g自相關(guān)分析作為空間分析的一個(gè)分支，用于探索空間變量的相關(guān)性和相關(guān)程度，和Tobler的地理學(xué)第一定律密切相關(guān)[28]。生態(tài)系統(tǒng)服務(wù)價(jià)值和各種地理實(shí)體一樣，存在著空間關(guān)聯(lián)性，但之前的大部分研究都只是靜態(tài)的快照式研究，忽視了它的復(fù)雜性、動(dòng)態(tài)性以及非線性相關(guān)性。對生態(tài)服務(wù)價(jià)值空間分布的聚集模式、關(guān)聯(lián)特征進(jìn)行定量研究，有利于理解研究區(qū)生態(tài)系統(tǒng)的空間格局特征，揭示生態(tài)服務(wù)價(jià)值變化的影響因素和驅(qū)動(dòng)機(jī)制。

本文在借鑒已有研究成果的基礎(chǔ)上對?？谑泻０稁У纳鷳B(tài)系統(tǒng)服務(wù)價(jià)值進(jìn)行研究，利用多期遙感影像分類數(shù)據(jù)及其它統(tǒng)計(jì)數(shù)據(jù)，結(jié)合當(dāng)?shù)貙?shí)際情況，構(gòu)建海口市海岸帶生態(tài)系統(tǒng)服務(wù)價(jià)值動(dòng)態(tài)評估模型，對各年份的生態(tài)系統(tǒng)服務(wù)價(jià)值的分布情況進(jìn)行了可視化展示，定量研究分析了生態(tài)系統(tǒng)服務(wù)的時(shí)空分布特征和規(guī)律，為區(qū)域可持續(xù)發(fā)展提供決策依據(jù)。

1研究區(qū)與數(shù)據(jù)準(zhǔn)備

?？谑械靥幬覈Ｄ鲜”辈垦睾＃Ｄ鲜∫彩俏覈诙蠛u海南島，?？谑呛Ｄ鲜〉氖?，分設(shè)秀英、龍華、瓊山、美蘭4個(gè)區(qū)。?？谑泻０毒€長約136.2km，其中自然岸線長101.7km，占全市海岸線的74.7%，人工岸線長34.5km，占全市海岸線的25.3%，主要用于旅游、港口碼頭、防潮堤和防波堤、工程用海、漁業(yè)養(yǎng)殖等。?？谑泻０稁郦M長，其東部的東寨港紅樹林自然保護(hù)區(qū)為國家級自然保護(hù)區(qū)，也是國際重要濕地保護(hù)區(qū)之一，對海岸帶起到重要的保護(hù)作用。海南島最長的河流南渡江從?？谑泻０稁е胁看┻^而入海，河流搬運(yùn)大量的泥沙堆積在入?？?，造成海岸淤積，對海岸線產(chǎn)生一定的影響。?？谑泻０稁Ы^大部分區(qū)域?qū)儆诤？谑兄行某菂^(qū)范圍，人口壓力大，道路密度高。此外，?？谑泻０稁碛胸S富的旅游資源，在旅游開發(fā)的同時(shí)，由旅游業(yè)帶動(dòng)的房地產(chǎn)行業(yè)也得到蓬勃發(fā)展，建設(shè)用地的擴(kuò)展趨勢非常強(qiáng)烈。海口市海岸帶的過度開發(fā)趨勢明顯，土地利用結(jié)構(gòu)變化劇烈，導(dǎo)致當(dāng)?shù)厣鷳B(tài)系統(tǒng)服務(wù)功能發(fā)生明顯的改變。

研究區(qū)海岸線東起澄邁灣，西至鋪前鎮(zhèn)，海岸帶研究范圍分為陸域部分和海域部分，陸域部分為自?？谑泻０毒€向陸延伸10km的濱海陸地，海域部分由沿海灘涂、紅樹林兩類濕地生態(tài)系統(tǒng)組成。

圖1　?？谑泻０稁а芯繀^(qū)示意圖Fig.1　Location of Haikou city coastal zone

序號Sequencenumber數(shù)據(jù)名稱Dataname數(shù)據(jù)格式Dataformat數(shù)據(jù)時(shí)間Datatime數(shù)據(jù)來源Datasource1遙感影像Landsat條帶號成像日期123/46,124/461988-06-08,1991-10-30123/46,124/461994-08-12,1995-09-23123/46,124/462004-02-13,2004-12-20123/46,124/462012-06-29,2012-10-26美國馬里蘭大學(xué)遙感數(shù)據(jù)下載網(wǎng)站2海南省糧食產(chǎn)量統(tǒng)計(jì)數(shù)據(jù)1990,1995,2004,2012海南省統(tǒng)計(jì)年鑒3海南省耕地面積統(tǒng)計(jì)數(shù)據(jù)1990,1995,2004,2012海南省統(tǒng)計(jì)年鑒4海南省糧食價(jià)格統(tǒng)計(jì)數(shù)據(jù)2001海南省統(tǒng)計(jì)年鑒5海南省GDP統(tǒng)計(jì)數(shù)據(jù)1990,1995,2004,2012海南省統(tǒng)計(jì)年鑒

本文基于研究區(qū)的4期(1990年、1995年、2004年、2012年)Landsat TM遙感影像數(shù)據(jù)，同時(shí)結(jié)合土地利用現(xiàn)狀數(shù)據(jù)進(jìn)行遙感影像分類，將研究區(qū)的土地利用分為耕地、林地、建設(shè)用地、沿海灘涂、紅樹林、水域和未利用土地7類。為研究生態(tài)系統(tǒng)服務(wù)價(jià)值的需要，結(jié)合?？谑泻０稁ё匀坏乩憝h(huán)境，將沿海灘涂和紅樹林劃為濕地生態(tài)系統(tǒng)，林地劃為森林生態(tài)系統(tǒng)，耕地劃為農(nóng)田生態(tài)系統(tǒng)，水域劃為水體生態(tài)系統(tǒng)，未利用土地劃為未利用地生態(tài)系統(tǒng)，另外，考慮到建設(shè)用地對生態(tài)環(huán)境產(chǎn)生的負(fù)效應(yīng)遠(yuǎn)遠(yuǎn)大于正效應(yīng)，故本文不考慮其生態(tài)系統(tǒng)服務(wù)價(jià)值[23]。在綜合考慮研究區(qū)范圍、數(shù)據(jù)的疏密性以及工作效率的基礎(chǔ)之上，將土地利用數(shù)據(jù)的格網(wǎng)大小定為500m×500m。除了遙感數(shù)據(jù)之外，其他數(shù)據(jù)還包括海南省的糧食產(chǎn)量、糧食價(jià)格等數(shù)據(jù)。

2研究方法

2.1生態(tài)系統(tǒng)服務(wù)價(jià)值動(dòng)態(tài)評估模型構(gòu)建

本文采用基于當(dāng)量因子的生態(tài)系統(tǒng)服務(wù)價(jià)值評價(jià)方法，構(gòu)建?？谑泻０稁鷳B(tài)系統(tǒng)服務(wù)價(jià)值動(dòng)態(tài)評估模型。根據(jù)謝高地等[27]的定義，生態(tài)系統(tǒng)生態(tài)服務(wù)價(jià)值1個(gè)當(dāng)量因子為1hm2全國平均產(chǎn)量的農(nóng)田每年自然糧食產(chǎn)量的經(jīng)濟(jì)價(jià)值，在生態(tài)系統(tǒng)服務(wù)體系中對應(yīng)為農(nóng)田食物生產(chǎn)服務(wù)，1個(gè)生態(tài)服務(wù)價(jià)值當(dāng)量因子的經(jīng)濟(jì)價(jià)值量為當(dāng)年全國平均糧食單產(chǎn)市場價(jià)值的1/7。

我國各地的自然地理環(huán)境差異顯著，這決定著各地的糧食產(chǎn)量及生態(tài)系統(tǒng)的生物量相差較大，所以需要對文獻(xiàn)[27]制定的代表全國平均水平的生態(tài)服務(wù)價(jià)值當(dāng)量因子表進(jìn)行改進(jìn)，使其符合區(qū)域特征。海口市海岸帶生態(tài)系統(tǒng)服務(wù)價(jià)值動(dòng)態(tài)評估模型的構(gòu)建如圖2所示，包括3個(gè)部分：(1)影像數(shù)據(jù)預(yù)處理，對獲得的遙感影像進(jìn)行分類、類型合并，得到各生態(tài)系統(tǒng)的面積；(2)生物量修訂，結(jié)合?？谑泻０稁У膶?shí)際生物量，得到各生態(tài)系統(tǒng)的生物量因子調(diào)整系數(shù)；(3)地區(qū)糧食產(chǎn)量、價(jià)格調(diào)整，計(jì)算一個(gè)當(dāng)量因子的價(jià)值量值。

圖2　?？谑泻０稁SV動(dòng)態(tài)評估模型構(gòu)建流程圖Fig.2　The flowchart of building dynamic evaluation model of coastal Ecosystem Services Values in Haikou city

關(guān)于模型中的生物量修訂，一般來說，生物量越大，生態(tài)服務(wù)功能越強(qiáng)[29]，于是為研究方便起見，假設(shè)生物量與生態(tài)系統(tǒng)服務(wù)價(jià)值呈線性關(guān)系，并在此基礎(chǔ)上利用植被凈第一生產(chǎn)力(net primary productivity, NPP)的比較分析進(jìn)行生物量因子的調(diào)整[30]。植物的凈第一性生產(chǎn)力(NPP)反映了植物群落在自然環(huán)境條件下的生產(chǎn)能力，因而能夠反映出各土地類型的生態(tài)服務(wù)的高低。另外，關(guān)于地區(qū)糧食價(jià)格的調(diào)整，本文剔除價(jià)格變化的影響，直接將2001年?？谑懈鞣N糧食作物的平均價(jià)格定為糧食單產(chǎn)值的計(jì)算價(jià)格。

模型中所涉及到的計(jì)算公式及其變量含義見表2。

2.2局部空間自相關(guān)分析

為定量研究?？谑泻０稁鷳B(tài)系統(tǒng)服務(wù)的空間分布差異和變化規(guī)律，本文在利用構(gòu)建好的價(jià)值評估模型計(jì)算出各年份的生態(tài)系統(tǒng)服務(wù)總價(jià)值的基礎(chǔ)上，采用Moran散點(diǎn)圖和空間聯(lián)系的局部指標(biāo)(LISA)相結(jié)合的方法，進(jìn)行局部空間自相關(guān)分析，探索局部空間關(guān)聯(lián)模式，反映空間關(guān)聯(lián)模式的顯著性水平。具體流程見圖3，分為A、B和C三部分。其中，A部分為數(shù)據(jù)預(yù)處理，得到ESV格網(wǎng)數(shù)據(jù)；B部分為探索局部空間關(guān)聯(lián)模式，通過計(jì)算區(qū)域單元的空間權(quán)重矩陣，得到Moran散點(diǎn)圖；C部分為反映空間關(guān)聯(lián)模式的顯著性水平，計(jì)算局部Moran′s I指數(shù)，并進(jìn)行顯著性水平檢驗(yàn)，結(jié)合Moran散點(diǎn)圖，得到LISA聚類地圖和LISA顯著性水平圖。

表2　評估模型中的公式及其變量

圖3　海口市海岸帶ESV時(shí)空變異分析流程圖Fig.3　The flowchart of spatial-temporal variation analysis of coastal Ecosystem Services Values in Haikou city

Moran散點(diǎn)圖的橫軸是各單元標(biāo)準(zhǔn)化處理后的屬性值，縱軸為該單元與周圍單元屬性值的加權(quán)平均，即空間滯后因子。Moran散點(diǎn)圖的4個(gè)象限表達(dá)了空間單元與鄰域的4種空間聯(lián)系類型，分別為“高高(HH)”、“高低(HL)”、“低高(LH)”和“低低(LL)”。HH類型和LL類型為正的空間自相關(guān)關(guān)系，表示空間單元和周邊單元的空間差異程度小，而HL類型和LH類型為負(fù)的空間自相關(guān)關(guān)系，表示空間單元和周邊單元的空間差異程度大。

空間聯(lián)系局部指標(biāo)值由局部Moran′s I指數(shù)計(jì)算得到，計(jì)算公式為：

(1)

在給定顯著性水平α的情況下，可用標(biāo)準(zhǔn)化統(tǒng)計(jì)量對局部Moran′s I指數(shù)的顯著性進(jìn)行檢驗(yàn)，如公式(2)所示。

(2)

式中，Z(Ii)表示Ii的標(biāo)準(zhǔn)化，E(Ii)表示Ii的數(shù)學(xué)期望，VAR(Ii)表示方差。

本文利用Geoda095i軟件計(jì)算空間權(quán)重矩陣，繪制Moran散點(diǎn)圖，計(jì)算局部Moran′sI指數(shù)，并在顯著性水平P=0.05的情況下，利用99次置換方法對局部Moran′sI指數(shù)進(jìn)行顯著性檢驗(yàn)。

3結(jié)果與分析

3.1數(shù)據(jù)假設(shè)

本文基于前人的研究成果[31- 32]，計(jì)算了海口市海岸帶生物量因子調(diào)整系數(shù)?；跀?shù)據(jù)的可獲得性，除農(nóng)田、森林之外的其他生態(tài)系統(tǒng)的生物量水平未做調(diào)整，具體情況見表3。

表3　生態(tài)系統(tǒng)生物量因子調(diào)整系數(shù)及其基礎(chǔ)數(shù)據(jù)

關(guān)于?？谑泻０稁ЪZ食產(chǎn)量的調(diào)整，考慮到本研究的重點(diǎn)以及數(shù)據(jù)的可得性，將海南省的糧食產(chǎn)量作為?？谑泻０稁У募Z食產(chǎn)量來參與計(jì)算。表4是通過查閱海南省的各年份統(tǒng)計(jì)年鑒，計(jì)算得到的海南省各年份單位耕地面積的糧食生產(chǎn)情況，以此為基礎(chǔ)得到海口市海岸帶各年份的1個(gè)生態(tài)服務(wù)價(jià)值當(dāng)量因子的經(jīng)濟(jì)價(jià)值量En。

表4　海南省各年份單位耕地面積糧食生產(chǎn)情況

3.2生態(tài)系統(tǒng)服務(wù)價(jià)值評估結(jié)果分析

表5和圖4顯示了?？谑泻０稁鷳B(tài)系統(tǒng)的面積變化情況和分布變化情況。從表5可以看出，從1990至2012年，?？谑猩鷳B(tài)系統(tǒng)的面積發(fā)生了很大的變化，變化最明顯的是未利用地生態(tài)系統(tǒng)和水體生態(tài)系統(tǒng)。首先，未利用地生態(tài)系統(tǒng)減少了89.27%，這體現(xiàn)了人類對未利用土地的開發(fā)力度。其次，水體生態(tài)系統(tǒng)增加了63.72%，這和水產(chǎn)養(yǎng)殖面積的增大密切相關(guān)。

表5　?？谑泻０稁Ц髂攴莞魃鷳B(tài)系統(tǒng)面積

圖4　?？谑泻０稁鷳B(tài)系統(tǒng)分布圖Fig.4　The distribution diagram of ecosystem in Haikou coastal zone

表6和圖5顯示的是?？谑泻０稁Ц髂攴萆鷳B(tài)系統(tǒng)服務(wù)價(jià)值情況?？傮w看來，對生態(tài)系統(tǒng)服務(wù)總價(jià)值貢獻(xiàn)最大的是森林生態(tài)系統(tǒng)，其次是農(nóng)田生態(tài)系統(tǒng)和水體生態(tài)系統(tǒng)。從時(shí)間角度來看，?？谑泻０稁鷳B(tài)系統(tǒng)服務(wù)總價(jià)值在1990至2004年之間處于增長趨勢，在2004至2012年之間出現(xiàn)下降的情況。對于單個(gè)生態(tài)系統(tǒng)而言，水體生態(tài)系統(tǒng)的生態(tài)系統(tǒng)服務(wù)價(jià)值持續(xù)增加，未利用地生態(tài)系統(tǒng)的生態(tài)系統(tǒng)服務(wù)價(jià)值持續(xù)減少，而森林生態(tài)系統(tǒng)、農(nóng)田生態(tài)系統(tǒng)和濕地生態(tài)系統(tǒng)的生態(tài)系統(tǒng)服務(wù)價(jià)值都處于有增有減的狀態(tài)。

表6　海口市海岸帶各年份生態(tài)系統(tǒng)服務(wù)價(jià)值表/(104元)

圖5　海口市海岸帶生態(tài)系統(tǒng)服務(wù)價(jià)值變化情況Fig.5　Changes in ecosystem services values in Haikou coastal zone

3.3生態(tài)系統(tǒng)服務(wù)價(jià)值時(shí)空變化總體特征

本研究利用GIS技術(shù)將生態(tài)系統(tǒng)服務(wù)總價(jià)值量化到每個(gè)格網(wǎng)，來分析生態(tài)系統(tǒng)服務(wù)價(jià)值的時(shí)空演變。圖6為各年度生態(tài)系統(tǒng)服務(wù)價(jià)值的空間分布情況圖。從空間角度來看，由于?？谥行某菂^(qū)西面的森林覆蓋率比較高，生態(tài)完整性好，所以生態(tài)服務(wù)價(jià)值高；東寨港周圍分布有國家級的紅樹林保護(hù)區(qū)，人類活動(dòng)干擾程度低，所以生態(tài)服務(wù)價(jià)值也較高；海口市中心城區(qū)人口密度高，城市化水平高，對自然生態(tài)系統(tǒng)的干擾程度大，所以生態(tài)服務(wù)價(jià)值低。從時(shí)間角度來看，1990—2012年間，生態(tài)系統(tǒng)服務(wù)高值區(qū)域在不斷縮小，且離散度變高，而低值區(qū)域以海口市中心城區(qū)為中心向周圍不斷擴(kuò)張。

圖6　?？谑泻０稁鷳B(tài)系統(tǒng)服務(wù)價(jià)值空間分布Fig.6　The spatial distribution of ecosystem services values in Haikou coastal zone

3.4生態(tài)系統(tǒng)服務(wù)價(jià)值局部空間自相關(guān)分析

Moran散點(diǎn)圖(圖7)展現(xiàn)了區(qū)域單元和周圍單元的空間聯(lián)系形式，從圖中可以看出，生態(tài)系統(tǒng)服務(wù)價(jià)值散點(diǎn)主要分布在第一象限(HH)和第三象限(LL)，第二象限(LH)和第四象限(HL)散點(diǎn)分布相對較少，說明生態(tài)系統(tǒng)服務(wù)價(jià)值具有較大的空間正相關(guān)性。從1990至2012年，沿趨勢線分布的散點(diǎn)越來越多，散點(diǎn)與趨勢線匹配的效果越來越好，這體現(xiàn)了自相關(guān)程度的增強(qiáng)。

LISA聚類圖(圖8)是對Moran散點(diǎn)圖中通過了顯著性檢驗(yàn)的區(qū)域單元的地理表達(dá)。可以看出，生態(tài)系統(tǒng)服務(wù)價(jià)值之間的空間集聚主要為空間相似值之間的集聚，空間奇異值極少，這印證了Moran散點(diǎn)圖所揭示的空間正相關(guān)性。和生態(tài)系統(tǒng)服務(wù)價(jià)值的分布情況相一致，高值集聚區(qū)主要分布在中心城區(qū)西面和東寨港周圍，低值集聚區(qū)主要是以中心城區(qū)為中心分布。

LISA顯著性水平圖(圖8)顯示了和周圍單元具有顯著的空間聯(lián)系的區(qū)域單元。從空間上看，大部分區(qū)域的相關(guān)性不顯著，主要在高值集聚區(qū)和低值集聚區(qū)表現(xiàn)出較高的顯著性水平，顯著性水平達(dá)到0.01的區(qū)域位于顯著性水平達(dá)到0.05區(qū)域的中間位置。從時(shí)間上看，較高顯著水平(P=0.01)的區(qū)域面積呈增長趨勢，主要發(fā)生在東寨港周圍以及?？谑兄行某菂^(qū)周圍，表明這部分區(qū)域的生態(tài)系統(tǒng)服務(wù)價(jià)值在空間分布上的差異性變小，空間趨同性增強(qiáng)。

圖8　?？谑泻０稁?個(gè)時(shí)期生態(tài)系統(tǒng)服務(wù)價(jià)值顯著性檢驗(yàn)Fig.8　The significance testing of the ecosystem services values in Haikou coastal zone

4結(jié)論與討論

本文在借鑒已有的研究成果基礎(chǔ)之上，采用基于當(dāng)量因子的生態(tài)系統(tǒng)服務(wù)價(jià)值評價(jià)方法，結(jié)合?？谑泻０稁У膶?shí)際自然地理和社會經(jīng)濟(jì)情況，構(gòu)建了?？谑泻０稁鷳B(tài)系統(tǒng)服務(wù)價(jià)值動(dòng)態(tài)評估模型，并利用局部空間自相關(guān)模型，對生態(tài)系統(tǒng)服務(wù)價(jià)值的時(shí)空變異性進(jìn)行研究。研究表明：

(1)?？谑泻０稁г?990、1995、2004和2012這4a的生態(tài)系統(tǒng)服務(wù)價(jià)值(ESV)分別為37.05億元、41.46億元、42.19億元和39.24億元。整體變化不大，但近年來ESV有所降低，生態(tài)環(huán)境出現(xiàn)退化趨勢。

(2)?？谑泻０稁鷳B(tài)系統(tǒng)服務(wù)價(jià)值有顯著的正空間相關(guān)性，空間集聚程度較高，高值區(qū)主要集中于中心城區(qū)西部和東寨港周圍，低值區(qū)主要集中于中心城區(qū)；1990—2012年，高值區(qū)以中心城區(qū)為中心向四周輻射擴(kuò)張，低值區(qū)在東寨港周圍有所增加，而在中心城區(qū)西部有所減少。

本文在對生態(tài)系統(tǒng)服務(wù)價(jià)值的研究過程中不僅考慮到區(qū)域差異，而且實(shí)現(xiàn)了動(dòng)態(tài)評估。在空間上，通過生物量指標(biāo)對全國平均水平的生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量因子表進(jìn)行修訂；在時(shí)間上，通過各年份的海南省糧食產(chǎn)量調(diào)整，實(shí)現(xiàn)生態(tài)系統(tǒng)服務(wù)價(jià)值的動(dòng)態(tài)評估。此外，在對生態(tài)系統(tǒng)服務(wù)的時(shí)空變異研究方面，不同于大多數(shù)對生態(tài)系統(tǒng)服務(wù)價(jià)值分布格局的空間展示，本研究利用GIS空間分析的方法，對生態(tài)系統(tǒng)服務(wù)價(jià)值空間聚集規(guī)律和空間關(guān)聯(lián)模式的動(dòng)態(tài)變化進(jìn)行定量研究。

本研究也存在著一些不足之處。首先，本文是從土地利用角度來估算生態(tài)系統(tǒng)服務(wù)價(jià)值，所以對土地利用類型的解譯精度要求較高。其次，由于中國陸地生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量表是以全國民眾對各項(xiàng)服務(wù)的生態(tài)服務(wù)的可能的支付意愿和支付能力為基礎(chǔ)的，所以在構(gòu)建地區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值評估的模型中，需要考慮到當(dāng)?shù)貙Ω黜?xiàng)生態(tài)系統(tǒng)服務(wù)的支付意愿和支付能力。受限于數(shù)據(jù)來源的局限，本研究對這些因素考慮不足，有待進(jìn)一步發(fā)展和改進(jìn)。

參考文獻(xiàn)(References):

[1]Daily G C. Nature′s Services: Social Dependence on Nature Ecosystems. Washington: Island Press, 1997.

[2]Costanza R, de Groot R, Sutton P, van der Ploeg S, Anderson S J, Kubiszewski I, Farber S, Turner R K. Changes in the global value of ecosystem services. Global Environmental Change, 2014, 26: 152- 158.

[3]陳偉琪, 王萱. 圍填海造成的海岸帶生態(tài)系統(tǒng)服務(wù)損耗的貨幣化評估技術(shù)探討. 海洋環(huán)境科學(xué), 2009, 28(6): 749- 754.

[4]索安寧, 于永海, 韓富偉. 環(huán)渤海海岸帶生態(tài)服務(wù)價(jià)功能評價(jià). 海洋開發(fā)與管理, 2011, 28(7): 67- 73.

[5]Costanza R, Arge R, Groot R, Farber S, Grasso M, Hannon B, Limburg K, Naeem S, O′Neill R V, Paruelo J, Raskin R G, Sutton P, van den Belt M. The value of the world′s ecosystem services and natural capital[J]. Nature, 1997, 387: 253- 260.

[6]Bateman I J, Harwood A R, Mace G M, Watson R T, Abson D J, Andrews B, Binner A, Crowe A, Day B H, Dugdale S, Fezzi C, Foden J, Hadley D, Haines-Young R, Hulme M, Kontoleon A, Lovett A A, Munday P, Pascual U, Paterson J, Perino G, Sen A, Siriwardena G, van Soest D, Termansen M. Bringing ecosystem services into economic decision-making: land use in the United Kingdom. Science, 2013, 341(6141): 45- 50.

[7]Boyd J, Banzhaf S. What are ecosystem services? The need for standardized environmental accounting units. Ecological Economics, 2007, 63(2/3): 616- 626.

[8]Fisher B, Polasky S, Sterner T. Conservation and human welfare: economic analysis of ecosystem services. Environmental and Resource Economics, 2011, 48(2): 151- 159.

[9]Nelson E, Mendoza G, Regetz J, Polasky S, Tallis H, Cameron D R, Chan K M, Daily G C, Goldstein J, Kareiva P M, Lonsdorf E, Naidoo R, Ricketts T H, Shaw M R. Modeling multiple ecosystem services, biodiversity conservation, commodity production and tradeoffs at landscape scales. Frontiers in Ecology and the Environment, 2009, 7(1): 4- 11.

[10]Polasky S, Nelson E, Pennington D,Johnson K A. The impact of land-use change on ecosystem services, biodiversity and returns to landowners: a case study in the state of Minnesota. Environmental and Resource Economics, 2011, 48(2): 219- 242.

[11]Bennett E M, Peterson G D, Gordon L J. Understanding relationships among multiple ecosystem services. Ecology letters, 2009, 12(12): 1394- 1404.

[12]Raudsepp-Hearne C, Peterson G D, Teng? M, Bennett E M, Holland T, Benessaiah K, MacDonald G K, Pfeifer L. Untangling the environmentalist′s paradox:why is human well-being increasing as ecosystem services degrade?. BioScience, 2010, 60(8): 576- 589.

[13]石龍宇, 崔勝輝, 尹鍇, 劉江. 廈門市土地利用/覆被變化對生態(tài)系統(tǒng)服務(wù)的影響. 地理學(xué)報(bào), 2010, 65(6): 708- 714.

[14]呂一河, 張立偉, 王江磊. 生態(tài)系統(tǒng)及其服務(wù)保護(hù)評估: 指標(biāo)與方法. 應(yīng)用生態(tài)學(xué)報(bào), 2013, 24(5): 1237- 1243.

[15]竇玥, 戴爾阜, 吳紹洪. 區(qū)域土地利用變化對生態(tài)系統(tǒng)脆弱性影響評估——以廣州市花都區(qū)為例. 地理研究, 2012, 31(2): 311- 322.

[16]吳霜, 延曉冬, 張麗娟. 中國森林生態(tài)系統(tǒng)能值與服務(wù)功能價(jià)值的關(guān)系. 地理學(xué)報(bào), 2014, 69(3): 334- 342.

[17]石洪華, 鄭偉, 丁德文, 呂吉斌, 張學(xué)雷. 典型海洋生態(tài)系統(tǒng)服務(wù)功能及價(jià)值評估——以桑溝灣為例. 海洋環(huán)境科學(xué), 2008, 27(2): 101- 104.

[18]陳春陽, 陶澤興, 王煥炯, 戴君虎. 三江源地區(qū)草地生態(tài)系統(tǒng)服務(wù)價(jià)值評估. 地理科學(xué)進(jìn)展, 2012, 31(7): 978- 984.

[19]陳為國, 許文杰. 湖泊生態(tài)系統(tǒng)服務(wù)功能影響因子分析與評價(jià)研究. 節(jié)水灌溉, 2010(12): 35- 37.

[20]李景文, 范樹陽. 達(dá)里諾爾國家級自然保護(hù)區(qū)生物多樣性保護(hù)功能價(jià)值評估研究. 內(nèi)蒙古環(huán)境科學(xué), 2009, 21(6): 110- 115.

[21]饒恩明, 肖燚, 歐陽志云, 鄭華. 海南島生態(tài)系統(tǒng)土壤保持功能空間特征及影響因素. 生態(tài)學(xué)報(bào), 2013, 33(3): 746- 755.

[22]肖滋民, 王立華, 郝亮, 吳麗云. 濰坊市城市綠地生態(tài)系統(tǒng)環(huán)境凈化服務(wù)價(jià)值研究. 湖北農(nóng)業(yè)科學(xué), 2011, 50(19): 3929- 3933.

[23]黃云鳳, 崔勝輝, 石龍宇. 半城市化地區(qū)生態(tài)系統(tǒng)服務(wù)對土地利用/覆被變化的響應(yīng). 地理科學(xué)進(jìn)展, 2012, 31(5): 551- 560.

[24]索安寧, 于永海, 韓富偉. 遼河三角洲盤錦濕地景觀格局變化的生態(tài)系統(tǒng)服務(wù)價(jià)值響應(yīng). 生態(tài)經(jīng)濟(jì), 2011(6): 147- 151.

[25]呂一河, 馬志敏, 傅伯杰, 高光耀. 生態(tài)系統(tǒng)服務(wù)多樣性與景觀多功能性——從科學(xué)理念到綜合評估. 生態(tài)學(xué)報(bào), 2013, 33(4): 1153- 1159.

[26]Fu B, Wang S, Su C H, Forsius M. Linking ecosystem processes and ecosystem services. Current Opinion in Environmental Sustainability, 2013, 5(1): 4- 10.

[27]謝高地, 魯春霞, 冷允法, 鄭度, 李雙成. 青藏高原生態(tài)資產(chǎn)的價(jià)值評估. 自然資源學(xué)報(bào), 2003, 18(2): 189- 196.

[28]陳彥光. 基于Moran 統(tǒng)計(jì)量的空間自相關(guān)理論發(fā)展和方法改進(jìn). 地理研究, 2009, 28(6): 1449- 1463.

[29]謝高地, 甄霖, 魯春霞, 陳操. 一個(gè)基于專家知識的生態(tài)系統(tǒng)服務(wù)價(jià)值化方法. 自然資源學(xué)報(bào), 2008, 23(5): 911- 919.

[30]艾志華. 潛江市土地利用變化對生態(tài)系統(tǒng)服務(wù)價(jià)值的影響研究[D]. 武漢: 華中師范大學(xué), 2013.

[31]謝高地, 肖玉, 甄霖, 魯春霞. 我國糧食生產(chǎn)的生態(tài)服務(wù)價(jià)值研究. 中國生態(tài)農(nóng)業(yè)學(xué)報(bào), 2005, 13(3): 10- 13.

[32]趙同謙, 歐陽志云, 鄭華, 王效科, 苗鴻. 中國森林生態(tài)系統(tǒng)服務(wù)功能及其價(jià)值評價(jià). 自然資源學(xué)報(bào), 2004, 19(4): 480- 491.

Coastal ecosystem services and their spatial-temporal variation in Haikou, China

YU Lulu1, ZHANG Xiaoxiang1, LI Yangfang2,*, ZHAO Xingxing1, MA Shengquan3, Huhua Cao4, CAO Junkuo3

1Instituteofgeographicinformationscienceandengineeringresearch,Hohaiuniversity,Nanjing210098 ,China2CoastalandOceanManagementInstitute,Xiamenuniversity,Xiamen361102 ,China3DepartmentofComputerScience&Technology,HainanNormalUniversity,Haikou571158 ,China4DepartmentofGeography,SchoolofInternationalDevelopmentandGlobalStudies,UniversityofOttawa,OttawaK1N6N5,Canada

Abstract：The coastal zone is in the land-ocean interaction area, which has great ecological importance and the ecological environment of this zone is often vulnerable. In the present research, the coastal zone in Haikou, China was selected as a case study to explore the coastal ecosystem services and their spatial-temporal variation. Firstly, four Landsat TM images collected in 1990, 1995, 2004, and 2012 are interpreted to obtain the land use/cover data. Secondly, based on the equivalent factor method of Ecosystem Services Values (ESV) evaluation and the local conditions, a dynamic assessment model of ecosystem services values in the Haikou coastal zone was constructed. Finally, using GIS-based spatial analysis technology and the local spatial autocorrelation model, the spatial-temporal variation characteristics of Ecosystem Services Values were quantitatively studied. Analysis of the results indicated: (i ) Change in ESV in the past twenty years has been insignificant, with ESV decreasing in recent years; thus the ecological environment has shown a trend of degradation; (ii) the spatial-temporal variation of ESV in the Haikou coastal zone is clear: spatial correlation is significant and the degree of spatial agglomeration is relatively high. High values are mainly distributed in the western part of the downtown area, and surrounding Dongzhaigang, whereas low values are mainly distributed in the downtown area. Recently, under the dual-strategy of economic development and environmental protection, both the low value downtown area and the high value Dongzhaigang mangrove nature reserve area have experienced continuous radiation expansion.

Key Words:ecosystem services; GIS; spatial correlation; Haikou coastal zone

基金項(xiàng)目:國家國際科技合作項(xiàng)目(2012DFA11270); 國家自然科學(xué)基金項(xiàng)目(41201394); 海南省國際科技合作專項(xiàng)(KJHZ2014-25)

收稿日期:2014- 10- 10; 網(wǎng)絡(luò)出版日期:2015- 08- 18

*通訊作者

Corresponding author.E-mail: yangf@xmu.edu.cn

DOI:10.5846/stxb201410101994

喻露露, 張曉祥, 李楊帆, 趙杏杏, 馬生全, Huhua Cao,曹均闊.?？谑泻０稁鷳B(tài)系統(tǒng)服務(wù)及其時(shí)空變異.生態(tài)學(xué)報(bào),2016,36(8):2431- 2441.

Yu L L, Zhang X X, Li Y F, Zhao X X, Ma S Q, Cao H, Cao J K.Coastal ecosystem services and their spatial-temporal variation in Haikou, China.Acta Ecologica Sinica,2016,36(8):2431- 2441.