Evaluation of landscape pattern evolution and spatial configuration of urban park green space in Haikou City based on SEDMs model

海口市(19°31′~20°04′N，110°07′~110°42′E)地处海南岛最北端，南接定安，北邻琼州海峡，包含秀英区、龙华区、琼山区和美兰区4个区，面积为3 126.82 km²。年均降水量为1 669 mm，年均日照时数为2 000 h以上，年均气温为23.8 ℃，属于热带海洋性季风气候。地形较为平缓，以滨海平原和台地为主。

海口市城市总体范围为椰海大道以北，粤海大道以东和南渡江以西连接北边海岸线，形成闭环式的建成区边界(图1)，包含综合公园、专题公园、游园等266个公园绿地，城市公园绿地面积为3 233.63 hm²，其中，建成区的公园绿地总面积为645.90 hm²。北部以世纪公园、白沙门公园、万绿园、五源河湿地公园等与南渡江沿岸周边绿地主体，南部以美舍河凤翔湿地公园、红城湖公园等构成建成区南部绿地景观空间。

Figure 1.  Boundary range of municipal district in the built-up area of Haikou

遥感数据包括QuickBird影像(2008年11月、2013年11月)、Worldview影像(2013年12月、2017年12月)和Landsat 8 OLI_影像(2020年1月)。在ENVI处理平台完成遥感数据的辐射定标和大气校正，并提取公园绿地。同时收集了海口市矢量路网、小区、公园和路口及2020年的兴趣点(point of interest，POI)数据。

在相关研究^[18]的基础上，引入斑块密度(PD)、景观破碎度(SPLIT)、景观蔓延度(CONTAG)、景观分离度(DIVISION)等景观格局指数，定量描述2008—2020年海口市建成区公园绿地景观格局时空演变特征。

基于ArcGIS网络分析模块，对建成区居民到达最邻近目的地的最小距离以及公园服务范围进行定量测算，从而进行可达性配置评价^[19−20]。利用海口市POI矢量数据模拟海口市建成区环境，综合考虑研究区域基础设施、交通、环境和阻力因素等，评估海口市民到达城市公园的便捷程度，探究演变过程中各城市公园之间的空间连通性。

由于海口市高温多雨的气候条件和错综复杂的热带植被状况，本研究选取11月至翌年3月数据计算NDVI^[21]，以此时计算的NDVI^[22]均值作为全年的NDVI均值，以此评估市民在一定范围内绿地数量的供给水平。利用NDVI对植被区域敏感的特征来反映市民的绿地可获得性水平。

SEDMs是用于监测景观格局变化的土地利用模型^[23−24]。该模型综合了景观模式时空动态分析模型(Spatiotemporal Dynamic Analysis of Landscape Patterns)的特点和景观模式变化类型，可以更准确地评估景观格局过程变化结果和土地利用变化过程。

斑块数(L₁)、斑块周长(L₂)、斑块面积(L₃)和平均斑块分形维度(L₄)是描述景观格局时空变化的量化指标，通过识别斑块空间形态差异能有效获取景观破碎化、景观聚合等景观格局演变特征的反馈信息^[24]。因此本研究选取以上4个景观指数作为评判指标模拟海口市城市公园绿地的景观格局变化。

景观状态变化类型分为景观破碎、景观扩张和景观聚合3种形式。SEDMs模型将这3种形式分为6种景观状态作为景观结构演变的评判依据。其中，景观破碎化包括分离型(D)、收缩型(S)和穿孔型(P) 3种景观状态，表示在自然或人工驱动下斑块被分割成更小的、独立的斑块；景观扩张包括创建型(C)和扩张型(E) 2种景观状态，表示景观格局中不同斑块类型斑块面积的增加；景观聚合仅有聚合型(A) 1种景观状态，表示同种类型斑块团聚，斑块连通性和面积增加。将统计的评判指标生成数据集，根据相邻调查时间段4个景观格局指数的差值(∆L_i)评判斑块类型。

式(1)~(3)中：a、b表示所选取的任意调查时间段，i表示第i个评判指标，ai、bi表示在a或b调查时间段第i个指标，L_ai、L_bi表示在a或b时间段对应的指标值，X_a、X_b则表示从X_a到X_b随机2个时间段斑块数、斑块周长、斑块面积和平均斑块分形维度4个景观格局指数的集合，T_S为相邻的调查时间段，S、P、D、E、A、C分别表示6种景观破碎化状态，N表示调查时间段内没有产生格局演变。

从表1可见：2008—2020年，7种公园绿地类型的斑块数、斑块面积、斑块周长、平均斑块分形维度的最佳像元为52 m×52 m。根据评判指标的变化表明：①城市综合公园斑块周长在研究期间增长了19.99 m，但平均面积下降了3.73 hm²，其格局呈现破碎化趋势。②社区公园斑块周长从801.59 m骤降至576.93 m，其他评判指标无明显变化，说明其格局逐渐简单规则化，且至少发生过1次景观格局演变。③历史名园斑块平均面积从1.81 hm²降至1.49 hm²，再增长到1.91 hm²；分维数和斑块周长在2013年达到峰值，随后有所降低，斑块形状趋于规则化。④其他专类公园斑块数年均增加8.33个，斑块分形维数数值则从1.06降至1.04，表明海口市建成区专类公园类型趋于单一、简单化。⑤遗址公园斑块面积由3.05 hm²降至1.64 hm²，周长增长1.25倍，表明景观破碎化现象明显。⑥游乐公园与遗址公园相似，呈现景观破碎化趋势。⑦游园的斑块数从35个增加到130个，但斑块面积从1.71 hm²降至0.84 hm²，景观破碎化程度逐渐加深。

2008—2020年海口市建成区各类型城市公园的景观格局总体上呈破碎化趋势，大部分城市公园处于穿孔型(图2)，部分城市公园呈现空间聚集化或景观扩张趋势。城市综合公园的景观格局处于分离型(∆L₁＞0，∆L₂＜0，∆L₄=0)，其格局演变过程可能受城市交通或者园林道路建设的影响；社区公园格局发生了多次演变，由分离型转变为扩张型(∆L₁=0，∆L₂＞0，∆L₄＜0)，最后演变为穿孔型(∆L₁=0，∆L₂＜0，∆L₃＞0，∆L₄≥0)，很大程度上归因于外界的人为干扰因素，演变过程较为复杂，致使其他性质的用地类型逐渐转变成社会公园建设；历史名园景观格局由原来的穿孔型转变为扩张型，表明海口市园林建设一定程度上促进了历史名园绿地面积的增长；其他专类公园空间格局连通性良好，绿地面积持续增长，景观格局出现了明显的景观聚合现象，处于聚合型(∆L₁＜0，∆L₂≥0，∆L₄＜0)；遗址公园、游乐公园、游园3种公园类型12 a间均未发生明显的景观格局变化，属于穿孔型的景观格局。

Figure 2.  Evolution characteristics of landscape pattern in the built-up area of Haikou from 2008 to 2020

由表2可知：海口市建成区公园绿地景观破碎度指数和景观蔓延度指数都呈增长趋势，2008—2020年分别增长了1.990和2.143，其中2013年蔓延度最低，而2020年蔓延度最高。表明公园绿地景观异质性在一定程度上增加，但优势斑块(其他专类公园、综合公园)之间具有良好的空间连接度，总体趋势向高集聚方向发展。

从表3可见：综合公园的景观分离度指数呈先升后降的趋势，2017年为0.991，而2020年下降为0.988。同时，其他专类公园在2013年达0.947，随后有所降低。总体来看，各类型公园的整体集成度较低，空间联系差。由于设计规划策略和政策的影响，在2013年前后，根据《海口市城市绿线管理办法》《海口市城市绿地绿线控制规划(第一批)》等相关规定和园林绿地建设的持续推进，海口市公园绿地建设获得“井喷式”发展，但整体空间分布规律分散程度较高，绿地空间结构单元复杂化，各类型绿地斑块碎片化。2017年，为响应“创建全国文明城市和国家卫生城市”攻坚战，提出了“城市更新”“生态城市”等工作战略部署，园林空间景观品质获得显著成效，绿地空间整体分布由离散式向连续集中式布局发展。2020年综合公园的规范化管理，有效保证了各综合公园之间的空间连续性，但部分专类公园景观分离度仍较高，需进一步加强公园管理措施。

将2008—2020年NDVI平均值、园林绿地面积、绿化覆盖率、公园数、公园面积与人均公园绿地面积等进行纵向比较，量化海口市公园绿地的可获得性水平。从表4和表5可见：2020年，建成区的绿化覆盖率和人均公园绿地均高于园林城市标准值(绿化覆盖率≥36.0%，人均公园绿地面积≥7.5 m²·人⁻¹)，城市园林绿地面积从3 558 hm²增加到8 008 hm²，呈线性增长趋势，近几年海口市建成区人均公园绿地面积为12.3~12.5 m²·人⁻¹，但绿化覆盖率在2015年达到峰值(42.9%)后下降趋势明显。7种类型公园的面积占比发生变化，综合公园、专类公园和游园3种类型公园面积占所有公园面积的比例从2008年的96.55%增加到2020年的97.04%，海口建成区公园建设的重心正在朝这3类公园偏移。同时，海口市建成区NDVI在12 a间发生了剧烈变化。虽然建成区公园面积由628 hm²增加到了1 062 hm²，但海口市建成区绿地植被面积在2013年绿化程度仍处于较低水平，在2017—2020年回升并保持稳步增长，总体趋势逐渐向高绿化水平发展。其中，NDVI高的区域主要分布在万绿园、世纪公园、白沙门公园等处于近郊或者绿地面积较广的区域，NDVI低的区域分布在人民公园、海瑞墓等处于市中心或者城市下垫面硬质面积较多的区域。

图3为Bigemap全能版地图提供的海口市POI位置。

Figure 3.  City profile simulated by POI location in the built-up area of Haikou in 2020

如图4所示：以距离为阻抗属性构建服务区，分析居民生活圈与公园设施点之间的网络拓扑关系，对服务区分析求解后获取公园所覆盖的服务范围，并探究小区居民到达附近公园设施点的便捷程度。建成区的公园服务范围主要覆盖了海口市秀英区西海岸、美兰区海甸区域、龙华区北部区域、琼山区东部区域。对于位于其他区域小区位置的居民来说，至少需要行驶1 km以上或步行时长超过15 min才能到达附近的公园。

Figure 4.  Analysis of green space service area of Haikou City Park in 2020　　　　　　　

运用ArcGIS以距离为阻抗属性构建最近设施点分析POI数据(图5)，综合考虑学校区域、路口、单行道等不利因素的影响，分析公园(事件点)与小区(设施点)之间的拓扑关系，对设施点分析求解后获取小区到公园的最佳路径。建成区东部尤其是东北部的小区到公园的最佳路径较短，且公园类型多、数量多，代表公园有世纪公园、万绿园、金牛岭公园等。因此它们的辐射范围更广，景观与人居环境空间连通性强，便捷程度较高。建成区西部的小区到公园的最佳路径较长，且公园数量少，主要位于近城区或者处于湿地保护区，景观空间布局简单分散，便捷程度较低。

Figure 5.  Analysis of the closest facility in the green space of Haikou City Park in 2020

海口市建成区城市公园形成景观破碎化整体空间格局，在局部区域产生集聚式、组团式格局演变，这一结果与李杨^[25]、雷金睿等^[26]的研究结果一致。本研究推测，道路建设和区域规划的进一步推进是导致海口市建成区景观破碎化的主要原因之一。2017年，《海口市总体规划(空间类2015—2030)》的编制和公园建设，在推进划定生态保护红线和城市开发边界基础上，完善中心城区空间布局的调整，使建成区局部区域保持一定的空间连续性。“孤岛式”的公园体系被逐步淘汰，“空间整体性连续性”是未来城市公园绿地发展的重要目标。因此，要在建设绿地生态体系的同时确保各类型公园绿地均衡发展，平衡好景观价值、服务价值和历史价值之间的关系。

受研究区域尺度制约，本研究用网格划分景观单元分析海口市建成区整体区域，更有利把控整体环境，也避免小型绿地对整体趋势的干扰。然而从宏观角度分析，仅使用4个年份的遥感数据对小尺度景观格局的趋势分析会存在一定误差，使用长时间多时相高精度的遥感图像，景观格局的演变趋势会更具科学性和严谨性。此外，研究过程中只发现了4种景观状态类型(分离、扩张、聚合、穿孔)和2种演变(分离→扩张→穿孔、穿孔→扩张)，可能是其他景观状态在研究区域内极少发生，也可能是其他景观状态类型对公园绿地演变特征不敏感。

基于此，海口市城市绿地规划应有效针对城市公园绿地景观单元分散和空间连接度差等问题寻求空间布局最优解。首先，对于龙华区和美兰区可达性好的区域，园林建设应突破传统的公园设计范式，保证在城市绿地系统规划和国土空间规划下充分利用好土地资源和公园的景观价值，进行景观提升和空间优化吸引周边居民，丰富居民的人居环境和绿地供给水平；对于琼山区和秀英区中公园可达性较差的区域，应在充分考虑游客和市民短途旅行或漫步休憩需求的基础上，在建成区辖区范围内的适宜区位增设新的城市公园，以保证公园可获得性水平和建成区绿地空间公平。

基于SEDMs模型和遥感影像研究分析，2008—2020年海口市各类园林绿化指标虽然总体保持线性增长，但各类型绿地之间发展不平衡，总体上呈现景观破碎化的趋势。海口市建成区公园绿地景观异质性明显增加，整体空间布局由全局离散复杂化趋于景观破碎化和局部斑块规则化，绿地景观空间结构由离散非连续布局向局部中心集聚式布局发展，总体趋势向高集聚空间模式发展，绿地植被总体趋势向高质量绿化模式发展。