【特稿】
穆兴民1，2，胡春宏3，高 鹏1，2，王 飞1，2，赵广举1，2
（1.西北农林科技大学黄土高原土壤侵蚀与旱地农业国家重点实验室，陕西杨凌712100；2.中国科学院水利部水土保持研究所，陕西杨凌712100；3.中国水利水电科学研究院，北京100038）
摘 要：近十几年黄河输沙量的突兀性减少是社会各界人士普遍未曾预料到的。黄河输沙量减少原因和未来趋势亟待研究，也需结合现状进一步提出新的适应策略。基于对黄土高原水土保持与黄河水沙变化等的长期研究，提出黄河输沙量研究需要关注的几个关键问题和认识：考虑到学科特点及发展水平，黄河输沙量以及人类活动对输沙量变化影响的研究“宜粗不宜细”；黄河输沙量变化研究可以年代为时段开展对比分析，1979年之前可作为基准期；在流域及区域侵蚀产沙性降雨的相关研究中，宜把12 mm日降雨量作为临界侵蚀降雨量指标；加强淤地坝淤积记录反演研究，揭示土壤侵蚀强度变化，以解析坡面土壤侵蚀减少与沟道坝库工程拦蓄对黄河输沙量减少的贡献，明确当前黄土高原水土流失治理方向；强化黄河未来输沙量变化趋势及治理对策研究。
关键词：输沙量；降雨；人类活动；黄土高原；土壤侵蚀；黄河
黄土高原的独特性质决定了黄河泥沙研究是一个永恒的课题。进入21世纪以来，黄河输沙量突兀性减少，输沙量和含沙量之低是社会各界普遍未曾预料到的，黄河泥沙减少这一重大问题再一次摆在科技工作者面前。黄河泥沙变化原因及未来状态关系到黄河及黄土高原的治理对策和决策，受到广大科技工作者和管理部门的高度关注。目前，有关高校及科研单位正在组织各方力量开展深入研究，笔者认为，当务之急是从宏观方面研究和解决如下几个关键问题。
黄河输沙量的年际变化过程具有显著的阶段性特征，其阶段可以采用累计距平法、双累积曲线法、有序聚类方法等［1-3］清晰辨识。不同阶段的输沙量均值、变差系数等存在显著差异。就黄河中游吴堡、龙门以及潼关水文站而言，输沙量突变点均发生在1979年；中游各支流水文站年输沙量发生突变的年份虽有所差异，但主要集中在20世纪70年代中后期［3-5］。不同学者采用不同方法辨识黄河干流或支流输沙量发生突变的更精确年份，结果往往大同小异。
目前，在水土保持措施等各种人类活动对产沙影响研究的定量评估中，参考期（基准期）及变化期各时段具体年份的划分“五花八门”。通过绘制面平均降雨量与输沙量的双累积曲线，分析发现黄河河口镇至潼关区间（简称河潼区间）输沙量存在三个显著的阶段（见图1），即20世纪70年代之前、80—90年代以及21世纪以来三个时段。
图1 黄河河潼区间降雨量与输沙量双累积曲线
根据我国使用年代的习惯，特别是在各种研究和政府文件中广泛使用年代的习惯，以及水土保持学科本身并非能达到很高精度，更出于便于比较同类研究成果，因此笔者建议：采用年代制，把基准期（参照期）临界年份确定为1979年，重点研究1979年之前、1980—1999年以及2000年以来三个阶段黄河及其支流输沙量变化情况。对个别支流，2000年（临界年份）输沙量突变并不显著，可以采用两段制进行分析。河潼区间输沙量阶段特征见表1。
表1 河潼区间年输沙量阶段性变化特征亿t
时段平均最大（年份）最小（年份）1958—1979 13.48 27.83（1958）3.72（1965）1980—1999 7.16 13.32（1988）3.18（1987）2000—2013 2.40 5.90（2003）0.65（2009）
河流输沙量受气候、地质地貌和人类活动等因素影响，理论上讲，可表示为气候、地质地貌和人类活动这3个因素的函数［6］：
式中：S表示输沙量；C表示气候因素；H表示人类活动因素（包括各种引起下垫面条件改变的因素，如水土保持工程、植被变化、开发建设项目等）；G表示土壤及地质地貌因素。
对于黄河流域或其他区域，在不足百年时间尺度的条件下，流域的土壤及地质地貌条件（G）可认为相对不变，即河流输沙量变化主要取决于气候和人类活动两大因素，因此式（1）可简化为
黄土高原土壤侵蚀以水蚀为主，气候因素常用降雨特征来表征，如降雨量、降雨日数以及降雨强度等。在坡面土壤侵蚀研究中，周佩华等［7］、王万忠［8］分别根据径流小区等小尺度上的试验结果，提出了小尺度侵蚀性雨量标准和雨强标准，其中土壤侵蚀雨量标准大致为1 h雨量10～13 mm、6 h雨量25～30 mm、24 h雨量50～60 mm。
对于大中尺度流域的河流输沙问题，基于小尺度的降雨指标进行研究显然不妥。另外，基于小尺度的试验结果，在中大尺度开展研究时，其资料的获取非常困难，流域尺度上的各站点降雨指标整合更是一个难以实现的问题。之前有关研究者提出用年降雨量、汛期降雨量或月降雨量作为降雨指标，但对降雨特征信息的表征都不够全面和准确。根据气象资料共享平台，日降雨资料可以直接获得。基于此，提出“流域侵蚀产沙性降雨”概念及流域侵蚀产沙雨量标准［8］。基于流域日降雨量，根据参考期的面平均日降雨量—输沙量累积百分关系曲线，按照流域侵蚀产沙性降雨所引起的输沙量占总输沙量的90%，确定流域侵蚀产沙性降雨的日雨量阈值。以皇甫川流域研究为例，考虑到资料的易获得性和降雨侵蚀力计算需要，确定日降雨12 mm为流域侵蚀产沙性降雨标准［9］。以1979年以前为参照期，其年径流量和输沙量与降雨指标的决定系数见表2，分析表明用大于等于12 mm的日降雨量可以较好地表征降雨对侵蚀产沙的作用。
表2 皇甫川流域径流量、输沙量与降雨指标的决定系数
注：∗∗∗表示达到0.01的显著水平，∗∗表示达到0.05的显著水平
降雨指标年径流量年输沙量年侵蚀产沙性降雨量0.71∗∗∗0.74∗∗∗年侵蚀产沙性降雨日数0.62∗∗∗0.65∗∗∗年降雨量0.60∗∗0.68∗∗年降雨日数-0.20-0.08
降雨径流是黄土高原土壤侵蚀的核心驱动因素，河流输沙量与流域降雨量关系密切，采用双累积曲线法可以解析降雨和人类活动（下垫面变化）对输沙量变化的作用程度（见图2），因此建议采用双累积曲线法区分降雨和人类活动的作用，尽管各种人类活动的作用是一笔复合在一起的“混账”。黄河中游不同时段降雨及人类活动对输沙量减少的作用分析结果见表3。由表3可知：20世纪80—90年代人类活动对黄河输沙量减少的贡献率接近80%，进入21世纪后人类活动的作用已超过90%。20世纪80—90年代输沙量的减少与黄土高原梯田、淤地坝的规模化建设相对应，而进入21世纪以来黄河输沙量的减少则主要与退耕还林（草）植被恢复过程相对应。
图2 1957—2013年黄河中游降雨量与输沙量双累积曲线
表3 黄河中游不同时段降雨及人类活动对输沙量减少的作用
时段实测年输沙量／亿t计算年输沙量／亿t实测输沙量变化情况降雨影响人类活动影响减少量／亿t减幅／%减少量／亿t占比／%减少量／亿t占比／% 1979年以前11.25 11.16 1980—1999年6.60 10.30 4.65 41.3 0.95 20.4 3.70 79.6 2000—2013年0.61 10.53 10.64 94.6 0.72 6.8 9.92 93.2
人类活动对流域侵蚀产沙量的影响非常复杂，既有减少土壤侵蚀量的作用，也有增加土壤侵蚀量的作用。就黄河输沙量减少的原因，如果能够明确是坡面土壤侵蚀量减少还是河道坝库拦蓄，那么就能明确黄土高原水土保持的基本方向。理论上，就水土保持措施而言，相同数量某项措施因其空间配置不同、水文连通性也不同，故其减沙贡献度会产生显著差异。如淤地坝会因其在沟道上下游位置的不同而产生差异，但这种差异目前仍然难以测度和验证。因此，鉴于目前水土保持措施等人为活动的作用机制尚不十分清楚以及受学科发展水平的限制，人类活动对土壤侵蚀产沙作用贡献率的研究“宜粗不宜细”，不必纠缠于各项措施作用贡献率的“精确”辨识。如果能够明确人类活动的作用在减少坡沟系统（即沟沿线以上的坡面及沟沿线以下的沟坡）土壤侵蚀和河道坝库淤积的程度，那么就可明确黄土高原水土流失治理的方向或对象。
一般来说，流域土壤侵蚀产沙量平衡方程可以表示为
式中：SES表示流域土壤侵蚀产沙量；OS表示水文站测定输沙量；BS表示坝库拦蓄量；SS表示河道冲淤变化量；PS表示灌溉等引沙量。
流域土壤侵蚀产沙量（SE S）实际上也是不可测度的，传统的径流小区所测定的侵蚀产沙量并不能代表一个流域或区域的土壤侵蚀产沙量。如果说可以测度的话，那么就可以在支毛小沟上设置卡口站进行测定，若干个支毛小沟测定的均值可以代表一个区域或流域的土壤侵蚀产沙水平。
对于黄土高原中小流域，式（3）中的水文站测定输沙量（OS）可采用实测值，河道冲淤变化量（SS）对多数流域可忽略不计，灌溉等引沙量（P S）可以相对准确测算，因此估算人类活动的作用就可以简化为测定河道坝库（其上游无坝库影响）的拦蓄量。
通过淤地坝（坝地）的泥沙淤积信息研究小流域土壤侵蚀产沙，是目前研究流域或区域土壤侵蚀的重要方法。通过对淤地坝沉积旋回的研究，结合库容曲线获取淤地坝内泥沙淤积量，反演控制区内的土壤侵蚀产沙强度［10-11］（见图3，图中从左至右分别为淤积剖面的137Cs活度、泥沙颗粒组成、沉积旋回层剖面、采样剖面），进而利用泥沙理化性质识别其主要来源。通过空间上若干个淤地坝淤积过程的研究并结合淤地坝建设情况，进而推算流域或区域的坝库拦蓄泥沙量。如果可以调查和统计淤地坝及水库泥沙淤积数量，进一步明确黄河泥沙减少的原因主要是坡沟系统土壤侵蚀产沙量减少还是河道坝库工程的拦蓄，那么就可以明确治理黄河泥沙的重点是加强坡沟系统的水土流失治理还是沟道（河道）坝库工程建设。
图3 皇甫川小石拉他坝泥沙淤积信息
解决泥沙问题是黄河治理的关键，而黄河适宜输沙量及未来输沙状态是黄河治理开发规划的重要基础。黄河年输沙量从1919—1959年的年均16亿t减少至近年来最小不足1亿t，这种现象或事实已经颠覆了人们对黄河泥沙现象的传统认识。在黄河输沙量减少原因分析的基础上，有关学者关注并预测了维持黄河健康的输沙量、未来黄河输沙的状况，但所得结果迥异。有的学者认为未来黄河的年输沙量可能在6亿t左右，也有学者认为在8亿t左右（修建三门峡水库时，预估全面实施水土保持后黄河输沙量减少一半），还有学者认为黄河的年输沙量未来将稳定在3亿t左右［12］。现实情况是，在降雨量没有显著减少的情况下，2000—2015年黄河潼关水文站年输沙量平均为2.56亿t，其中最大值出现在2003年（6.08亿t），其余各年均不超过5亿t，而最小仅为0.55亿t（2015年），如图4所示。
图4 黄河潼关站年输沙量变化情况
30 a前，几乎没有人预测到黄河输沙量会降低到目前的状态。黄河输沙量变化已经打破人们的传统思维，受气候、政策及下垫面条件变化不确定性的影响，黄河未来输沙量变化存在更大的不确定性，应该加强对未来治理和开发规划下、不同降雨频率下黄河输沙量状态的研究。
系统分析黄河输沙量变化原因可知，人类活动是黄河输沙量减少的主要原因，政策是降低土壤侵蚀、减少入黄泥沙的根本驱动力，各项水土保持措施是表现形式。在气候变化与人类活动共同作用下，近百年来黄河水沙呈现出了前所未有的剧烈变化。进入21世纪以来，黄河输沙量呈现出新的状态。黄河输沙量的大幅减少是黄河中游黄土高原各项水利水土保持措施综合作用的结果。针对生态环境脆弱的黄土高原地区，考虑到人类活动与自然和谐的特殊性，必须长期坚持保育水土，维系其生态和谐，才能维持该地区的自然生态环境稳定。实践经验证明，尊重科学、尊重自然、科学认识自然规律，是人类改造自然过程中必须遵循的真理。在黄土高原生态环境建设与黄河治理中，必须在科学研究的基础上做好相关的调研、规划，进行科技创新与制度创新，慎重对待入黄泥沙锐减引起的一系列新问题，保障区域水资源安全、防止洪涝灾害，确保为国家和区域社会经济发展与粮食安全提供支撑。
针对目前黄河输沙量大幅减少的现状，黄土高原治理措施格局应进行相应调整。在新的水沙条件下，黄土高原生态环境建设与黄河治理应注重以下两点：首先，黄河水沙异源、水少沙多、水沙关系不协调的问题仍未解决。在考虑黄河不同区间影响水沙关键要素的基础上，合理配置调控措施，有效控制水土流失，在尽可能减少入黄泥沙的同时较小影响径流量，缓解黄河干支流水资源短缺问题。其次，黄河流域是我国人口集中、经济社会快速发展的重要区域。在新的水沙条件下，要在继续大力开展流域水土保持和生态建设的基础上，提高自然资源利用率，建立区域经济可持续发展新模式，促进区域生态功能与社会经济协调可持续发展。
参考文献：
［1］ 穆兴民，张秀勤，高鹏，等.双累积曲线方法理论及在水文气象领域应用中应注意的问题［J］.水文，2010，30（4）：47-51.
［2］ 穆兴民，王万忠，高鹏，等.黄河泥沙变化研究现状与问题［J］.人民黄河，2014，36（12）：1-7.
［3］ 岳晓丽.黄河中游径流及输沙格局变化与影响因素研究［D］.杨凌：西北农林科技大学，2016：18-32.
［4］ 何毅.黄河河口镇至潼关区间降雨变化及其水沙效应［D］.杨凌：西北农林科技大学，2016：122.
［5］ 李二辉.黄河中游皇甫川水沙变化及其对气候和人类活动的响应［D］.杨凌：西北农林科技大学，2016：14-23.
［6］ 穆兴民.黄土高原水土保持对河川径流及土壤水文的影响［D］.杨凌：西北农林科技大学，2002：22-59.
［7］ 周佩华，王占礼.黄土高原土壤侵蚀暴雨标准［J］.水土保持通报，1987，7（1）：38-44.
［8］ 王万忠.黄土地区降雨特性与土壤流失关系的研究Ⅲ——关于侵蚀性降雨的标准问题［J］.水土保持通报，1984，4（2）：58-63.
［9］ 毕彩霞.黄河中游皇甫川流域产沙性降雨及其对径流输沙的影响［D］.杨凌：中国科学院教育部水土保持与生态环境研究中心，2013：29-34.
［10］ ZHAO G J，KLIK A，MU X M，etal.Sediment Yield Estimation in a Small Watershed on the Northern Loess Plateau，China［J］.Geomorphology，2015，241：343-352.
［11］ 穆兴民，李朋飞，高鹏，等.土壤侵蚀模型在黄土高原的应用述评［J］.人民黄河，2016，38（10）：100-110，114.
［12］ 胡春宏.黄河水沙变化与治理方略研究［J］.水力发电学报，2016，35（10）：1-11.
【责任编辑 张智民】
Key Issues and Reflections of Research on Sediment flux of the Yellow River
MU Xingmin1，2，HU Chunhong3，GAO Peng1，2，WANG Fei1，2，ZHAO Guangju1，2
（1.State Key Laboratory of Soil Erosion and Dryland Farming on the Loess Plateau，Northwest A＆F University，Yangling 712100，China；2.Institute of Soil and Water Conservation，CAS＆MWR，Yangling 712100，China；3.China Institute of Water Resources and Hydropower Research，Beijing 100038，China）
Abstract：It is unexpected that the Yellow River has experienced such a significant reduction in sediment load during the past decades.The causes and future changing trends need to be greatly investigated，as well as the strategy for adapting to the new condition of the Yellow River basin.According to our long-term study，this manuscript proposed some key issues associated with the sediment load changes of the Yellow River.Firstly，itis better to identify the spatial-temporalvariation ofsedimentload and their causes in a general way atlarge scale.Secondly，for the long-term investigation，the transition year is 1979，which is widely used to divide it into reference period（with relative limited human activities influence）and changing period（affected by intensive human activities）.Thirdly，a daily rainfall threshold of 12 mm is proposed to studies the soil erosion and sediment yield atcatchmentscale.In addition，itis necessary to clarify whether the reduction ofthe sediment load is due to decreased soil erosion on the hillslopes or trapping effects of dams and reservoirs on the Loess Plateau.It is also emphasized that further study are needed to obtain real sediment yields through dams′sedimentation，to make reasonable strategies for soil erosion control and to predict future changing trends of sediment load in the Yellow River basin.
Key words：sediment load；precipitation；human activities；Loess plateau；soil erosion；Yellow River
中图分类号：TV141；TV882.1
文献标志码：：A
doi：10.3969／j.issn.1000-1379.2017.08.001
收稿日期：2017-05-31
基金项目：国家重点研发计划项目（2016YFC0402401，2016YFC0501707）。
作者简介：穆兴民（1961—），男，陕西华阴人，研究员，博士生导师，长期从事水土保持、生态水文研究工作。
E-mail：xmmu＠ms.iswc.ac.cn
通信作者：胡春宏（1962—），男，浙江慈溪人，中国工程院院士，高级工程师（教授级），博士生导师，长期从事泥沙运动力学、河床演变与河道整治等领域的理论与应用研究工作。
E-mail：huch＠iwhr.com
联系客服
