城市与城郊木本冠层对反复干旱的光谱响应：一项基于匹配三城市的先导研究

《Plants》：Spectral Responses of Urban and Peri-Urban Woody Canopy to Recurrent Drought: A Matched Three-City Pilot

　　

城市条件可能加剧木本植被的干旱胁迫，尽管管理和冠层组成可以改变城市与城郊之间的差异。本研究利用2013至2024年间的协调型Landsat与Sentinel-2观测数据（Harmonized Landsat–Sentinel-2, HLS），比较了柏林、马德里和维也纳三个城市中持续存在的木本冠层的光谱异常。研究人员共匹配了70个城市-城郊配对，这些配对在海拔、长期气温、降水量和保留的木本面积上保持一致。在预先定义的26个候选反复干旱事件中，有12个满足了80%的观测支持标准。研究结果以配对内城市减去城郊异常的中位数呈现，指标包括双波段增强型植被指数（two-band enhanced vegetation index, EVI2）和归一化差值水分指数（normalized difference moisture index, NDMI）。两个指数均产生了七次正对比和五次负对比；EVI2的中位数为0.00047，NDMI的中位数为0.00065。针对事件水平的模型在调整了城市和当前干旱严重程度后，检验了前期干旱严重程度和事件间间隔的影响。标准化预测因子的系数范围在指数单位上为-0.00225至0.00252。所有四个异方差一致性95%置信区间均包含零。同样的模式也出现在配对内重采样、空间缓冲以及基于生长季严重程度的模型中，每个区间的置信区间均跨越零值。由于主要回归模型包含12个事件和五个拟合项，仅剩七个残差自由度，导致估计不精确且统计功效较低。所有26个候选事件的气候前期条件均可识别，而是否纳入分析取决于植被观测支持情况。在所保留的案例中，研究人员既未发现一致的城市光谱劣势，也未发现对比度与干旱历史之间存在明确关系。该匹配设计是可复制的，但其提供的证据来自一个基于观测可用性筛选的小样本，并且仅涵盖在整个研究期间持续存在的密集冠层。因此，这些估计不应推广至整个城市森林，也不应被解释为生理响应、死亡率或冠层损失的度量。

论文解读：城市与城郊木本冠层对反复干旱的光谱响应——一项基于匹配三城市的先导研究

一、研究背景与科学问题

全球气候变化背景下，干旱事件的频率和强度不断增加，对城市生态系统构成严峻挑战。城市木本植被不仅承受区域性的水文气候异常，还同时面临城市热岛效应、根系生长空间受限、土壤性质改变、生境破碎化以及不均匀的人工管理等多重压力。理论上，这些城市特有的环境条件可能加剧植被对干旱的敏感性。然而，另一方面，灌溉、物种选择、局部微气候调节或养护措施也可能抵消甚至逆转预期的城市-城郊差异。因此，城市与周边自然或半自然环境中的木本植被在面对干旱时的实际表现差异，是一个尚未被充分解答的科学问题。

以往的研究多集中于单一地点或特定树种，缺乏跨城市、可比性强的系统性评估。此外，多数研究关注单次干旱事件后的恢复力或抵抗力，而对反复干旱事件序列如何影响植被响应，特别是前期干旱事件与当前事件之间的相互作用机制，探讨尚不充分。为此，本研究选取欧洲三个气候与环境背景迥异的城市——柏林、马德里和维也纳作为比较案例，旨在回答以下核心问题：（1）在反复干旱事件中，城市与城郊持续存在的木本冠层光谱异常如何变化？（2）城市减去城郊的光谱异常是否与前期干旱严重程度或事件间隔时间存在一致关联？（3）这种关联在不同时间聚合方式、事件间隔定义、测量重采样、空间缓冲以及生长季对齐的严重程度下是否保持稳定？

二、研究方法与技术路线

本研究采用了一套严谨的多阶段技术方案，主要包括：

1. 1.
   研究设计与匹配策略：基于全球人类住区层（Global Human Settlement Layer, GHS-SMOD）1990年和2020年的数据，将研究单元划分为城市（稠密城市集群或城市中心）和城郊（郊区或城郊网格单元）。在每个功能城市区域内，以半径1公里的圆形区域为候选单元，采用混合整数线性规划进行匹配。匹配协变量包括海拔、1981-2010年平均气温、1981-2010年年降水量以及保留的木本面积，要求配对间各协变量的绝对标准化均值差低于0.10，以确保城市与城郊单元在非城市化因素上的可比性。
2. 2.
   持续木本冠层队列构建：利用Hansen全球森林变化数据集和ESA WorldCover 2020土地覆盖产品，筛选出Hansen树覆盖率≥70%、WorldCover树木分数≥0.70且在2001-2024年间无记录损失的高密度木本像素，形成持续存在的冠层队列。每个1公里单元至少需保留1.8公顷的木本冠层。
3. 3.
   干旱事件识别：基于标准化降水蒸散指数（Standardized Precipitation Evapotranspiration Index, SPEI），采用六个月累积窗口（SPEI-6）从1961年至2024年的气候记录中识别干旱事件。一个合格事件需至少连续两个月SPEI严格低于-1，且至少有一个月落在设定的生长季内。事件严重程度定义为所有SPEI<-1月份中(-1 - SPEI_t)的总和。前期事件指在当前事件之前发生的最近一次合格干旱事件。
4. 4.
   遥感数据处理与光谱异常计算：整合2013-2024年的协调型Landsat与Sentinel-2（HLS）观测数据，计算双波段增强型植被指数（EVI2）和归一化差值水分指数（NDMI）。为每个事件、单元、像素和指数分别建立参考模型（采用谐波模型，无趋势或含线性趋势），通过剔除当前事件年份及可能的后续年份数据进行训练。每日异常值为观测值减去参考模型预测值。年度异常值通过对生长季内各月异常值等权平均获得。
5. 5.
   统计分析：采用普通最小二乘回归模型，以事件水平的中位数城市-城郊光谱对比度为因变量，以前期干旱严重程度（或事件间隔时间的对数）和当前干旱严重程度为核心预测变量，并以城市为固定效应。通过异方差一致性标准误（HC3）计算95%置信区间，并进行多种稳健性分析，包括事件删除、城市删除、配对内自举重采样、30米内部缓冲区分析、生长季对齐严重程度、固定参考模型以及传感器组成调整等。

三、主要研究结果

2.1 分析样本与匹配结果

在所有26个候选当前干旱事件中，每个事件均在完整气候记录中找到了紧邻的前期事件。最终，14个事件因未能满足80%的观测支持标准而被排除，剩余12个分析事件：柏林2个、马德里5个、维也纳5个。这三个城市共贡献了140个匹配单元、70个匹配对、16,721个持续木本像素和约1,399.96公顷的精确裁剪木本面积。所有四个匹配协变量在匹配后的绝对标准化均值差均低于预设的0.10阈值，表明匹配效果良好。

2.2 城市与城郊光谱异常

在等月聚合条件下，EVI2显示出七次正对比和五次负对比（中位数0.00047；范围-0.00955至0.01133）。NDMI同样显示了七次正对比和五次负对比（中位数0.00065；范围-0.01054至0.01660）。城市和城郊各自的异常值往往相似，尽管它们之间的差值在不同事件间大小和方向均有变化。这表明在城市与城郊之间并未观察到一致的、方向固定的光谱差异。

2.3 干旱序列关联

所有四个主要模型的估计值相对于其不确定性区间都非常小。对于前期干旱严重程度，一个标准差（SD）的增加对应于EVI2对比度变化-0.00225，NDMI对比度变化0.00180。对于事件间隔的对数，一个SD的增加对应于EVI2变化0.00071，NDMI变化0.00252。所有四个HC3 95%置信区间均包含零。NDMI的前期严重程度估计在所有删除拟合中保持了相同的符号，而其他三个估计在至少一次删除拟合中改变了符号。这些结果表明，在现有数据和模型框架下，无法确定前期干旱特征与当前城市-城郊光谱对比度之间存在可靠的关联。

2.4 时间定义与事件间隔的敏感性

采用不同时间聚合方式（等有效天数、采集权重、受影响月份）的分析，以及要求事件间隔至少六个月的分析，所得结果均与全样本分析的方向一致，但所有HC3置信区间依然包含零。随着事件数量的减少，估计精度和统计功效进一步下降。

2.5 扩展稳健性与敏感性分析

配对内自举重采样、30米内部缓冲区分析以及生长季对齐严重程度分析的结果均显示，所有焦点系数的置信区间都围绕零分布，未发现稳定的方向性证据。个别事件在Cook's距离诊断中被标记为具有较大影响，但没有事件被移除。

2.6 设计限制与审计结果

匹配后的配对质心距离在各城市间有所差异。固定参考模型和传感器组成调整的检查增加了不确定性，而非提供了共同方向的证据。所有相关的HC3区间均跨越零。持续冠层代表了基线木本面积的56.7%至93.4%，表明研究结论仅适用于幸存下来的高密度木本斑块。

四、讨论与结论

本研究的核心发现是，在来自柏林、马德里和维也纳的12个受观测支持的反复干旱事件中，城市与城郊持续木本冠层的光谱对比度既有正值也有负值，其中位数接近零。因此，研究人员未能在所保留的案例中发现城市冠层存在一致的光谱劣势。但这并不意味着城市与城郊的干旱敏感性相等。由于主要模型仅有七个残差自由度，估计精度和统计功效都很低，合理的效果可能在正负两个方向上延伸。

研究强调了几个重要的方法论和解释性限制。首先，匹配设计虽减少了区域气候差异，但无法控制物种组成、灌溉、管理、土壤条件等未被测量的因素。其次，研究仅分析了持续存在的密集冠层，不涉及树木死亡率或冠层损失。第三，干旱严重程度和事件间隔的前期效应并未与观测到的对比度表现出一致关联。第四，研究包含了大量测量值（数千个像素），但独立的暴露历史数量很少（12个城市事件），这限制了推断的普遍性。

研究结论明确指出：这项研究提供了一个可重复的匹配设计和透明的推断限制。然而，它并未确立城市与城郊的等价性、生理性干旱记忆、干旱序列效应的缺失或整个城市森林的健康与死亡动态。更广泛的推论需要更多城市和城市事件历史、对依赖关系和基于观测的选择进行明确处理，以及关于物种、管理、土壤、局部热量、根区水分和冠层损失的直接信息。