内蒙古草原流域降水-径流模拟试验及水文演变驱动因素辨析

2021-08-07张阿龙于庆峰房丽晶高瑞忠

中国农村水利水电 2021年7期

张阿龙，于庆峰，2，于婵，房丽晶，高瑞忠

（1.内蒙古农业大学水利与土木建筑工程学院，呼和浩特010018；2.内蒙古农业大学职业技术学院，内蒙古包头014109；3.内蒙古自治区水文总局，呼和浩特010020）

近年来，气候波动和人类活动导致极端降水事件频繁、径流锐减［1-3］，流域产汇流过程发生变化［4-6］，进而草原流域草地退化、土地沙化和沙尘暴等生态环境问题日益突出。国内外学者对于流域产汇流过程从物理模拟、数学及数值模拟进行研究分析，管晓祥［7］对VIC 模型、新安江模型、WBM 模型和GR4J 模型在黄河典型子流域径流过程模拟效果进行对比；Hamed［8］和Martinez［9］发现，在相同的降水历时内，径流量和产沙量对降水强度存在正向响应；姚冲［10］通过模拟降水对南方典型黏土坡面进行研究，查明降水量、降水强度和坡度对坡面产流产沙量和土壤侵蚀过程的影响，并建立坡面土壤侵蚀模型；SWAT模型在域水文研究中应用广泛，Aronld［11］、Van Liew［12］运用SWAT 模型分别对美国不同流域径流和小沃希托河流域进行水文模拟，均取得了满意的模拟效果，Demirel［13］同时建立人工神经网络模型和SWAT模型，得出SWAT 模型显示的均方误差值更小，模拟效果更好。

综上所述，人工物理模拟降水-径流试验对地表产汇流模拟具有科学性，SWAT 模型对水文模拟效果具有优越性，因此，本文将人工物理模拟降水-径流试验和SWAT 模型共同引入径流模拟及变化归因分析中，研究草原流域产汇流特征，定量分析驱动径流变化的气候波动与人类活动的贡献，旨在探讨环境变化对草原天然径流演变的影响，进而为草原流域植被建设、生态环境保护和生态水文研究提供参考。

1 材料与方法

1.1 研究区概况

巴拉格尔流域（117°32'～118°30'E，43°26'～44°39'N）位于内蒙古锡林郭勒盟（图1），是以针茅和羊草为主的典型草原内陆河流域，其流域面积2 703 km2，冬季寒冷、夏季干热，平均气温2.39 ℃、年均降水量330.3 mm；流域主要经济产业为畜牧业和工业［2，5］。

1.2 数据来源与处理

土壤植被数据通过野外调查获取。结合流域多年降水资料、植被特征、地形地貌和草地管理方式等设计试验方案，共布设取样点15个（图1），分别测得0～30 cm 土壤干容重、饱和渗透系数、饱和含水率、毛管上升含水率和田间持水率；选取50 cm×50 cm 的植被样方，将鲜草剪去装入保鲜袋，并在该区挖取深30 cm 土层，筛洗出根系，装入保鲜袋，在实验室进行烘干称重测取地表和地下植被干物质量。

图1 研究区位置及土壤植被采样点分布Fig.1 Location of balaguer river basin and distribution of sampling points

降水-径流模型由供水、降水和径流装置3 部分构成（图2），供水装置由水箱，水泵，供水槽组成，由雨强调节器控制水量大小，人工模拟雨强均匀的降水，降水器支架为角钢焊接成的矩形框架，降水量大小利用雨量筒测定，经过多次重复试验后在雨强调节器上标定。降水-径流模拟试验方案设计为在典型区随机采集100（长）cm×50（宽）cm×30（深）cm 的低扰动土样块，进行不同降水强度下草原产汇流及水土侵蚀过程的试验模拟，分别设置不同降水强度、持续时间和土地坡度等多组合模拟场景，即降水强度为60、90 和120 mm/h，降水持续时间为1、2、3、4、5、6、7、8、9、10、11、12、13、14、15 min，坡度为10°、16°、28°等。

图2 降水-径流模拟装置示意图（单位：cm）Fig.2 Schematic diagram of simulated precipitation-runoff device

研究区气候数据包括降水、气温、风速、辐射和相对湿度等，来源于研究区内气象站实时监测和中国气象数据网（http：//data.cma.cn），DEM 数据来源于地理空间数据云（http：//www.gscloud.cn），土壤数据来源于寒区旱区科学数据中心（http：//wastdc.westgis.ac.cn）（图3），土地利用数据来源于美国地质勘查局（United States Geological Survey，USGS）数据共享网站（http：//glovis.usgs.gov/）提供的Landsat TM/OLI 影像（时间分辨率16 d，空间分辨率为30 m）解译获取（图4）。

图3 土壤类型分布Fig.3 Distribution of the soil type

图4 土地利用和植被覆盖分布Fig.4 Distribution of the land use and vegetation cover

统计分析采用Excel@2010软件完成，土壤-植被-根系的空间分布图通过ArcGIS@10.3软件制作。

1.3 研究方法

Mann-Kendall（M-K）检验是一种广泛用于时间序列趋势和突变分析，该法计算简便，可以识别突变开始时间及指出突变区域［14］。

克里金插值法是常用的地质统计格网化方法，其考虑了变量空间相关特性，广泛应用于地质、水文、气象等领域的空间分析［15］。

利用降水-径流模拟试验建立不同坡度降水强度与产流系数的数理回归模型，由降水强度经验频率与降水强度、不同坡度空间分布面积占比的变化范围分别划分区间，计算各区间经验频率，公式如下：

式中：Ri为坡度在i区间的径流系数；Qi为坡度在i区间的面积占比；Pj为降水强度在j区间的经验频率；rij为坡度在i区间，经验频率在j区间的径流系数；i为坡度，0°～10°、10°～16°、16°～28°；j为降水强度，0～10、10～60、60～90、90～120 mm/h；R为模拟径流系数。

SWAT 模型进行水文演变驱动因素分析，以Nash-Suttcliffe效率系数（Ens）和决定系数（R2）进行模型参数率定和模型验证评价及检验模拟拟合过程是否合理［17］。

以流域径流突变前的多年平均径流量为天然时期流域径流的基准值，则天然径流量与径流突变后期的多年实测平均径流量之间的差值就是受气候变化与人类活动的影响部分，因此，气候变化与人类活动的影响分割公式如下［18］：

式中：QH为突变后的实测径流；QN为突变前径实测流均值；QA为突变前实测径流均值；QB为突变前模拟径流均值；QC为突变后模拟径流均值；QCI为修正突变后模拟径流均值；I为修正系数；ΔQC为气候变化引起的均值变化度为；ΔQH为人类活动引起的均值变化度；WC为气候变化引起的贡献率；WH为人类活动引起的贡献率。

2 结果与分析

2.1 气象水文突变检验

通过M-K 检验法得到气象水文要素的突变年份（表1），突变主要集中在20 世纪90年代到20 世纪末，1982年最低气温突变最早，2001年降水量突变最晚，1988年径流深突变，突变前后变化率为-35.97%，而径流系数突变年份为1994年，突变前后变化率为-38.57%，径流深与径流系数突变前后变化率较为接近。

表1 气象水文要素突变前后变化Tab.1 Distinction of the hydrometeorological factors before and after the shift year

2.2 土壤物理及水力特性

土壤物理参数（表2）中饱和渗透系数介于0.19～14.23 m/d，但均值仅为2.87 m/d，表明土壤物理性质空间异质性显著，干容重最大值、最小值与均值接近，毛管上升含水率极差最大，饱和含水率次之。

表2 土壤物理特征统计Tab.2 Statistical parameters of the soil physical characteristics

2.3 土壤类型及植被分布

流域下垫面特征与土壤质地、植被分布和植被根系密切相关。研究区土壤类型、草干重、草根干重空间分布特征见图5，土壤以砂质壤土为主，粉砂质壤土相间的多砾沙质壤土和少砾砂质壤土呈岛状分布的分布特征，草干重最大值80～90 g 主要分布在上游山谷区，分布面积较小，呈现出中部最少，由南向北，由东向西递减的趋势，草根重以200～300 g 分布区域为主，最大值介于400～500 g，分布于上游东北部，100～200 g呈星点状分布，草根重于草种分布存在较大差异。

图5 土壤类型、草干重、草根干重空间分布图Fig.5 Spatial Distribution of soil type，grass dry weight and grass root dry weight

2.4 降水-径流模拟与分析

研究区多年降水频率分布曲线表明流域降水强度主要集中在0～10 mm/h，约占80%，地面坡度空间分布图（图6）中流域坡度主要集中在0～10°，约占流域面积的86%。

图6 流域坡度空间分布图Fig.6 Distribution of the basin slope

根据降水经验频率特征和不同坡度空间分布面积占比，分别将降水-径流模拟试验分为4组和3组（表3）进行。

表3 流域坡度面积与降水强度经验频率占比Tab.3 Proportion of different slope area and empirical frequency of precipitation Intensity

降水径流模拟试验中每场降水前土壤平均含水率保持在10%～15%之间，产汇流特征参数见表4，随着降水强度和坡度的增加，径流总量和径流系数均增大，地面坡度为10°时，降水强度为120 mm/h 的产流总量是降水强度为90 mm/h 的2.19 倍，是降水强度为60 mm/h 的4.69 倍；而地面坡度为16°时，降水强度为120 mm/h 的产流总量是90 mm/h 的2.07 倍，分别是60 mm/h 的4.71 倍；而坡度为28°时，降水强度120 mm/h 的产流总量是90 mm/h的1.92倍，分别是60 mm/h的4.66倍，综上，总产流量不仅与降水强度有关，而且受坡度影响较大。当坡度为10°时，3种降水强度产流机制主要以蓄满产流为主，随着坡度的增加，产流量不断加大，产流机制由蓄满产流为主逐渐向超渗产流过度，径流系数加大。

表4 不同坡度和降水强度下产流特征参数Tab.4 Statistical parameters of runoff characteristics

根据模拟降水试验结果，建立不同坡度降水强度与产流系数的多项式回归方程（图7），各回归方程R2值均大于99%，拟合效果较好。

图7 不同坡度、降水强度与径流系数关系Fig.7 Relationship between slope，rainfall intensity and runoff coefficient

不同区间降水强度值以中值代替，不同区间坡度值以最大值代替，计算不同坡度与降水强度对应的径流系数，结果如表5。

表5 不同坡度和降水强度下的径流系数Tab.5 Runoff coefficient of the different slope and rainfall intensity

根据公式（1）、（2）计算得到天然径流系数R值为0.424（表6），可看出该区坡度为16°～28°区间径流系数相对较大，而坡度为0°～10°区间径流系数较小，表明流域坡度对径流系数影响较大。

表6 流域天然径流系数计算Tab.6 Calculation of the basin runoff coefficient

2.5 流域径流变化驱动因素

基于人工降水模拟试验建立降水强度-坡度-径流系数的数学模型结果，可得突变前人类活动对巴拉格尔流域径流系数变化的贡献率为24%，气候变化贡献率为76%，突变后人类活动对巴拉格尔流域径流系数变化的贡献率为53%，而气候变化贡献率为47%（表7）。

表7 流域径流系数变化归因分析Tab.7 Attribution analysis of runoff coefficient

基于水文气象要素实测数据，利用SWAT 模型研究流域径流变化驱动因素，率定期模拟值和实测值的纳什系数Ens为0.88，决定系数R2为0.87，验证期模拟值和实测值的Ens为0.78，R2为0.81，由SWAT 的评价标准可知，SWAT 模型在该流域模拟效果较好，验证期的模拟值总体略大于实测值，但径流极大值实测值均大于模拟值（图8）。

图8 流域年径流模拟值与实际值对比Fig.8 Comparison between the observed and simulated runoff

由公式（3）～（9）分析流域气候和人类活动驱动影响径流的贡献，气候变化贡献率为79%，人类活动贡献率为21%，径流变化主要受气候变化影响，而受到的人类活动影响较小（表8）。

表8 流域径流变化归因分析Tab.8 Attribution analysis of basin runoff

3 讨论

巴拉格尔流域水文气象要素突变主要发生在20 世纪90年代到20 世纪末，突变前后水文气象要素差异较大，主要表现为气温升高，降水量、径流量减少，生态环境恶化加剧，导致生态环境与气象水文要素互馈陷入恶性循环，人类活动（对地下水吸夺，过度放牧，工业兴起，城市扩张等）加重了恶化趋势，造成巴拉格尔河丰水期断流、草场退化等生态水文现象凸显。Zhang等［2］通过对巴拉格尔流域气象水文要素突变时序研究，发现平均最低气温、平均气温、平均最高气温、径流量、降水对环境敏感性由强到弱，由气温的突然升高导致了巴拉格尔河流域春季径流量的急剧下降，因春季径流占全年的40%以上，春季径流的突变减少和温度的突变增加，导致降水的突然减少，突变时序与白勇［19］、高瑞忠等［20］对巴拉格尔流域水文气象突变年份一致；土壤类型和植被的空间分布对产汇流过程影响较大，不同的土壤类型渗透性差异较大，不同的植被对降水的截留和产汇流速率的影响较大，由于研究区土壤类型的质地和土壤水力参数相近，植被分布主要以草地为主，其他地物类型星点分布，虽然对径流系数有一定的影响，但相对流域整体而言可忽略不计。鲁克新等［21］对黄土区降水产流产沙过程研究认为裸坡坡面下，认为雨强是影响径流的主要因素，坡度是影响产沙的主要因素，忽略了模拟降水过程中蒸散发作用的影响和流域下垫面的空间异质性等因素的干扰，径流系数可以间接反映流域内自然地理要素对降水-径流关系的影响，由于人工模拟降水-径流关系可排除人类活动造成的干扰，且各要素均可人为控制，根据研究区地形地貌多年水文气象资料计算出各要素的分布频率，再通过人工降水模拟得理想条件较接近天然状态的径流系数，虽然该方法对天然径流系数的刻画存在一定的误差，但从模拟结果看，通过人工降水模拟试验建立降水强度-坡度-径流系数的数学模型得出天然径流系数值为0.424，模拟值与突变前的径流系数接近，较突变前多年平均实测径流系数大，因巴拉格尔河流域产流机制主要为超渗产流，多年平均径流系数的改变主要受下垫面条件控制，下垫面的变化主要受人类活动的影响；试验结果表明突变前人类活动对巴拉格尔流域径流系数变化的贡献率为20%左右，突变后为50%左右。利用分布式水文模型SWAT对径流变化归因分析发现巴拉格尔流域上游气候变化贡献率为79%，人类活动贡献率为21%，径流变化主要受气候变化影响，而受到的人类活动影响较小，与人工降水模拟试验建立降水强度-坡度-径流系数的数学模型得出的贡献率存在较大差异，由于SWAT 模型假设突变前的多年平均径流量为天然时期流域径流的基准值，缩小了人类活动贡献率，进而导致突变后人类活动对径流的影响的贡献率偏小，但SWAT 模型模拟结果与苏辉东等［22］对黄河、黑河、长江、雅鲁藏布江流域四个寒区流域得出气候变化对径流的影响贡献率达到78%以上是径流演变的主导作用的结论相近，而与高瑞忠［20］对巴拉格尔流域和王威娜［23］、于婵［24］等对锡林河流域的研究相反，可能是使用的水文数据序列和径流归因分析的方法不同造成。

变化环境下径流量的减少逐年加剧，由于突变后经济飞速发展，人类活动范围的不断扩大（过度放牧对草场的破坏、地下水开采对径流量的吸夺、矿产资源开发对地下水的输干和对下垫面条件的破坏等）和全球持续增温等因素的影响，导致草原生态系统退化，地下水位递减，草场面积萎缩等环境问题突出。气温突变增高导致径流减少，径流突变减少和气温的突变增高共同导致了降水的突变减少，促使巴拉格尔流域水文气象极端事件增多，干旱加剧，水文气象要素与生态环境互馈机制改变，但气候变化也受人类活动的干扰。张树磊等［25］认为在气候较为干燥的地区，径流对气候和下垫面变化都更为敏感，且区域差异性明显；刘昌明等［26］和张晓明等［27］认为森林植被增加不仅能显著减少丰水、平水及枯水期年的径流，也改变了地下水与地表水的分配。

学界对高原内陆河草原流域降水强度-坡度-径流系数之间关系的深入研究相对较少，需进一步探讨。对不同降水模式下流域产汇流试验研究中由于缺乏不同植被盖度的对比研究，降水-径流模拟设备还需进一步改进完善，需要加入测量土壤实时含水率、径流无间断测定等设备，考虑壤中流与地下径流的水文过程，以及各径流中的化学成分变化等，径流系数模型建立未考虑蒸散发，壤中流和地下径流等，由于降水强度取值较少，对降水强度与径流系数的拟合曲线可能存在一定的影响，对径流系数数学模型还需进一步研究，对径流归因分析方法间的差异仍需要进一步探究。