基于土地利用模拟预测模型分析的城市绿色空间发展多情景模拟及建设时序研究——以湛江市中心城区为例

赵霁雨,崔 柳,*,王 佳,陈 思

1 北京林业大学园林学院, 北京 100083

2 北京林业大学林学院, 北京 100083

3 德克萨斯农工大学地理系, 德克萨斯州 77843

国土空间是由自然地理环境和人类社会活动两个子系统相互作用构成的人地耦合系统,是自然地理演化过程和经济社会发展过程的综合体[1]。建立国土空间规划体系,实现“多规合一”是党中央和国务院在新时代作出的重大部署,其核心思想是生态文明建设和空间治理现代化[2]。明确了资源空间“三区三线”的区域划定及城市未来空间发展的三类基本型:城镇空间、农业空间、生态空间,而国土空间规划即对该三类空间进行资源体系管理与城市机制体系整合的重要举措[3—4]。三生空间格局优化有利于提升生态环境及社会经济效益,有助于明确区域发展重点与方向,三生空间合理布局是城市空间规划的发展目标[5]。

城市绿色空间 (urban green spaces,UGS) 是新时代背景下国土空间总体格局优化战略中的重要项目,对“三区”空间内容均有涉略,是城市社会-经济-自然复合生态系统中的重要组成[6—7]。近年来,城市绿色空间发展潜力备受研究者青睐,相关学者从城市绿色空间功能[8—9]、城市绿色空间结构与指标[10—11]、城市绿色空间优化策略[12—13]等方面入手开展了大量研究。表明城市绿色空间在提高城市环境质量、增强地域性绿色空间弹性、促进可持续生活方式等方面具有潜力、在稳定城市生态与环境、建构城市形态与秩序、供给高密度城区游憩服务等方面具有不可替代的价值[14]。这些价值可概括为城市生态系统服务 (ecosystem services,ES),包括供应服务、调节服务和支持服务、文化服务。城市绿色空间结构及发展模式与城市生态系统服务价值体现休戚相关[15]。生态系统服务研究的最终目标是辅助决策者优化生态保护规划与管理措施,以促进人类社会与自然环境的可持续发展[16]。湛江市中心城区城市绿色空间在发展前一阶段中与建设用地冲突明显,被挤占严重,支持服务价值和供给服务价值较低,同时《湛江市国土空间总体规划(2020—2035年)》对中心城区绿色空间保护及建设目标较其他地区更为严格,故本文选取湛江市中心城区绿色空间为研究区域。

城市生态服务系统之间的权衡或协同关系可通过不同土地利用方式和利用强度表征显现[17—18]。探究如何保持合理高效的土地利用结构并发挥其效益最大化[19],分析主观政策指导下发展得到的土地利用结构是否可以达到预期,对后续建设重点区域及时序问题提供理论支持具有重要意义[20]。目前,我国针对土地利用变化的研究多样化,研究集中于土地利用变化过程、驱动因素分析[21—22]、土地利用模拟预测[23—24]、耦合模型评估生态特征[25]等方面。土地利用变化研究常用的模型主要有Markov模型、CA模型、系统动力学模型、GeoSOS模型、CLUE-S模型等,现有研究表明,Markov模型可以进行长期预测但不能描述空间尺度上的变化;系统动力学模型可以反映复杂系统之间的相互作用但很难处理空间信息;传统CA模型没有限制元胞状态变化的模块,只能模拟一种土地利用类型的变化[26];CLUE-S模型的非空间模块需要借助独立的数学模型进行计算,忽略了非优势地类转化的可能性[27];GeoSOS模型在识别非线性转换规则的空间变异性方面存在局限性[28]。

本研究使用高性能空间计算智能实验室(HPSCIL)[29]开发的PLUS模型(Patch-Generating Land Use Simulation Model),其耦合土地利用扩张分析策略(LEAS)和多类型随机斑块种子的CA模型(CARS),对土地利用变化机制的理解力提升,对模拟自然土地利用类型及其与建设用地冲突边界的土地利用动态变化模拟具有优势[30—31],且随机森林算法具有时间属性,同时集成了规划交通更新机制和规划开发区内的随机种子机制,填补了已有研究只能考虑规划的约束作用(保护区、禁建区)、无法考虑规划政策的驱动和引导作用的问题[32]。以湛江市中心城区两期(2010年和2020年)土地利用数据为基础数据,模拟研究2035年自然发展情景与国土空间规划政策情景下研究区土地利用动态变化,对比各情景生态系统服务价值(ESV),并针对研究区域ESV冷热点分空间聚集度分析,合理评估实际效益,以此为依据进行城市绿色空间核心区域有效识别及空间层级规划,为湛江中心城区绿色空间开发时序提出建议,以希为规划政策管理及资源配置提供数据支持与策略依据。

1 研究区域概况与数据来源

1.1 研究区域

湛江市位于中国大陆南端、广东省西南部,介于东经109°40′—110°58′,北纬20°13′—21°57′之间,位于粤、桂、琼三省(区)邻接处,又处于东亚、东南亚、我国大西南等三大地区的交汇地带,是大小双三角区位的核心[33]。湛江市三面环海,是全国少有的内海湾型城市,陆地大部分由半岛和岛屿组成,地势多为平原和台地,呈现南北高而中间低,中轴高而东西低的整体趋势。

《湛江市国土空间总体规划(2020—2035年)》包括市域和中心城区两个规划层次,中心城区包括霞山区、赤坎区、开发区乐华街道和泉庄街道、麻章区麻章镇和湖光镇,以及坡头区部分地区,陆域国土空间628.39 km2。以促进一体化生态格局为目标,国土空间规划对中心城区绿色空间保护及建设目标较其他地区更为严格,《湛江市碧道建设总体规划(2020—2035年)》对空间总体格局规划如图1所示,空间类型面积排序为城镇空间>生态空间>农业空间。本研究关注湛江市中心城区范围内的绿色空间发展潜力,通过生态系统服务价值评估及人类活动强度计算,耦合研究区域空间热点分析,为城市绿色空间建设的空间落位及时序安排提出建议,以期通过空间优化最大限度提升城市空间生态系统服务价值。

图1 研究区域国土空间总体格局Fig.1 Study of the overall regional territorial space pattern

1.2 数据来源

本研究所采用的数据中,2010年及2020年土地利用数据来源于遥感信息处理研究所[34],土地利用依据《土地利用现状分类》(GB/T 21010—2017)分为6类:耕地、林地、草地、水域、建设用地、未利用地。空间分辨率为30m×30m;DEM高程数据来源于地理空间数据云平台(http://www.gscloud.cn),空间分辨率为30m×30m;平均温度数据来源于WorldClim数据网站[35],降水量数据来源于国家科技基础条件平台—国家地球系统科学数据中心(http://www.geodata.cn),空间分辨率为1km×1km;现状道路、铁路、河流水域等矢量数据来源于OpenStreetMap;人口数据来源于美国宇航局地球观测系统数据和信息系统(EOSDIS)中的社会经济数据应用中心(SEDAC)(https://sedac.ciesin.columbia.edu);经济数据来源于中国科学院资源环境科学与数据中心[36],空间分辨率为1km;其他社会经济数据来自《湛江统计年鉴2020》、《2020全国农产品成本收益资料汇编》。以上数据运用ArcGIS进行投影变换、裁剪、重采样等工具进行数据预处理转化为30m×30m的栅格数据。

2 研究方法

研究过程如图2,提取规划政策条件,分析规划政策情景下城市绿色空间的发展趋势及结果,进行生态系统服务价值量化及空间落位与评估,计算研究区内人类活动强度在空间上的聚集程度,进而对生态系统服务价值空间分布及人类活动强度空间分布叠置分析,得出规划政策下湛江市中心城区绿色空间发展建设优先级分区,为规划政策提供空间政策的量化数据基础,为其可行性实施、时序安排及预期结果提供数据支撑与建议。

图2 技术路线Fig.2 Technology route

2.1 生态系统服务价值量化

当量因子法是目前最常用的生态系统服务量化方法,最早由Costanza等[37]提出。本文以谢高地等建立的“单位面积生态系统服务价值当量表”[38]为基础,建立研究区的原始生态系统服务价值当量表(表1),具体方法为:耕地取用旱地和水田的平均值;林地取83.9%阔叶、16.1%针叶[39]、针阔混交及灌木的平均值;草地取草原、灌草丛、草甸的平均值;水域取湿地和水系的平均值;未利用地用荒漠代表;建设用地计为0。通过公式(1)和公式(2)估算2020年实际情景、2035年自然发展情景、2035年规划政策情景3期生态系统服务价值。

表1 研究区原始单位面积生态系统服务价值当量Table 1 Original ecosystem service value equivalent per unit area in the study area

(1)

(2)

式中,ESV为生态系统服务价值,单位为元/年;i为土地利用类型;Ai为第i类土地利用类型的面积,单位为hm2;VCi为第i类土地利用类型单位面积生态系统服务价值,单位为元hm-2a-1;j为生态系统服务类型;ECj为某类土地利用类型第j项生态系统服务价值当量;k为生态系统服务类型数量;Ea为1单位生态系统服务价值的经济价值,单位为元hm-2a-1。

结合实地情况,运用公式(3)修正单位生态系统服务的经济价值[40],计算出湛江市地均生态系统服务价值及修正后的粮食产量的经济价值约为1884.53元 hm-2a-1,

(3)

式中,Ea为1单位生态系统服务的经济价值,单位为元hm-2a-1;i为粮食作物种类;m为第i种粮食作物在研究区内的平均价格,单位为元/kg;pi为第i种粮食作物的单产,单位为kg/hm2;qi为第i种粮食作物的种植面积,单位为hm2;M为粮食作物的总种植面积,单位为hm2。

2.2 土地利用格局情景模拟

2.2.1PLUS模型

由高性能空间计算智能实验室(HPSCIL)开发的斑块生成土地利用变化模拟模型PLUS(Patch-generating Land Use Simulation Model)包括用地扩张分析策略(LEAS)和基于多类随机斑块种子的CA模型(CARS)两部分,描述通过经济、生态和社会系统之间空间动态的相互作用而产生的土地利用变化状态,可以更好地支持规划政策以实现可持续发展[29]。

LEAS自两期土地利用变化中提取各类用地扩张部分,并从中采样,采用随机森林算法逐一挖掘各类土地利用扩张和驱动力的因素,获取各土地利用类型发展概率以及各驱动因子对该时段各土地利用类型扩张的贡献。CARS模块是基于多类随机斑块种子元胞自动机(CA)模型,并对其基于随机种子生成和阈值递减机制进行了改进[41—42]。

2.2.2驱动因子的选取

本研究以2期(2010年与2020年)土地利用数据为基础,综合已有研究成果[43—44],从研究区实际情况、数据可获取性、驱动因子可量化性三方面进行考虑,选取驱动土地利用变化的因素共计15项(图3)作驱动因子输入LEAS模块,过程计算各驱动因素对每种土地利用类型扩张的贡献程度,得到各土地利用类型在各像元上的发展概率(图4),经精度检验验证其可作为后续载入CARS模块的前期计算结果。

图3 土地利用变化驱动因素Fig.3 Driving factors of land use change

图4 各土地利用类型的发展概率Fig.4 Growth probability for each land use type

2.2.3土地利用模拟情景预设

研究通过多情景分析探究不同可能、对比不同路径的发展状态,为决策提供依据。模型预测通过输入不同驱动因素、空间政策的限制、设置土地利用需求、转移矩阵和邻域权重等训练参数来设置两种发展情景,模拟未来土地利用变化。一是延续2010年至2020年发展规律的自然发展情景;二是将国土空间规划、地区政策作为情景原则与限制条件输入模拟之中的规划政策情景(表2)。

表2 湛江市中心城区未来土地利用发展情景模式及原则Table 2 Development scenario and explanation of driving factors of land use change in the central urban area of Zhanjiang

2.3 人类活动强度计算

人类活动强度 (HAILS, Human Activity Intensity of Land Surface)从土地利用与人类活动关系角度客观反映人类活动对陆地表层影响和作用程度[45—46]。湛江市中心城区人类活动强度计算采用徐勇[47]等提出的人类活动强度计算模型,区域内建设用地当量面积占区域土地总面积的百分比即为该区域的人类活动强度。

(4)

(5)

式中,HAILS为陆地表层人类活动强度;SCLE为建设用地当量面积;S为区域总面积;SLi为第i种土地利用面积;CIi为第i种土地利用/覆被类型的建设用地当量折算系数;n为区域内土地利用数量。

2.4 热点分析

热点分析 (Getis-Ord Gi*)用于识别具有显著性的高值(热点)和低值(冷点)的空间聚类[48],热点分析可以直观反映湛江市中心城区生态系统服务价值是否存在高值集聚和低值集聚[49],以及确定热点区、次热点区、过渡区、次冷点区、冷点区在空间上发生聚类的位置[50],为后续判断城市绿色空间建设开发时序提出建议。

3 结果与分析

3.1 湛江市中心城区未来土地利用变化模拟结果

3.1.1模型模拟精度检验

模型验证是判断模型参数设置是否合理以及在研究区是否适用的标准,Kappa系数是一种衡量精度的指标,它可以定量反应模拟效果,也可以用来评价两幅图的一致性[51]。计算公式为:

(6)

式中:P0表示模拟结果与实际之间的相对一致性;Pe表示偶然一致性的假设概率。通常Kappa系数取值范围为[0,1],数值越大表示模拟的精度越高,模拟结果与实际的土地利用分布越接近。

以2010年土地利用数据为基础进行2020年土地利用变化模拟,将模拟结果与实际土地利用状况进行精度检验,Kappa系数为0.815,模拟预测结果与实际土地利用类型分布较为吻合,说明可用PLUS模型及相关参数设置预测未来湛江市中心城区土地利用变化。

3.1.2多情景下湛江市中心城区未来土地利用变化模拟结果分析

以2020年土地利用数据为基础,进行自然发展情景和规划政策情景下2035年土地利用变化模拟(图5)。从空间分布来看,自然发展情景下2035年研究区内建设用地范围延续此前发展规律向周边扩张,面积为122.00km2,与2020年相比增加了11.78km2,相应地侵占了耕地及水域面积,耕地面积小幅减少,相比2020年减少了2.71%。西南部局部水域减少明显,相比2020年减少了23.72%;而在规划政策情景下2035年研究区建设用地多围绕海域建设,面积小幅减少为106.16 km2。中心城区中部范围内耕地、林地增多,整体水域得到有效保护,面积减少甚微。

图5 不同情景下2035年湛江市中心城区土地利用变化模拟Fig.5 Simulation of land use change in the central urban area of Zhanjiang in 2035 under different scenarios

对模拟结果分别进行转移矩阵分析(表3、4),发现其中耕地、林地、水域及建设用地转移面积差异较大。自然发展情景延续前一阶段土地利用变化规律和发展模式,在2020—2035阶段,建设用地面积大幅增加,多为耕地及水域向建设用地转换,导致耕地及水域面积减少较多;规划政策情景下,因加入了划定生态保护红线及水域保护、保护耕地等规划限制条件,模拟结果建设用地呈减少趋势,耕地、林地增多,整体水域得到有效保护。由于规划交通发展,建设用地向东南方向转移。

表3 自然发展情景下2020—2035年湛江市中心城区土地利用类型转移矩阵/km2Table 3 Transfer matrix of land use types in the central urban area of Zhanjiang from 2020 to 2035 under the natural development scenario

表4 规划政策情景下2020—2035年湛江市中心城区土地利用类型转移矩阵/km2Table 4 Transfer matrix of land use types in the central urban area of Zhanjiang from 2020 to 2035 under the planning policy scenario

3.2 湛江市中心城区生态系统服务价值特征分析

3.2.1多情景下湛江市中心城区生态系统服务价值变化特征分析

由表5可知,湛江市中心城区的生态系统服务价值主要来自水域,贡献率达55%以上,其次是耕地及林地,分别占比30%左右和15%左右,草地和未利用地贡献率较小。自然发展情景下除林地外,ESV整体呈大幅下降趋势,从11.88亿元降低至10.53亿元,总体生态系统服务价值明显下降,降低率为11.41%,主要原因是对生态系统服务价值贡献率较高的水域、耕地与林地的减少以及不利于生态系统服务价值增加的建设用地的扩张;规划政策情景下2035年湛江市中心城区生态系统服务价值呈上升趋势,从11.88亿元升至12.22亿元。相对于自然发展情景,水域和耕地生态系统服务价值高了1.75亿元,总量增加了1.69亿元,这主要与《湛江市国土空间总体规划(2020—2035年)》中规划政策协调了各类用地面积,划定生态保护红线、水域保护等约束条件有关,保护了研究区的生态安全。但由于林地面积减少,林地生态系统服务价值降低了0.06亿元,故还应对林地予以进一步约束,控制林地向耕地转换。

表5 模拟2035年多情景下各类土地生态系统服务价值/万元Table 5 Simulate the service value of various land ecosystems under the multi scenario in 2035

由图6可得,自然发展情景下2035年生态系统服务价值提升区域较少,分布分散,下降区域较多,相对集中于麻章区内部、赤坎区与霞山区中城市建设用地范围内;规划政策情景下2035年生态系统服务价值总量提升且提升区域面积相较于自然发展情景显著增加,赤坎区出现集中成片提升区域。

图6 不同情景下2035年与2020年生态系统服务价值变化图Fig.6 Change in the value of ecosystem services in 2035 and 2020 under different scenarios

3.2.2多情景下湛江市中心城区生态系统服务价值冷热点分析

创建1km×1km渔网单元,利用ArcGIS软件中热点分析工具分析得出2020年和多情景下2035年研究区的生态系统服务价值的冷热点区域分布图,采用自然断点法将研究区生态系统服务价值空间聚集程度从低到高划分为5个等级:冷点区、次冷点区、过渡区、次热点区及热点区。(图7)

图7 不同情景下2035年湛江市中心城区ESV冷热点分布Fig.7 Distribution of ESV cold and hot spots in the central urban area of Zhanjiang in 2035 under different scenarios

研究区ESV值的冷点区域分布较广且集中成片,热点区域分布较为分散。各情景下的热点及次热点空间分布范围:2035年规划政策情景>2020年实际情景>2035年自然发展情景。自然发展情景下2035年研究区生态系统服务价值高值在空间上的聚集情况有所减弱,西部、南部及东北部冷点区域范围扩大;规划政策情景下2035年研究区南侧热点区域连片聚集,中部自然保护地生态系统服务价值得以延续,北部建设用地范围内大量冷点区域向过渡区及热点区转换,东北部近海域过渡区向原冷点区扩张。这是由于相应规划政策要求,新增热点区的形成是由于水域保护、自然保护地保护等措施的实施。

3.3 城市绿色空间发展优先级分区分析

3.3.1人类活动强度冷热点分析

通过公式(4、5)计算得出规划政策情景下2035年人类活动强度指数,将人类活动强度数据连接至前文所建立1km×1km的渔网单元中,应用ArcGIS中热点分析工具,利用自然断点法将人类活动强度指数从低到高分为冷点区、次冷点区、过渡区、次热点区、热点区五个等级,得到规划政策情景下2035年人类活动强度冷热点分布图(图8)。湛江市中心城区人类活动强度空间分布呈现东西低,中部高的特征,以建设用地为中心发散减弱,近海域区域人类活动强度最弱。

图8 规划政策情景下2035年人类活动强度冷热点分布Fig.8 Distribution of cold hotspots of human activity intensity in 2035 under planning policy scenarios

3.3.2城市绿色空间发展优先级分区结果

统筹考虑自然与人文,以平衡生态系统服务价值及发展建设间的关系为出发点,对湛江市中心城区城市绿色空间发展优先级进行分区规划。通过对生态系统服务价值与人类活动强度的空间分布进行叠置分析,将ESV及HAILS空间聚集程度由冷点区至热点区分别赋值为0、1、2、3、4,根据矩阵(表6)可得到计算分值为0—8的25类特征空间,将25类空间根据0—2、3—5、6—8进行三级分区,分为先序发展区、中序发展区和后序发展区,得到发展优先级分区结果(图9)。其中先序发展区指生态系统服务价值高且人类活动强度较大的区域,但现状研究区内并没有分值为8的网格单元;中序发展区指生态系统服务价值及人类活动强度综合处于过渡阶段的区域;后序发展区指生态系统服务价值及人类活动强度均较弱的区域,在研究区中面积占比最大。

表6 各网格单元ESV及HAILS空间分布矩阵分析Table 6 Spatial distribution matrix analysis of ESV and HAILS for each grid cell

图9 城市绿色空间发展优先级分区结果Fig.9 Results of urban green space development priority zoning

先序发展区聚集于赤坎区中部,此区域生态系统服务价值提升潜力最大,并处于城市建设用地范围内,人类活动强度大,城市绿色空间承载压力较大,应重点修复原有绿色空间基底,控制建设用地在此区域内的扩张程度,强调其综合效益;中序发展区为原有绿色基底较好同时人类活动强度较低的区域与人类活动强度较高但生态服务价值较弱的区域,主要包含赤坎区、霞山区及坡头区的建设用地范围及水域范围,较集中于湛江市中心城区中部,应强调水域保护,城区中部耕地保护,围绕先序发展区周边逐步展开建设,强化先序发展区生态效益的同时向周边扩散至后序发展区,形成发展区过渡;后序发展区面积最大,人类活动强度相对较小,用地类型以耕地为主, 主要分布于湛江市中心城区的边缘范围内。通过城市绿色空间发展优先级分区,将规划政策进行空间响应,基于地理单元考虑生态系统的完整性、资源环境的综合承载能力、人地关系耦合的整体性,逐步形成环境友好,自然人文统筹发展的城市绿色空间。

4 讨论与结论

4.1 讨论

本文通过城市绿色空间的分布变化和利用强度表征显现城市生态系统之间的权衡或协同关系,对推进中国城市绿色空间发展与生态文明建设具有借鉴意义。在模型选择方面优化于此前研究常用的Markov模型、CLUE-S模型[21—22,52]等,应用PLUS模型对各情景下土地利用状况中城市绿色空间的地理分布及价值量化进行了分析,综合考虑自然条件变化与人类活动的影响,生成了真实性更高的土地利用格局。耦合人类活动强度及生态系统服务价值进行综合分析及空间响应,更直观地反映了规划政策情景下发展结果的实际效益。进而研究得到基于生态系统服务价值的城市发展模式建议:自然发展情景下2035年研究结果表明在经济快速发展情景下,生态平衡遭到了破坏,而湛江市中心城区建设用地面积较大,且位于陆地与海域过渡段,对市域生态结构具有重要意义,属于湛江市中部的生态安全屏障,因此湛江市中心城区发展要避免蔓延扩张形式,注意林地保护,限制建设用地扩张;规划政策情景下2035年生态系统的调节服务价值升高显著,主要原因在于赤坎区林地显著增多、建设用地面积小幅减少至106.16km2,耕地增多,水域得到有效保护,通过成片建设用地局部转换至未利用地、建设用地与耕地、林地统筹发展来促进经济、生态共赢的合理化模式。为有效优化城市绿色空间结构,为空间规划决策和可持续的城市发展与生态平衡提供重要的指导意义。

本文主要有以下局限性:(1)采用当量因子法进行ESV评估,虽然对ESV系数进行了基于粮食价格的修订,但研究以假设2020年价值量不变进行预测,存在高估或低估未来ESV变化的可能。(2)为表征人类活动强度及ESV空间聚集程度,选择1km网格的方式进行研究,尺度较为单一,未来可采用多尺度对比分析,探索尺度效应对ESV聚集程度演变的影响。(3)现状城市绿色空间发展滞后于建设空间,成为了制约城市可持续发展的主要障碍[53—54],研究分析虽得出了城市绿色空间发展时序建议,但对其质量及数量是否可以满足人类生活需求这一方面没有落实到指标上。湛江市中心城区作为湛江市本轮国土空间规划发展的重中之重,有必要继续探讨在高质量发展背景下如何将国土空间格局优化措施精细化,统筹考虑以人文及生态服务两方面的可持续联系[55]。

4.2 结论

(1)PLUS模型适用于湛江市中心城区,对未来城市土地利用变化模拟预测结果精度较高,Kappa系数为0.815,表明未来与本区域相关研究可继续考虑应用该模型。在对湛江市中心城区进行不同情景模拟后发现,自然发展情景下2035年研究区建设用地面积为122km2,相比2020年增加了11.78km2,向周边扩张明显,水域及耕地面积减少;规划政策情景下建设用地面积小幅减少至106.16km2,耕地、林地增多,水域得到有效保护。规划政策限制林地转换建设用地,未限制林地向耕地转换导致林地面积减少,应对林地予以进一步约束。

(2)2035年规划政策情景ESV>2020年实际情景ESV>2035年自然发展情景ESV,分别约为12.22亿元、11.89亿元、10.53亿元,规划政策情景下生态系统服务价值总量较自然发展情景下提升约1.69亿元。证明国土空间规划的实施将作为一种精细化治理服务的领域化手段,引导湛江市中心城区生态系统得到有效保护。赤坎区、霞山区ESV升高较显著,而麻章区东南部、经开区及坡头区沿海域部分有所下降,两极变化趋势需要加以重视。规划政策情景下2035年ESV总体提升,其中支持服务和调节服务提升最为显著,文化服务略有下降,在未来建设过程中,切勿忽视景观美学的营造。

(3)从ESV热点分析结果来看其在空间上的聚集程度,热点区域沿海域聚集性较强;对生态系统服务及人类活动强度的空间分布进行叠置分析,进行绿色空间发展优先级分区,建议以赤坎区中部、沿海域建设用地及水域范围、其他区域顺序进行分区建设。