这篇论文聚焦一个过去常被忽略、但对城市热环境研究非常关键的问题：公园不仅能给周边降温，公园内部本身也存在明显的温度梯度。以往大量研究更关注“公园比周边城区冷多少”，也就是公园对外部建成区的冷却溢出效应，但对于从公园边界走向公园内部时，地表温度究竟如何变化、变化有多强、什么样的公园更容易形成稳定冷核，缺乏全球尺度、统一方法的系统定量。论文正是从这一空白切入，提出了一个新的全球分析框架，试图回答“公园内部降温是否普遍存在”以及“哪些景观因子决定了这种内部降温效应”这两个基础问题。
研究对象的尺度非常大：作者基于全球709个城市、58,957个城市公园的数据，覆盖不同大陆和不同柯本-盖格气候带，用Landsat反演地表温度，并结合NDVI、土地覆盖、再分析气候数据，构建了一个“内部缓冲区—函数拟合”的分析框架，专门刻画公园边界到公园内部的温度梯度特征。论文最有价值的地方，不只是证明了全球多数公园都具有“冷核”特征，更进一步识别出关键阈值：当公园面积超过4.65 ha、园内平均NDVI超过0.60、树木覆盖率超过0.59时，公园内部降温会明显增强。这使论文从“现象描述”迈向了“规划可用”的层面。
从研究意义上看，这篇论文的贡献不在于简单重复“绿色空间有降温作用”这一共识，而在于把“公园内部热环境异质性”系统拉到了全球尺度，并将其转化为可量化、可比较、可解释的指标体系。也就是说，它不再只把公园看成一个整体冷斑，而是把公园视为具有内部结构、边界效应和景观阈值响应的复杂热环境单元。这对城市公园设计、树种配置、空间布局优化，以及气候韧性城市建设都具有直接参考价值。
在全球变暖和快速城市化叠加的背景下，城市热环境恶化已成为全球城市面临的共同挑战。城市公园作为典型的自然基础解决方案，被普遍认为能够通过植被遮阴、蒸散和地表能量再分配缓解热压力，因此在城市降温研究中长期占据重要位置。现有研究已经较充分地揭示了公园相对于周边建成区的冷却效应，也表明这种冷却存在显著的区域差异、季节差异和气候带差异。
但问题在于，过去很多研究都把公园当成一个“整体对象”来处理，更多比较“园内平均温度”和“园外平均温度”的差别，而忽略了公园内部本身并不是均一的。实际上，公园内部不同功能分区、不同植被结构、不同硬化程度，会造成明显的微气候差异。已有少数局地研究已经发现，从公园边缘向内部移动时，温度通常会逐步下降；在高树冠覆盖、植被连续性更好的区域，温度往往更低，而广场、建筑和道路等不透水区域则可能形成局部热斑。然而，这些成果主要停留在个别公园、个别城市层面，研究方法也不统一，导致目前仍缺乏一个全球尺度上可比较的认知框架。
因此，这篇论文提出的核心背景逻辑非常清楚：如果我们真正想理解公园是如何发挥降温作用的，就不能只看公园对外部的影响，还必须看公园内部温度是如何组织起来的。 换句话说，公园并不是一个简单的“冷块体”，而是一个存在从边缘到内部热梯度变化的空间系统。只有识别这种内部温度梯度的空间格局及其驱动机制，才能真正为公园设计提供更精确的依据。
🌍研究目标与研究区
📌研究目标
这项研究围绕两个核心问题展开。第一，如何在全球尺度上系统量化公园内部的降温效应，并揭示其空间分布格局？ 第二，哪些因素主导了这种公园内部降温效应，这些因素是否存在明确的阈值响应？ 这两个问题一个解决“怎么看”，一个解决“为什么”，共同构成了论文的主线。
📌研究区
在研究区设置上，作者选择了全球709个城市中的58,957个城市公园作为样本，覆盖热带、干旱、温带和寒带四大主要气候区，并分布于亚洲、非洲、欧洲、北美洲、南美洲和大洋洲六大洲。城市选择基于2020年城市建成区面积大于200 km²的标准，公园面积控制在2–100 ha之间，以兼顾遥感分辨率适配性、样本规模和结果的可比性。这个尺度设计很关键，因为它既避免了过小公园在30 m热红外分辨率下难以稳定识别内部梯度的问题，也避免了超大公园因功能过于复杂而带来的额外偏差。
从研究对象特征上看，这篇论文并不是研究“全球公园总体是否降温”这么宽泛的问题，而是更具体地锁定在城市内部、公园边界到内部这一段空间过程。因此，它的研究区不是以国家或城市为主要分析单元，而是以公园本身为核心单元，再进一步聚合到城市、气候区和洲际尺度进行比较。这种多层级设计使得论文既能识别全球总体规律，也能保留不同地理背景下的差异性。
🖼️ 图 1. 全球709 个城市及城市公园的全球格局和统计数据。城市及其对应气候带的空间分布 (a)；以中国北京为例，展示城区和城市公园 (b)，并附局部放大图 (c)；各气候带和各大洲的城市数量 (d) 和城市公园数量 (e)；不同面积区间（2–100 公顷）的城市公园数量，左侧 y 轴表示面积为 2–20 公顷的公园数量，右侧 y 轴表示面积为 20–100 公顷的公园数量 (f)。
🌍 研究方法
· 数据与方法体系
在方法上，这篇论文最核心的创新，是提出了一个“buffer-function fitting”框架，也就是“内部缓冲区—函数拟合”方法，用来描述公园内部从边界到内部的LST梯度。作者首先从公园边界向内，以10 m间隔连续生成多个内部缓冲带，计算每个缓冲带的平均地表温度；然后用线性函数和指数函数分别对“距离—温度”关系进行拟合，选择拟合效果更好的函数作为该公园的最终温度梯度模型。最后，将最外层边界与最内层缓冲区的拟合温差定义为IPCIs（intra-park surface cooling intensity），即公园内部地表降温强度。若IPCIs大于0，说明公园内部比边界更冷，属于cool-core park；反之则属于warm-core park。
数据方面，研究主要使用了2018–2022年夏季Landsat-8 LST与NDVI数据，并结合GLDAS空气温度和风速、ERA5-Land降水数据来表征背景气候；同时引入Global Urban Boundary界定城市边界，使用OpenStreetMap提取公园边界，并基于ESA WorldCover 2020获取树木、灌丛、草地和建成地比例等土地覆盖信息。也就是说，作者把影响公园内部降温的因素划分为四大类：公园形态、公园内部景观、公园周边缓冲区景观、背景气候。其中共纳入了15个潜在驱动因子，如公园面积、形状指数、园内NDVI、树木覆盖率、周边建成地比例、气温、降水和风速等。
在驱动机制分析上，作者进一步结合LightGBM + SHAP开展解释性机器学习分析。一方面，LightGBM用于处理全球大样本、多变量、非线性关系，识别各因子的相对贡献；另一方面，SHAP用来揭示每个变量在不同取值区间内对IPCIs的促进或抑制作用，并识别其“拐点”或“阈值”。为了检验稳健性，作者还使用XGBoost进行了对照验证。这使得整篇论文不仅停留在统计相关层面，而是能够更清晰地呈现“哪个因子重要、重要到什么程度、超过什么值后作用会显著增强”的完整逻辑链。
🖼️ 图 2. 全球城市内部公园表面冷却效应两步法的分析框架。
🌍 研究结果
论文首先从“是否存在冷核”这一层面出发，发现全球约74%的城市公园属于cool-core parks，其共同特征是地表温度从公园边界向内部逐步下降，说明“公园内部比边缘更冷”是一个相当普遍的全球现象；相对地，warm-core parks约占26%。而且这种冷核现象在不同气候区都有体现，其中干旱区和寒带城市中cool-core park比例更高，分别达到79%和77%，说明在干冷背景下，公园内部热梯度更容易形成。
进一步分析cool-core parks的降温强度后，作者发现全球城市尺度上的平均IPCIs达到1.15 ± 0.42 K。这意味着，对绝大多数公园而言，从边界走向最内部区域，地表温度平均会下降约1 K以上，而且这种内部降温在全球尺度上具有较高稳健性：96%的城市平均IPCIs超过0.5 K的不确定性阈值，约75%的单个公园也超过这一阈值。分气候区看，寒带城市公园的内部降温最强，其次是干旱区，温带和热带相对较弱。也就是说，背景气候确实会调节公园内部冷却能力，但这种作用并未改变“公园内部普遍存在降温梯度”的总体事实。
在温度梯度形态上，作者发现68%的cool-core parks呈指数衰减型，即从边界进入公园后，温度最初下降得较快，随后逐步趋于平缓。这表明公园内部降温并不是线性均匀发生的，而是具有明显的边界—内部过渡特征：边界附近受外部建成环境干扰较强，而越进入内部，微气候越稳定，但降温增益也会逐渐减弱。这个结果很重要，因为它说明公园内部降温并不是简单“越往里越一样冷”，而是存在一个具有边际递减特征的空间过程。
在驱动因素方面，研究表明，公园形态和公园内部景观是决定IPCIs的两大主导维度，其贡献率分别为33%和29%；相比之下，周边缓冲区景观和背景气候的贡献相对较小，均为19%。如果看单个变量，公园面积是最重要的因子，占30%，其次是园内平均NDVI（15%）和树木覆盖率（14%）。这说明真正决定公园内部是否形成强冷核的，不只是“有没有绿”，而是空间尺度、植被整体状态和树冠结构是否足够支撑稳定冷却过程。
更进一步，论文识别出了三个尤其关键的阈值：当公园面积超过4.65 ha时，其对IPCIs的作用由限制转为促进；当园内NDVI超过0.60、树木覆盖率超过0.59时，内部降温明显增强。与此同时，公园面积与IPCIs之间呈现出明显的先快速增长、后逐步饱和的关系，在约20 ha附近接近潜在最大值，说明面积扩张存在边际收益递减；而形状越不规则，内部降温表现反而越弱。对于周边环境而言，较高的周边建成地比例会放大园内外热对比，从而增强IPCIs；反之，若周边植被也非常丰富，则园内外温差缩小，内部降温强度会有所减弱。
🖼️ 图 3. 709个全球城市中冷核公园和暖核公园的比例。城市中冷核公园 (a) 和暖核公园 (c) 比例的空间分布；箱线图显示了全球不同气候带和大陆的城市中两种公园类型比例的变化 (b 和 d)。
🖼️ 图 4. 全球709 个城市冷核心公园的IPCI值和模式。全球城市层面IPCI值的空间分布(a)；箱线图显示全球不同气候带和大陆城市间IPCI值的变化(b)；内部地表温度梯度呈指数型 (R² > 0.8)、线性型 (R² > 0.8) 和其他模式（指数型和线性型 R² < 0.8）的公园比例 (c)；具有不同数量内部缓冲区（颜色越深，缓冲区数量越多）的公园在指数型 (d) 和线性型 (e) 中的平均地表温度梯度拟合曲线；这两种模式公园的IPCI值分布 (f)。
🖼️ 图 5. 基于中国北京城市公园的园内地表降温效应的空间尺度敏感性分析。小型公园（2–3 公顷）与所有公园（2–100 公顷）分别展示。比较基于 10 米和 30 米缓冲区得到的园内地表降温指数（IPCIs）（a）；两种缓冲区尺度下，具有指数型（R² > 0.8）、线性型（R² > 0.8）和其他降温模式（指数型和线性型 R² < 0.8）的冷核公园比例（b）；基于两种缓冲区尺度得到的大型公园（c）和小型公园（d）的代表性园内地表温度梯度曲线。
🖼️ 图 6. 全球冷核公园关键驱动因素与 IPCI 之间的响应关系。蓝色散点表示每个因素区间内的平均IPCI值，红色曲线表示拟合的函数关系。
🖼️ 图 7. 影响各气候带IPCI的关键因素的 SHAP 依赖性图。蓝色散点表示因素值及其对应的 SHAP 值，红色曲线表示使用 GAM 模型拟合的平滑 SHAP 值。临界点表示 SHAP 值在负值和正值之间转换的因素值，对应于 GAM 曲线上的 SHAP 值为零。
🖼️ 图 8. 展示了 709 个全球城市（涵盖不同气候带和大陆）中四种公园热效应类型（冷核降温、冷核升温、暖核升温、暖核降温）的比例变化。降温标签代表公园周边净热效应，定义为公园的平均地表温度低于（降温）还是高于（升温）等面积周边缓冲区。该缓冲区采用与驱动因素分析相同的方法构建。
🧩 讨论部分
这篇论文最重要的讨论价值，在于它重新定义了我们理解城市公园降温效应的方式。过去很多研究默认“公园是一个整体冷岛”，而本文表明，更准确的理解应该是：公园内部其实存在从边缘到核心的连续降温梯度，且这种梯度在全球大多数公园中都真实存在。 这意味着公园降温并不是一个单一均值问题，而是一个与边界效应、内部结构和植被配置密切相关的空间过程。换句话说，公园能否真正发挥热缓解作用，不只是看“有没有公园”，更要看这个公园是否形成了稳定的内部冷核。
从机制解释上看，研究强调了面积与植被结构的协同作用。面积越大，越有可能削弱边界热干扰，为内部形成稳定微气候提供空间基础；但面积的作用并不是无限增强的，当公园扩大到一定程度后，新增面积带来的温度收益会逐渐减弱。因此，论文并没有简单得出“公园越大越好”的结论，而是指出在城市土地资源紧张的现实下，适度扩大规模 + 提高植被质量与树冠覆盖，往往比单纯追求大面积更有效。尤其是NDVI和树木覆盖率在不同尺度分组中都保持较强作用，这表明高质量绿化结构比单纯面积扩张更具持续稳定的降温意义。
论文还给出了很强的规划启示。它所识别的4.65 ha、NDVI 0.60、树木覆盖率0.59并不是绝对意义上的统一设计标准，但它们为全球城市公园设计提供了一个非常有价值的“经验阈值区间”。尤其对于高密度建成区的小微公园建设而言，这些阈值提示我们：若空间受限，至少要通过提高植被连续性、树冠配置和内部绿量质量，尽可能推动公园越过“有效冷核形成”的门槛。也就是说，这项研究真正提供的不是抽象结论，而是具有可操作性的设计导向。
当然，这项研究也存在一定局限。首先，Landsat热红外的有效空间分辨率仍然有限，虽然作者采用10 m内部缓冲区进行曲线采样，但底层LST仍来自30 m像元，因此在较小公园中，相邻缓冲区之间可能存在像元重叠与重采样影响。其次，研究依赖OSM和全球土地覆盖产品，不同地区数据精度并不完全一致；此外，背景气候变量被重采样到30 m仅用于配准，并不代表真实具备30 m尺度细节。换句话说，这篇论文已经尽可能在全球尺度上做到统一和稳健，但其结果仍应理解为一种全球宏观规律总结，而不是对每一座具体公园微气候过程的完全替代。
🧩 研究结论
本研究采用结合内部缓冲区和函数拟合的方法，量化了全球709个城市超过5万个城市公园的园内地表降温效应，并进一步考察了四个主要因素的相对贡献和阈值效应：公园形态、公园内部景观、周边缓冲区景观和背景气候。研究结果总结如下：
首先，全球74%的城市公园属于冷核公园，几乎是暖核公园（26%）的三倍。全球冷核公园的平均IPCI值达到1.15 ± 0.42 K，其中干旱和寒冷气候区的冷核公园分布最广，降温强度也最强。其次，全球68%的冷核公园呈现非线性指数衰减模式，其特征是边界附近地表温度迅速下降，然后向内部逐渐降低。14%的公园呈现线性衰减模式，其余公园则表现出不均匀的温度梯度。第三，全球冷核和暖核公园的形成主要受公园内部景观（占40%）和周边缓冲区景观（占37%）的影响。对于冷核公园而言，IPCI s主要受公园形态（33%）和内部景观（29%）驱动，其中最主要的因素是公园面积（30%）、平均 NDVI（15%）和树木覆盖率（14%）。当这些因素分别超过 4.65 公顷、0.60 和 0.59 的阈值时，会增强降温效果，但不同规模的公园和气候带之间存在轻微差异。
这些发现强调了系统研究公园内部降温梯度的重要性，并填补了公园降温效应研究中的一个重要知识空白。本研究提出的定量框架和响应阈值为城市公园设计和气候适应型城市规划提供了切实可行的指导。
🛰️🌆GeoAI 星球 · 洞见时刻
如果把这篇论文放到更大的GeoAI与城市气候研究语境下看，它的启发其实非常明显。过去我们谈城市绿色基础设施时，常常停留在“面积更多、分布更广、总量更高”的层面，但这篇文章提醒我们：热环境响应本质上是空间组织问题，而不只是面积统计问题。 公园内部边界—核心的梯度结构、植被配置质量以及周边建成背景，共同决定了公园降温是否真正有效。也就是说，未来的城市热环境治理，需要从“有多少绿”进一步转向“这些绿是如何组织起来的”。
