驱动因素对黄河中游泥沙负荷的解耦效应

《CATENA》:Decoupling effects of driving factors on sediment load in the middle reaches of the Yellow River

【字体: 时间:2026年01月07日 来源:CATENA 5.7

编辑推荐:

  黄河中游河口-龙门区泥沙负荷时空演变及驱动机制研究。基于1960-2020年连续水文数据,采用PLS-SEM模型解析植被、气候与径流对泥沙负荷的协同作用,揭示植被-径流耦合主导机制及空间异质性。结果表明:1980-1999年气候主导(总效应0.37),2000年后植被主导(总效应-0.33);径流对泥沙变率影响极显著(P<0.01);Tou-Wu区气候驱动(-0.85~-0.52),Wu-Long区植被驱动(-0.41~0.57),植被通过径流的中介效应(-0.55~0.36)显著超越直接效应(-0.10~0.10)。为黄土高原泥沙治理提供时空靶向管理依据。

  
马超辉|郭晓雪|徐瑞瑞|高鹏|周佳慧|张玉轩|穆兴民
中国西北农林科技大学水土保持与荒漠化控制国家重点实验室,水土保持研究所,陕西省杨陵市712100

摘要

广泛的水土保持措施深刻改变了黄河中游的泥沙输送情况。科学量化这些变化背后的驱动因素对于有效的水域管理至关重要。以往的研究主要关注气候和人类活动对整个流域泥沙负荷的直接影响,常常忽略了径流介导的途径和空间异质性。本研究利用1960年至2020年的连续水文记录,分析了黄河中游及其子流域泥沙负荷的时空变化。通过偏最小二乘结构方程建模(PLS Structural Equation Modeling),阐明了植被和气候对泥沙负荷的直接影响,以及驱动因素的空间异质性。结果显示,1960年至2020年间泥沙负荷显著下降(P < 0.01)。自2000年以来,吴龙地区已成为主要的泥沙来源。研究区域的泥沙负荷驱动因素表现出时间和空间的变异性。时间上,1980年至1999年气候起主导作用(总效应 = 0.37),而2000年至2020年植被成为主导因素(总效应 = ?0.33)。空间上,头吴地区主要受气候影响(?0.85 ≤ 总效应 ≤ ?0.52),而吴龙地区则以植被为主(?0.41 ≤ 总效应 ≤ 0.57)。径流对泥沙负荷变化有显著影响(P < 0.01)。在头吴和吴龙地区,植被通过径流的间接效应(分别为?0.55至0.36和?0.59至0.32)超过了直接效应(分别为?0.10至0.10和0.18至0.25),证实了径流是关键的中介因素。本研究揭示了河口-龙门地区泥沙负荷驱动因素的空间异质性,为黄河中游的精准管理和可持续生态提供了科学依据。

引言

河流泥沙对地表水文过程至关重要,因为它促进了物质循环和能量转移。它与生物地球化学循环、河口栖息地的发展、泛滥平原农业土壤中的养分动态以及水质密切相关,这些因素共同决定了生态系统的健康和社会经济的可持续性(Best, 2019; Wang and Wang, 2023; Yin et al., 2023)。近年来,气候的加剧变化和人类活动的频繁干预导致全球许多主要河流流域的泥沙负荷发生了显著变化(Best, 2019; Cohen et al., 2014; Song et al., 2020; Walling and Fang, 2003)。大约40%的流域(Li et al., 2021),尤其是亚洲的主要河流(Das et al., 2025; Wang et al., 2011),经历了泥沙负荷的显著减少。这些变化引发了一系列新问题,包括河道不稳定、河口栖息地退化、水库淤积以及传统水土保持措施的效力下降(Li et al., 2024; Liu et al., 2018b; Lu et al., 2011; Ma et al., 2022)。
黄河是世界第五长的河流,输送的泥沙量最大,是全球水文循环的重要组成部分。其中游流经黄土高原,该地区以严重的土壤侵蚀和水资源枯竭为特征,每年产生超过8000吨/平方公里的泥沙(Zhao et al., 2015)。这一区域占黄河总泥沙流入量的约90%(Gao et al., 2017b)。在该区域内,河口-龙门地区(从头道沟站到龙门站,简称头吴)的土壤侵蚀最为严重。尽管该地区仅占黄河流域总面积的14.8%,但它产生了河流中约73%的粗粒泥沙(Gao et al., 2017b),使其成为该流域泥沙相关问题的主要来源(Zhang et al., 2022a)。多年来,该地区的泥沙输送率一直接近1(Zheng, 2017),加剧了下游河道的淤积,导致河床抬高(Ma et al., 2022; Xu et al., 2022)。此外,气候变化加剧了黄河流域和华北平原的洪水频率和强度(Hessel et al., 2003),对公共安全构成威胁,并威胁区域社会经济稳定。近年来,河口-龙门地区实施了多种水土保持措施,包括“ Grain for Green”计划、拦沙坝、梯田和水库,以减少土壤侵蚀。这些干预措施有效控制了侵蚀,显著减少了河流中游的泥沙产量(Wang et al., 2015)。2010年至2017年间,该段的年均泥沙负荷降至1.1294 × 10^6吨,仅相当于20世纪70年代前的8%(Kong et al., 2023)。然而,泥沙负荷的急剧减少带来了新的挑战。泥沙不足的流动增强了冲刷能力,导致河床切割和黄河三角洲的严重退化(Li et al., 2024; Qiao et al., 2024; Yi et al., 2022; Yin et al., 2023)。此外,许多水库和拦沙坝由于基础设施老化,现在面临淤积和储水能力不足的问题(Ran et al., 2013, Ran et al., 2021)。当前的情况表明,传统的水土保持措施已不再适应现代流域动态。
以往的研究通常采用三种主要方法来解释泥沙负荷的变化:经验关系法、半经验公式法和水文模型(Devia et al., 2015)。经验关系法简单、高效且相对准确,因此被广泛使用(Miao et al., 2011; Zhang et al., 2022a)。然而,这种方法仅依赖统计分析来建立因变量和自变量之间的数学关系,忽略了潜在的物理机制,仅关注输入和输出之间的统计相关性。在泥沙归因分析中,这种方法无法捕捉泥沙变化的动态过程和特征,从而降低了准确性(Zhang et al., 2019)。相比之下,水文建模结合了降雨、入渗、泥沙输送和径流等物理过程,在研究区域的空间单元中模拟不同情景下的泥沙产量(Sun et al., 2020b; Wang et al., 2022b; Yang et al., 2024)。Liu et al.(2024)在汾河流域使用土壤和水资源评估工具(SWAT)发现,保土坝、土地利用和气候变化分别对泥沙输送贡献了33-35%、34-36%和19-22%。尽管水文模型提供了高精度估计,但它们需要大量数据且实施复杂。半经验公式法则结合了物理机制和统计原理,构建了简化的水和泥沙过程模型。这种方法便于数据收集,支持更合理的影响因素选择,并提供了比纯经验方法更强的解释力(Shi et al., 2019; Yao et al., 2016; Zheng et al., 2018)。例如,在黄土高原,Zhang et al.(2019)结合了Budyko假设和分形理论来研究泥沙输送变化的驱动因素,明确考虑了人类活动引起的径流和泥沙浓度变化。他们的结果表明,气候、人类活动和径流分别对泥沙负荷变化的贡献为44%-6%和43%-76%,其中46%和16%的变化可归因于径流和泥沙浓度的波动。类似地,Das(2021)和Das et al.(2023)应用偏最小二乘回归(PLSR)方法,确定了印度半岛十二条主要河流及其最大流域戈达瓦里河流域泥沙输送的控制因素。后一项研究进一步量化了地貌、岩性和气候在调节泥沙负荷动态中的相对重要性,分别解释了十二条河流和戈达瓦里河流域87%的观测方差。与传统多变量统计分析相比,PLS–SEM模型基于物理水文过程选择适当的参数,提供了对多个因素之间复杂因果关系的更深入理解(Jing et al., 2015; Ren et al., 2021; Schweiger et al., 2016)。
总之,许多定量研究探讨了影响黄土高原泥沙减少的因素。然而,大多数研究直接将泥沙减少归因于人类活动和气候变化,而忽视了它们通过径流变化对泥沙动态的间接影响。此外,许多以往的研究使用离散水文站的泥沙负荷数据推断整个流域的泥沙趋势。这种方法在空间尺度转换过程中引入了不确定性,难以全面准确地描述整个流域的水和泥沙分布的空间变异性(Gao et al., 2017a; Ni et al., 2023)。为了解决以往研究的局限性,本研究的主要目标是:(1)将研究区域划分为两个子流域,利用相邻主流水文站之间的泥沙负荷差异以及支流控制站的泥沙负荷数据,准确全面地描述河口-龙门地区泥沙输送的时空变化;(2)基于气候、植被和径流数据建立PLS–SEM模型,阐明链式驱动机制,并识别单个潜在变量对泥沙负荷变化的影响;(3)识别主要控制泥沙负荷变化的区域特定因素。本研究旨在阐明气候因素和植被覆盖对河口-龙门地区泥沙负荷变化的直接影响,以及通过径流介导的间接影响,从而为流域内的水和泥沙的精准管理提供科学依据。

研究区域

河口-龙门地区位于黄河的上游和中部,地理坐标为北纬35°25′02″–40°39′22″,东经107°13′48″–112°53′49″,海拔高度在1,100至1,500米之间。主河道长度为723公里,平均坡度为0.84%,流域面积为111,586平方公里(图1)。该区域是黄河泥沙产量和粗粒泥沙的主要来源,也是土壤侵蚀最严重的区域。

年际变化

根据1960年至2020年头吴及其两个子流域的泥沙负荷变化观测(图2a),主流域、头吴和吴龙子流域的多年平均泥沙负荷分别为5.21、2.69和2.53 × 10^8吨。这三个区域的泥沙负荷均显著下降(p < 0.01),下降率分别为?0.16、?0.10和?0.06 × 10^8吨/年。根据BGHS算法结果,1980年和1999年被确定为转折点。

水文要素对泥沙变化的影响

降水是流域水文循环的基本要素,也是土壤侵蚀的主要驱动因素。其时间和空间变化直接影响侵蚀和径流过程,最终影响泥沙产量(Wang et al., 2022b; Zhao et al., 2014)。由于厚厚的、高度易侵蚀的黄土沉积物和频繁的季风暴雨,河口-龙门地区极易发生水蚀(Ping et al., 2013)。长期统计分析表明……

结论

本研究分析了1960年至2020年河口-龙门地区的泥沙负荷变化,该地区是黄土高原粗粒泥沙的主要来源。应用PLS–SEM模型量化了植被、气候和径流对泥沙负荷的影响。结果揭示了显著的环境变化,包括温度升高、土地利用和植被覆盖的变化,以及广泛的水土保持措施。这些因素共同导致了……

作者贡献声明

马超辉:写作——审稿与编辑、初稿撰写、可视化、软件使用、方法论研究、概念构思。郭晓雪:验证、监督、方法论指导。徐瑞瑞:监督、概念构思。高鹏:方法论研究、资金筹集。周佳慧:监督、数据管理。张玉轩:监督、数据管理。穆兴民:监督、项目管理。

利益冲突声明

作者声明没有已知的财务利益或个人关系可能影响本文的研究结果。

致谢

本项工作得到了中国国家自然科学基金(U2243211)的资助。感谢“国家地球系统科学数据中心”的数据支持(https://www.geodata.cn)。作者还要感谢黄河水利委员会的水文局提供宝贵的气候和水文数据。同时感谢审稿人的意见,这些意见大大提高了论文的质量。
相关新闻
生物通微信公众号
微信
新浪微博
  • 搜索
  • 国际
  • 国内
  • 人物
  • 产业
  • 热点
  • 科普
  • 急聘职位
  • 高薪职位

知名企业招聘

热点排行

    今日动态 | 人才市场 | 新技术专栏 | 中国科学人 | 云展台 | BioHot | 云讲堂直播 | 会展中心 | 特价专栏 | 技术快讯 | 免费试用

    版权所有 生物通

    Copyright© eBiotrade.com, All Rights Reserved

    联系信箱:

    粤ICP备09063491号