《Remote Sensing》：Long-Term Elevational Redistribution of Vegetation Greenness Along the Dushanzi–Kuqa Transect on the Southern Flank of the Tianshan Mountains
编辑推荐：
山区植被变化通常被概括为平均绿化或褐化趋势，或是离散边界的位移；这两种方式都可能遗漏植被绿度海拔剖面形状的变化。本研究利用沿独库公路（G217国道）布设的Landsat样带，测试了1989年至2024年间天山南坡绿度剖面是否发生了重组。研究人员在固定的1000
山区植被变化通常被概括为平均绿化或褐化趋势，或是离散边界的位移；这两种方式都可能遗漏植被绿度海拔剖面形状的变化。本研究利用沿独库公路（G217国道）布设的Landsat样带，测试了1989年至2024年间天山南坡绿度剖面是否发生了重组。研究人员在固定的1000–4100米海拔域内，仅使用在所有比较时段均有效的采样位置构建归一化差值植被指数（NDVI）和改进型双波段增强型植被指数（EVI2）剖面；对剖面进行均值中心化处理，并通过时空窗口、样带宽度、海拔分箱宽度和指数四种参数的36种组合下的空间块置换检验评估早期与近期剖面的差异。所有36种组合均显示显著的剖面形状变化（p ≤ 0.005），且具有跨情景一致性，该对比在校正跨传感器差异后依然成立。对比呈现一致的三带结构：1100–1950米处相对下降，2550–3300米处上升，3550–4100米处再次下降。绝对绿度几乎全域增加（25334个样本中超过97%），但中海拔带的增幅是低海拔带和高海拔带的四到五倍，沿样带方向变化幅度达八倍，且主要由海拔和沿样带位置决定（空间交叉验证随机森林R2 ≈ 0.66）。该植被对比与十年尺度温度、土壤水分或积雪对比的方向不一致。因此，沿该样带的绿度变化是整个剖面的海拔依赖性再分配，而非均匀绿化或单一边界位移，且在当前的十年尺度分析框架下未检测到其与十年水文气候对比之间存在稳健的方向对应关系。
一、研究背景与科学问题
山区海拔梯度在短距离内压缩了广泛的环境条件谱系，长期以来被视为生态学中的天然实验场。数十年的卫星观测记录了山区植被绿度的广泛变化。然而，现有研究通常将这种变化概括为两种模式之一：要么是向绿化或褐化的平均趋势，要么是林线或植被等值线等离散特征的位移。这两种方式都将景观的海拔维度简化为一个标量或单一的位移估计值，从而可能遗漏一种结构上截然不同的变化模式——即绿度在不同海拔带之间发生再分配（某些带相对于剖面均值增加而另一些带减少），而不表现为统一的方向性或连贯的边界位移。
已有证据表明，山地绿度变化受海拔组织这一现象在多个山系中均有积累。例如，喜马拉雅和兴都库什地区的Landsat记录显示了植被向亚冰雪带的扩张、高山绿度的海拔偏移及其强烈的经向变异，以及在干旱的西喜马拉雅地区无法用温度、降水和积雪指标解释的绿化趋势。欧洲阿尔卑斯山的积雪减少与高海拔生产力增加相关联；安第斯山脉的森林样地重访调查显示出广泛但异质性强的、由海拔组织的群落响应；而在干旱的中亚山地-盆地系统中，绿度-气候关系的符号随海拔发生逆转。区域汇总的绿化还可能掩盖随时间变化的同期褐化现象。在这些遥感研究中，变化主要被归纳为特定条带趋势、绿化幅度、植被扩张或边界位移，而整个绿度剖面的形状变化尚未得到系统检验。
离散边界作为长期变化检测的目标尤其具有挑战性。生态交错带很少是清晰的界线，它们通常占据混合组成和逐渐过渡的地带，因此其制图位置对分类不确定性高度敏感。在遥感分析中，估计出的边界不仅是生态特征，更是由空间尺度、测量类型、空间支撑和检测方法共同塑造的统计构造。当过渡渐变、局部对比弱或中等分辨率像元整合了异质性植被结构时，检测的不确定性最大。多指标的一致性可以缩小过渡范围，但每个指标特定的边界仍然依赖于估计方法。因此，林线研究已从单纯追踪边界位移转向对交错带宽度、结构和动态的多维描述。剖面水平的分析遵循同样的逻辑，它将海拔-绿度关系作为一个整体进行评估，无需预先假定某个离散边界是稳定的、唯一定义的或在生态上有明确标签。
中亚的天山山脉为此类分析提供了一个严苛的试验场，因为其海拔梯度上的植被变化在空间上异质性强，且受到干旱大陆性气候的强烈约束。基于Landsat的重建表明，1989–2015年间各垂直植被带和子区域并未均匀响应。区域NDVI分析显示，表观绿化在时间上非平稳，部分区域在1998年后因土壤水分限制加剧而转向褐化，这与更广泛的中亚地区植被生长日益受水分限制的证据一致，且绿化-褐化趋势方向和干旱敏感性随分析时段而变化。MODIS时代的记录还显示了海拔依赖但空间变异的返青日期变化。单一的位移度量难以捕捉此类变化。相比之下，固定支撑的剖面分析可以在不预设特定边界移动或机制的前提下，检验植被绿度的海拔组织是否发生了变化。
将该框架应用于Landsat存档需要超越历史上主导变化检测的双时相比较。表观剖面变化可能源于观测伪影，包括不同海拔数据可用性的变化、特定时期的云和云阴影掩膜以及不匹配的海拔域；特别是高山绿化趋势会因夏季观测密度随时间增加而被夸大。Landsat变化估计还取决于指数选择、预处理和算法决策。冠层结构的复杂性可能进一步使基于指数的绿度与植被功能脱钩，而不断扩大的指数家族使得依赖单一指数越来越难以辩护。因此，一项严格的检验必须：（i）保持各时期采样总体和海拔域恒定；（ii）将剖面形状变化与均匀的平均绿度变化分离；（iii）在多种合理的时相、空间、分箱和指数规格下保持稳定。这种结合空间匹配采样、形状聚焦检验和跨规格确认的方法此前很少应用于干旱山区的海拔植被变化。对于天山而言，已有的海拔分辨分析涉及垂直带位置、区域绿度趋势和返青时间，但据研究人员所知，在这些控制条件下整个海拔-绿度剖面形状的长期变化及其与十年水文气候对比的方向对应关系尚未得到检验。
然而，仅检测到变化本身并不足以确定其含义。一条走廊尺度的剖面对比可能掩盖沿样带的强异质性，因此在剖面模式可被解释之前，必须建立对比的空间分布及其由地形结构（如坡度和坡向）的组织方式。独立地，在水分受限的山地系统中，海拔依赖的绿度变化自然引发水文气候解释，但这种解释应被检验而非假设，例如通过评估植被对比是否与十年尺度的温度、土壤水分和积雪条件对比方向一致。将检测、空间组织和环境一致性作为三个独立的、顺序排列的问题来处理，可确保归因仅在检测和表征之后进行。
本研究沿穿越中国天山的独库公路（G217国道），使用1989至2024年的Landsat观测数据，聚焦于山脉南坡，依次处理上述三个问题。该序列反映了逻辑上的依赖关系：首先必须在控制观测伪影的条件下证明剖面形状变化是真实的；只有在此时，表征变化发生在何处以及如何受地形组织才具有意义；并且只有在现象及其空间范围被确定之后，才能对环境一致性进行检验，以避免对伪影或空间模糊的信号进行归因。该公路仅用于定义一条连续、可重复的地理采样线，不作为暴露因素处理，且分析中不引入任何道路距离变量。研究人员（1）测试了南坡植被绿度（NDVI和EVI2；以下简称绿度剖面）的海拔剖面在早期和近期Landsat时期之间是否发生了重组，采用了结果盲法空间匹配采样、均值中心化剖面、空间块置换检验以及跨越36种规格组合的交集规则；（2）通过结合样本水平变化图、沿样带汇总以及空间交叉验证随机森林分析，量化了变化发生的位置及其受地形结构组织的强度；（3）在一个单独冻结的设计下，测试了植被对比是否与十年水文气候对比方向一致。
二、关键技术方法
本研究采用的主要技术方法包括：基于Landsat Collection 2 Tier 1 Level-2地表反射率数据（TM、ETM+、OLI、OLI-2）的生长季（6月1日至9月30日）年度中位数合成与跨年medoid选取规则；在UTM Zone 45N投影下建立固定500米采样网格，并基于Copernicus DEM GLO-30提取高程、坡度等地形参数；实施三时期空间匹配采样策略以确保时间对比不受采样位置变化影响；采用均值中心化处理将剖面形状变化与均匀绿度变化分离，并使用5公里空间块级别的标签置换（4999次）生成零分布进行显著性检验；通过36种参数组合（两种时窗×三种走廊宽度×三种高程分箱宽度×两种植被指数）的交集规则和余弦相似度一致性检验确认结果的稳健性；运用跨传感器桥接检验评估传感器链对观测对比的可能贡献；采用分层对比分析和空间交叉验证（十折）随机森林回归（500棵树）量化地形关联性；最后，基于ERA5-Land再分析数据和CHIRPS降水数据进行十年水文气候一致性检验，使用符号内积和Spearman秩相关系数进行评估。
三、研究结果
3.1 南坡绿度剖面在早期与近期之间发生重组
在南坡走廊的主设计（10年时窗、20公里走廊、50米分箱）中，均值中心化后的NDVI和EVI2剖面均显示出清晰的早期至近期形状变化，近期剖面并非早期剖面的简单垂直平移。全剖面L2统计量NDVI为2.92、EVI2为1.82（均为p < 0.001）。尽管平均绿度在该时期对比中有所增加（NDVI +0.070，EVI2 +0.040），但由于每个剖面在差分前均已中心化，这一均匀增加本身不能产生形状对比。两种指数呈现出相同方向的对比，且全部36种规格组合均通过了置换检验和一致性标准，表明该重组在时窗、走廊宽度、分箱宽度和指数选择上是稳健的。中期剖面（2002–2011年，不参与主检验）落在早期与近期剖面之间（NDVI 63个分箱中的54个、EVI2 63个中的57个），且早期至中期与早期至近期的中心化对比相关性达r = 0.94–0.97，表明变化呈渐进式而非阶跃式。
3.2 一致的三带海拔结构
重组集中在特定海拔带而非均匀分布在1000–4100米范围内。NDVI方面，超出置换包络线的分箱集中在1100–1950米（负对比）和3550–4100米（负对比），正偏离则间断出现在2550–3250米；EVI2呈现相同结构，负偏离位于1150–1950米和3550–4100米，正偏离位于2800–3300米。由于剖面经过中心化处理，这些符号描述的是域内的相对再分配而非绝对增益或损失。中海拔带独自贡献了平方L2幅度的45.0%（NDVI）和46.8%（EVI2），但没有任何单个50米分箱贡献超过17.9%，两个域边缘分箱合计最多贡献20.1%。
3.3 走廊范围内普遍绿化伴随强烈中海拔放大效应
绝对样本水平变化将中心化对比置于物理情境中：25334个南坡样本中97.7%的NDVI增加（面积加权样本均值ΔNDVI = +0.065，ΔEVI2 = +0.040），295个5公里块中292个（99.0%）具有正块均值，前20%贡献块仅占总变化的46–48%。然而，变化幅度高度结构化：中海拔带（2550–3300米）平均NDVI增益为+0.127，是低海拔带（1100–1950米，+0.027）和上海拔带（3550–4100米，+0.030）的四到五倍。相对于早期基线，中海拔带增长约34–37%，而低海拔和上海拔带分别增长42–46%和56–70%。沿样带方向，10公里段ΔNDVI中位数变化幅度达八倍（0.017–0.139），20个段中有18个排除了走廊中位数。南坡走廊末端（520–530公里，约1100米）出现明显的次级增长，对应于库车绿洲农田，但通过ESA WorldCover 2021掩膜实验证实，低海拔带相对下降是自然和半自然植被的属性而非农田变化所致。
3.4 变化的地形关联性
三带结构在所有地形类别中保持稳定：在每个有数据的坡向和坡度层内，分层中心化对比的符号与全样本结构一致（与全样本剖面的Spearman相关系数0.912–0.988），地形类别仅调节对比的振幅和峰值位置而非其组织方式。空间交叉验证随机森林模型解释了约三分之二的出样空间方差（ΔNDVI出样R2 = 0.656，ΔEVI2 = 0.653；入样R2 = 0.87），表明大部分变化可由地形结构预测。共线性感知的排列重要性分析显示，海拔-路线里程联合排列导致R2下降1.09（NDVI）和1.13（EVI2），远大于其余地形属性联合排列的效果（0.24和0.20）。条件排列分析将联合重要性的主要部分归因于海拔：海拔在里程层内保留了0.57（NDVI）和0.63（EVI2）的重要性，而里程在海拨层内仅保留约0.10。偏依赖图再现了三带海拔形式和里程约430公里处的阶梯状转变。模型残差在整个走廊中较小（中位数≈0，5–95百分位?0.03至+0.03），但在库车绿洲农田内大且符号混杂，证实低海拔农业变化不以与山地信号相同的方式受地形组织。
3.5 跨规格稳健性与传感器链检验
重组并非主参数设置所特有：全部36个南坡L2置换检验均通过（其中30个p < 0.001，其余p ≤ 0.005），所有36个一致性检查也通过（中位数余弦相似度0.90–0.98，有效复制比例≥ 0.987）。四种额外光谱指数（kNDVI、NDMI、SAVI、MSAVI）均恢复了相同的三带结构（全部p < 0.001）。提高早期最小有效观测阈值至5、8或10同样保持了三带结构、上带相对下降和显著性。三项针对时间观测不平衡的直接检验——将近期时期重新合成为早期观测密度水平、用多年中位数合成替代跨年medoid、以及应用早期最小阈值重新合成——均复现了三带结构。交叉传感器检验表明，估计的传感器对比无法解释该重组：传感器对比投射到观测对比上的分数NDVI为0.249（95%区间0.168–0.325）、EVI2为0.058（?0.026至0.139），传感器校正后的对比保留了观测结构（余弦相似度分别为0.877和0.900）。
3.6 与十年水文气候对比无方向一致性
在单独冻结的十年设计下，三个ERA5-Land变量组均与南坡植被对比无方向一致性，均被归类为无十年关联，该分类在所有六种时窗×走廊宽度情景中保持一致。温度方面，符号内积与零无显著差异（NDVI: IP = 0.303, p = 0.466）。土壤水分方面，内积显著为负（NDVI: ?0.618, p < 0.001），即水分对比与植被对比方向相反。积雪方面，内积为负但与零兼容（NDVI: ?0.089, p = 0.072）。独立的水分代理CHIRPS降水对比则显著为正（5–8月NDVI IP = 538.7, p < 0.001），两个水分代理指向相反方向，在本设计中无法优选任一，因此无法做出稳健的水分方向推断。总之，在当前十年尺度、海拔分辨的设计下，未检测到已确认的南坡重组与十年水文气候对比之间存在稳健的方向对应关系。
四、讨论与结论
讨论部分总结指出，南坡走廊经历了植被绿度海拔剖面的稳健重组。绝对绿度几乎全域增加，但中海拔带增幅是低海拔和上海拔带的四到五倍，这种差异放大重塑了整个剖面的形状。低海拔和上海拔带的负中心化对比是相对于剖面均值的相对下降，而非绝对褐化。这种海拔放大的绿化在其他高山系统中也有报道。同时，已确认的对比与十年温度、土壤水分或积雪对比无方向一致性，且两个水分代理相互矛盾，因此该重组不能被简单解读为十年水文气候均值的响应。几种非互斥的解释仍然开放，包括相对于十年平均状态的滞后或非线性响应、未测试变量的控制作用（如生长季积雪动态、放牧或土地利用变化），以及海拔分箱气候聚合值与结构化山地植被过程之间的尺度不匹配。研究人员的剖面水平主张与边界风格分析形成鲜明对照：后者在探索阶段未能通过预定的确认链，说明边界位移估计可能依赖于估算器选择和观测支持。空间分析以三种方式改变了剖面结果的解读方式：普遍性排除了局部伪影的可能性；地形组织的强度界定了需寻求解释的范围；位置-海拔梯度对坡度和坡向的主导地位支持将该模式解读为海拔现象而非纯粹的地形现象。
研究结论部分原文翻译如下：本研究将植被绿度的海拔剖面本身作为分析对象，探究了1989至2024年间中国天山南坡独库公路沿线其形状是否发生了重组。结果表明确实发生了重组：相对于剖面均值，绿度在1100–1950米和3550–4100米处下降，同时在2550–3300米处上升——这一三带结构对NDVI和EVI2均成立，通过了36种规格组合的交集规则，并在跨传感器证伪检验中幸存。由于绝对绿度在几乎整个走廊范围内增加（25334个样本中97.7%），这些下降是相对于剖面均值的相对再分配，而非植被损失：该剖面是由一个近乎普遍的增幅被大约四到五倍的差异放大所重组的。由此得出两个对山地植被监测的意义。第一，剖面形状捕捉了一种平均趋势和边界位移度量都可能遗漏的变化模式，两者不应混淆：同一影像的边界风格分析未能通过预定的确认链，说明了位移声明的估算器依赖性及其与数据实际支持的偏差。第二，空间分析界定了需寻求解释的范围。由于约三分之二的变化可在出样评估下由地形结构预测，且符号混杂的残差局限于库车绿洲，走廊研究中跨越山地-绿洲梯度的低海拔农业段必须与山地信号分开处理。研究结果的范围和解读应保持在设计之内。所有推论均来自采样的独库公路南坡走廊，本身不能确立全范围模式；负中心化对比描述的是相对再分配而非褐化；在当前十年尺度、海拔分辨的设计下，未检测到与十年温度、土壤水分或积雪对比之间的稳健方向对应关系——这是未检测到的对应关系，而非水文气候不起作用的证据。解释该重组将需要土地利用、积雪动态和植被组成的独立约束。将海拔绿度剖面作为分析对象——配合空间匹配采样、均值中心化和跨规格确认——提供了一种可转移的方法来检测山地植被再分配，而不会将其误读为边界迁移或简单的气候响应。
