山坡地形曲率分布特征及其水文效应分析——真实流域的野外实验及相关分析研究

研究了山坡地形曲率的空间分布特征,发现流域略呈现出微凹、发散的形状特征,平均的剖面和水平曲率分别为-4.62×10-4和3.49×10-4。通过划分源头型和边坡型山坡,发现源头山坡多是收敛的,边坡型山坡则多为发散的;在山坡内部收敛、发散、凹和凸等地形地貌类型是同时存在的,即山坡一般由相互组合的各种坡型所组成。基于野外采样观测的成果,分析了和睦桥子流域90个采样点的土壤含水量及其对应曲率的相互关系。结果显示:一般来说,凹形山坡、收敛形山坡对应的土壤含水量更高,而凸形和发散形山坡对应的土壤含水量较低;不论是收敛的凹坡,还是发散的凸坡,其采样点土壤含水量都有随高程下降而升高的趋势,即愈靠近沟谷土壤含水量愈大。     

Quantifying the effects of subsurface heterogeneity on hillslope runoff using a stochastic approach

The role of heterogeneity and uncertainty in hydraulic conductivity on hillslope runoff production was evaluated using the fully integrated hydrologic model ParFlow. Simulations were generated using idealized high-resolution hillslopes configured both with a deep water table and a water table equal to the outlet to isolate surface and subsurface flow, respectively. Heterogeneous, correlated random fields were used to create spatial variability in the hydraulic conductivity. Ensembles, generated by multiple realizations of hydraulic conductivity, were used to evaluate how this uncertainty propagates to runoff. Ensemble averages were used to determine the effective runoff for a given hillslope as a function of rainfall rate and degree of subsurface heterogeneity. Cases where the water table is initialized at the outlet show runoff behavior with little sensitivity to variance in hydraulic conductivity. A technique is presented that explicitly interrogates individual realizations at every simulation timestep to partition overland and subsurface flow contributions. This hydrograph separation technique shows that the degree of heterogeneity can play a role in determining proportions of surface and subsurface flow, even when effective hillslope outflow is seen. This method is also used to evaluate current hydrograph separation techniques and demonstrates that recursive filters can accurately proportion overland and base-flow for certain cases.     

天然实验流域降雨径流现象发生的悖论

作为滁州水文实验系统一部分的1号天然实验流域,面积7897m2,以安山岩为基底,上覆平均厚度2.46m的第四系沉积物。实测了包括地面径流和地面下径流的各种降雨径流响应,后者包括来自非饱和带的壤中流和饱和带的地下水径流。表明这些径流成分有着复杂的组合类型,主要是以地面径流为主的SR型和以地面下径流为主的SSR型,以及中间的和演化的类型。SR型实例中的地面径流量可占总径流的65%,而SSR实例中的地面下径流量可有90%。主要降水的7月,地面下贡献占54.5%,其中地下水径流即占33%。大部分地面径流与降水的18O组成有不小的差异,比较了同时进行测验的3个实验流域,在1400min的降水径流过程中,降水的平均δ18O为-1.210%,而同期地面径流的平均δ18O,2号水文山流域(512m2)为-1.132%,1号南大洼天然流域为-1.065%,3号只有薄层风化碎屑的牵牛花流域(4573m2)为-0.801%。这质疑了现行同位素流量过程线划分方法8个假定中的两个:地面径流的同位素组成不同于降水,天然流域尤其如此;各种水源在汇集过程中的同位素分馏影响并非都可忽略不计。实验流域因降水而产生径流,但所产生的径流却有着非本次降水的组分,这一降水径流悖论发生于各个实验流域。SR型和SSR型总的非本次降水组分分别达16%和64%。     

山坡表层关键带结构与水文连通性研究进展

山丘区是洪水的"策源地",山丘区坡地、沟谷及间歇性河道为洪水的形成提供了通道,同时也是水文连通时空变化最为强烈的地带。然而,对流域表层关键带结构特征及其水文连通机制等的认识尚存不足,限制了产汇流理论及模型方法的发展和应用。通过对比国内外山坡水文实验,发现山坡物理结构连通性控制并深刻影响着水流的连通过程,现有水文连通实验侧重孔隙等微观尺度的规律研究,与水文模型理论存在尺度上的巨大偏差。为此,提出水文连通性应侧重揭示水流在山坡地表、地下的宏观表象通道及分布特征,探索径流连通的动力学机制,即山坡水文连通性研究重在剖析其结构特征的水文累积效应,应保持关键带结构特征合理概化与产汇流理论适度复杂之间的平衡。     

Modelling of hillslope runoff processes

 The present study is aimed at modelling hillslope flows with emphasis on subsurface stormflows that involve macropores. The physical processes connected with the runoff process on a hillslope are identified. The components which are considered in modelling the hillslope flow are the nature of flows in the macropore and micropore domains, the spatial and temporal characteristics of the macropore network, the interaction between the domains, and the initiation of flow in the macropores. Both Horton and Dunne's variable source area generation mechanisms are explicitly incorporated in the model. The dominant physical processes governing hillslope runoff are conceptualized in terms of parameters which are derived from the physical properties of the soil, the nature of macropores, and hillslope geometry. The conceptualization of the model is then used to examine infiltration and runoff production. This helps to compute the development of the groundwater table, runoff hydrograph, and soil moisture profile. Received: 5 October 1996 · Accepted: 25 June 1997     

Effect of using multiple stream gauging stations on calibration of hydrologic parameters and estimation of hydrograph of ungauged neighboring basin

In this study, Hydrologic Engineering Center-Hydrologic Modeling System is used to simulate hydrologic processes in a watershed in Western Black Sea Region that frequently experiences flooding. The region is mountainous with steep hill slopes and receives high precipitation throughout the year. There are three stream gauging stations in the basin whose data are available for calibration and validation of hydrologic parameters. Simulations are performed for different scenarios to investigate the effect of using multiple stream gauging stations’ data on catchment wide calibration and validation of hydrologic parameters. Furthermore, performance of using calibrated internal stream gauging stations’ flow data in the estimation of hydrologic parameters in an assumed neighboring ungauged basin was assessed. It is found that using data of multiple stream gauging stations for calibration and validation gives satisfactory results for direct runoff hydrograph but the peak discharge predictions are not improved. The study results suggest that using data of internal stream gauging stations enables improved understanding of internal dynamics and transport in the basin and better predicted direct runoff hydrograph for the assumed neighboring ungauged basin.     

Wildfire impacts on the processes that generate debris flows in burned watersheds

Every year, and in many countries worldwide, wildfires cause significant damage and economic losses due to both the direct effects of the fires and the subsequent accelerated runoff, erosion, and debris flow. Wildfires can have profound effects on the hydrologic response of watersheds by changing the infiltration characteristics and erodibility of the soil, which leads to decreased rainfall infiltration, significantly increased overland flow and runoff in channels, and movement of soil. Debris-flow activity is among the most destructive consequences of these changes, often causing extensive damage to human infrastructure. Data from the Mediterranean area and Western United States of America help identify the primary processes that result in debris flows in recently burned areas. Two primary processes for the initiation of fire-related debris flows have been so far identified: (1) runoff-dominated erosion by surface overland flow; and (2) infiltration-triggered failure and mobilization of a discrete landslide mass. The first process is frequently documented immediately post-fire and leads to the generation of debris flows through progressive bulking of storm runoff with sediment eroded from the hillslopes and channels. As sediment is incorporated into water, runoff can convert to debris flow. The conversion to debris flow may be observed at a position within a drainage network that appears to be controlled by threshold values of upslope contributing area and its gradient. At these locations, sufficient eroded material has been incorporated, relative to the volume of contributing surface runoff, to generate debris flows. Debris flows have also been generated from burned basins in response to increased runoff by water cascading over a steep, bedrock cliff, and incorporating material from readily erodible colluvium or channel bed. Post-fire debris flows have also been generated by infiltration-triggered landslide failures which then mobilize into debris flows. However, only 12% of documented cases exhibited this process. When they do occur, the landslide failures range in thickness from a few tens of centimeters to more than 6 m, and generally involve the soil and colluvium-mantled hillslopes. Surficial landslide failures in burned areas most frequently occur in response to prolonged periods of storm rainfall, or prolonged rainfall in combination with rapid snowmelt or rain-on-snow events.     

Limits of sediment transfer in an alpine debris-flow catchment,Illgraben, Switzerland

The 9.5 km² Illgraben catchment, located in the Rhône valley in the Central Alps of Switzerland, is one of the most active debris flow torrents in the Alps. In this paper we present sediment yield data collected in 2006 for segments where hillslopes and channels form a fully connected network and contrast these with sediment yields measured for disconnected hillslopes. The data reveal that sediment yields are 1–2 orders of magnitude larger in segments where hillslopes are connected with the channel network than on disconnected hillslopes. Support for this conclusion is provided by observations made on 1959, 1999 and 2004 aerial photographs that the vegetation cover in the disconnected segments is still intact, whereas denudation rates of several centimeters per year in the connected segments have inhibited the establishment of a stable vegetation cover. Furthermore, sediment supplied from hillslopes during the past 40 years has temporarily accumulated along the Illgraben channel, indicating that the channel aggraded over this period and has not yet recovered. An implication of this observation is that initiation of debris flows in the Illgraben catchment is limited more by the availability of intense rainfall than sediment. In contrast, on disconnected hillslopes, sediment flux does not appear to be driven by precipitation.The petrographic composition of the Illgraben fan deposits indicates two distinct sediment sources, one related to rockfall and the other to landslides and debris flows. The presence of clasts from both sources implies multiple processes of erosion, deposition, mixing and re-entrainment in the catchment before the material is exported to the Illgraben fan and to the Rhône River. In addition, delivery of large amounts of coarse-grained sediment to the Rhône causes a modification of the flow pattern from meandering or anastomosing upstream to braided downstream. Hence, the direct connectivity between hillslope and channelized processes in the Illgraben catchment causes not only rapid topographic modifications in the catchment, but also morphologic adjustment in the Rhône valley downstream.     

推求地面径流单位线的Free Search算法

地面径流单位线在水文计算和水文预报工作中应用十分普遍,但当净雨时段较多时,传统的单位线推求方法存在工作量大、效率不高、过分依赖专业经验、计算过程复杂、模拟结果可能出现负值等问题。为此,本文引入源于动物群体迁移行为的自由搜索(Free Search,FS)算法,将FS单个动物每步探查行走的位置向量作为一条初始单位线,在FS算法迭代过程中再将所有的初始单位线转化为符合定义的单位线,从而利用动物群体的迁移行为推求出最优单位线。该方法的优点是概念清晰,操作简便,能够直接给出单位线的各时段流量值。三个典型实例模拟结果表明:利用FS推求单位线时,寻优一次即可得到满意结果,所得单位线不需要修正,并且能够有效克服锯齿现象,对地面径流量的模拟精度也较高。     

集水区降雨径流晌应的环境同位素实验研究

利用环境同位素氚和氧-18对实验集水区进行降雨和径流响应的研究表明:(1)地面径流必源于本次降雨的概念不确,其中往往有非本次降雨的水量.经对1986年～1989年各次降雨估算表明,非本次降雨贡献最高可达50.5%.(2)非饱和带壤中流和饱和带地下水径流中必有非本次降雨的水量,并与地面径流一样,在次降雨径流过程中有时程变化.(3)对不同径流组成的流量过程,非本次降雨所占的比重不同.通过分析,可知降雨径流相关关系中的一一对应假定不确切,认为:传统的降雨径流经验关系和单位线概念需重新考虑;传统的过程线经验划分方法和现行同位素划分方法的有关基本假定不完全符合实际.     

一个用于瞬时单位线参数更新的递归公式的推导及应用

在应用瞬时单位线(IUH)进行流域汇流模拟时,需要先率定其2个参数:线性水库的个数n及调蓄系数K.研究表明,用传统的矩法率定得到的IUH,对地面径流的模拟精度较差,需要对参数进行更新.此前袁作新给出了一个用于参数更新的递归公式.本文推导了一个新的递归公式,并与袁作新公式的应用效果进行了对比.新递归方程的推导原理是针对地面径流过程线的2个特征值:洪峰流量和峰现时间,按模拟值与观测值之间的比来调整参数n和K,以使调整后得到的地面径流尽量接近观测值,调整的目标函数选取确定性系数(DC)和误差的均方值(MSE).将新递归公式和袁作新公式应用于2个案例.结果表明,新递归公式在提高地面径流模拟精度上要优于袁作新公式,对于案例1,新递归公式能提高模拟精度,而袁作新公式则不能;时于案例2,新递归公式提高模拟精度的幅度要高于袁作新公式,DC值提高了1.1%,MSE值降低了14.2%.分析原因认为较大的参数调整幅度可能跨越最优点,而新递归公式对参数的调整幅度小于袁作新公式,从而使其能够找到最优点.应用本文的递归公式时,不宜调整次数过多,否则会降低地面径流模拟精度,对2个案例的最优调整次数均为2次,产生此现象的原因是增加调整次数虽然可以提高时洪峰流量和峰现时间的模拟精度,但不一定能提高对整个地面径流过程的模拟精度.     

Hydrochemical and hydrological processes in the different landscape zones of alpine cold region in China

Investigation of water sources and flow pathways is crucial to understand and evaluate the characteristics of surface water and groundwater systems. This article aims to identify the hydrochemical and hydrological processes in different landscape zones based on hydrochemical analyses of various samples, including samples from glacier, snow, frozen soil meltwater, surface water, groundwater, and precipitation, in the alpine cold region of China. Hydrochemical tracers indicated that chemical compositions are characterized by the Ca-HCO₃ type in the glacier-snow zone; the Mg-Ca-SO₄ type in the alpine cold desert zone; the Ca-HCO₃-SO₄ type in the marsh meadow zone; the Ca-Mg-HCO₃ type in the alpine shrub zone; and the Ca-Na-SO₄ type in the mountain grassland zone. An end-member mixing model was used for hydrograph separation. The results showed that the Mafengou River in the wet season was recharged by groundwater in the alpine cold desert and alpine shrub zones (67%), surface runoff in the glacier-snow zone (11%), surface runoff in the alpine cold desert zone (8%), thawed water from frozen soil in the marsh meadow and mountain grassland zones (9%), and direct precipitation on the river channel (5%). This study suggests that precipitation from the whole catchment yielded little direct surface runoff; precipitation was mostly transformed into groundwater or interflow and was then concentrated into the river channel. This study provides a scientific basis for evaluation and management of water resources in the basin.     

用BATS模型模拟径流的个例研究

为了解陆气间水热交换在径表成中的作用,本文采用ＢＡＴＳ模型模拟了淮河流域山区和平原在１９９１年汛期５０天的暴雨洪水过程,计算了径流,土壤温度和感热,通量,并用常规的汇流计算方法得到了流域出口断量过程线。     

考虑地表径流与地下渗流耦合的斜坡降雨入渗研究

为简化分析，在模拟斜坡降雨入渗暂态渗流时，通常没有考虑入渗和产流的耦合过程。笔者提出了一种新的考虑地表径流与地下渗流耦合的斜坡降雨入渗分析方法，它较好地模拟了入渗和产流的过程，并通过算例，研究了考虑与不考虑该耦合作用对斜坡孔隙水压力分布的影响。     

利用SWMM模型模拟岩溶峰丛洼地系统降雨径流过程——以桂林丫吉试验场为例

文章以桂林丫吉试验场为例,验证是否可以利用SWMM模型模拟以管道为主的岩溶峰丛洼地系统降雨径流过程。以洼地为单元,把研究区概化为由管道相连的6个次级汇水流域,选择Green-Amp入渗计算方法,同时考虑包气带裂隙水对管道的补给,运行SWMM模型计算出研究区管道总出口S31泉的流量曲线。结果显示：模拟流量变化过程与实测流量变化过程基本一致。说明该模型可以用来模拟岩溶峰丛洼地地区降雨径流过程。经验证,模拟时段内S31泉总量相对误差为19.1%。     

植被分布对小流域产流影响的数值实验

基于二维运动波模型,建立了一个适用于小流域场次降雨产汇流过程的动力学模型。通过对典型小流域内不同区域(坡顶、坡中及坡底区域及距离流域出口的远近)种植植被时产流过程的数值模拟,分析了小流域内植被分布对产流过程的影响。结果表明:植被分布及其特性对小流域场次降雨产流有较明显的影响;下游区域种植植被的减水效果和延滞洪峰作用优于上游区域,陡坡区域减水效果优于缓坡区域,且郁闭度越大,这种差别越明显;在该研究条件下,下游区的减水效果可达到上游区的3倍;30%郁闭度条件下减水效果可达10%和20%郁闭度条件下的3倍和1.4倍。     

Surface runoff in a torrent catchment area in Middle Europe and its prevention

The Schesa, a sinister contributory torrent to the Ill river near Bludenz (federal province of Vorarlberg) is the largest basin-shaped gully of Middle Europe and endangers the underlying villages by torrential debris flow and gigantic mass movements. The catchment is characterized by a complex geological situation, high annual precipitation and torrential rains from spring to early autumn, which cause enormous amounts of surface runoff. Based on field investigations comprising rain simulation experiments on representative plots, investigations on land-use, vegetation cover, soil physical characteristics, geology, hydrogeology and other features of the catchment area, surface runoff coefficient maps were developed. They formed the basis for assessment of runoff potential for different scenarios in vegetation cover and land-use intensity. Calculation of runoff for the recurrent design event by use of an improved run-time method showed the urgent necessity of runoff reduction measures in large parts of the catchment area above the gully. Based on the modelling results a concept for reduction of both, surface runoff and amount of deep percolating water has been elaborated.     

滦河流域土地利用/覆被变化的水文响应

以滦河流域为研究区,利用1985和2000年土地利用数据,结合SWAT分布式水文模型定量评价了流域土地利用/覆被变化的水文效应,并分析了流域地表径流变化与主要景观类型的响应关系。结果表明:SWAT（Soil and Water Assessment Tool）模型可以较好地模拟滦河流域的月流量过程,在研究区具有较好的适用性;1985—2000年流域林地向草地和耕地的转变导致流域年均地表径流和总径流量分别增加了12.6%和5.1%;并使得流域年均地表径流变化空间差异显著,整体呈增加趋势,且主要受到林地变化的影响,而在三道河子以上集水区地表径流的变化则主要受到耕地景观的影响。合理规划土地利用格局,对于流域水资源可持续利用具有重要意义。     

Test of a distributed modelling approach to predict flood flows in the karst Suoimuoi catchment in Vietnam

The major obstacles for modelling flood processes in karst areas are a lack of understanding and model representations of the distinctive features and processes associated with runoff generation and often a paucity of field data. In this study, a distributed flood-modelling approach, WetSpa, is modified and applied to simulate the hydrological features and processes in the karst Suoimuoi catchment in northwest Vietnam. With input of topography, land use and soil types in a GIS format, the model is calibrated based on 15 months of hourly meteorological and hydrological data, and is used to simulate both fast surface and conduit flows, and groundwater discharges from karst and non-karst aquifers. Considerable variability in the simulation accuracy is found among storm events and within the catchment. The simulation results show that the model is able to represent reasonably well the stormflows generated by rainfall events in the study catchment.     

Study of cloudburst and flash floods around Leh, India, during August 4–6, 2010

Leh and surrounding region of the Ladakh mountain range in the trans-Himalaya experienced multiple cloudbursts and associated flash floods during August 4–6, 2010. However, 12.8 mm/day rainfall recorded at the nearest meteorological station at Leh did not corroborate with the flood severity. For better understanding of this event, hydrological analysis and atmospheric modeling are carried out in tandem. Two small catchments (<3 km²) were studied along the stream continuum to assess the flood characteristics to identify the cloudburst impact zones. Peak flood discharges were estimated close to the head wall region and at the catchment outlet of the Leh town and the Sabu eastern tributary catchments. Storm runoff depth is estimated by developing a triangular hydrograph by using the known time base of the flood hydrograph. This triangular hydrographs have been transformed further into storm hydrographs to gain a better understanding of the storm duration by using the dimensionless hydrograph method at selected cross sections. Storm duration is estimated by using the relationship between time to peak and time of concentration of the catchment. The peak flood estimates ranged from 122(±35 %) m³/s for Leh town catchment (2.393 km²), 545(±35 %) m³/s for Sabu eastern tributary catchment (2.831 km²) to 1,070(±35 %) m³/sec for Sabu catchment (64.95 km²). To assess the atmospheric processes associated with this event, a triple nest simulation (27, 9 and 3 km) is performed using Advanced Research Weather Research and Forecasting (WRF) modeling system. The simulation does show the evolution of the event from August 4 to 6, 2010. Observation constraints, orographic responses, etc. make such analysis complex at such scale. Independent estimate by the atmospheric process model and the hydrological method shows the storm depth of 70 mm and 91.8(±35 %) mm, respectively, in catchment scale. Hydrological evaluation further refined the spatial and temporal extents of the cloudbursts in the respective catchments with an estimated storm depth of 209(±35 %) mm in 11.9 min and 320(±35 %) in 8.8 min occurring in an area of 0.842–1.601 km², respectively. This study shows that the insight developed on the cloudburst phenomena by the atmospheric and the hydrological modeling is hugely constrained by the spatial and temporal scales of data used for the analysis. Apart from this, study also highlighted the regular occurrence of cloudburst events over this region in the recent past. Most of such events go unreported due to lack of monitoring mechanisms in the region and weaken our ability to understand these events in complete perspective.