《International Soil and Water Conservation Research》：Assessment of vegetation sediment reduction potential based on future climate system models, Loess Plateau, China
编辑推荐：
摘要专业翻译
近几十年来，黄土高原发生了显著变化。水土保持措施改变了侵蚀、储存与输沙之间的平衡。为检验植被措施的 climate robustness（气候稳健性），研究人员使用了第六次耦合模式比较计划（CMIP6）中五个全球气候模式（GCMs：CanESM5
摘要专业翻译
近几十年来，黄土高原发生了显著变化。水土保持措施改变了侵蚀、储存与输沙之间的平衡。为检验植被措施的 climate robustness（气候稳健性），研究人员使用了第六次耦合模式比较计划（CMIP6）中五个全球气候模式（GCMs：CanESM5、EC-Earth3、IPSL-CM6A-LR、MIROC6 和 MRI-ESM2-0）在四种共享社会经济路径（SSP）情景（SSP1-2.6、SSP2-4.5、SSP3-7.0 和 SSP5-8.5）下的预估数据。研究人员采用 delta method（增量法）对这些数据进行降尺度处理，并将降尺度结果与站点观测进行了验证。随后，利用中国土壤流失方程（CSLE）在流域尺度上模拟了植被对产沙量的影响。该方案解析了未来侵蚀的时空格局，并追踪了植被减蚀能力随情景的变化。增量法有效降低了偏差，尤其是对温度和降水的偏差，且多模式集合与观测吻合最佳。两个变量在高原上均预计呈上升趋势。在 SSP1-2.6、SSP2-4.5、SSP3-7.0 和 SSP5-8.5 情景下，年平均降水量预计分别为 503.8、539.9、513.6 和 538.4 mm，其中 SSP2-4.5 情景既呈现最高的时段平均降水量，也表现出最陡峭的时间递增趋势。这些气候变化改变了与植被相关的侵蚀控制效应，并产生了显著的空间异质性。年均侵蚀强度排序为 SSP2-4.5 > SSP5-8.5 > SSP3-7.0 > SSP1-2.6。许多区域预计将超过植被减蚀阈值，尤其是在陕北、山西和高原南部。这些结果阐明了未来侵蚀机制，并为流域尺度保护规划提供了支撑。
论文解读文章
一、研究背景与科学问题
土壤侵蚀是全球范围内导致土地退化的主要驱动因素，也是影响生态系统稳定性、农业生产力和可持续土地利用的关键环境问题。通过剥离表土和加速养分流失，土壤侵蚀降低了土壤肥力、扰乱了水文过程，并引发更广泛的生态退化。植被作为陆地生态系统的重要组成部分和气候变化的敏感指示器，在抑制土壤侵蚀方面发挥着关键作用。
黄土高原是世界上侵蚀最严重的地区之一，长期以来一直是黄河泥沙的主要来源。然而，近几十年来，由于重大生态保护与修复工程的实施，黄土高原已成为中国生态恢复最显著的范例之一。过去一个世纪，黄土高原的泥沙输入量减少了近85%，从1919—1959年间的16亿吨/年降至2000—2018年间的2.4亿吨/年。这些显著的侵蚀控制成效被广泛归因于退耕还林工程下的大规模植被恢复。
尽管生态恢复已大幅减少了土壤侵蚀和产沙量，但未来气候变化可能通过改变降水格局和降雨侵蚀力，为这些成果引入新的不确定性。由于降雨侵蚀力是侵蚀强度的主要驱动因子，气候的变化可能影响未来土壤侵蚀的幅度和空间分布，即使在植被覆盖已改善的区域也是如此。因此，评估未来侵蚀风险需要同时考虑植被恢复的持续效益和气候变化对侵蚀驱动力的潜在影响。
目前尚不清楚在不同气候变化情景下，植被能在多大程度上持续抵消未来的侵蚀风险。为此，本研究将CMIP6在四种SSP情景下的五个GCM预估数据应用于黄土高原，采用增量法进行降尺度处理并与站点观测进行验证，基于CSLE模型推演2020、2030、2040和2050年的侵蚀格局，并检验不同SSP情景下植被对减沙的贡献，进一步评估植被控沙能力能够维持多久，为黄河流域水沙趋势预测提供科学依据。
二、主要技术方法
研究人员采用了以下关键技术方法：一是统计降尺度方法，利用黄土高原295个气象站点的日观测数据，采用增量法对CMIP6五个GCM模式（CanESM5、EC-Earth3、IPSL-CM6A-LR、MIROC6、MRI-ESM2-0）的日降水、温度、短波辐射和相对湿度数据进行偏差校正和空间降尺度，并通过两折时间交叉验证评估降尺度性能；二是饱和水汽压差（VPD）计算，基于温度和相对湿度估算VPD，作为植被预估中的气候预测因子；三是中国土壤流失方程（CSLE）模型，综合降雨侵蚀力（R）、土壤可蚀性（K）、坡长坡度因子（LS）、生物措施因子（B）、工程措施因子（E）和耕作措施因子（T）估算土壤侵蚀模数；四是因子贡献分解方法，将未来侵蚀变化分解为降雨侵蚀力变化贡献和植被覆盖变化贡献；五是植被减沙潜力评估方法，以当前植被减蚀模数与未来阈值之差定义植被减沙潜力。土壤属性数据来源于协调世界土壤数据库（HWSD），地形高程数据来源于NASA的SRTM DEM，所有栅格统一重采样至500米分辨率。
三、研究结果
3.1 增量降尺度模型评估
通过两折时间交叉验证，降尺度后的温度与观测高度吻合，五个模式的R2值分别为0.944、0.964、0.951、0.948和0.958，回归斜率介于0.982至1.029之间，其中EC-Earth3表现最优。降水的表现则更具变异性，R2值介于0.460至0.681之间，MIROC6表现最佳，MRI-ESM2-0吻合度最弱。多模式集合在温度方面R2达0.983，降水方面R2达0.811，显著改善了与观测的一致性。校正后的模拟在主要雨季（6—9月）仍存在对降水的低估。
3.2 预测期气象要素的时空动态
以2011—2020年为基线，2030—2050年间四种情景下的多年平均温度分别为9.38°C、9.48°C、9.42°C和9.87°C，SSP5-8.5情景升温最为显著。降水整体呈增加趋势，增长率分别为2.422、4.319、1.053和1.203 mm/年，SSP2-4.5情景递增趋势最陡。相对于基线期，各情景降水变化幅度介于-12%至21%之间，最大增幅出现在SSP5-8.5情景。空间上，黄土高原呈现西北—东南降水梯度。
3.3 基线期土壤侵蚀强度的时空变化特征
2020年基线期，黄土高原降雨侵蚀力自西北向东南增强，平均值为1326.49，61.49%的区域低于1000。土壤可蚀性因子K在中心丘陵沟壑区最大，LS因子在西部较高。生物措施因子B平均值为0.18，高值区主要分布在西北部植被茂密区域，低值区集中在东南部。2020年平均侵蚀模数为23.86 t hm?2 a?1，近88%的区域处于中度以下水平，剧烈侵蚀（>80 t hm?2 a?1）集中在丘陵沟壑区和黄土塬区。
3.4 未来黄土高原土壤侵蚀强度的时空变异性
在保持K、LS、E、T因子不变的条件下，仅改变R因子和B因子进行预估。2030、2040和2050年四种情景的平均降雨侵蚀力约为1198、1438、1318和1340，SSP2-4.5情景下2050年R值最大，达1675.87。B因子总体稳定，介于0.164（2050年SSP2-4.5）至0.194（2050年SSP1-2.6）之间。侵蚀响应呈现明显的场景间差异：从2030年到2040年，SSP1-2.6、SSP2-4.5和SSP5-8.5情景下平均侵蚀模数预计下降，而SSP3-7.0情景下预计上升。到2050年，多个情景下土壤流失趋于再次增加，尤其是SSP2-4.5和SSP5-8.5情景。空间上，侵蚀呈现西北向东南加剧的格局，西部丘陵沟壑区和黄土高原沟壑区为侵蚀热点区域。
3.5 未来情景下黄土高原植被减沙潜力
在SSP1-2.6情景下，预测期植被减蚀阈值为31.4 t hm?2 a?1，2030、2040和2050年的减沙潜力分别为-15.6、-10.4和-15.3 t hm?2 a?1，超过阈值的面积占比分别为52.2%、53.9%和49.8%。在SSP2-4.5情景下，阈值为35.2 t hm?2 a?1，潜力分别为-18.3、-14.1和-23.6 t hm?2 a?1，超阈值面积占比54.0%、51.1%和55.6%。在SSP3-7.0情景下，阈值为33.5 t hm?2 a?1，潜力分别为-15.3、-19.4和-15.7 t hm?2 a?1，超阈值面积占比51.0%、53.9%和50.5%。在SSP5-8.5情景下，阈值为34.0 t hm?2 a?1，潜力分别为-16.6、-16.4和-18.4 t hm?2 a?1，超阈值面积占比50.8%—51.3%。超阈值区域主要分布在陕北、山西、河南及内蒙古北部等地。
四、讨论与结论
讨论部分指出，黄土高原降水极端事件增加伴随着降雨侵蚀力因子的上升，SSP2-4.5情景下降水在2050年达到最高水平，相应降雨侵蚀力也最大，表明气候对土壤侵蚀具有更强的强迫作用。退耕还林工程通过增加冠层截留、地表覆盖和入渗，有效减少了土壤剥离和径流，显著降低了产沙量。预估显示土地利用总体变化不大，但林草地面积和覆盖度略有增加，意味着侵蚀压力将继续缓解。然而，随着气候变化日益复杂，植被能否维持当前的减沙效益是一个亟待回答的问题。
研究存在以下局限性：月降水数据难以完全表征短历时强降雨事件的频率和强度；受数据限制，淤地坝和水库等工程措施未能显式建模；所选GCM在降水预估方面存在不确定性；研究仅关注年际变化，未评估季节和月际变化。
研究结论如下：（1）验证表明增量降尺度具有较强时空拟合能力，多模式集合均值误差最低，最适合黄土高原。预测期各情景温度上升速率分别为0.014、0.041、0.056和0.077 °C/年，降水上升速率分别为2.422、4.319、1.053和1.203 mm/年。（2）四种SSP情景下预测的土壤侵蚀强度呈现非单调时间格局，大多数情景下侵蚀预计在2040年下降、2050年再次上升。SSP2-4.5情景的平均侵蚀强度最大。空间上以西北—东南梯度为主导格局，西部丘陵沟壑区和黄土高原沟壑区持续为侵蚀热点。（3）与2020年相比，2030年四种情景的植被减蚀潜力分别为15.6、18.3、15.3和16.6 t hm?2 a?1，2040年为10.4、14.1、19.4和16.4 t hm?2 a?1，2050年为15.3、23.6、15.7和18.4 t hm?2 a?1。超过阈值的面积占比在2030年为52.2%、54.0%、51.0%和50.8%，2040年为53.9%、51.1%、53.9%和51.3%，2050年为49.8%、55.6%、50.5%和51.3%。这些区域主要集中在陕北、山西和黄土高原的河南部分，表明在未来气候变化下，植被预计将在这些地区持续提供显著的减蚀效益。
