气候变化与土地利用变化情景下的土壤侵蚀：青尼罗河上游流域的一种集成建模方法

《Journal of Hydrology: Regional Studies》：Soil erosion under climate and land-use change scenarios: An integrated modelling approach for the Upper Blue Nile Basin

　　

研究区域：本研究在埃塞俄比亚青尼罗河上游流域的Gumara–Maksegnit流域进行，该区域是一个易侵蚀景观，排入塔纳湖（Lake Tana），一个水文和生态上至关重要的淡水系统。该地区以陡峭的悬崖、高度可蚀性土壤和快速扩张的农业为特征，使其成为流域内沉积物产生和土地退化的热点区域。 研究焦点：本研究通过整合RUSLE（Revised Universal Soil Loss Equation，修正通用土壤流失方程）模型与基于人工神经网络（ANN）的土地利用/土地覆盖（LULC）模拟、五个CORDEX（Coordinated Regional Climate Downscaling Experiment，协调区域气候降尺度实验）区域气候模式（RCM）的贝叶斯模型平均（BMA）以及偏差校正降雨预测，评估了过去（2003年）、现在（2023年）和未来（2083年）的土壤侵蚀。土地利用分析显示，农业用地从53.6%（2003年）扩大到68.7%（2023年），预计到2083年将达到72.1%，森林和放牧区相应减少。降雨侵蚀力（R因子）从565–647 MJ mm ha-1h-1yr-1增加到613–689 MJ mm ha-1h-1yr-1，预计到2083年将上升到756–853（RCP4.5）和809–907（RCP8.5）。平均年土壤流失量从2003年的30.1 t ha-1yr-1增加到2023年的39.2 t ha-1yr-1，预计仅土地利用变化下将达到42 t ha-1yr-1，在RCP 4.5和RCP 8.5气候情景下分别为43.4和45.9 t ha-1yr-1，当两种驱动因素结合时分别为46.3和49.2 t ha-1yr-1。 区域新水文见解：未来土壤流失在LUC+RCP8.5综合情景下将增加高达63.5%，北部和东北部陡坡仍然是持续的侵蚀热点。地形（LS因子：39.6%）是主导侵蚀驱动因素，而耕作坡地上的高C（覆盖管理因子）和P（支持实践因子）值加剧了沉积物产生。结果表明，气候变化比土地利用变化具有更强的边际影响，并且它们的相互作用非线性地放大了侵蚀，这强调了需要采取以坡度为焦点、适应气候的土壤和水资源保护措施，以保护塔纳湖的水文稳定性。

论文解读：《气候变化与土地利用变化情景下的土壤侵蚀：青尼罗河上游流域的一种集成建模方法》

研究背景与意义

土壤侵蚀是全球土地退化的主要驱动力，降低土壤肥力、农业生产力和生态系统稳定性。在埃塞俄比亚高地，由于陡峭地形、易蚀土壤、集约化耕作和植被覆盖下降，侵蚀热点持续存在。Gumara–Maksegnit流域是青尼罗河上游流域的重要支流，也是快速退化的景观，面临加速的气候和土地利用压力。气候变化通过温度升高、降雨增强和极端事件改变了径流生成和沉积物通量；同时，农业扩张、森林丧失和放牧地减少等土地利用/土地覆盖（LULC）转变继续降低植被覆盖并放大侵蚀风险。社会经济限制和不可持续的耕作方式进一步加剧了这些压力。尽管埃塞俄比亚有广泛的侵蚀研究，但很少有研究通过集成的多模型框架联合考察历史、当前和未来的气候-LUC相互作用，也没有研究将此方法应用于Gumara–Maksegnit流域。因此，研究人员通过开展这项研究，旨在填补这一空白，为流域管理提供科学依据。该论文发表在《Journal of Hydrology: Regional Studies》。

主要技术方法

研究人员采用了几项关键技术方法。首先，利用贝叶斯模型平均（BMA）框架对五个CORDEX区域气候模式（RCM）的降雨预测进行偏差校正和加权集成，生成历史时期、当前和未来时期的降雨侵蚀力（R因子）。其次，基于Google Earth Engine平台，使用随机森林（RF）分类器对2003、2013和2023年的Landsat影像进行分类，提取LULC信息；并利用MOLUSCE插件中的人工神经网络（ANN）和多层感知器（MLP）架构模拟2003至2083年的未来LULC变化。第三，应用修正通用土壤流失方程（RUSLE）模型，整合R、K（土壤可蚀性因子）、LS（坡度长度-陡度因子）、C（覆盖管理因子）和P（支持实践因子）六个因子，计算年土壤流失量。第四，采用皮尔逊相关分析、多元线性回归（MLR）和随机森林回归来评估各RUSLE因子对空间土壤侵蚀的相对影响。

研究结果

3.1 降雨和降雨侵蚀力趋势

研究发现，历史降雨范围为1021.1至1165.1毫米，平均1092.6毫米；当前降雨增至1107至1240.4毫米，相应的R因子从565.7–646.9 MJ mm ha^-1h^-1yr^-1升至613–689.2 MJ mm ha^-1h^-1yr^-1。未来RCP4.5和RCP8.5情景下，降雨将进一步增加，R因子预计分别达到756.6–852.6和809.4–906.5 MJ mm ha^-1h^-1yr^-1。这表明降雨侵蚀力持续上升，RCP8.5情景下增幅最大，意味着更强烈的侵蚀性降雨事件。

3.2 土壤可蚀性（K）

K因子在流域内介于0.1至0.38 t h MJ^-1mm^-1之间，较高值集中在南部和西南部耕地主导区，较低值在北部和东北部植被较好区域。这反映了土地利用和土壤管理的影响：连续耕作降低有机质和结构稳定性，增加可蚀性；森林土壤因有机质高而更稳定。

3.3 LS因子

LS因子变化巨大，从0到742.6，平均为6.8。高值出现在北部、东部和东北部陡坡和切割地形区，低值在靠近塔纳湖的平缓南部。LS量化了坡长和坡度对径流和土壤输送的增强作用，是山地景观中侵蚀的主要驱动因素。

3.4 土地利用/土地覆盖动态和C因子

过去二十年，农业用地从2003年的208.4平方公里（53.6%）扩展到2023年的267.1平方公里（68.7%），预计2083年达280.1平方公里（72.1%）。同期，放牧地从130.6平方公里（33.6%）降至80.2平方公里（20.6%），森林从48.8平方公里（12.6%）减至40.2平方公里（10.3%）。LULC分类精度很高，总体精度从85%升至93.4%。C因子（覆盖管理因子）随之变化：农业用地C值较高（0.15），森林最低（0.01），表明植被保护能力下降。

3.5 P因子

P因子（支持实践因子）反映土壤和水资源保护（SWC）措施的效果。农业用地根据坡度分配较低P值（0.1–0.33），非农业用地为1。随着农业向陡坡扩张，P值呈上升趋势，因为许多保护结构（如石坎、梯田）退化或缺乏维护，降低了保护效能。

3.6 基准期土壤侵蚀

2003年平均土壤流失量为30.1 t ha^-1yr^-1，2023年升至39.2 t ha^-1yr^-1，增长30.2%。极高值出现在深切河谷和陡坡，但面积很小（小于0.1%）。北部和东北部高地因陡坡和耕地扩张成为持续热点。

3.7 LUC和气候变化对土壤侵蚀的影响

情景分析显示，仅LUC情景下，2083年平均土壤流失为42 t ha^-1yr^-1（比2023年增6.7%）；仅气候驱动下，RCP4.5为43.4，RCP8.5为45.9 t ha^-1yr^-1；LUC+RCP4.5为46.3，LUC+RCP8.5为49.2 t ha^-1yr^-1，较历史条件增加63.5%。气候变化比LUC单独作用具有更大的边际影响，两者结合产生协同放大效应。

3.8 LUC+气候强迫下的侵蚀相互作用

皮尔逊相关和随机森林重要性分析表明，LS因子是最强驱动（贡献39.6%），其次是P（18.5%）、C（16.7%）、K（13.9%）和R（11.3%）。空间上，高降雨区往往位于缓坡，导致R和K与土壤流失呈负相关，凸显地形的主导作用。MLR系数强调P和C的关键作用，反映保护结构退化和植被减少。

3.9 局限性和管理启示

研究承认一些不确定性：RCM输出的固有变率、LULC模拟假设过去趋势延续、P因子基于文献可能低估实际退化情况，以及RUSLE未模拟沟壑和沉积过程。管理上，优先在陡坡和高LS区域实施石坎、梯田、检查坝等措施；恢复退化SWC结构；在农业区推广农林复合、草带和覆盖作物以降低C；适应气候变化的等高耕作和覆盖。

结论翻译

本研究首次对Gumara–Maksegnit流域的历史、当前和未来土壤侵蚀进行了综合评估，捕捉了土地利用转变和气候驱动的降雨增强的综合效应。通过整合RUSLE与基于ANN-MOLUSCE的LULC预测、BMA优化的CORDEX降雨情景以及随机森林因子诊断，分析揭示了过去20年及未来土壤退化长期增加的明确轨迹。土地利用分析显示了强烈且持续的农业扩张转变——从2003年的53.6%到2023年的68.7%——预计到2083年将上升到72.1%，主要以牺牲森林和放牧地为代价。这种转变增加了C和P值，并使陡峭景观暴露于加剧的侵蚀中。降雨侵蚀力也显著增加，从过去的565–647 MJ mm ha^-1h^-1yr^-1到当前的613–689 MJ mm ha^-1h^-1yr^-1，预计到2083年将超过850（RCP4.5）和900（RCP8.5）MJ mm ha^-1h^-1yr^-1。相应地，平均土壤流失量从2003年的30.1 t ha^-1yr^-1增加到2023年的39.2 t ha^-1yr^-1。未来预测表明，仅LUC下平均土壤流失为42 t ha^-1yr^-1，RCP4.5和RCP8.5下分别为43.4和45.9 t ha^-1yr^-1，LUC+RCP4.5和LUC+RCP8.5下分别为46.3和49.2 t ha^-1yr^-1。后者代表最严重的结果，与历史条件相比增加了63.5%。地形被确定为主要侵蚀驱动因素：LS占随机森林重要性的39.6%，并与土壤流失显示出最强的相关性。然而，土地覆盖（C）和保护有效性（P）在陡坡上强烈调节敏感性。这些相互作用解释了北部和东北部高地持续的热点，那里农田扩张与高LS和不断上升的降雨侵蚀力重叠。结果还表明，气候变化对未来侵蚀的边际影响大于单独的LUC，它们的综合效应产生了超越各自贡献的协同影响。管理启示很明确：必须优先考虑通过梯田、石坎和检查坝稳定坡面；修复和维护现有SWC结构；在高LS区域重新造林和农林复合；以及采用残茬保留和等高耕作等适应气候的农业做法。在确定的热点实施这些措施将最大程度地减少土壤侵蚀和流入塔纳湖的沉积物。总之，该研究强调迫切需要协调的土地利用规划和气候响应型保护，以在日益增加的环境变化中保护土壤生产力和流域韧性。