Impacts of water conservancy and soil conservation measures on annual runoff in the Chaohe River Basin during 1961-2005

National Key Technology RD Program [2006BAD03A02]

辽河流域降雨侵蚀力的时空变化分析

降雨侵蚀力是反映流域降雨侵蚀能力的综合指标之一。根据辽河流域10个气象站的日降雨量资料,利用日降雨侵蚀力模型估算辽河流域的降雨侵蚀力。结果表明:辽河流域降雨侵蚀力的空间变异与降雨量的空间分布趋势基本一致,由东南向西北递减,变化于1 000~3 800 MJ.mm/(hm2.h.a)之间;降雨侵蚀力年内集中度高,6—8月3个月约占全年的80%;降雨侵蚀力年际变化大,年际变率Cv在0.367~0.649之间,采用时序系列的Mann-Kendall检验表明,降雨侵蚀力并无显著变化趋势;特别是在流域水土流失严重的西辽河地区,年降雨侵蚀力较小,但年内集中程度大,年际变化更突出。

基于日雨量的延安地区降雨侵蚀力动态特征分析

延安地区是黄土高原水土流失最为严重的地区之一。利用延安气象站1951-2005年的日降雨量数据,采用日雨量侵蚀力模型估算延安地区降雨侵蚀力,结果表明:该地区降雨侵蚀力主要集中在6-9月,占到了全年的85.6%。年降雨侵蚀力的平均值为1765.73MJ.mm/(hm2.h),55年间,年降雨侵蚀力变异程度适中,从整体上看,趋势保持平稳,其离差系数Cv和变异趋势系数r分别为0.41和-0.071。

施加PAM的坡地降雨入渗过程及其模型研究

作为一种新型高效的土壤结构改良剂 PAM,在一定条件下能显著地提高土壤入渗能力 ,减小坡面径流。采用室内人工模拟降雨试验 ,研究了不同 PAM覆盖度下降雨产流随时间的动态变化及其与雨强、坡度的关系 ,以及入渗率随时间的变化规律 ,分析了 PAM、坡度、雨强对入渗率的影响。结果表明 ,地表施加 PAM后土壤的入渗率及稳定入渗率都比未施 PAM显著提高。通过对比施加 PAM后 Kostiakov入渗模型与 Horton入渗模型的显著性 ,表明 Horton入渗公式的适用性更好

Modeling on Runoff Concentration Caused by Rainfall on Hillslopes and Application in Maoping Slope

Based on the fact that the concentration flowlines of overland flow depend on the surface landform of hillslope, a kinematic wave model was developed for simulating runoff generation and flow concentration caused by rainfall on hillslopes. The model-simulated results agree well with experimental observations. Applying the model to the practical case of Maoping slope, we obtained the characteristics of runoff generation and infiltration on the slope. Especially, the simulated results adequately reflected the confluent pattern of surface runoff, which offers a scientific foundation for designing the drainage engineering on the Maoping slope.

Modeling of Hillslope Runoff and Soil Erosion at Rainfall Events Using Cellular Automata Approach

National Science Fund for Distinguished Young Scholars of China [40225004]; National Natural Science Foundation of China [40471048]

Modeling on Runoff and Sediment Yield in Small Watersheds

The prediction and estimate of water and soil loss is fundamental important for understanding the effect of the spatial heterogeneity of underlying surfaces and preventing ecological environment deterioration. In this paper, a dynamic model of runoff and sediment yield in small watersheds is established. The proposed model includes three components: runoff generation caused by rainfall, soil erosion on hillslopes by overland flow, and runoff concentration and sediment transport on watersheds. Applying the proposed model, the runoff and sediment yield processes in a typical catchment on the loess plateau was estimated, which exhibited a good agreement between predicted results and observation.

分布式水文模型的误差分析

分布式水文模型以其具有明确物理意义的参数结构和对空间分异性的全面反映,能够准确详尽地描述和模拟流域内真实的降水径流过程而被广泛需求和关注。在模拟土地利用、土地覆盖、水土流失等各种变化过程的水文响应,面源污染、陆面过程、气候变化影响评价等诸多领域都有广泛的应用。模型的预报精度和误差至关重要,决定了模型的应用和推广。在分析分布式水文模型建立和验证过程的基础上,提出了模型的4类误差来源:被排除在外的因素引起的误差,实测历史记录资料的随机或系统误差,参数误差和模型结构误差,讨论了各类误差的分析与计算方法,为模型的发展和成长提供了依据。

Effects of soil type on leaching and runoff transport of rare earth elements and phosphorous in laboratory experiments

Through leaching experiments and simulated rainfall experiments, characteristics of vertical leaching of exogenous rare earth elements (REEs) and phosphorus (P) and their losses with surface runoff during simulated rainfall in different types of soils (terra nera soil, cinnamon soil, red soil, loess soil, and purple soil) were investigated. Results of the leaching experiments showed that vertical transports of REEs and P were relatively low, with transport depths less than 6 cm. The vertical leaching rates of REEs and P in the different soils followed the order of purple soil > terra nera soil > red soil > cinnamon soil > loess soil. Results of the simulated rainfall experiments (83 mm h(-1)) revealed that more than 92% of REEs and P transported with soil particles in runoff. The loss rates of REEs and P in surface runoff in the different soil types were in the order of loess soil > terra nera soil > cinnamon soil > red soil > purple soil. The total amounts of losses of REEs and P in runoff were significantly correlated.

Numerical Analysis of Rainfall Infiltration in the Slope with a Fracture

With the finite volume method, a 2D numerical model for seepage in unsaturated soil has been established to study the rainfall infiltration in the fractured slope.The result shows that more rain may infiltrate into the slope due to existing fracture and then the pore pressure rises correspondingly. Very probably, it is one of the crucial factors accounting for slope failure. Furthermore a preliminary study has been conducted to investigate the influence of various fracture and rainfall factors such as the depth, width and location of a crack, surface condition, rainfall intensity and duration. Pore pressure and water volumetric content during the transient seepage are carefully examined to reveal the intrinsic mechanism.

Numerical Study on Rainfall Infiltration in Rock-Soil Blocks

A mathematical model for the rain infiltration in the rock-soil slop has been established and solved by using the finite element method. The unsteady water infiltrating process has been simulated to get water content both in the homogeneous and heterogeneous media. The simulated results show that the rock blocks in the rock-soil slop can cause the wetting front moving fast. If the rain intensity is increased, the saturated region will be formed quickly while other conditions are the same. If the rain intensity keeps a constant, it is possible to accelerate the generation of the saturated region by properly increasing the vertical filtration rate of the rock-soil slop. However, if the vertical filtration rate is so far greater than the rain intensity, it will be difficult to form the saturated region in the rock-soil slop. The numerical method was verified by comparing the calculation results with the field test data.

A Continuation Method of Parameter Inversion for Non-Equilibrium Convection-Dispersion Equation

Based on the homotopy mapping, a globally convergent method of parameter inversion for non-equilibrium convection-dispersion equations (CDEs) is developed. Moreover, in order to further improve the computational efficiency of the algorithm, a properly smooth function, which is derived from the sigmoid function, is employed to update the homotopy parameter during iteration. Numerical results show the feature of global convergence and high performance of this method. In addition, even the measurement quantities are heavily contaminated by noises, and a good solution can be found.

Three-Dimensional Network Model And Its Parameter Calibration

A material model, whose framework is parallel spring-bundles oriented in 3-D space, is proposed. Based on a discussion of the discrete schemes and optimum discretization of the solid angles, a 3-D network cell consisted of one-dimensional components is developed with its geometrical and physical parameters calibrated. It is proved that the 3-D network model is able to exactly simulate materials with arbitrary Poisson ratio from 0 to 1/2, breaking through the limit that the previous models in the literature are only suitable for materials with Poisson ratio from 0 to 1/3. A simplified model is also proposed to realize high computation accuracy within low computation cost. Examples demonstrate that the 3-D network model has particular superiority in the simulation of short-fiber reinforced composites.

Comparison of Various Adhesion Contact Theories and the Influence of Dimensionless Load Parameter

Three models, JKR (Johnson, Kendall and Roberts), DMT (Derjaguin, Muller, and Toporov) andMD (Maugis-Dugdale),are compared with the Hertz model in dealing with nano-contact problems. It has been shown that both the dimensionless load parameter, P D P=.1/4

Plot Erosion Model Using Gray Relational Analysis Method

The main factors affecting interrill erosion-including runoff discharge, rainfall intensity, mean flow velocity, and slope gradient-were analyzed by using a gray relational analysis. An equation for interrill erosion was derived by coupling this analysis with dimensional and regression analyses. The values of erosion rates predicted by this equation were in good agreement with experimental observations.