哈萨克斯坦南部气候变化生态敏感性评估：基于NDVI–气候复合指数（Composite NDVI–Climate Index）方法（2010–2025）

《Diversity》：Assessment of Ecological Sensitivity to Climate Change in Southern Kazakhstan: A Composite NDVI–Climate Index Approach (2010–2025)

　　

摘要：气候变化威胁中亚干旱区的生态系统稳定性，但在土地利用政策可操作尺度上，区域植被对气候变化的响应仍缺乏明确解析。研究人员整合了哈萨克斯坦水文气象局（Kazhydromet）2000–2024年长期气象记录，以及哈萨克斯坦南部五个行政区（约710,000 km2）四个时期（2010、2015、2020、2025）的Landsat地表反射率影像。经Landsat 5 TM与Landsat 8 OLI跨传感器一致性校正后，发现高密度植被覆盖（归一化植被指数NDVI > 0.4）在各区域均小幅增加，累计面积从909万公顷增至960万公顷（+5.6%），其中2020年的暂时性最低值与2018–2020年干旱有关。由于四期采样（n = 4）限制，各区普通最小二乘法（OLS）趋势斜率在 p < 0.05 水平上无统计显著性。构建的复合生物多样性–气候敏感性指数（Biodiversity–Climate Sensitivity Index, BCSI）由四个归一化分量组成——温度趋势、降水亏缺、NDVI趋势，以及作为生物多样性–脆弱性替代指标的密植被覆盖变异系数（Coefficient of Variation, CV），识别出锡尔河下游洪泛平原及前咸海沿岸为最敏感区域，北天山为最具韧性区域。该框架为契合可持续发展目标（SDG 13与SDG 15）的气候适应性保护提供了可操作循证基础。

论文解读：哈萨克斯坦南部基于NDVI–气候复合指数的生态系统气候变化敏感性评估

一、研究背景与立项依据

气候变化被广泛认为是全球生物多样性丧失与生态系统退化的主要驱动力之一，气温升高、降水格局改变及极端气候事件频发直接影响物种分布、生态系统生产力与生态稳定性。上述影响在干旱与半干旱区尤为显著，因其生态系统运行已接近生理与水文极限。哈萨克斯坦南部涵盖沙漠、半沙漠、草原、山麓及天山西脉与准噶尔阿拉套山脉生态系统，拥有众多特有种与孑遗种，近年观测到明显增温与干旱化趋势，对生物多样性保护与可持续土地利用构成挑战。尽管中亚干旱区已有变暖记录，但植被对气候变化的响应并不单向——部分地区出现褐变（browning），部分地区因CO₂施肥效应、灌溉扩张及放牧减少出现绿化（greening）。加之哈萨克斯坦南部特有的人类活动压力（如什姆肯特等次级城市固体废物堆积）叠加，使植被变化归因复杂化。目前在土地利用政策可操作尺度（即行政区划）上，尚缺乏分区、多时相、多指标的解耦分析。为响应联合国可持续发展目标SDG 13（气候行动）与SDG 15（陆地生物），有必要在区域尺度量化气候–植被交互作用，识别气候适应性保护的优先区，因此研究人员开展了此项研究。该论文发表于《Diversity》期刊。

二、主要关键技术方法概述

研究人员选取哈萨克斯坦南部五行政区（克孜勒奥尔达Kyzylorda、突厥斯坦Turkistan、江布尔Zhambyl、阿拉木图Almaty、杰特苏州Zhetysu，总面积约710,000 km^{2）。气候数据源自哈萨克斯坦国家水文气象局（Kazhydromet）2000–2024年气象站年均温与年降水量记录，采用普通最小二乘法（Ordinary Least Squares, OLS）回归及Mann–Kendall检验与Theil–Sen斜率估算分析温度线性趋势，降水变率用变异系数（Coefficient of Variation, CV_P）表征。植被数据采用USGS EarthExplorer获取的Landsat Collection 2地表反射率产品（Landsat 5 TM对应2010年，Landsat 8 OLI对应2015、2020、2025年），生成7–8月生长季峰值合成影像，利用CFMask进行云掩膜，按标准公式择波段计算NDVI（Landsat 5 TM用Band 3红、Band 4近红外；Landsat 8 OLI用Band 4红、Band 5近红外），并按Roy等(2016)报道的OLS系数对OLI反射率做TM等效转换以实现跨传感器一致性校正。以NDVI > 0.4划定密植被覆盖，在ArcGIS Pro 3.2中以WGS 1984 UTM Zone 42N投影进行分区统计。构建复合生物多样性–气候敏感性指数（Biodiversity–Climate Sensitivity Index, BCSI），含四归一化分量——2000–2024年均温趋势（热胁迫代理X₁）、年降水量趋势取反（水文胁迫代理X₂）、密植被覆盖趋势取反（植被响应代理X₃）、2010–2025密植被覆盖CV（生态系统脆弱性代理X₄），经min–max标准化后等权（各0.25）加和，并做三种替代权重方案敏感性测试，BCSI连续表面按自然断点法划为五等级敏感性类别。整体工作流包括数据预处理与跨传感器校准、气候与NDVI时序趋势分析（OLS并行Mann–Kendall）、分区密植被统计量计算（均值、标准差、CV、OLS斜率、R2、干旱响应与恢复）、BCSI计算与权重敏感性测试、五区及三大生态系统类型（沙漠、草原、山地）比较评估。}

三、研究结果

3.1. Climatic Context, 2000–2024（2000–2024年气候背景）

通过对Kazhydromet长期记录分析，研究人员发现哈萨克斯坦南部呈全区一致变暖趋势，每十年升温速率约+0.28至+0.42 ℃，增温最强出现在克孜勒奥尔达与突厥斯坦南部低地干旱区，阿拉木图与杰特苏州山地生态系统升温较温和但仍显著。降水变率增大：干旱与草原区降水稳定或略降，部分山区略增。2012年及2018–2020年出现明显干旱脉冲，与区域报道一致。

3.2. Spatial NDVI Patterns and Dense-Vegetation Dynamics（NDVI空间格局与密植被动态）

四期 harmon化Landsat NDVI空间分布显示，高NDVI值集中于阿拉木图与杰特苏州山麓山地及突厥斯坦中部灌区，克孜勒奥尔达沙漠低地与部分江布尔区NDVI低（<0.2）；2020年各区间中等NDVI（0.2–0.4）缩减反映2018–2020干旱，2025年有所恢复。密植被覆盖（NDVI > 0.4）2010–2025年间各区净变化均为正：阿拉木图由29.50%升至31.13%（+1.63个百分点），杰特苏州由25.40%升至26.74%（+1.34个百分点），突厥斯坦由8.57%升至8.96%（+0.39个百分点），江布尔由6.75%升至7.13%（+0.38个百分点），克孜勒奥尔达由4.51%升至4.80%（+0.29个百分点）。哈萨克斯坦南部密植被总面积由2010年909.2万公顷增至2025年960.2万公顷（+5.6%），2020年出现暂时性低谷（938.5万公顷）对应干旱影响。四期OLS趋势斜率各区在 p < 0.05 水平无统计显著性（n = 4），OLS拟合R²为0.508（突厥斯坦）至0.691（阿拉木图），表明线性时间分量可解释相当比例时期间方差但样本量限制无法确证持续绿化轨迹。各区平均覆盖度与CV呈反比——阿拉木图与杰特苏州山地（均值≈26–31%）CV最低（1.97–2.05%），克孜勒奥尔达沙漠（均值≈4.7%）CV最高（2.38%），支持将CV纳入BCSI生态敏感性组分。2020年干旱响应以干旱低地为最强（?2.0至?2.5%），山地最弱（?0.8%）；2020→2025恢复幅度亦呈相同梯度（克孜勒奥尔达+4.4%，阿拉木图+1.6%）。

3.3. Composite Biodiversity–Climate Sensitivity Index (BCSI) Results（复合生物多样性–气候敏感性指数BCSI结果）

BCSI空间分布显示明显异质性。克孜勒奥尔达大部为中–高BCSI，锡尔河下游洪泛平原及前咸海沿岸广泛分布极高敏感区；突厥斯坦呈梯度——中部灌区为低–中BCSI，克孜勒库姆与莫因库姆沙漠边缘敏感性升高；江布尔以中等BCSI为主，吉尔吉斯山脉山麓及塔拉斯河与楚河走廊见局部低敏感区；阿拉木图区内北天山山地东南 sector 为极低–低BCSI（最具韧性），低地草原及卡普恰盖水库周边为中–高BCSI；杰特苏州呈准噶尔阿拉套山麓低敏感向巴尔喀什流域中BCSI过渡。山地湿润山麓生态系统BCSI持续偏低，干旱低地、水文改造区及沙漠–草原过渡带为高敏感优先保护区。BCSI空间格局在三种替代权重方案下极高与极低敏感类空间位置保持不变。

四、讨论与结论翻译

讨论部分指出：观测到的每十年0.28–0.42 ℃增温符合中亚更广气候型态，低地蒸散增强与水分胁迫加剧，2018–2020干旱脉冲致2020年植被暂时缩减而后全面恢复，干旱响应与恢复幅度自干旱低地向山地递减。密植被面积净增5.6%与部分中亚干旱区退化叙事相悖，但四期采样无力确证趋势，方向一致支持区域韧性而非确立持续绿化；观测到的绿化信号可能叠加入渗CO₂施肥、后苏联时代放牧压力未恢复至1990年代峰值致草原与半沙漠区部分再生、以及灌区植被人为绿化贡献，未来可用全球灌溉农田掩膜（如GFSAD30）分解。BCSI揭示韧性分布不均——锡尔河下游洪泛平原与前咸海沿岸最敏感（高/极高BCSI），北天山山地最韧性（极低BCSI），区内山地–低地双峰反差强调亚区分辨率在保护规划中的重要性；表2中高CV、低均值覆盖区对应高BCSI，支持以密植被覆盖CV作生态系统脆弱性替代指标。结果与近期中亚区域评估相符，BCSI在行政可操作尺度集成多指标是本研究独特贡献。作者也阐明局限：五年度四期NDVI采样统计效力不足，建议未来用Google Earth Engine构建Landsat全时序年度序列并做Mann–Kendall与Sen斜率分析；NDVI > 0.4阈值可能低估稀疏草原生产力，宜改用生长季积分NDVI；BCSI生物多样性组分（X₄）现为密植被覆盖CV遥感替代而非物种丰富度或Shannon指数，未来应接入哈萨克斯坦国家生物多样性信息系统或GBIF数据；未显式分离气候驱动与人类活动（放牧强度、灌溉取水、土地利用变化、固体废物污染），建议叠加人为压力层及灌溉农田掩膜；BCSI等权虽经多方案敏感性验证空间高/低敏感类稳定，未来可用层次分析法（Analytic Hierarchy Process, AHP）由专家赋权并结合物种级生物多样性数据优化。

结论翻译如下：

本研究对2010–2025年哈萨克斯坦南部生态系统气候–植被交互作用进行了分区综合评估，融合Kazhydromet公开气候数据、跨传感器一致性校正多时相Landsat NDVI影像及复合生物多样性–气候敏感性指数（BCSI），得出三项核心发现。第一，多时相NDVI分析揭示生态系统特异性格局：阿拉木图与杰特苏州山地生态系统密植被覆盖最高（均值分别为30.54%与26.27%）且时期间波动最低（CV ≈ 2.0%）；江布尔与突厥斯坦草原及山麓景观覆盖中等且较稳定；克孜勒奥尔达干旱生态系统虽绝对植被覆盖最低（4.66%），但未出现统计显著下降，表现出意外稳定性。2020年干旱响应干旱低地最强（?2.0至?2.5%）、山地最弱（?0.8%），2020–2025年恢复幅度遵循相同梯度。第二，哈萨克斯坦南部密植被覆盖总面积从2010年909.2万公顷增至2025年960.2万公顷——15年间净增+5.6%，2020年出现暂时性最低值938.5万公顷；每一行政区均录得正净变化（克孜勒奥尔达+0.29至阿拉木图+1.63个百分点）。重要的是，四期采样设计下各区OLS趋势在 p < 0.05 均无统计显著性；观测变化方向与区域韧性一致，但尚不能解释为持续绿化轨迹的确凿证据，需用连续年度时序确认，此为本工作最高优先级延展方向。第三，BCSI将气候与植被指标综合为决策支持工具，有效刻画行政区操作尺度内的空间敏感性——克孜勒奥尔达锡尔河下游洪泛平原及前咸海沿岸为哈萨克斯坦南部最脆弱区，凸显针对性适应管理必要性；反之阿拉木图周边北天山山地具韧性，提示保护现有生态条件的重要性。BCSI优先级分类在不同权重方案下保持一致，强化其保护决策实用价值。框架依赖公开数据源确保可复现性与可移植至面临相似气候挑战的其他干旱半干旱区，支撑契合SDG 13与SDG 15及哈萨克斯坦国家气候战略的循证环境管理。与全球其他干旱半干旱区比较进一步表明，哈萨克斯坦南部遵循气候驱动植被转型的更广泛全球趋势，包括局部绿化、干旱诱发不稳定性及变暖下生态系统异质性增强并存过程。