《Ecosphere》：Integrating multi-decadal satellite imagery and species distribution models for conservation planning in oak woodlands
编辑推荐：
橡树林是加州保护最薄弱但生态最重要的生态系统之一，为多种动植物提供关键栖息地。这些基础植物群落日益受到干旱、野火及其他气候驱动干扰的威胁，这些干扰已导致更新和成年树木存活率的显著下降。理解橡树林韧性与持续存在的时空动态，对于评估生态系统在不断变化条件下的响应至
橡树林是加州保护最薄弱但生态最重要的生态系统之一，为多种动植物提供关键栖息地。这些基础植物群落日益受到干旱、野火及其他气候驱动干扰的威胁，这些干扰已导致更新和成年树木存活率的显著下降。理解橡树林韧性与持续存在的时空动态，对于评估生态系统在不断变化条件下的响应至关重要。本研究提出了一种新方法，将基于多年代卫星时间序列的时间混合残差分析得出的植被韧性指标，与预测高排放情景下21世纪中叶蓝橡（*Quercus douglasii*）栖息地适宜性的物种分布模型（Species Distribution Models，SDMs）相结合。研究人员将该方法应用于加州大自然保护协会（The Nature Conservancy）管理的两个研究区域——兰德尔保护区（Randall Preserve）和戴伊溪保护区（Dye Creek Preserve），生成了结合历史植被格局与预测气候适宜性的空间显式地图。研究发现，历史上具有韧性且预计到21世纪中叶气候适宜的区域通常位于海拔较高、温度较低且冬季降水较多的地区，而低适宜性区域则集中在较温暖的低海拔地区。为测试该方法在广泛使用的卫星平台上的适用性，研究人员比较了30米Landsat和10米Sentinel-2数据的结果。这些比较显示出大致相似的模式，其中Sentinel-2捕捉到了更大的时间变异性。在兰德尔保护区18个蓝橡更新样地进行的实地验证证实，高更新位点与未来栖息地适宜性高且历史上植被普遍具有韧性的区域相对应。该方法提供了一种可扩展、可转移的方式来描述历史和预测的橡树林动态，为气候变化下的土地管理、恢复和土地保护规划提供支持。
**整合多年代卫星影像与物种分布模型的橡树林保护规划：研究解读**
**一、研究背景与科学问题**
加州植物区系省（California Floristic Province，CFP）是全球公认的生物多样性热点区域，拥有超过6000种维管植物。以橡树为主的林地（oak-dominated woodlands）和稀树草原（savannas）覆盖约400万公顷，是该区域的基础生态系统（foundational ecosystems），为大量动植物提供栖息地。然而，这些生态系统正面临干旱、气温升高、火情格局改变和土地转换等多重威胁，导致橡树死亡率上升、更新失败（recruitment failure）和栖息地丧失。尽管具有重要的生态价值，橡树林却是加州保护最薄弱的生态系统之一，预计到2040年将有超过30万公顷面临城市发展风险。
在此背景下，理解橡树林如何响应过去的环境压力，以及哪些区域未来可能免受气候变化影响，对于适应性保护规划至关重要。生态学中，韧性（resilience）指系统在干扰后重组并恢复其结构、过程和功能的能力，而抵抗力（resistance）则指系统在环境扰动后保持不变的程度。卫星遥感技术的进步使得研究人员能够在广阔的空间和时间尺度上监测植被动态，但传统分析方法多依赖影像对或少量影像，难以揭示长期动态。时间混合分析（temporal mixture analysis）方法利用完整的卫星观测档案，可刻画跨数十年的植被动态，识别偏离稳定年际格局的位置，从而作为历史植被韧性的指示指标。
另一方面，相关物种分布模型（Species Distribution Models，SDMs）通过关联物种出现记录与环境变量（如气候、地形、土壤），生成栖息地适宜性的空间显式地图，可用于预测未来栖息地动态。尽管SDMs存在气候不平衡、采样偏差和缺乏扩散动态等局限性，它仍是保护研究和规划中广泛使用的工具。然而，历史植被动态评估与未来栖息地适宜性预测通常被独立应用，缺乏综合考虑历史韧性和未来抵抗力的整合框架。这一空白对于识别那些既在过去表现稳定、又预计在未来气候条件下保持适宜的区域至关重要，这类区域正是气候避难所（climate refugia）的候选地。
**二、研究内容与意义**
为填补这一空白，研究人员开发了一种新颖的整合方法，将基于多年代卫星时间序列的时间混合残差分析得出的植被韧性指标，与预测蓝橡（*Quercus douglasii*）栖息地适宜性的SDMs相结合。蓝橡是加州特有的冬季落叶橡树物种，是加州最广泛橡树林植物群落的代表。研究选取了加州大自然保护协会管理的两个保护区——位于蓝橡分布区北端的格雷·戴维斯·戴伊溪保护区（Dye Creek）和南端的弗兰克与琼·兰德尔保护区（Randall），这两个保护区代表了蓝橡栖息地的南北极端，增强了结果向其他橡树林景观推广的普适性。
研究人员应用此方法生成了空间显式地图，将历史植被格局与预测的气候适宜性相结合，识别出（1）历史上植被韧性高的区域，（2）预计未来气候下保持蓝橡适宜的（高抵抗力）区域，以及（3）兼具历史植被稳定性和预测栖息地稳定性的区域。研究还分析了这些区域的气候、地形和土壤条件变化。通过比较30米Landsat和10米Sentinel-2两种卫星平台的数据，检验了该方法在不同空间分辨率和辐射特性平台上的稳健性。野外实地验证利用兰德尔保护区18个蓝橡更新监测样地的数据，评估了景观分类与实地更新的对应关系。这项研究的意义在于提供了一个可扩展、可转移的方法，将历史地球观测与SDM气候预测相结合，为气候变化下的土地管理、恢复和土地保护规划提供科学支持，特别适用于像蓝橡这样繁殖周期缓慢、自然扩散能力有限的物种。
**三、关键技术与方法**
研究人员使用了以下关键方法：**（1）时间混合残差分析（Temporal Mixture Residual，TMR）**——首先利用线性光谱混合分析从Landsat（1982–2025）和Sentinel-2（2015–2025）影像中提取亚像元植被覆盖度时间序列，然后构建线性时间混合模型，用常绿多年生、落叶多年生和一年生草本三种时间端元表示稳定的植被物候轨迹，计算残差以量化偏离预期物候的植被不稳定性，最后对残差时间序列进行基于协方差的主成分分析（PCA），提取主导的植被变化模式（TMR 3）。**（2）物种分布模型**——利用已有的蓝橡SDMs（基于盆地特征模型v8的水文气候变量、土壤属性、地形因子构建），预测RCP 8.5排放情景下2040–2069年的栖息地适宜性，并平均两个通用环流模型（CNRM-CM5和HadGEM2-ES）的结果以降低不确定性。**（3）k-means聚类**——将TMR 3图层与未来栖息地适宜性图层叠加，对每个保护区和卫星数据源分别进行k-means无监督分类，生成五个可解释的景观类别。**（4）环境特征分析**——将水文气候、土壤和地形预测因子进行PCA，表征各类别的基线环境条件。**（5）野外数据验证**——利用兰德尔保护区2022–2023年间18个35米半径样地的蓝橡幼苗和幼树调查数据，通过泊松广义线性模型评估景观类别与实地更新的关系。
**四、主要研究结果**
**时间混合模型**：对两个保护区分别分析了Landsat和Sentinel-2时间序列。PCA显示前三个主成分解释了兰德尔>80%、戴伊溪>85%的方差。时间混合模型拟合良好，均方根误差值<0.05。Landsat衍生的TMR特征向量捕捉了过去四十年四次全州性干旱的不同程度响应，而Sentinel-2的TMR 3在2020–2022年干旱期间表现出高值。
**栖息地适宜性预测**：共识蓝橡SDM表现良好（AUC=0.95；TSS=0.81；Sorensen指数=0.84）。两个保护区未来蓝橡栖息地适宜性总体呈现高海拔地区高于低海拔山麓的空间格局。戴伊溪中部和东部的河岸区域（以柳树、西梧桐、谷橡等物种为主）当前和未来适宜性均较低。
**历史植被覆盖的时空格局**：TMR 3代表从落叶植被衰退或不稳定（低值）到植被恢复或覆盖度增加（高值）的梯度。在兰德尔，高TMR 3值区域自2010年起植被增加，与TNC在2008、2010和2011年保护首批土地地块的时间吻合，主要分布在高海拔地区；低值区域在2012–2016年干旱后急剧下降且仅部分恢复。在戴伊溪，高TMR 3值指示长期稳定覆盖加上近期增加，低值区域集中在河岸走廊，2020–2022年干旱和2024年公园大火（Park Fire）后出现植被损失。
**蓝橡林与稀树草原的景观分类**：通过k-means聚类将历史植被动态（TMR 3）与预测的21世纪中叶栖息地适宜性结合，在每个保护区识别出五个景观类别。兰德尔分类中，高适宜性区域分为植被覆盖增加（类别1）和受干旱影响但恢复中（类别2）两类；中等适宜性区域分为稳定（类别3）和干旱影响后恢复（类别4）；类别5为低适宜性区域，历史条件混合。戴伊溪的分类中，类别1为高适宜性、年际变异且近期增加，类别2为低-中等适宜性但同样变动和增加，类别3为中等适宜性且植被稳定，类别4为混合适宜性、受干旱影响部分恢复，类别5为低适宜性、历史条件混合。环境PCA显示，类别1位于高海拔、冬季降水较多的区域，类别5占据较温暖的低海拔区域。
**与野外蓝橡更新数据的比较**：兰德尔保护区18个监测样地的蓝橡更新数量（幼苗和幼树数）从0到10株不等，表现出显著的空间变异。更新数量在中等至低未来栖息地适宜性的景观类别（类别3–5）中显著下降。高更新样地更集中于Sentinel-2分类的类别1，该类别对应高未来栖息地适宜性和近期植被恢复或增长区域。100%拥有超过3株幼年蓝橡的样地位于预测21世纪中叶高栖息地适宜性的区域。
**五、讨论与结论**
**讨论部分总结**：该研究推进了先前在橡树林动力学监测和树冠覆盖度变化表征方面的工作，首次将卫星时间序列得出的长期植被轨迹与气候预测SDM输出相结合，作为蓝橡生态系统生态韧性的指标。Landsat和Sentinel-2分析产生了大致相似的空间格局，但两者提供互补信息：Landsat捕捉了四十年来的重大干旱影响，Sentinel-2突出了近期植被动态。环境梯度（海拔、温度和冬季降水）构建了景观分类，高海拔、凉爽、冬季降水丰富的区域被归类为未来气候避难所。火情历史和管理背景对观测到的植被覆盖格局有重要贡献。研究强调不应仅优先保护“安全”区域，而应保持景观镶嵌体的复杂性。局限性包括缺乏戴伊溪的野外更新数据，以及韧性/抵抗力指标可能仅部分捕捉这些复杂过程，但研究展示了该方法作为灵活、空间显式工具的价值，可与景观基因组学结合用于恢复优先级排序。
**结论翻译**：橡树林是加州保护区网络中保护最薄弱的生态系统之一，也是加州30×30保护区倡议（30×30 protected land initiative）应优先考虑的硬木景观之一。需要以稳健、一致的空间技术为指导的针对性保护工作，以确保橡树的持续存在，并支持加州各地的生态保护和管理。本文提出的生态分类框架支持一种灵活的、基于景观的方法，使土地管理者能够优先考虑一系列蓝橡保护和恢复地点，包括那些在变化条件下具有恢复或适应潜力的地点。最终，将历史生态动态与前瞻性环境变化预测相结合的工具，对于加州橡树林和其他脆弱生态系统的气候适应性保护战略至关重要。
