摘要

以海口市為研究對象，分析 2010—2020 年土地利用/土地覆蓋（LULC）數(shù)據(jù)和生態(tài)系統(tǒng)服務價值的變化，采用PLUS模型和當量因子法模擬2035年自然發(fā)展（BAU）、生態(tài)保護（EC）和經(jīng)濟發(fā)展（ED）3種情景下生態(tài)系統(tǒng)服務價值的時空變化特征，探究不同土地利用模式對生態(tài)系統(tǒng)服務價值的影響。結(jié)果表明：?？谑型恋乩妙愋鸵粤值睾秃档卣急茸罡?，生態(tài)系統(tǒng)服務結(jié)構(gòu)相對穩(wěn)定；受建設用地增加影響，水文調(diào)節(jié)和水資源供給服務功能下降幅度最大，表明LULC變化與生態(tài)系統(tǒng)服務價值變化密切相關；EC情景下，各項生態(tài)系統(tǒng)服務均得到改善，LULC的結(jié)構(gòu)調(diào)整尤其是提高濕地占比，顯著改善區(qū)域整體生態(tài)系統(tǒng)服務。

Abstract

Using Haikou City as the research object， this paper analyzed the changes in land use and ecosystem service value from 2010 to 2020， and employed the PLUS model and the equivalent factor method to simulate the spatial and temporal characteristics of ecosystem service value under three scenarios of natural development （BAU）， ecological conservation （EC）， and economic development （ED） in 2035， and explored the impacts of different land use patterns on ecosystem service value. The results showed that land use types in Haikou City had the highest proportion of forested land and dry land， and the structure of ecosystem services was relatively stable. Hydrological regulation and water resource provisioning services declined the most under the influence of the increase of construction land， indicating that the change of LULC was closely related to the change of ecosystem service value. All ecosystem services improved under the EC scenario， and the structural adjustments of LULC， especially the increase of wetland proportion， significantly improved the overall ecosystem services in the region.

中圖分類號：TU986 文獻標志碼： A DOI：10.12233/j.gdyl.2024.02.010

Spatial-temporal Analysis and Muti-scenario Simulation of Ecosystem Service Value in Haikou Based on PLUS Model

FU Hui，CHEN Jie，LI Shuling

基金項目：海南省自然科學基金項目（編號：722QN288）；海南省高等學校教育教學改革研究資助項目（編號：Hnjg2024-10）；海南大學教育教學改革研究項目（編號：hdjy2420）；海南大學課程思政教育專項（編號：hdsz2205）；海南大學校級完全學分制改革研究項目（編號：hdxfz2204）

文章亮點

1）利用多情景模擬探究未來生態(tài)系統(tǒng)服務價值；2）區(qū)域LULC的結(jié)構(gòu)調(diào)整尤其是提高濕地占比，顯著改善區(qū)域整體生態(tài)系統(tǒng)服務；3）生態(tài)保護情景對未來城市生態(tài)規(guī)劃具有重要參考意義。

文章編號：1671-2641（2024）02-0065-08

關鍵詞

土地利用模式；生態(tài)系統(tǒng)服務價值；多情景模擬；PLUS模型

Keywords

Land use pattern；Ecosystem service value; Multiple scenarios simulation; PLUS model

收稿日期：2023-10-19

修回日期：2024-02-05

生態(tài)系統(tǒng)服務評估是關于生態(tài)系統(tǒng)結(jié)構(gòu)、功能與人類福祉之間關聯(lián)的研究領域。由于Robert Costanza等學者的引領作用[1～2]，國內(nèi)外開展大量對生態(tài)系統(tǒng)服務價值（Ecosystem service value）相關的評估工作，早期主要是針對全球和國家等大的空間尺度[3～4]。隨著研究的推進，對濕地、森林、草原和湖泊等中小尺度的研究也在逐漸增

多[5～6]。近年來，學者們從土地利用變化角度對生態(tài)系統(tǒng)服務開展評估和動態(tài)分析，包括基于土地利用的生態(tài)系統(tǒng)服務價值評估、對土地利用與生態(tài)系統(tǒng)服務的復雜相互影響的分析[7～10]，由此來探討不同生態(tài)系統(tǒng)服務的形成機理、驅(qū)動機制及土地利用優(yōu)化配置的策略。

科學政策的制定必須建立在對生態(tài)系統(tǒng)服務價值變化深層次原因的分析之上，不僅要厘清土地利用變化對生態(tài)系統(tǒng)服務的影響機制，更要分析土地利用變化的驅(qū)動機制，即這些變化是如何受到復雜的自然、經(jīng)濟和社會等因素的影響。從而揭示影響生態(tài)系統(tǒng)服務價值的關鍵指標，建立從土地利用變化驅(qū)動因素到生態(tài)系統(tǒng)服務價值的關聯(lián)系統(tǒng)，把生態(tài)系統(tǒng)服務價值的評估切實納入土地利用策略、經(jīng)濟、社會與環(huán)境方面的政策方案的制定中來。

?？谑惺呛Ｄ鲜〉氖?，也是政治、經(jīng)濟、文化和生態(tài)文明建設的前沿陣地。自2009年國際旅游島建設開始，海口進入了城市發(fā)展的高速階段，短期內(nèi)城市的急速擴張加劇了土地結(jié)構(gòu)的變化，同時也影響了城市的生態(tài)系統(tǒng)結(jié)構(gòu)與功能。研究基于2010、2015及2020年海口市土地利用/土地覆蓋（LULC）數(shù)據(jù)，識別?？谑蠰ULC及生態(tài)系統(tǒng)服務價值的時空變化特征，并運用PLUS模型模擬2035年3種情景下?？谑猩鷳B(tài)系統(tǒng)服務的變化發(fā)展，分析不同LULC變化對生態(tài)系統(tǒng)服務價值的影響，以期為?？谑谐擎?zhèn)建設開發(fā)和生態(tài)保護協(xié)調(diào)發(fā)展提供方法論指導和實踐依據(jù)。

1 材料與方法

1.1 材料

1.1.1 研究區(qū)概況

?？谑形挥诤Ｄ鲜”辈垦睾^(qū)域，國土面積2 296.82 km2， 屬熱帶季風氣候。?？趽碛胸S富的地表水和地下水資源，水資源總量達到11.08億m3。截至2021年末，?？谑邢螺犘阌?、龍華、瓊山和美蘭4個區(qū)，下設21個街道、22個鎮(zhèn)[11]（圖1）。

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

本文LULC數(shù)據(jù)采用中國科學院資源環(huán)境科學數(shù)據(jù)中心（https：//www.resdc.cn）2010、2015和2020年全國土地利用類型遙感監(jiān)測空間分布數(shù)據(jù)，空間分辨率為

30 m。參考2021年付杰的海南省生態(tài)系統(tǒng)類型分類[12]，將用地劃分為水田、旱地、林地、灌木林、草地、水域、濕地、建設用地和未利用地9類。數(shù)字高程模型（DEM）數(shù)據(jù)來自地理空間數(shù)據(jù)云（http：/www.gscloud.cn/）；氣溫和降雨數(shù)據(jù)來自國家氣象科學數(shù)據(jù)中心（http：//data.cma.cn/wa）；城市市政和公共服務設施分布來自百度在線地圖API接口獲取的興趣點（POI）數(shù)據(jù)；路網(wǎng)數(shù)據(jù)來源于OpenStreetMap（OSM）數(shù)據(jù)集；人口密度和國內(nèi)生產(chǎn)總值（GDP）數(shù)據(jù)由海南省各市縣統(tǒng)計公報、統(tǒng)計年鑒及中華人民共和國國民經(jīng)濟和社會發(fā)展第十四個五年規(guī)劃和社會經(jīng)濟數(shù)據(jù)等匯總得出；生態(tài)系統(tǒng)服務價值當量計算所需的主要糧食平均價格、播種面積和平均糧食產(chǎn)量數(shù)據(jù)來自2010—2020年的《中國農(nóng)產(chǎn)品價格調(diào)查年鑒》《海南省統(tǒng)計年鑒》及《中國統(tǒng)計年鑒》。

1.2 研究方法

1.2.1 土地利用單一動態(tài)度分析

土地利用單一動態(tài)度分析指不同土地利用類型的年均變化速度，即土地利用年變化率，該指標能夠反映土地利用變化的整體強度。公式如下：

（1）

式中，K表示土地利用年變化率，Ua和Ub分別表示某土地利用類型在研究階段初期和末期的面積，T表示初期和末期間的年份數(shù)。

1.2.2 PLUS 模型

PLUS模型是綜合了土地利用擴展分析策略和多類型隨機斑塊種子的新型土地利用模擬模型，通過結(jié)合隨機種子生成和閾值降低機制，實現(xiàn)對土地利用的精確模擬[13]。首先，根據(jù)?？谑袑嶋H情況和數(shù)據(jù)的可獲取性，選取社會和自然2方面的10個影響因子，包括高程、坡度、坡向、降雨、溫度、GDP、人口密度、城市道路、公共設施和市政設施，并進行標準化處理。其次，通過用地擴張分析（LEAS）模塊整合土地利用數(shù)據(jù)和影響因子數(shù)據(jù)，獲得各土地利用類型的發(fā)展概率?；?015—2020年的土地利用類型轉(zhuǎn)移概率矩陣，通過Markov-chain對2035年的不同情景下的土地需求進行數(shù)量預測。最后，將數(shù)量預測結(jié)果導入PLUS模型進行預測模擬。

1）轉(zhuǎn)換成本矩陣設置

轉(zhuǎn)移成本矩陣表示2種地類間相互轉(zhuǎn)換的可能性，1表示可以相互轉(zhuǎn)化，0則表示不能轉(zhuǎn)化?；诤？谑袑嶋H情況和相關研究經(jīng)驗[7， 9]，當前的社會生產(chǎn)力基本可以實現(xiàn)任意地類間的轉(zhuǎn)換，但考慮到各地類尤其是建設用地向水域的轉(zhuǎn)變成本較高且極少發(fā)生，因此主要對建設用地向水域和濕地的轉(zhuǎn)入進行限制，其余各地類在空間上均可進行交換（表1）。

2）領域權重設置

領域權重表示不同用地類型的擴張強度，取值范圍為0～1，越接近1，擴張能力越強，反之越弱。以2015年和2020年數(shù)據(jù)為基礎，對各土地利用的擴張能力進行計算。根據(jù)下式對擴張強度進行無量綱處理，所得參數(shù)如表2所示。

（2）

式中，表示權重值，為土地利用類型面積，、，分別代表面積最大和最小的土地利用類型的數(shù)量。

3）模型精度驗證

為驗證PLUS模型模擬的準確性，以2010年和2015年土地利用數(shù)據(jù)為基礎對2020年土地利用進行模擬，將模擬結(jié)果與2020年實際土地利用數(shù)據(jù)進行精度校驗。結(jié)果顯示， Kappa 系數(shù)為0.87，總體精度為0.92，表明模擬結(jié)果科學可靠。

1.2.3 生態(tài)系統(tǒng)服務價值計算方法

本文采用改進的中國陸地單位面積生態(tài)服務基礎當量表進行生態(tài)系統(tǒng)服務價值評估，并根據(jù)研究區(qū)實際情況進行調(diào)整。1）參照過往研究經(jīng)驗[12， 14～15]，對建設用地當量進行補充；2）依據(jù)2010—2020年?？谑衅骄Z食產(chǎn)量

（4 110.73 kg/hm2）、中國平均糧食產(chǎn)量（5 430.74 kg/hm2）

修正區(qū)域糧食產(chǎn)量。研究表明，1個標準當量的價值量等于該區(qū)域單位面積糧食市場價值的1/7。?？谑?010—2020年平均糧食收購價格為2.65元/kg，由此計算得出?？谑袉挝划斄恳蜃拥慕?jīng)濟價值為1 183.11元，并由此得到?？谑袉挝幻娣e生態(tài)系統(tǒng)服務價值系數(shù)表（表3）。?？谑懈饔玫仡愋偷纳鷳B(tài)系統(tǒng)服務計算基于Robert Costanza的生態(tài)系統(tǒng)服務價值模型[1～2]進行計算，公式如下：

（3）

（4）

式中，ESV為生態(tài)系統(tǒng)服務的總價值；為研究區(qū)第j項生態(tài)系統(tǒng)服務功能價值；為研究區(qū)第i種類型的土地面積；為第i種類型土地的單位ESV系數(shù)；為第i種土地利用類型的第j項生態(tài)系統(tǒng)服務功能價值系數(shù)[16]。

為了探究海口市生態(tài)系統(tǒng)服務價值在空間上的分異特性，采用Arcgis漁網(wǎng)工具建立500 m×500 m的方形網(wǎng)格，計算每個網(wǎng)格內(nèi)的ESV，按照ArcGIS標準分類中的幾何分隔法將?？谑蠩SV劃分為7個等級，Ⅰ～Ⅶ級取值區(qū)間依次為[-37.0，7.7）、[7.7，29.5）、[29.5，40.1）、[40.1，61.8）、[61.8，106.5）、[106.5，198.2）、[198.2，386.5]。

2 結(jié)果分析

2.1 ?？谑型恋乩妙愋蜁r空變化特征

從?？谑?010—2020年各地類的面積變化及空間分布情況（表4及圖2）可知，?？谑杏玫仡愋椭饕獮榱值?，廣泛分布于?？谑械臇|部和南部區(qū)域，其次為旱地，多呈塊狀分散分布。2010—2020年海口市土地利用變化特征主要表現(xiàn)為建設用地擴張導致農(nóng)業(yè)和生態(tài)用地大幅縮減。10年間建設用地面積占比快速上升，由8.26%增長至13.22%，增長主要發(fā)生在北部城市建成區(qū)南側(cè)，2015年后建設用地逐漸向東西兩側(cè)延伸。在年變化率方面，2015年后建設用地擴張得到一定控制，增速明顯放緩。水田和旱地受嚴格的基本農(nóng)田政策影響，在2015年后減速開始下降，但也使林地和灌木林成為建設用地擴張的主要來源，年減少速率有所加快。此外，濕地受到主城區(qū)西側(cè)沿海地帶港口建設的影響，年減少速率開始上升。

2.2多情景下的?？谑型恋乩妙愋蜁r空變化特征

結(jié)合相關研究經(jīng)驗[6， 17]，根據(jù)?？谑形磥戆l(fā)展規(guī)劃和土地變化特征，設定自然發(fā)展情景（BAU）、生態(tài)保護情景（EC）和經(jīng)濟發(fā)展情景（ED）。

1）自然發(fā)展情景：不考慮任何規(guī)劃政策對土地利用變化的約束性影響，利用Markov模型以2015—2020年土地利用轉(zhuǎn)移概率為基礎，預測2035年的土地利用數(shù)量需求。2）生態(tài)保護情景：依據(jù)《?？谑猩鷳B(tài)保護第十四個五年規(guī)劃》，嚴格控制城市擴張和人為開發(fā)活動對自然生態(tài)的影響，因此降低林地、灌木林、草地、水域和濕地轉(zhuǎn)化為建設用地和耕地的概率，并通過退耕還林、還湖的方式恢復城市主要生態(tài)用地，提高水田和旱地轉(zhuǎn)林地、灌木林和濕地的概率。同時，根據(jù)?？谑猩鷳B(tài)廊道規(guī)劃，將區(qū)域內(nèi)的水域及濕地設置為不可轉(zhuǎn)化區(qū)域。3）經(jīng)濟發(fā)展情景：在平衡經(jīng)濟與生態(tài)過程中，優(yōu)先滿足城市居民的生產(chǎn)、住房和娛樂需求。海口市以第三產(chǎn)業(yè)作為主要經(jīng)濟支柱，與城市建設用地開發(fā)密切相關，因此根據(jù)往年轉(zhuǎn)換情況，提高作為建設用地主要來源的水田、旱地、林地、灌木林、濕地和未利用地轉(zhuǎn)為建設用地的概率。

結(jié)果表明（表5），BAU情景延續(xù)了當前城市發(fā)展的趨勢，建設用地持續(xù)擴張，面積增加10 953.09 hm2，但年增速進一步下降。該情景下農(nóng)田生態(tài)系統(tǒng)得到較好的保護，旱地以0.25%的年增速實現(xiàn)正向增長，與當前《海南自由貿(mào)易港建設總體方案》提出的補充耕地指標儲備庫和保障糧食安全的耕地保護政策相一致。其余地類面積均保持不同程度的下降，年減速介于EC和ED之間，且較2015—2020年有所下降。

EC情景下建設用地增速得到控制，面積增量最小，僅為9 366.75 hm2，同時林地、灌木林和水域等重要生態(tài)用地的減少幅度明顯小于其余2個情景，并且濕地實現(xiàn)正向增長，漲幅達5.26%。該情景下，除由于退耕還林等生態(tài)措施推進而面積減少明顯的旱地外，其余地類雖然面積仍有不同程度的減少，但退化速度明顯放緩，遠低于其余2種情景。

ED情景下建設用地擴張規(guī)模遠高于其余2種情景，面積增加l14 271.93 hm2，但年增速仍低于2010—2020年水平，這表明未來城市空間開發(fā)將趨于飽和。同時，該情景下更大規(guī)模的建設用地擴張也導致各地類的萎縮程度高于其余2個情景，其中林地、旱地和灌木林的減量最高，未利用地、濕地和草地則由于基數(shù)較小而損失速率更高，尤其是未利用地以每年3.4%的速度消失，相關生態(tài)系統(tǒng)將受到嚴重破壞。

3種情景下的海口市生態(tài)系統(tǒng)類型空間的變化模式整體相同（圖3），建設用地的擴張均集中于西北沿海地帶，符合“以港口為中心，沿海岸線進行建設”的濱海城市發(fā)展特征。但區(qū)別于2010—2020年城市空間主要由沿海向內(nèi)陸擴張的模式，至2035年，海口市建設用地開發(fā)呈現(xiàn)出向東西兩翼延展的趨勢，其中主城區(qū)西側(cè)的長流片區(qū)和東側(cè)的江東新區(qū)最為明顯，因此未來這些區(qū)域周邊的土地利用結(jié)構(gòu)變化也更為劇烈。

2.3海口市生態(tài)系統(tǒng)服務價值時空變化特征

2.3.1生態(tài)系統(tǒng)服務價值數(shù)量結(jié)構(gòu)變化

?？谑械纳鷳B(tài)系統(tǒng)服務結(jié)構(gòu)相對穩(wěn)定，其中水文調(diào)節(jié)和氣候調(diào)節(jié)因?？谑写竺娣e的林地和濕地資源，成為最主要的生態(tài)系統(tǒng)服務類型，二者占比超過價值總量的70%。從總量變化上看，生態(tài)系統(tǒng)服務的損耗逐年加重，但2015年后生態(tài)系統(tǒng)服務的損耗速度明顯放緩，這得益于2016年開始落實的生態(tài)文明試驗區(qū)政策，對生態(tài)環(huán)境指標進行了細化和嚴格管控，降低了城市建設用地對其他地類的轉(zhuǎn)換。各類生態(tài)系統(tǒng)服務中以水資源供給的損耗最為劇烈，其次是凈化環(huán)境和氣體調(diào)節(jié)，三者的年損失率遠高于其余生態(tài)系統(tǒng)服務，表現(xiàn)出更高的敏感性和不穩(wěn)定性（圖4）。

預測數(shù)據(jù)顯示（圖5），BAU情景下各項ESV均持續(xù)下降，水文調(diào)節(jié)和水資源供給損失最高，約占損失總量的59%。除糧食生產(chǎn)、水資源供給、氣體調(diào)節(jié)和維持養(yǎng)分循環(huán)的ESV減速低于2010—2020年，其余ESV都表現(xiàn)出進一步下降的趨勢。

EC情景下，各項生態(tài)系統(tǒng)服務均得到改善，以水文調(diào)節(jié)和水資源供給改善最為明顯。水文調(diào)節(jié)實現(xiàn)正向增長，而水資源供給的損失則大幅下降，年損失率由2015—2020年的5.5%縮小至2%，這是該情景ESV高于其他情景的主要原因。由于退耕還林、還湖等生態(tài)修復工程主要依賴耕地資源，因此該情景下旱地的損失量高于其余2種情景，但仍高于?？谑?035年總體規(guī)劃的耕地保有量指標。

ED情景下，除糧食生產(chǎn)外，各項ESV較其余2種情景大幅下降，其中水文調(diào)節(jié)減少占總損失量的41%。除糧食生產(chǎn)、水資源供給和維持養(yǎng)分循環(huán)的ESV減速低于2010—2020年外，其余損失速率大幅上升。一方面是由于建設用地的擴張強化了建設用地的溢出效應，另一方面是削弱了自然生態(tài)系統(tǒng)尤其是草地和濕地生態(tài)系統(tǒng)，導致其提供的服務下降。

2.3.2生態(tài)系統(tǒng)服務價值空間分布變化

在空間分布上，海口市ESV以中部南渡江為分界，大致呈現(xiàn)出西低東高的格局（圖6）。低值區(qū)在城市建成區(qū)，高值區(qū)在北部東寨港自然保護區(qū)以及東南部的農(nóng)林生產(chǎn)保障區(qū)。研究周期內(nèi)低值區(qū)不斷擴大，其變化范圍與建設用地高度吻合，主要集中在新海港臨港片區(qū)和江東新區(qū)。區(qū)別于BAU和ED情景，EC情景下ESV的增長主要集中在南渡江以東、大坡鎮(zhèn)以北的區(qū)域。該區(qū)域森林和濕地增加，因此區(qū)域生態(tài)服務改善顯著。

3 結(jié)論

近年來“國際旅游島”和“海南自貿(mào)港”政策的相繼發(fā)布，帶動了?？谑薪?jīng)濟進入高速發(fā)展時期，同時也加劇了生態(tài)環(huán)境和經(jīng)濟發(fā)展的矛盾沖突。本文對2010—2020年?？谑蠰ULC和ESV的時空變化特征進行研究，并運用PLUS模型與當量因子法模擬2035年3種情景下海口市生態(tài)系統(tǒng)服務的變化發(fā)展，分析了ESV與LULC變化的關系，得到結(jié)果如下：

1）?？谑型恋乩妙愋鸵陨鷳B(tài)用地為主，林地和旱地占比最高。2010—2020年，海口市土地利用變化以建設用地增加和林地、旱地的減少最為顯著，2015年前建設用地增速最快，空間上集中于主城區(qū)周邊，2015年后有所放緩，向西北部長流片區(qū)擴張。2020年主城區(qū)東西兩側(cè)均形成城市擴張熱點，在建設用地擴張下生態(tài)用地類型的面積大幅下降。

2）?？谑猩鷳B(tài)系統(tǒng)服務結(jié)構(gòu)相對穩(wěn)定，水文調(diào)節(jié)和氣候調(diào)節(jié)占據(jù)主導地位，林地和濕地是調(diào)節(jié)服務的主要來源，對維持?？诘纳鷳B(tài)穩(wěn)定性具有重要意義。研究期間，?？谑械腅SV呈持續(xù)下降趨勢，但損失速率逐漸放緩，這與LULC變化速率一致。受建設用地增加影響，水文調(diào)節(jié)和水資源供給服務功能下降最高，表明LULC變化與ESV變化密切相關?？臻g上，ESV低值區(qū)擴散范圍也與建設用地擴張范圍有很高的一致性。

3）在BAU情景下，各項ESV均持續(xù)下降。在ED情景下，水文調(diào)節(jié)大幅下降，進一步表明林地和濕地轉(zhuǎn)變?yōu)榻ㄔO用地有很強的負面生態(tài)系統(tǒng)服務效應。在EC情景下，各項生態(tài)系統(tǒng)服務均得到改善，以水文調(diào)節(jié)和水資源供給改善最為明顯；增值區(qū)主要集中在南渡江以東、大坡鎮(zhèn)以北的農(nóng)林生產(chǎn)保障區(qū)域。進一步表明區(qū)域LULC的結(jié)構(gòu)調(diào)整，尤其是提高濕地占比，能顯著改善區(qū)域整體生態(tài)系統(tǒng)服務。

本文中結(jié)合PLUS模型與當量因子法進行量化和模擬的研究方法，避免了單一模型的缺點，充分發(fā)揮2種方法在數(shù)量預測和空間模擬方面的優(yōu)勢[18]。

沿海城市的沿海地區(qū)通常是土地利用變化最頻繁的地區(qū)，也是生態(tài)系統(tǒng)服務下降最嚴重的地區(qū)[19～20]。?？谑凶鳛榈湫偷难睾！皢魏恕背鞘校磥沓鞘虚_發(fā)將集中在西北部的沿海岸線，形成西北區(qū)域的生態(tài)系統(tǒng)服務價值洼地。研究表明，推動城市由單核向多核轉(zhuǎn)變能夠轉(zhuǎn)移部分產(chǎn)業(yè)和人口，以疏解大城市功能，減輕生態(tài)壓力，有效緩解大城市空間功能的集中分布帶來的生態(tài)系統(tǒng)結(jié)構(gòu)單一等壓力[21]。2021年10月公示的《?？谑袊量臻g總體規(guī)劃（2020—2035）（公眾版）》提出，在國土空間保護開發(fā)策略中按照“東進、西提、南育、北聯(lián)、中優(yōu)”的原則，扭轉(zhuǎn)海口現(xiàn)狀單中心的空間格局，向多組團、網(wǎng)格化的空間結(jié)構(gòu)轉(zhuǎn)變。這樣不僅可以平衡生態(tài)質(zhì)量和區(qū)域經(jīng)濟增長的需求，還可為周邊城鎮(zhèn)提供發(fā)展空間，縮小城鄉(xiāng)發(fā)展差距[21～22]。這與本研究的結(jié)果一致。

對比3種模擬情景，EC、BAU和ED的ESV差異顯著，分別為447 559.2萬元、399 556.3萬元和386 073.6萬元。BAU情景延續(xù)了當前土地政策的發(fā)展趨勢，但ESV損失與ED情景相近，遠高于EC情景，這表明在沒有任何干預的情況下，海口市目前的土地和空間開發(fā)模式仍然不足以平衡經(jīng)濟發(fā)展和生態(tài)保護，需要重新配置土地資源，以促進生態(tài)和經(jīng)濟的協(xié)調(diào)發(fā)展。EC、BAU和ED的建設用地的年變化率分別為1.94%、2.27%和2.96%，EC情景能夠在保持城市擴張速率與BAU情景接近的同時，維持較低的ESV損耗，更符合未來城市綠色發(fā)展的需求。未來ESV損失區(qū)域?qū)⒓性诤？谑形鞅辈康某鞘薪ǔ蓞^(qū)周邊，而改善區(qū)域集中在?？谑袞|南部的農(nóng)林生產(chǎn)保障區(qū)?！逗？谑袊量臻g總體規(guī)劃（2020—2035）（公眾版）》提出，構(gòu)建“兩核三帶、八廊多園”的市域生態(tài)保護格局，其中北部濱海生態(tài)帶和南部生態(tài)綠帶可分別改善城市建成區(qū)北部沿海區(qū)域和西部建成區(qū)的生態(tài)系統(tǒng)服務功能。2035年ESV低值擴大區(qū)主要位于新海港臨港片區(qū)和江東新區(qū)，因此未來還需關注“西提”，即西海岸和美安科技新城建設片區(qū)、“東進”江東新區(qū)，在新城規(guī)劃時充分留足生態(tài)用地空間，采取最優(yōu)的土地利用模式規(guī)避生態(tài)損失。

注：圖1來自海南測繪地理信息局，審圖號為瓊S（2021）023號；其余均由作者自繪。

參考文獻：

[1]COSTANZA R，D'ARGE R，DE GROOT R，et al. The value of the world's ecosystem services and natural capital[J]. Nature，1997，387：253-260.

[2]COSTANZA R，DE GROOT R，BRAAT L，et al. Twenty years of ecosystem services： How far have we come and how far do we still need to go？[J]. Ecosystem Services，2017，28：1-16.

[3]GRIMM N B，F(xiàn)AETH S H，GOLUBIEWSKI N E，et al. Global change and the ecology of cities[J]. Science，2008，319：756-760.

[4]SCHNEIDER A，LOGAN K E，KUCHARIK C J. Impacts of urbanization on ecosystem goods and services in the U.S Corn Belt[J]. Ecosystems，2012，15（4）：519-541.

[5]GAO J，YU Z，WANG L，et al. Suitability of regional development based on ecosystem service benefits and losses： A case study of the Yangtze River Delta urban agglomeration，China[J]. Ecological Indicators，2019，107：105579.

[6]黎佳君，廖秋林，沈守云，等. 長株潭城市群綠心地區(qū)生態(tài)系統(tǒng)服務價值時空變化研究[J]. 中國園林，2022，38（1）：100-105.

[7]張曉瑤，張瀟，李冬花，等. 城市土地利用變化對生態(tài)系統(tǒng)服務價值影響的多情景模擬——以深圳市為例[J]. 生態(tài)學報，2022（6）：1-12.

[8]高星，楊劉婉青，李晨曦，等. 模擬多情景下白洋淀流域土地利用變化及生態(tài)系統(tǒng)服務價值的空間響應[J]. 生態(tài)學報，2021，41（20）：7974-7988.

[9]歐陽曉，賀清云，朱翔. 多情景下模擬城市群土地利用變化對生態(tài)系統(tǒng)服務價值的影響——以長株潭城市群為例[J]. 經(jīng)濟地理，2020，40（1）：93-102.

[10]陳兵飛，廖鐵軍，張莉坤. 生態(tài)紅線約束下萬州區(qū)土地利用情景模擬及生態(tài)價值評估[J]. 水土保持研究，2020，27（5）：349-357，364.

[11]?？谑械胤绞分巨k公室. 2022年?？谀觇b[M]. ?？冢耗虾３霭嫔?，2022.

[12]FU J，ZHANG Q，WANG P，et al. Spatio-Temporal Changes in Ecosystem Service Value and Its Coordinated Development with Economy： A Case Study in Hainan Province， China[J]. Remote Sensing，2022，14（4）：970.

[13]KURNIA A A，RUSTIADI E，F(xiàn)AUZI A，et al. Understanding Industrial Land Development on Rural-Urban Land Transformation of Jakarta Megacity's Outer Suburb[J]. Land，2022，11（5）：670.

[14]YE Y，BRYAN B A，CONNOR J D，et al. Changes in land-use and ecosystem services in the Guangzhou-Foshan Metropolitan Area， China from 1990 to 2010： Implications for sustainability under rapid urbanization[J]. Ecological indicators，2018，93：930-941.

[15]DAI X，JOHNSON B A，LUO P，et al. Estimation of urban ecosystem services value： A case study of Chengdu， Southwestern China[J]. Remote Sensing，2021，13（2）：207.

[16]CHEN J，F(xiàn)U H，CHEN S. Multi-Scenario Simulation and Assessment of Ecosystem Service Value at the City Level from the Perspective of “Production–Living–Ecological” Spaces： A Case Study of Haikou， China[J]. Land，2023，12（5）：1021.

[17]王婷，韋仕川. 2009-2018年?？谑型恋乩棉D(zhuǎn)型特征及其影響因素[J]. 水土保持通報，2022，42（1）：344-352.

[18]SHI M，WU H，F(xiàn)AN X，et al. Trade-offs and synergies of multiple ecosystem services for different land use scenarios in the Yili river valley， China[J]. Sustainability，2021，13（3）：1577.

[19]O’HIGGINS T G，O’DWYER B. A vulnerability framework to protect coastal social ecological systems[J]. Anthropocene Coasts，2019，2（1）：279-292.

[20]POLITI E，PATERSON S K，SCARROTT R，et al. Earth observation applications for coastal sustainability： Potential and challenges for implementation[J]. Anthropocene Coasts，2019，2（1）：306-329.

[21]肖金成. “多中心組團”模式應對“城市病”[N]. 中國社會科學報，2012-04-09（A07）.

[22]史乃聚，段培君. 城市病治理的戰(zhàn)略哲學處方[J]. 智庫理論與實踐，2021，6（3）：45-54.

作者簡介：

付暉/1985年生/女/山東肥城人/博士/海南大學（?？?570228）/副教授/研究方向為景觀生態(tài)規(guī)劃

陳杰/1996年生/男/江蘇南京人/碩士/海南大學（海口 570228）/專業(yè)方向為生態(tài)系統(tǒng)服務評估

黎淑翎/1986年生/女/廣東佛山人/博士/海南大學（?？?570228）/講師/研究方向為中外園林史、城鄉(xiāng)規(guī)劃與土地利用管理、風景園林規(guī)劃管理與法規(guī)