植被恢复对黄土丘陵沟壑区产流模式的影响研究.pdf

收稿日期：基金项目：国家自然科学基金资助项目（）；国家自然科学基金联合基金项目（）作者简介：穆兴民（），男，陕西华阴人，研究员，博士生导师，主要从事水土保持、生态水文、河流水沙研究工作：【“黄河流域生态保护和高质量发展”专题】植被恢复对黄土丘陵沟壑区产流模式的影响研究穆兴民，刘卓昕，王双银，温永福，高 鹏，赵广举，孙文义，马雪燕（中国科学院 水利部 水土保持研究所，陕西 杨凌；西北农林科技大学，陕西 杨凌；吕梁市水利局，山西 吕梁；南京水利科学研究院 农村水利研究所，江苏 南京；新疆水利水电科学研究院，新疆 乌鲁木齐）摘 要：探究植被恢复对产流模式的影响是认识变化环境下流域水文及水土流失过程的基础。通过分层原位双环土壤入渗试验和植被恢复前后流域实测场次洪水过程解析，研究植被恢复对黄土丘陵沟壑区典型小流域产流模式的影响。研究结果表明：植被恢复显著增强了土壤入渗性能，不同深度土层的入渗速率均有不同程度的增大。植被恢复增大了表层土和 土层入渗性能的差异。土壤入渗性能随土壤密度减小和孔隙度增大而显著增强，植被恢复通过降低土壤密度和增大孔隙度影响土壤入渗和产流过程。植被恢复后小流域场次洪水的洪峰流量、径流深、径流系数以及地表径流占比均减小，壤中流增加，产流模式从超渗产流向蓄满产流转变，形成黄土高原特有的浅层蓄满产流模式。植被恢复影响了流域的产流模式，为深入研究黄土高原土壤水文过程提供了科学依据。关键词：产流模式；土壤物理特性；土壤入渗；植被恢复；黄土高原中图分类号：；文献标志码：引用格式：穆兴民，刘卓昕，王双银，等植被恢复对黄土丘陵沟壑区产流模式的影响研究人民黄河，（）：，（，；，；，；，；，）：），（），：；长期的水土流失治理特别是 年退耕还林（草）工程实施以来，黄土高原植被的生产能力及覆盖度大幅提升、结构和功能明显改善。植被恢复对降低坡面土壤侵蚀、减少入河水沙具有不可替代的重要作用，黄土高原生态恢复也改变了河流水文情势。一般认为黄土区包气带厚，土壤入渗能力差，雨强大，植被覆盖度低，产流模式以超渗产流为主。随着黄土高原土地利用调整和植被结构及功能的恢复，流域产流模式也在发生改变。年 月延安市发生长历时强降雨，浅层土壤含水量出现大于田间持水量的现象，部分区域出现饱和地表径流，延河第 卷第 期 人 民 黄 河 ，年 月 ，流域并未发生较大洪水，表明通常所认为的超渗产流并没有对流域产流过程产生主导作用。等采用水文模型模拟了渭河流域产流模式，发现渭河上游的主要产流模式可能从超渗产流转变为蓄满产流。李楠选择具有代表性的黄土高原上的汾河上游静乐站控制流域和祖厉河流域，基于次降雨洪水特征辨识了流域产流模式变化，发现汾河上游静乐站控制流域植被恢复后流域超渗产流次数占比减小，混合和蓄满产流的次数占比增加，而下垫面变化未导致祖厉河流域产汇流机制发生较大变化。胡彩虹等对佳芦河流域的洪水过程数据进行分析，发现随植被覆盖度提高，年均流量大幅度减小，以超渗产流为主的洪水场次减少，以混合及蓄满产流为主的洪水场次增加，该流域发生蓄满产流的概率增大。顾朝军和穆兴民等研究指出，植被恢复显著改善了土壤水文物理性质，提高了土壤入渗能力，进而改变了产流模式，延河一级支流西川河的产流模式偏向混合产流和蓄满产流，流域在部分降雨过程中发生蓄满产流情况。洪艳艳利用 水文模型从物理层面探讨了植被恢复对流域产流模式的影响，发现近地表层的土壤饱和导水率在植被恢复后显著提高，土层间导水率差异性增强，这使得土壤形成上层滞水，从而导致黄土高原饱和地表径流产流量增加。黄土高原的植被恢复为什么能以及如何改变流域产流模式？场次洪水径流组成成分的划分从水文学上说明了产流模式是否改变这一问题，但不足以回答为什么会改变的问题，也就难以从机理层面揭示植被恢复对流域产流模式的影响，更难以预测未来产流模式如何变化。影响产流的主要因素有气候因素（降雨特征、蒸散发特性）和下垫面因素（地形地貌、植被状况、土壤特性），而对于特定区域而言，植被状况是影响产流机制的主要因素。随植被覆盖度及生产力的提升，植被会通过冠层、近地表枯枝落叶层降低雨滴动能、增加地表糙率，通过地下根系等增强土壤孔隙度、提高土壤水分入渗能力、影响坡面的产（汇）流机理，进而影响流域的产流模式。目前研究多从场次洪水过程线研判产流模式的变化，缺乏产流模式与土壤性质、土壤入渗之间关系的研究，难以揭示黄土高原产流模式变化的内因。在百年尺度内，流域基本地形地貌及降雨特征并未发生本质性变化。因此，本研究采用传统水文方法揭示受植被恢复影响的流域产流模式变化原因，同时分析土壤性质、土壤入渗变化与产流之间的关系，以期为科学认识黄土高原植被恢复驱动下产流模式的变化及其水文模型构建提供理论依据。研究区与研究方法 研究区概况试验研究在陕西省绥德县的桥沟流域进行。桥沟位于黄土高原黄土丘陵沟壑区第一副区，是裴家峁流域下游右岸的一级支沟，流域面积 ，主沟道长，流域平均宽 ，沟壑密度 ，流域内右岸有一支沟、二支沟两条较大支沟。桥沟一支沟是桥沟流域下游右岸的一个支沟（靠近桥沟流域出口），流域面积 ，主沟道长 ，沟道比降，流域内植被类型以草地为主。桥沟二支沟是桥沟流域中游右岸的一个支沟（远离桥沟流域出口），流域面积 ，主沟道长 ，沟道比降，流域内植被类型以草地为主（有少许柠条和刺槐）。土壤类型主要为黄绵土。气候类型为半干旱大陆性季风气候，多年平均气温约为 ，年平均降水量约为 ，主要集中在 月，且多以暴雨形式出现。与黄河水利委员会绥德水土保持科学试验站资深研究者座谈得知，年以前桥沟流域内大部分土地为坡耕地，尤其是流域梁峁坡和沟底约八成坡耕农田，沟坡及沟道多为荒地，亦因长期过度放牧，流域植被覆盖度不足。年以来，随着国家退耕还林（草）工程的实施，梁峁坡和部分沟床地种植乔灌木林和人工草，植被覆盖度显著增大，目前流域植被类型以草本植物为主，主要有艾蒿、狗尾草、胡枝子、猪毛蒿等草本植物，没有大片林地，植被覆盖度达。研究方法 土壤分层入渗试验试验于 年 月进行。在桥沟流域的一支沟和二支沟中选择植被生长典型且较为平坦的样地进行不同深度土层的双环入渗试验，土层深度分别为（表层土）、，每个样地均进行 次重复试验。入渗试验开始前先去除地表植被和枯落物，并将地表铲平，用铁锤将内环直径 、外环直径 的双环垂直打入土壤表面以下 。试验开始后，同时向内外环迅速注水，当内外环水深达 时开始计时，采用定水头注水法测定入渗速率。每次双环入渗试验进行 ，以、为计数时刻，将前 的平均入渗速率视为初始入渗速率，最后 的平均入渗速率视为稳定入渗速率，的平均入渗速率为平均入渗速率。同时在入渗样地周围使用 环刀采集原状土样用于测定土壤密度与孔隙度，分析植被恢复后不同土层土壤密度和孔隙度变化对土壤入渗性能的影响。场次降水和洪水数据来源选取桥沟流域雨量站和水文站的场次降水和洪水人 民 黄 河 年第 期观测数据，数据序列为 年（植被恢复前）以及 年（植被恢复后）。监测数据已经过多次人工校核，确保数据准确、可靠。产（汇）流特征值的确定水文观测站径流成分可以分为 个主要部分：地表径流、壤中流和地下水。本研究采用退水曲线法对径流进行分割，采用蓄泄关系法计算总径流深。洪水历时、洪峰流量、洪峰滞时等汇流特征值均可由实测径流过程推求（见图）。洪水历时是指洪水起涨点至消退点的时长；洪峰流量是指一次洪水事件中流域内达到的最大流量；洪峰滞时是指净雨质心到洪峰出现的时间。图 流域产（汇）流特征值示意 结果与分析 植被恢复对不同深度土层土壤入渗能力的影响桥沟一支沟、二支沟和农地土壤入渗规律总体相似，均表现为土壤入渗速率在前 内迅速下降，在 后基本趋于稳定，且随着土层深度增加，入渗速率降低（见图）。植被恢复对土壤的入渗能力有显著影响（），可以明显看出，在整个入渗过程中，农地不同深度土层的入渗速率整体上小于一支沟和二支沟，且随着土层的加深，农地和植被恢复地入渗性能的差距逐渐增大（见图）。植被恢复后，一支沟和二支沟不同深度土层的土壤入渗能力均有不同程度的提升，且平均入渗速率均显著大于农地（）（见图，图中不同大写字母表示同一植被类型下不同土层的土壤特征指标差异显著，不同小写字母表示相同深度下不同植被类型差异显著）。相较于农地，一支沟 各土层的初始入渗速率依次增大、，稳定入渗速率依次增大、；二支沟 各土层的初始入渗速率依次增大、，稳定入渗速率依次增大、。图 不同深度土层土壤入渗过程图 不同深度土层土壤入渗特征值人 民 黄 河 年第 期 随着土层的加深，所有样地的土壤入渗能力均显著下降（），各项土壤入渗特征值均随着土层的加深显著降低。由入渗过程及特征值可知，一支沟和二支沟的土壤分层入渗速率差异明显，尤其是表层土与 土层之间差异性显著（）。桥沟一支沟、二支沟和农地 土层的初始入渗速率相较于表层土分别降低、和，稳定入渗速率分别降低、和。不同深度土层土壤入渗影响因素分析土壤密度和孔隙结构是影响土壤入渗能力的主导因素。由各样地不同深度土层的土壤密度和孔隙度（见表）可知，土层深度对土壤密度和孔隙度有显著影响（），随着土层的加深，土壤密度增大，孔隙度降低。相较于农地，受植被恢复影响的草地土壤孔隙结构差异显著（），尤其在 以上一支沟和二支沟的土壤孔隙度显著高于农地（）。表 各样地不同深度土层土壤密度和孔隙度土层深度样地密度（）总孔隙度毛管孔隙度 非毛管孔隙度一支沟（）（）（）（）二支沟（）（）（）（）农地（）（）（）（）一支沟（）（）（）（）二支沟（）（）（）（）农地（）（）（）（）一支沟（）（）（）（）二支沟（）（）（）（）农地（）（）（）（）一支沟（）（）（）（）二支沟（）（）（）（）农地（）（）（）（）一支沟（）（）（）（）二支沟（）（）（）（）农地（）（）（）（）注：“（）”外为要素平均值，“（）”内为要素均值的标准偏差，其值越大代表离散程度越高。不同大写字母表示同一植被类型下不同土层的土壤特征指标差异显著（）；不同小写字母表示相同深度下不同植被类型差异显著（）。为进一步明确植被恢复后土壤密度和孔隙结构对土壤入渗能力的影响程度，对土壤入渗能力与土壤密度和孔隙结构进行了相关性分析和回归分析。土壤密度和各类孔隙度对土壤入渗性能有显著影响（见表）。初始入渗速率、稳定入渗速率和平均入渗速率与土壤密度均显著负相关（），与总孔隙度、毛管孔隙度和非毛管孔隙度显著正相关（）。土壤入渗特征值与密度、总孔隙度和毛管孔隙度的拟合程度均较高，而与非毛管孔隙度拟合程度较低（见图、图）。其中，密度、总孔隙度和毛管孔隙度对初始入渗速率解释程度分别为、和（分别为、和）；对稳定入渗速率解释程度分别为、和（分别为、和）；对平均入渗速率解释程度分别为、和（分别为、和）。而非毛管孔隙度对初始入渗速率、稳定入渗速率和平均入渗速率解释程度仅为、和（分别为、和）。表 土壤入渗参数与土壤特性指标的相关系数入渗参数相关系数密度总孔隙度毛管孔隙度 非毛管孔隙度初始入渗速率稳定入渗速率平均入渗速率注：表示土壤入渗参数与土壤特性指标相关性显著（）。图 初始入渗速率与土壤特征参数的关系图 稳定入渗和平均入渗速率与土壤特征参数的关系人 民 黄 河 年第 期 植被恢复下的流域产流模式桥沟流域植被恢复前、后 场洪水主要特征值见表、表。植被恢复前 场和植被恢复后 场降雨的平均降雨历时分别为 和 ，平均降雨量分别为 和 ，平均降雨强度分别为 和 ，场次平均最大降雨强度为 和 。植被恢复后平均降雨历时、降雨量、降雨强度和最大降雨强度分别增大、和，平均降雨量显著增大。洪水历时平均值分别为 和 ，洪峰滞时均值分别为 和 ，多年平均洪峰流量为 和 ，平均径流深分别为 和 ，径流均以地表径流为主，平均占比分别为和，植被恢复后洪水历时和洪峰滞时分别增加 和，多年平均洪峰流量和径流深分别减小 和，洪峰滞时显著增大。综上，植被恢复后多年产洪降雨的平均降雨历时、平均降雨量、平均降雨强度和最大降雨强度均大于植被恢复前，但是洪峰流量、径流深、径流系数以及地表径流占比均减小，地表径流占比的下降说明降水入渗量增大，基流量增加。表 桥沟流域 年洪水主要特征值统计序号洪号降雨历时 雨量 平均雨强（）最大雨强（）洪水历时 洪峰流量（）洪峰滞时 峰型径流深地表径流占比 径流系数 双峰单峰双峰单峰单峰双峰双峰单峰单峰单峰双峰单峰双峰单峰单峰多峰单峰表 桥沟流域 年洪水主要特征值统计序号洪号降雨历时 雨量 平均雨强（）最大雨强（）洪水历时 洪峰流量（）洪峰滞时 峰型径流深地表径流占比 径流系数 双峰单峰单峰单峰单峰多峰单峰多峰双峰单峰单峰 植被恢复后，陡涨缓落型的洪水场次占比显著增大，陡涨陡落型占比显著减小，两个时期均无缓涨缓落型洪水（见表）。陡涨陡落型洪水比例减小，而陡涨缓落型洪水比例增大，表明流域产汇流速率和退水速率减小，流域对径流的调蓄作用增强。表 桥沟流域 场洪水不同涨落类型场次统计时间陡涨陡落陡涨缓落缓涨缓落合计 年 年合计人 民 黄 河 年第 期 为了进一步说明流域产流模式的改变是由植被恢复引起的，从 场洪水中筛选出两组植被恢复前后的场次洪水进行对比分析。同一组中，植被恢复后产洪降雨的雨量和雨强大于植被恢复前，或与植被恢复前接近。若植被条件不变，按照传统水文学方法推求可知，植被恢复后场次洪水的洪峰流量、地表径流占比相较于植被恢复前场次洪水均应增大，但对比实测数据（见表、表）发现在植被恢复条件下组 洪水（和）的洪峰流量、径流深、地表径流占比和径流系数依次降低、和，组 （和）的 各 参 数 依 次 降 低、和，说明植被恢复改变了桥沟流域产流模式。对比场次洪水过程线（见图）发现，植被恢复前洪水过程线陡涨陡落，洪峰流量大，植被恢复后洪水过程线转变为陡涨缓落，洪峰流量显著减小，表明植被恢复影响了流域产流机制，改变了流域产流模式。桥沟流域无常流水，因此无地下径流补充。植被恢复后地表径流比例降低，洪水退水速率变小原因可能是土壤蓄积水分，对流域洪水径流起到一定调蓄作用。图 洪水过程线对比 讨论 植被恢复后流域下垫面变化已有研究表明，植被恢复过程中，土壤有机质的增加和地下生物活动的加剧导致土壤密度降低，孔隙度增大，土壤渗透性增强，深刻改变了流域下垫面条件。本研究发现，植被恢复可显著降低 土层土壤密度，改变土壤孔隙结构，提高土壤入渗能力。顾朝军虽也取得相似规律，但探究的土层深度较浅（），缺乏植被恢复对深层土壤入渗能力的影响。本研究还发现，植被恢复后，不同深度土层的入渗能力均有不同程度的增大，由于浅层土壤中具有发达的植物根系，因此植被恢复导致表层土和 土层入渗能力差异增大，这使得降雨时入渗量变大，而深层土壤入渗能力较差，容易产生“表层蓄满产流”。土壤密度和孔隙特征是影响土壤入渗能力的重要因素。本研究发现土壤密度与初始入渗率、稳定入渗率和平均入渗率显著负相关（），土壤总孔隙度和通气孔隙度与初始入渗率、稳定入渗率和平均入渗率显著正相关（），进一步可以得出植被恢复改善了土壤物理性质，进而提高了土壤入渗能力。同时，植物根系和土壤生物在土壤中产生的大孔隙使得土壤水入渗时出现优先流，这也会导致土壤入渗速率变大。下垫面变化对流域产流模式的影响影响产流模式的因素可概括为气候因素和下垫面因素两类。气候因素主要指降雨、蒸发对产流的影响。植被恢复后多场产洪降水的平均降雨量和降雨强度均大于植被恢复前，且有研究表明 年桥沟降雨量整体呈缓慢增加趋势，因此可以排除降雨特征变化对流域产流模式造成的影响。下垫面因素主要包括地形地貌、地表覆盖、土壤特性，地形地貌和土壤特性在特定区域内随时间变化甚微，在本研究所涉及的时间及空间范围内可以忽略，因此可以认为以植被恢复为主的地表覆盖变化是引起产流模式变化的主要因素。桥沟流域自 年开始进行植被恢复后，植被覆盖度显著提高，土壤孔隙、密度及入渗性能变化受到植被恢复影响显著。本研究结果表明，植被对土壤的影响主要表现在较浅的地表及根系所达土层的土壤中，尤其对 土层影响显著，且植被恢复增强了表层土和 以下土层入渗性能的差异性，这种变化亦造成壤中流的产生，促使土壤蓄积更多水分。在自然界，由于包气带结构的复杂性和降雨情况的多变性，因此存在多种产流模式组合的产流类型。典型的黄土区，年均 的降雨量显然不能使传统的包气带土壤完全饱和。但在黄土的形成过程中，常常存在相对较浅的不透水层，随植被恢复及其对土壤孔隙和入渗性能的改善，上层土壤透水性好，降水后径流成分主要为饱和地面径流和壤中流，若降雨强度超过表层土壤的入渗速率，径流成分则多为超渗地面径流和壤中流，而相对不透水层之下土层的水分仍未人 民 黄 河 年第 期饱和，从而在黄土高原形成特有的浅层蓄满产流模式。本研究中的场次洪水特征以及产流过程分析表明，植被恢复后桥沟流域产流过程由陡涨陡落向陡涨缓落的方向转变，其原因是浅层土壤水文性质显著变化，包气带中浅层相对不透水层促使形成了壤中流，从而导致洪水退水速率变小，产流模式发生变化。黄土高原在植被恢复驱动下，下垫面中的土地利用及地表覆盖条件发生深刻变化，导致浅层土壤水文性能改变，从而引起流域产流模式的变化。顾朝军、李楠、胡彩虹等、马雪燕采用次降雨洪水特征辨识流域产流模式的方法，分别对黄土高原的西川河、祖厉河、佳芦河、杨家沟和董庄沟进行研究，均发现植被恢复以来，陡涨陡落型洪水比例有所减小，而陡涨缓落型洪水比例有所增大，指出黄土高原产流模式有从传统认为的超渗产流向蓄满产流转变的趋势。但这里所谓的蓄满产流模式并非包气带土层全部饱和蓄满的产流模式，实际上是浅层蓄满的产流模式。结论针对黄土高原植被恢复条件下土壤入渗变化导致的流域产流模式转变问题，以绥德桥沟小流域为研究区，分析了植被恢复对土壤入渗性能的影响，探讨了植被恢复对桥沟小流域产流模式的影响，得到主要结论如下。）植被恢复显著增强了土壤入渗性能（），不同深度土层的入渗速率均有不同程度的增大，增大幅度随土层深度增加而减小。其中，表层土和 土层入渗性能的差异最大。）密度、总孔隙度和通气孔隙度与土壤入渗性能显著相关（），其中，密度与土壤初始入渗率、稳定入渗率和平均入渗率负相关，总孔隙度和通气孔隙度与土壤初始入渗率、稳定入渗率和平均入渗率正相关。）相较于植被恢复前，植被恢复后多年场次洪水的洪峰流量、径流深、径流系数以及地表径流占比均减小，壤中流增加。植被恢复前后陡涨陡落型洪水和陡涨缓落型洪水占比分别为、和、，场次洪水过程线及特征值的变化表明流域的产流机制发生改变，形成黄土高原特有的浅层蓄满产流模式。参考文献：穆兴民，李朋飞，刘斌涛，等 年黄土高原土壤侵蚀格局演变及其驱动机制人民黄河，（）：，：，（）：，（）：穆兴民，贾昊，高鹏，等水土保持生态建设对黄河中游水文情势影响研究人民黄河，（）：李欢玮，赵广举，穆兴民，等渭河水文情势变化评估及归因分析水土保持研究，（）：，刘昌明，洪宝鑫，曾明煊，等黄土高原暴雨径流预报关系初步实验研究科学通报，（）：赵人俊流域水文模拟：新安江模型与陕北模型北京：水利水电出版社，：顾朝军，穆兴民，孙文义，等极端暴雨洪水及侵蚀产沙对延河流域植被恢复响应的比较研究自然资源学报，（）：，（）：李楠下垫面变化条件下的黄土高原产汇流机制辨析研究郑州：郑州大学，：胡彩虹，冉广，荐圣淇佳芦河流域林草覆被变化对产流机制的影响水土保持学报，（）：顾朝军黄土区土壤水文物理特性及流域产汇流机制变化对植被恢复的响应杨凌：西北农林科技大学，：穆兴民，顾朝军，孙文义，等植被恢复改变黄土高原产流模式问题初探人民黄河，（）：洪艳艳黄土高原植被恢复对地表产流机制的影响研究杭州：浙江大学，：董玉婷，穆兴民，王双银，等产流及其研究进展华北水利水电大学学报（自然科学版），（）：包为民水文预报 版北京：中国水利水电出版社，：，：王国梁，刘国彬，周生路黄土丘陵沟壑区小流域植被恢复对土壤稳定入渗的影响自然资源学报，（）：（下转第 页）人 民 黄 河 年第 期统中的水库调蓄能力及城市水资源协调能力。参考文献：王利娜，苏静，郑晓东，等我国农业干旱风险研究进展简述水科学与工程技术，（）：赵丽，冯宝平，张书花国内外干旱及干旱指标研究进展江苏农业科学，（）：金菊良，宋占智，崔毅，等旱灾风险评估与调控关键技术研究进展水利学报，（）：张继权，刘兴明，严登华综合灾害风险管理导论北京：北京大学出版社，：金菊良，马强，崔毅，等基于三元链式传递结构的区域旱灾实际风险综合防范机制分析灾害学，（）：章国材气象灾害风险评估与区划方法北京：气象出版社，：倪长健论自然灾害风险评估的途径灾害学，（）：金菊良，郦建强，周玉良，等旱灾风险评估的初步理论框架灾害学，（）：赵克勤集对分析及其初步应用大自然探索，（）：王雅燕，周玉良，周平，等基于集对分析的安徽省农业旱灾风险评估水电能源科学，（）：黄崇福，刘新立，周国贤，等以历史灾情资料为依据的农业自然灾害风险评估方法自然灾害学报，（）：刘亚彬，刘黎明，许迪，等基于信息扩散理论的中国粮食主产区水旱灾害风险评估农业工程学报，（）：，（）：，：，（）：徐昆，朱秀芳，刘莹，等采用 作物生长模型研究中国玉米干旱脆弱性农业工程学报，（）：金菊良，孔令茹，崔毅，等基于五元半偏减法集对势的区域农业旱灾脆弱性评价农业机械学报，（）：金菊良，沈时兴，郦建强，等基于联系数的区域水资源承载力评价与诊断分析方法华北水利水电大学学报（自然科学版），（）：董涛，王振龙，金菊良，等基于风险矩阵和五元减法集对势的区域旱灾风险链式传递诊断评估方法灾害学，（）：汪洁，金菊良，周戎星，等基于连续半偏减法集对势的区域旱灾风险评估方法灾害学，（）：陈守煜水资源与防洪系统可变模糊集理论与方法大连：大连理工大学出版社，：金菊良，吴开亚，魏一鸣基于联系数的流域水安全评价模型水利学报，（）：赵克勤集对分析及其初步应用杭州：浙江科技出版社，：金菊良，陈鹏飞，陈梦璐，等水文水资源学家陈守煜先生学术研究的知识图谱分析水利学报，（）：金菊良，陈鹏飞，张浩宇，等五元减法集对势及其在水资源承载力趋势分析中的应用华北水利水电大学学报（自然科学版），（）：，金菊良，赵行玥，崔毅，等半偏减法集对势在区域旱灾风险评估中的应用水利水运工程学报，（）：李健，夏文菊，赵剑锋济南市水资源开发利用存在问题与建议山东水利，（）：，桑国庆区域干旱风险管理研究济南：山东大学，：【责任编辑 简 群】（上接第 页），（）：，：马雪燕植被恢复对黄土高塬沟壑区产流模式的驱动效应杨凌：西北农林科技大学，：刘卓昕植被恢复对黄土丘陵沟壑区土壤水文性质及产流模式的影响研究 杨凌：西北农林科技大学，：，：刘姗姗，李小兵，刘思君，等黄土丘陵沟壑区典型小流域不同空间尺度水沙变化研究泥沙研究，（）：李蔚，冉启华沟谷横断面类型对产流过程的影响水利学报，（增刊）：芮孝芳，宫兴龙，张超，等流域产流分析及计算水力发电学报，（）：【责任编辑 张 帅】人 民 黄 河 年第 期