长期气候变暖背景下中亚石竹科稀有植物Acanthophyllum gypsophiloides不同生境的土壤条件对比与NDVI动态分化

《Soil Systems》：Contrasting Soil Conditions and Divergent NDVI Dynamics Under Long-Term Climate Warming Across Acanthophyllum gypsophiloides (Caryophyllaceae) Habitats in Central Asia

　　

中亚稀有植物Acanthophyllum gypsophiloides（石竹科）分布于锡尔河卡拉套山脉环境差异显著的山地与山麓地带。研究人员评估了三个种群——Sayasu、Kuyik和Birlik——在共同区域气候背景下是否伴随剖面尺度的土壤条件对比，以及这些局部对比是否伴随着种群结构配置与当地植被绿度（NDVI）的差异。三个采样土壤剖面的物理化学特征差异显著。Sayasu和Kuyik以粗质地物质为主，细颗粒含量低，而Birlik剖面含有显著更高的<0.01 mm颗粒组分，但其腐殖质和大量营养元素总浓度最低，pH值最高（8.55–8.69），碳酸盐含量明显更高（CO2，4.91–7.04%）。WorldClim和CHELSA均显示1980–2019年期间存在长期变暖趋势，尽管两个产品估算的速率不同，而降水格局则更依赖于数据集选择。三个地点的种群密度和个体发育结构也存在差异。在种群位置提取的Sentinel-2 NDVI值显示，2025年Sayasu、Kuyik和Birlik分别比相应的2016年值下降了46.3%、59.5%和61.0%；这些观测代表的是像元尺度上的植被绿度，而非物种特异性的种群表现。总体而言，三个地点主要在共同区域变暖背景下因采样土壤剖面条件的差异而分化，而个体发育结构与当地植被绿度提供了补充性的生态证据。这些共现格局识别出了与保护监测相关的地点特异性环境对比，但并未建立土壤条件、气候、植被绿度与种群动态之间的因果关系。

研究背景、问题与研究意义

中亚地区属于强烈大陆性气候和水分受限区域，经历了显著的升温、水资源变化以及降水的空间异质性。山地生态系统中的植被响应呈现异质性，卫星与地面研究同时记录了绿化与褐化现象，其驱动因素包括气候、干旱、放牧、积雪状况及人类活动等多种因子。然而，关于物种水平的生态响应与适应过程的证据仍然分布不均。对于地理分布不连续的稀有植物而言，单一区域的普遍气候背景无法完全捕捉各个种群所经历的环境约束。因此，亟需在区域气候框架下，结合局地土壤条件、种群结构与遥感植被指标，对稀有植物生境进行多维度对比分析。

针对上述问题，研究人员选取了分布于南哈萨克斯坦锡尔河卡拉套地区的稀有植物Acanthophyllum gypsophiloides的三个自然种群——Sayasu（山地坡面）、Kuyik（山口高海拔）和Birlik（山麓平原）——开展研究。核心目标是评估这三个地点在共享区域气候背景下是否存在剖面尺度的土壤条件分化，以及这种分化是否伴随种群个体发育结构和当地植被绿度的差异。研究的意义在于揭示局地环境配置的差异性，为保护监测提供地点特异性的生态依据，而非建立因果机制。

主要技术方法

研究人员采用了以下几项关键技术方法：1）土壤调查与理化分析：于2022年春夏季在三个种群地点各挖掘一个土壤剖面，按遗传层次分层采样，采用卡钦斯基沉降法测定颗粒组成，丘林氧化法测定腐殖质含量，凯氏定氮法测定全氮，分光光度法测定全磷和全钾，电位法测定pH，气量法测定碳酸盐含量，水浸提法分析可溶性盐及离子组成。2）气候数据处理：使用WorldClim历史月气象数据和CHELSA两套网格化气候产品，提取1980–2019年月平均气温和降水量数据，经ArcGIS 10.4进行空间处理和时间聚合，计算年度、季节及年代际指标。3）NDVI遥感分析：基于Google Earth Engine平台处理Sentinel-2 MSI影像，选取2016、2020和2025年6月的无云影像合成月均值，计算归一化差值植被指数（NDVI）。4）统计建模：采用Spearman秩相关分析土壤变量间关系，主成分分析（PCA）综合气候变量，季节性自回归积分滑动平均（SARIMA）模型对气候时间序列进行外推预测。

研究结果

3.1 土壤理化性质与土壤分化

3.1.1 土壤剖面形态与质地

Sayasu和Kuyik剖面以粗质地、多碎石为特征，其中0.25–1.0 mm砂粒组分分别占66.58%和56.95%；Birlik剖面则以较细质地为主，0.01–0.05 mm粉粒占33.05%，<0.01 mm细颗粒组分高达39.14%。USDA等效质地分类进一步确认了这一差异。

3.1.2 碳酸盐-碱性条件、腐殖质含量与大量营养元素

Birlik剖面表现出显著不同的化学特征：腐殖质含量最低（0.83–0.96%），全氮、全磷、全钾浓度均低于其他两个剖面，pH值最高（8.55–8.69），碳酸盐含量（CO₂）达4.91–7.04%，远高于Sayasu和Kuyik。Sayasu和Kuyik化学性质较为接近，但Kuyik腐殖质略高，Sayasu全氮略高。

3.1.3 综合剖面尺度土壤对比

Spearman相关分析揭示了颗粒组成、碳酸盐-碱性特征、可溶性离子与大量营养元素之间的共变关系。砂粒与粉粒呈强负相关（r_s = -0.98），总可溶性盐与Ca²⁺呈正相关（r_s = 0.94），与全氮呈负相关（r_s = -0.93）。整体上区分出粗质地-有机养分组分与细质地-碳酸盐-碱性-盐分组分两大关联群。

3.2 区域气候背景：共同变暖与数据集依赖的降水格局

3.2.1 共同的长期变暖信号

WorldClim和CHELSA均一致显示1980–2019年三个地点存在长期升温趋势，但变暖幅度存在差异：WorldClim约为0.38–0.40°C/十年，CHELSA约为0.68–0.73°C/十年。Birlik为最温暖地点，Kuyik最凉爽，Sayasu居中。春季升温最为显著。

3.2.2 数据集依赖的降水格局

降水估计表现出更强的数据集依赖性。WorldClim显示Birlik降水最少，而CHELSA则给出Birlik较高的降水估计，尤其在后几十年。两个产品在Birlik处的降水差异最大。

3.2.3 跨数据集时间对应关系

温度指标的跨数据集Pearson相关性较高，而降水相关性较弱且变异较大，尤其在Birlik处最为明显。PCA前两个主成分解释了71.1%的气候变异，Birlik位于PC1负侧（最温暖），Kuyik位于正侧且分散更广。

3.3 探索性SARIMA统计外推

SARIMA模型将历史月气候序列外推至2020–2040年，结果显示所有地点-数据集组合的外推平均温度均高于2010–2019年历史均值。降水外推结果一致性较差，尤其在Birlik处数据集间分歧更大。

3.4 三个观测年份的生境水平NDVI变化

Sentinel-2 NDVI数据显示，所有地点2020年NDVI均低于2016年。Sayasu在2020–2025年间变化不大，Kuyik在2020年极低值后有所回升，Birlik则逐年下降。与2016年相比，2025年Sayasu、Kuyik和Birlik的NDVI分别下降了46.3%、59.5%和61.0%。

3.5 种群个体发育背景

已有调查数据显示，Sayasu种群密度最低（1.5±0.5株/m2），个体发育结构偏成熟；Kuyik密度较高（2.7±1.07株/m2），前期阶段个体占比最大；Birlik密度最高（2.9±1.09株/m2），以成熟生殖个体为主。

讨论与结论

讨论部分指出，三个Acanthophyllum gypsophiloides生境的核心差异在于局地上壤特征，而非区域气候本身。粗质地土壤（Sayasu和Kuyik）可能影响水分入渗与持留，而细质地且富含碳酸盐的Birlik土壤则面临较强的碱性与营养限制。但这些潜在约束并未被直接验证为实际的水分胁迫或营养缺乏。气候分析确认了区域变暖背景，但降水估计的高度数据集依赖性限制了对其生态意义的推断。NDVI变化反映了像元尺度植被绿度，不能等同于物种特异性种群表现。种群个体发育结构虽提供了补充信息，但不能单独用于判断种群稳定性或生存力。

研究结论翻译如下：本研究在长期气候变暖的区域背景下，识别出三个Acanthophyllum gypsophiloides地点在剖面尺度上的土壤条件对比。采样土壤在粗质地物理特性与碳酸盐-碱性化学特性方面表现出不同权重，指示了潜在环境约束的不同组合，而非单一的生境适宜性梯度。气候、遥感和个体发育分析提供了互补视角：WorldClim和CHELSA一致显示区域变暖，但降水估计更依赖于数据集选择；Sentinel-2观测记录了选定年份间生境水平植被绿度的差异；已有种群数据描述了不同的个体发育结构。这些观察无法建立土壤条件、气候变率、植被绿度与种群表现之间的直接因果关系。研究结果为地点特异性的保护监测提供了支持，建议整合土壤水力与化学测量与纵向的种群更新、存活和个体发育转变评估。未来需要通过空间重复的土壤采样、局地气象观测和时间分辨的植被监测，来确定已识别的环境对比是否影响A. gypsophiloides种群的长期存续。