《Sustainability》：Widespread Greening but Uneven Soil Carbon Change: Twenty-Year Evidence from a Semi-Arid Eurasian Steppe
编辑推荐：
准确评估长期土壤有机碳（SOC）变化需要重复的实地测量和对时间差异下土壤质量的明确处理。研究人员量化了位于哈萨克斯坦中部半干旱栗钙土带的20个永久监测站点在2004年至2024年期间的净SOC储量差值，使用了固定深度（FD）和等效土质量（ESM）计算方法，并结
准确评估长期土壤有机碳（SOC）变化需要重复的实地测量和对时间差异下土壤质量的明确处理。研究人员量化了位于哈萨克斯坦中部半干旱栗钙土带的20个永久监测站点在2004年至2024年期间的净SOC储量差值，使用了固定深度（FD）和等效土质量（ESM）计算方法，并结合Landsat归一化植被指数（NDVI）和水文气候摘要数据。在0–50 cm土层，基于FD计算的SOC储量从69.70 ± 17.43 Mg C ha?1增加到79.39 ± 25.00 Mg C ha?1，平均配对差值为+9.69 Mg C ha?1（95%置信区间: 2.78至16.60; p = 0.008）。在根据特定站点的2004年参考土质量进行ESM调整后，2024年的平均值为75.43 ± 25.24 Mg C ha?1，配对差值为+5.73 Mg C ha?1（95%置信区间: ?0.78至12.25; p = 0.081）。五个特定深度的FD比较在经过Holm校正后均不再显著。正的年复合NDVI斜率约占制图区域的94.7%，但站点层面的ΔNDVI与ESM校正后的ΔSOC之间的关联性较弱且不显著（皮尔逊相关系数 r = 0.259, R2 = 0.067, p = 0.270; 斯皮尔曼等级相关系数 ρ = 0.241, p = 0.307）。平均干旱指数从2004–2008年的0.345下降到2020–2024年的0.308，这被解释为水文气候背景，而非SOC变化的已证明驱动因素。由于大多数管理-土壤组合缺乏重复，因此仅将其作为探索性的站点层面比较进行了检验。这些结果表明，ESM校正会实质性地影响推断，而两期重采样设计支持对净终点差值的解释，而非连续的SOC变化。
论文解读文章
研究背景与科学问题
土壤有机碳（SOC）是陆地碳循环的关键组成部分，对维持土壤功能和缓解气候变化至关重要。在半干旱地区，土壤因生物量输入有限和高蒸发而脆弱，但其巨大的全球分布面积也意味着潜在的碳汇能力。然而，长期、基于实地观测的SOC动态研究在这些区域十分稀缺，尤其是在中亚哈萨克斯坦的广阔半干旱景观中。现有研究多基于短期测量或模型模拟，难以真实反映田间条件下的SOC储量变化。此外，多数研究集中于表层土壤（0–30 cm），而对更深层次SOC储量的关注不足。一个关键的方法论问题是：当两次采样间的土壤容重发生变化时，传统的固定深度（FD）估算可能会产生偏差，因为固定几何深度内所含的土壤质量并不恒定。因此，采用等效土质量（ESM）方法进行校正对于获得可靠的储量变化比较至关重要。同时，遥感指标如归一化植被指数（NDVI）虽能指示地上绿色度，但与地下碳输入和稳定过程的关联是间接的。为填补这一知识空白，本研究旨在利用重复采样的长期监测网络，量化半干旱欧亚草原区SOC的净变化，并评估不同计算方法及环境协变量对推断结果的影响。
研究方法概述
研究人员采用了一项为期20年（2004–2024）的重访式定点采样设计，在哈萨克斯坦中部的20个永久监测点进行土壤采样。关键技术方法包括：
- 1.
长期定位监测与土壤采样：于2004年建立20个代表不同土地利用历史和土壤条件的永久监测样地，并在2024年使用相同方案进行重访采样，获取0–50 cm剖面每10 cm一层的土壤样品。
- 2.
实验室分析与SOC转换：采用Tyurin湿氧化法测定土壤有机质（SOM），并通过两步换算（Van Bemmelen因子和湿燃烧校正系数）得到SOC浓度。使用环刀法测定土壤容重。
- 3.
SOC储量估算与ESM校正：分别使用传统固定深度（FD）方法和等效土质量（ESM）方法计算SOC储量。ESM方法以2004年各站点的0–50 cm累积土质量为参考，对2024年的SOC储量进行校正，以消除因容重变化导致的土质量差异带来的偏差。
- 4.
遥感和气候数据分析：利用Google Earth Engine（GEE）平台，基于Landsat影像提取2004–2024年生长季（5–9月）的年际NDVI数据，并进行像元尺度和站点尺度的趋势分析。同时，从CHIRPS和ERA5-Land数据集获取降水、温度等气候数据，并计算干旱指数（AI）。
研究结果
3.1. SOC储量与土壤物理性质的时间变化
- •
通过配对t检验发现，在0–50 cm剖面上，FD法估算的SOC储量在20年间平均增加了+9.69 Mg C ha?1（p = 0.008），统计上显著。
- •
然而，在进行ESM校正后，平均增加量降至+5.73 Mg C ha?1，且不再具有统计学显著性（p = 0.081）。这表明FD法的显著结果部分源于土壤质量的变化，而非纯粹的碳积累。
- •
对各10 cm分层进行的FD比较显示，所有层均为正值，但在经过Holm多重比较校正后，没有一个分层的差异达到显著水平，表明SOC变化在垂直方向上存在高度异质性，而非集中在某一特定深度。
- •
站点水平的ESM校正ΔSOC显示出巨大差异，13个站点为正，7个站点为负，最大值（+38.40 Mg C ha?1）和最小值（?19.33 Mg C ha?1）并存，强调了空间异质性。
3.2. 2004–2024年生长季NDVI模式
- •
基于年度生长季中值NDVI合成的空间趋势分析显示，约94.7%的研究区域呈现正向NDVI斜率，表明广泛的区域性绿化趋势。
- •
交叉传感器敏感性分析表明，尽管Landsat 8的引入可能高估了斜率的绝对值，但正向趋势的方向在站点层面是稳健的。
3.3. 植被动态与SOC终点差值的关系
- •
将站点按ESM校正ΔSOC的正负分组，两组站点的NDVI均呈上升趋势，表明绿化信号并非SOC增加所独有。
- •
站点层面的相关性分析显示，ΔNDVI与ESM校正的ΔSOC之间的关联非常弱且不显著（R2 = 0.067, p = 0.270），说明地上植被绿度的变化不能很好地预测地下SOC储量的变化。
3.4. 干旱指数变异性与水文气候背景
- •
研究期间干旱指数（AI）波动剧烈，但近期（2020–2024年）的平均AI（0.308）低于基线期（2004–2008年）（0.345），表明研究后期气候趋于干燥。此结果被视为解释SOC变化的背景信息。
3.5. 土地利用历史与土壤条件的探索性比较
- •
对站点按土地利用类型进行分组描述发现，活跃农田的两个站点表现出最大的平均SOC增加，但由于样本量极小，无法推广。
- •
管理-土壤条件矩阵图展示了巨大的站点间变异，即使在相同的管理类别内，SOC变化方向也可能相反，例如M2类别下的ΔSOC范围从?7.17到+38.40 Mg C ha?1。这表明管理、土壤和站点特异性效应相互混杂，无法分离。
讨论与结论
讨论部分总结：
讨论指出，本研究的核心发现是SOC变化的推断强烈依赖于计算方法。ESM校正揭示了FD法可能因土壤质量变化而产生误导性结论，强调了在长期重访研究中报告质量校正敏感性的重要性。NDVI与SOC的弱关联性表明，仅凭遥感绿化信号不足以推断土壤碳库的变化，因为SOC的储存受控于更复杂的根系输入、分解和矿物稳定过程。尽管观察到近期气候趋于干燥，但由于只有两个时间点的SOC数据，无法建立因果关系。土地利用和土壤类型的探索性分析虽然提供了有价值的背景，但因其非平衡的实验设计，无法区分管理和土壤的独立效应。
研究结论翻译：
本研究量化了哈萨克斯坦中部半干旱栗钙土带20个重访监测站点2004–2024年的净SOC差值。在0–50 cm土层，固定深度（FD）法显示平均配对增加量为+9.69 Mg C ha?1（p = 0.008），而等效土质量（ESM）校正后的平均差值为+5.73 Mg C ha?1（95% CI: ?0.78至12.25; p = 0.081）。因此，剖面尺度的推断对土壤质量的时间差异很敏感，且ESM校正后的均值在α = 0.05水平上与零无统计学差异。特定深度的FD差异在各站点间具有异质性。尽管原始检验提示0–10 cm和40–50 cm存在差异，但经Holm校正后，五个深度区间均不显著。因此，数据并未提供统计学证据表明正的SOC差异优先集中在20 cm以上或其他特定深度区间。年复合NDVI斜率在研究区域内占主导地位，但站点层面的ΔNDVI与ESM校正的ΔSOC之间的关联性较弱且不显著。近期较低的干旱指数提供了水文气候背景，但不能被解释为SOC变化的已证明时间驱动因素。同样，管理-土壤矩阵是探索性的，无法分离管理、土壤和站点效应。因此，本研究的主要贡献在于方法论和实证方面：重复站点的实地测量结合明确的ESM校正揭示了在一个异质性半干旱景观中，储量变化推断在很大程度上取决于土壤质量的处理方式。未来的监测应增加时间重复性，在必要时将采样深度延伸至足以覆盖参考土质量，并使用重复的管理和土壤分层来检验当前站点层面模式所暗示的机制。
