变暖引发的降水相态变化对青藏高原土壤侵蚀的影响

《Journal of Hydrology: Regional Studies》：Impacts of warming-induced precipitation phase shifts on soil erosion over the Qinghai-Tibet plateau

　　

中文翻译摘要 研究区域：研究覆盖了青藏高原（QTP），在该地区气候变暖正促使降水从固态向液态转变，可能改变土壤侵蚀风险。 研究焦点：研究人员利用修正通用土壤流失方程（RUSLE）估算侵蚀量，采用降雨比率（RRP）量化降水相态变化，以最大60分钟雨强（I60）表征事件尺度强度，并借助可解释的XGBoost-SHAP框架识别主导驱动因子。 区域水文学新见解：青藏高原大部分地区的RRP呈增加趋势，但其与降雨侵蚀力的关联表现出强烈的季节性依赖，秋季的耦合最为紧密（ρ = 0.683, p?<?0.001）。在年际尺度上，RRP与I60呈正相关（ρ = 0.710, p?<?0.001），而I60也与降雨侵蚀力呈正相关（ρ = 0.647, p = 0.002），表明降水相态变化、短历时降雨强度和侵蚀动力之间存在过程链式关联。降雨量（38.4%）和坡度（27.2%）仍是空间侵蚀异质性的首要控制因素，而RRP（7.8%）则作为次要调节因子。在海拔4500–5500米区间，侵蚀响应发生显著变化，此区域植被缓冲作用减弱，而对降水相态变化的敏感性增强。这些发现加深了对降水相态变化和事件尺度强度如何塑造侵蚀动力的理解，并为持续变暖背景下制定海拔特异性的水土保持措施提供了支持。

论文解读文章

研究背景与科学问题

土壤是陆地生态系统的基础资源，支撑着碳循环、水分调节和粮食生产等关键功能。水力侵蚀作为最普遍的土地退化形式之一，会剥离表土并输出土壤有机碳，从而损害土壤质量和生态系统服务。近几十年来，全球范围内侵蚀风险加剧，这反映了气候变暖和人类活动的综合影响。气候变暖强化了水文循环，导致频繁降雨和高径流事件，显著增加了动能和侵蚀动力。同时，土地利用变化、放牧和工程扰动减少了地表覆盖并破坏了土壤结构，降低了坡面的抗蚀性，加剧了土壤侵蚀。因此，在全球变化持续背景下，厘清侵蚀动态的主导控制因素及其空间异质性，对于制定针对性的水土保持和生态恢复策略至关重要。

青藏高原（QTP）是一个关键的高海拔、气候敏感区域，也是“亚洲水塔”的核心组成部分。冰川、永久冻土和季节性积雪构成了以冰冻圈为主的水文情势，控制着河源区的径流和泥沙动态，使该区域对温度变化高度敏感。长期的高山环境限制了高原上的成土过程，由于频繁的冻融循环，土壤通常薄而松散。一旦发生坡面剥蚀或沟蚀扩张，恢复过程往往极其漫长。近几十年来，高原经历了暖湿化趋势，这一变化加速了冰冻圈演变，改变了降水和径流的时机与能量，增加了侵蚀动力的变异性。在此背景下，即使植被覆盖持续改善，气候驱动的侵蚀力增强也可能抵消部分地区特定季节植被保护的效益，削弱生态恢复对侵蚀控制的边际效应。

尽管修正通用土壤流失方程（RUSLE）已被用于评估QTP上土壤侵蚀对降水总量、植被覆盖和人类活动的响应，但依赖于总量或平均状态可能会忽视侵蚀动力机制的关键变化。变暖引起的降雪向降雨的转变可能通过两条不同但相互作用的途径影响侵蚀。首先，用液态降雨替代降雪会增加土壤表面直接暴露于雨滴撞击的机会。与降雪相比，降雨将更大的动能传递到土壤表面，事件尺度强度的增加可以进一步增强土壤颗粒分离、溅蚀和降雨侵蚀力。其次，降水相态变化改变了水的储存和释放。降雪暂时储存降水并通常延迟径流，而降雨可以产生更快速的水文响应；雨雪混合事件还可能加速融雪并放大径流驱动的土壤分离和泥沙输移。与IPCC AR6关于寒冷地区降雪比例下降的证据一致，高海拔地区的极端降雨强度可能以每升温约15%的比例缩放，预估在高排放情景下QTP的雪转雨转变将尤为显著。这些变化表明，降水相态和事件尺度降雨强度可以为理解变暖如何改变高山区域的侵蚀动力提供一个有用的基于过程的视角。

基于上述物理机制的理解，本研究将降水相态转变视为分析高山区域水蚀动态空间异质性的关键切入点。具体而言，研究旨在解决三个相互关联的目标：（1）QTP上降雨比率（RRP）如何变化，以及在哪些季节RRP与最大60分钟雨强（I₆₀）和降雨侵蚀力（R因子）的关系最强；（2）降水量、地形、植被、土壤性质和降水相态结构如何共同作用于空间侵蚀异质性及非线性阈值响应；（3）降水相态变率、植被条件和土壤侵蚀之间的关系如何随海拔变化。通过整合降水相态变率、非线性归因和海拔依赖性响应，本研究提供了一个更具过程导向的理解，揭示了变暖相关的降水特征变化如何与高山区域的土壤侵蚀相关联。该论文发表于《Journal of Hydrology: Regional Studies》。

主要技术方法

为开展此项研究，研究人员采用了以下关键技术方法：1）利用ECMWF ERA5-Land再分析数据集获取气象强迫数据，并引入第三极高精度降水数据集（TPHiPr）进行对比验证，以确保降水数据的可靠性；2）采用动态阈值参数化方案（DTPS），以湿球温度（T_w）作为关键诊断变量来区分固态降水，并与ERA5-Land固态降水产品进行多尺度比较验证；3）运用修正通用土壤流失方程（RUSLE）模拟土壤侵蚀动态，计算降雨侵蚀力（R）、土壤可蚀性（K）、坡长坡度（LS）、覆盖管理（C）和支持措施（P）五个因子，量化2000–2020年的平均土壤流失量（A）；4）利用Spearman秩相关分析评估RRP、I₆₀和R因子之间的季节和年际关联，并通过95%偏差校正加速（BCa）自助法置信区间量化关键相关系数的不确定性；5）应用可解释的XGBoost-SHAP机器学习框架识别主导驱动因子并表征其非线性阈值响应；6）采用海拔带统计揭示侵蚀动态的海拔依赖性。所有输入层均统一至250米网格，投影至WGS_1984_UTM_Zone_47N坐标系，并在Google Earth Engine（GEE）平台上完成大规模预处理和因子推导。

研究结果

RRP的变化及其与I₆₀和R因子的关系**

通过分析QTP上RRP的空间分布和时序变化，研究发现RRP在高原大部分地区呈增加趋势，表明降水正从固态向液态转变。这种转变与降雨侵蚀力的关联表现出强烈的季节性依赖，其中秋季的耦合最为紧密（ρ = 0.683, p?<?0.001）。在年际尺度上，RRP与I₆₀呈显著正相关（ρ = 0.710, p?<?0.001），而I₆₀也与R因子呈显著正相关（ρ = 0.647, p = 0.002）。这一系列相关性证实了降水相态变化、短历时降雨强度和侵蚀动力之间存在过程链式关联：变暖导致的降雪减少和降雨增加，伴随着事件尺度降雨强度的增大，进而增强了降雨侵蚀力。

土壤侵蚀空间异质性的主导驱动因子及非线性阈值响应

利用XGBoost-SHAP框架对多维环境驱动因子进行归因分析，结果表明降雨量（贡献率38.4%）和坡度（贡献率27.2%）是控制QTP土壤侵蚀空间异质性的首要因素，而RRP（贡献率7.8%）则作为一个重要的次级调节因子发挥作用。此外，研究还识别出关键的非线性阈值响应。例如，当降雨量超过一定阈值后，其对侵蚀的促进作用趋于平缓；而在陡峭坡度条件下，侵蚀速率急剧上升。这些非线性关系揭示了不同环境因子在不同状态下对侵蚀影响的差异性。

侵蚀动态的海拔依赖性

通过对不同海拔带的统计分析，研究发现侵蚀响应在海拔4500–5500米区间发生了显著变化。在这一海拔范围内，植被的缓冲作用明显减弱，而土壤侵蚀对降水相态变化的敏感性却显著增强。这表明在中高海拔区域，随着植被覆盖度的降低和气温的进一步降低，降水相态转变（即降雪转为降雨）成为驱动侵蚀动态的一个更为关键的因素。低海拔地区植被覆盖良好，能有效拦截降雨能量；而极高海拔地区则以固态降水为主，侵蚀动力相对较弱。

总结讨论与研究结论

讨论部分强调了降水相态变化作为连接气候变暖与土壤侵蚀的关键中间过程的重要性。研究指出，传统的基于降水总量的侵蚀评估方法可能低估了变暖背景下侵蚀风险的加剧程度，因为降水相态转变和事件尺度强度的增加共同放大了侵蚀动力。研究还讨论了海拔梯度上植被与气候交互作用的复杂性，指出在特定海拔区间，气候变暖驱动的侵蚀力增强可能超过植被恢复带来的保护效益，这对当前的生态恢复实践提出了挑战。此外，研究结果支持了在高海拔地区实施海拔特异性水土保持措施的必要性。

研究结论总结如下：青藏高原大部分地区RRP呈增加趋势，且其与降雨侵蚀力的关联具有显著的季节性，秋季最强。年际尺度上，RRP、I₆₀和R因子之间存在显著的正相关链，揭示了降水相态变化通过影响短历时降雨强度进而调控侵蚀动力的过程机制。降雨量和坡度是控制QTP土壤侵蚀空间异质性的首要因素，而RRP作为次级调节因子，在特定条件下（如中高海拔区域）发挥着关键作用。侵蚀响应在海拔4500–5500米区间发生显著转变，表现为植被缓冲作用减弱和对降水相态变化敏感性的增强。这些发现不仅加深了对降水相态变化和事件尺度强度如何塑造高山区域侵蚀动力的理解，也为持续变暖背景下制定海拔特异性的水土保持策略提供了科学依据。