增强还是减弱？城市绿地降温效应的多年量化研究及其关键二维/三维驱动因素的比较

《International Journal of Applied Earth Observation and Geoinformation》：Strengthening or weakening? Multi-year quantification and comparison of urban green space cooling effects and their key 2D/3D drivers

本研究以北京为对象，利用Landsat影像和机器学习模型，量化了1985-2023年间城市绿地（UGS）冷却效果的时空演变，提出三维冷却强度指标（ICI、NCI、OCI）。发现NDVI、平均植被高度和聚合指数是主要驱动因素，传统冷却强度指标显示整体减弱，但新指标表明核心区域冷却增强。研究揭示了时空异质性和多尺度驱动机制，为城市热岛治理提供科学依据。

　　北京作为中国首都，其城市化进程在过去四十年中经历了显著的变化，这些变化不仅影响了城市景观结构，也对城市热岛效应（UHI）和居民的生活质量产生了深远的影响。本文的研究聚焦于城市绿地（UGS）的冷却效应，旨在揭示其在不同时间尺度上的变化趋势以及其背后的驱动因素。通过分析1985年至2023年的数据，研究不仅关注UGS的覆盖变化，还探索了其在空间分布上的演变，并提出了三种新的冷却强度（CI）指标，以更全面地理解UGS的冷却能力。

### 1. 研究背景与意义

随着城市化的加速和全球气候变化的加剧，城市热岛效应的强度和频率预计会显著增加。这种现象对城市居民，尤其是老年人和儿童等脆弱群体，带来了更严重的健康风险，如热浪、空气污染、心血管疾病和哮喘等。城市绿地因其在调节城市微气候、提供生态服务等方面的重要性，被广泛认为是缓解热岛效应的关键手段。然而，尽管已有大量研究探讨了UGS的冷却效应，仍存在一些不足之处，例如对冷却效应的定义不够明确，缺乏对长期变化的深入分析，以及不同研究方法之间的结论差异较大。

因此，本研究选择北京作为案例，该城市经历了快速的城市化过程，为长期分析UGS的冷却效应提供了理想的条件。研究目标包括：分析UGS及其城市热环境在1985至2023年间的时空变化，评估城市是否变得更加绿色或灰暗；提出具有空间特性的UGS冷却效应指标，以判断其冷却能力是否增强或减弱；比较传统逐步回归分析与机器学习模型（如XGBoost）在揭示冷却效应机制方面的表现，从而识别对冷却效率具有重要影响的UGS特征。

### 2. 研究方法与数据来源

本研究主要采用遥感技术，特别是Landsat系列卫星数据，来分析UGS的变化。Landsat数据具有较高的空间分辨率（30米）和良好的时间连续性，自1972年以来一直被用于城市景观和热环境的监测。此外，研究还利用了高分辨率植被结构数据和补充数据，如北京的经济社会统计数据、Google Earth Pro的高分辨率影像以及百度街景数据，以构建更全面的分析框架。

在方法上，研究首先通过计算归一化植被指数（NDVI）来估算UGS的覆盖情况，并利用FVC（植被覆盖度）作为衡量UGS演变的核心指标。FVC与冷却能力之间存在较强的正相关关系，因此被用于评估UGS的冷却效应。其次，通过分析不同时间段的LST（地表温度）数据，研究揭示了城市热环境的变化趋势。LST的获取基于辐射校正方程，并通过将原始数字数（DN）转换为大气顶层光谱辐射亮度（Radiance），再通过转换公式计算出亮温（Tb），最终估算出地表温度（LST）。

为了更深入地理解UGS的冷却效应，研究提出了三种基于网格的冷却指标：内部冷却强度（ICI）、邻近冷却强度（NCI）和外部冷却强度（OCI）。这些指标分别反映了绿地内部、相邻区域以及整个研究区域的冷却效应。此外，研究还采用逐步回归分析和XGBoost机器学习模型，以评估不同UGS特征对冷却强度的影响，并通过SHAP（Shapley Additive exPlanations）分析来解释模型预测结果，从而量化各特征对冷却效应的贡献。

### 3. 研究结果与分析

#### 3.1 UGS的变化趋势

研究结果显示，北京的UGS覆盖度经历了先下降后上升的过程。1985年时，大部分位于二环路内的网格的FVC低于20%，形成了连续的低覆盖区域。然而，到2023年，FVC低于20%的网格比例显著下降，仅占不到1%，表明UGS的分布更加均衡。总体来看，FVC在1985年为0.54，在2005年下降至0.37，但随后又恢复到0.52，呈现出明显的“V”形变化趋势。特别是在二环路内，FVC持续增长，从0.23上升至0.47，显示出城市中心区域的显著绿化。

#### 3.2 LST的空间分布

地表温度（LST）的空间分布与UGS的格局密切相关。在1985年至1995年间，LST主要集中在三环路至四环路之间，形成了连续的热岛区域。然而，自2005年起，LST的分布逐渐分散，形成了多个孤立的热点区域，这可能与UGS的整合以及城市结构的变化有关。总体而言，夏季平均LST呈现上升趋势，但城市热岛强度（SUHI）在1985年为2.94°C，到2023年下降至1.15°C，表明热岛效应有所缓解。然而，这种缓解可能是由于非绿地区域的温度变化，而非绿地区域的降温可能由遮阳效果和高大建筑物的增加所驱动。

#### 3.3 冷却效应的时空变化

传统的冷却强度（CI）指标显示，UGS与非绿地区域之间的温度差持续下降，从1985年的2.65°C降至2023年的1.83°C。这表明，尽管UGS的覆盖度有所恢复，但其与非绿地区域的温差却在缩小，可能意味着冷却能力的减弱。然而，研究通过提出三种新的冷却指标（ICI、NCI、OCI），发现冷却效应在空间分布上呈现出不同的趋势。例如，三环路内的网格显示了明显的冷却增强趋势，而四环路至五环路之间的东部和东南部区域则表现出冷却能力的下降。

研究进一步发现，虽然传统的CI指标显示出冷却能力的下降，但使用更精细的冷却指标（如OCI）时，发现有31%的研究区域表现出冷却能力的增强，而27%的研究区域则显示出冷却能力的减弱。这表明，冷却效应的评估需要结合不同的指标和空间尺度，不能仅凭单一指标得出结论。

#### 3.4 冷却效应的驱动因素

研究通过逐步回归分析和XGBoost模型，揭示了UGS特征对冷却效应的影响。结果显示，NDVI是影响冷却效应的最主要因素，解释了超过50%的CI变化。景观指标和三维特征分别贡献了38%和10%的解释力。结合所有二维和三维特征，可以解释超过70%的CI变化，其中NDVI、平均植被高度和聚集指数（AI）是影响最大的三个特征。

此外，研究发现，NDVI与UGS面积、植被密度和生长活力之间存在正相关关系。这表明，植被的丰富度和健康状况是冷却效应的重要驱动因素。同时，景观配置指标如斑块数量（NP）、斑块密度（PD）和形状指数（SHAPE）对冷却效应的影响也较为显著，显示出空间分布对冷却能力的重要性。

### 4. 研究讨论与启示

#### 4.1 冷却效应的增强或减弱

冷却效应的变化是一个复杂的过程，受到多种因素的影响。传统的CI指标显示冷却能力有所下降，但这可能只是单一指标的局限性所致。通过分析三种新的冷却指标（ICI、NCI、OCI），研究发现，虽然传统的CI指标显示了降温幅度的缩小，但空间尺度上的分析表明，冷却效应在某些区域显著增强。特别是三环路内的区域，显示出明显的冷却增强趋势，而四环路至五环路之间的东部和东南部区域则显示出冷却能力的下降。

因此，判断冷却效应是否增强或减弱，需要结合不同的指标和空间尺度。单一指标可能无法全面反映冷却效应的变化，而多尺度分析能够提供更准确的评估。

#### 4.2 冷却效应的主要驱动因素

研究发现，NDVI是影响冷却效应的最主要因素，这与其反映植被覆盖度、密度和健康状况的能力密切相关。此外，景观配置指标如聚集指数（AI）和形状指数（SHAPE）对冷却效应也有重要影响，显示出空间结构对冷却能力的调节作用。三维特征如植被高度对冷却效应的贡献相对较小，这可能是因为地表温度主要受二维植被特征的影响，而三维结构的变化可能在某些情况下起着辅助作用。

研究还指出，斑块数量（NP）和斑块密度（PD）与LST呈负相关，这意味着分散的绿地配置可能降低冷却效率。因此，在城市规划中，应注重绿地的连贯性和集中性，以最大化其冷却效应。

#### 4.3 对城市规划的启示

UGS的冷却效应对于调节城市微气候具有重要作用。研究结果表明，冷却效率不仅受UGS本身的特性影响，还与周围非绿地区域的特征密切相关。因此，在城市规划中，应综合考虑UGS的组成特征（如NDVI、FVC）、空间配置（如AI、PD）以及三维结构（如植被高度），以实现更有效的冷却效果。

此外，研究还指出，对于特定的UGS区域，其冷却能力可能受到周围高大建筑物的遮阳效果的影响。因此，在规划绿地时，应考虑其与高大建筑的相对位置，以最大化其冷却效应。同时，选择合适的UGS指标和CI指标对于不同空间尺度的分析至关重要，因为不同的指标可能得出不同的结论。

### 5. 研究的局限性与未来方向

尽管本研究取得了重要进展，但仍存在一些局限性。首先，Landsat数据的分辨率限制了对UGS细节的分析，尤其是在早期缺乏高分辨率数据的情况下。其次，LST只能反映特定时间点的地表温度，无法有效捕捉日内变化，特别是夜间变化，这需要未来进一步研究。此外，尽管研究探讨了多种方法，但对冷却效应背后的物理机制的深入分析仍显不足，缺乏模拟和预测研究。

未来的研究应结合多源数据（如高分辨率影像、LiDAR数据）进行长期、系统和精细的监测，以更全面地理解UGS的冷却效应。同时，应进一步探索冷却效应的物理机制，通过模拟和预测研究，为城市规划和气候适应策略提供更坚实的科学基础。

### 6. 研究结论

本研究通过对北京UGS冷却效应的时空变化进行分析，揭示了其在不同时间段和空间尺度上的变化趋势以及主要驱动因素。研究发现，UGS的覆盖度经历了先下降后上升的过程，尽管总体略有减少，但其空间分布更加均衡。三种新的冷却指标（ICI、NCI、OCI）显示，冷却效应在某些区域显著增强，而其他区域则有所减弱。NDVI作为最主要的驱动因素，对冷却效应的解释力超过50%，而景观指标和三维特征分别贡献了38%和10%。结合所有特征，可以解释超过70%的冷却效应变化。

这些发现为设计气候适应性城市和制定基于证据的绿化策略提供了科学依据，特别是在快速城市化地区，合理规划UGS对于缓解热岛效应和提升居民生活质量具有重要意义。