土质陡坡降雨侵蚀的数学模型及求解方法

建立了土质陡坡降雨侵蚀的数学模型,模型包括陡坡细沟间坡面流控制方程、变沟宽陡坡细沟流控制方程、陡坡细沟间坡面侵蚀泥沙连续方程及陡坡变沟宽细沟侵蚀泥沙连续方程。通过紊流冲击分布概率确定土壤剥蚀方程,确定了泥沙源/汇项的表达形式及细沟高程与宽度的表达式。通过有限元和有限差分相结合的方法求解该数学模型,利用有限单元法对数学模型的水动力方程、运动方程及泥沙连续方程进行空间上的离散,利用有限差分法解决时间域的问题,得出了顺序求解的数值计算公式及模型数值求解的具体方法和步骤。     

Laboratory experiments of sediment transport from bare soil with a rill

Abstract

Mathematical models developed for quantification of sediment transport in hydrological watersheds require data collected through field or laboratory experiments, but these are still very rare in the literature. This study aims to collect such data at the laboratory scale. To this end, a rainfall simulator equipped with nozzles to spray rainfall was constructed, together with an erosion flume that can be given longitudinal and lateral slopes. Eighty experiments were performed, considering microtopographical features by pre-forming a rill on the soil surface before the start of each experiment. Medium and fine sands were used as soil, and four rainfall intensities (45, 65, 85 and 105 mm h^-1) were applied in the experiments. Rainfall characteristics such as uniformity, granulometry, drop velocity and kinetic energy were evaluated; flow and sediment discharge data were collected and analysed. The analysis shows that the sediment transport rate is directly proportional to rainfall intensity and slope. In contrast, the volumetric sediment concentration stays constant and does not change with rainfall intensity unless the slope changes. These conclusions are restricted to the conditions of experiments performed under rainfall intensities between and 105 mm h^-1 for medium and fine sands in a 136-cm-wide, 650-cm-long and 17-cm-deep erosion flume with longitudinal and lateral slopes varying between 5 and 20%.

Editor Z.W. Kundzewicz; Associate editor G. Mahé

Citation Aksoy, H., Unal, N.E., Cokgor, S., Gedikli, A., Yoon, J., Koca, K., Inci, S.B., Eris, E., and Pak, G., 2013. Laboratory experiments of sediment transport from bare soil with a rill. Hydrological Sciences Journal, 58 (7), 1505–1518.     

紫色土坡耕地细沟发育的形态演变

研究细沟发育的形态演变规律是认识细沟发生机理的前提,本研究通过人工降雨与微地貌扫描相结合的方法,调查了紫色土坡耕地细沟发生过程的形态演变特征.采用三维激光扫描仪扫描不同降雨时段坡面,生成高精度坡面DEM图像,统计各时段不同坡段等高线偏离度,量算坡面凹陷区的表面积、周长、直径以及沟头下切深度;结果表明,等高线偏离度越大细沟发育越完全,有细沟发生坡面等高线偏离度是无细沟坡面的1倍以上;凹陷区的表面积、周长、直径相对于单位降雨量的扩展速率均在细沟贯穿时段出现最大值;跌坑位置与坡面水平起伏关系不明显,统计的3条沟头分别发育于水平方向起伏的最高点、中间点和低点;跌坑下切深度及其单位降雨量的扩展速率在跌坑贯穿形成细沟时段出现最大值,之后逐渐降低,说明跌坑的形成是细沟造床的关键.     

Effect of particle density and inflow concentration of suspended sediment on bedload transport in rill flow

Laboratory flume experiments were carried out to evaluate the effect of particle density on bedload transport of sand‐sized particles and the effect of a suspended load of clay particles (kaolinite) on bedload transport of sand‐sized particles in rill flow conditions. Three materials in the range 400–600 µm were selected to simulate bedload transport of primary particles and aggregates: sand (2650 kg/m³), crushed brick (2450 kg/m³) and anthracite (1300–1700 kg/m³). In the two first experiments, two different methods were applied to determine bedload transport capacity of coarse particles for various conditions of flow discharge (from 2 to 15 L/min) and slope (2.2, 3 and 4%). In the third experiment, clear water was replaced with kaolinite–water mixture and bedload transport capacity of crushed brick particles was determined for a 4% slope and different concentrations of kaolinite (0, 7, 41 and 84 g/L). The results showed that bedload transport increased significantly with the decrease in particle density but the effect of particle density on transport rates was much less important than flow discharge. Velocity measurements of clear flow, flow mixed with coarse particles and coarse particles confirmed the existence of a differentiation between suspended load and bedload. In these experimental conditions, suspended load of kaolinite did not affect bedload rates of crushed brick particles. Three transport capacity formulae were tested against observed bedload rates. A calibration of the Foster formula revealed that the shear stress exponent should be greater than 1.5. The Low and the Govers unit stream power (USP) equations were then evaluated. The Low equation was preferred for the prediction of bedload rates of primary particles but it was not recommended in the case of aggregates of low density because of the limited experimental conditions applied to derive this equation. Copyright © 2008 John Wiley & Sons, Ltd.     

Quantifying the effects of plant litter in the topsoil on the soil detachment process by overland flow in typical grasslands of the Loess Plateau,China

Plant litter can be incorporated into topsoil by a natural process, affecting the soil erosion process. This is a widespread phenomenon in erosion-prone areas. This study was conducted to investigate the effect of litter incorporation on the process of soil detachment on the Loess Plateau, China. Four common plant litters (Bothriochloa ischaemum L. Keng., Artemisia sacrorum Ledeb., Setaria viridis L. Beauv., and Artemisia capillaris Thunb.) were collected, then incorporated into the silt loam soil at five rates (0.1, 0.4, 0.7, 1.0, and 1.3 kg m⁻²) on the basis of our field investigation. Twenty litter–soil treatments and one bare soil control were prepared. After 50 days of natural stabilization, 30 soil samples of each treatment were collected. We used a flume test to scour the soil samples under six flow shear stress conditions (5.66, 8.31, 12.21, 15.55, 19.15, and 22.11 Pa). The results showed that the different incorporated litter masses and morphological characteristics, such as litter tissue density (ranging from 0.52 to 0.68 g cm⁻³), length density (2.34 to 91.00 km m⁻³), surface area density (LSAD; 27.9 to 674.2 m² m⁻³), and volume ratio (0.003 to 0.050 m³ m⁻³), caused varied soil detachment capacities (0.043 to 4.580 kg·m⁻²·s⁻¹), rill erodibilities (0.051 to 0.237 s m⁻¹), and critical shear stresses (2.02 to 6.83 Pa). The plant litter incorporated within the soil reduced the soil detachment capacities by 38%–59%, lowered the rill erodibilities by 32%–46%, and increased the critical shear stresses by 98%–193% compared with the bare soil control. The soil containing B. ischaemum (L.) Keng. litter was more resistant to erosion. By comparing different parameters, we found that the contact area between the litter and soil was the main factor affecting the soil detachment process. The soil erosion resistance increased with the increasing contact area between the soil and litter. Furthermore, the litter incorporation effect on rill erodibility can be comprehensively reflected by LSAD (R² = .93; Nash–Sutcliffe efficiency = 0.79), which could be used to adjust the rill erodibility parameter in physical process-based soil erosion models.     

稀土元素示踪法在坡面土壤侵蚀垂直分布研究中的应用

通过野外全坡长侵蚀小区试验,证明了稀土元素示踪法在土壤侵蚀垂直分布研究中的可应用性。结果表明该径流小区产沙强度坡长存在三种变化模式:增长型,增长 递减型和波动型,其形成机理可能与径流过程有关。利用稀土元素在坡面的再分析,初步建立了估算坡面沟蚀率的关系式。     

黄土坡面细沟侵蚀发生的临界条件

通过玻璃水槽试验和径流冲刷试验，研究了坡面径流的流速分布和坡面细沟侵蚀发生的临界条件。研究表明：径流在坡面上并非是以均匀流形式运动，而是以滚波的形式运动的，同时本文运用能量守恒原理分析了坡面土壤侵蚀率（Dr）与径流能耗（△E）之间的关系，建立了给定土壤条件下坡面土壤侵蚀率估算模型：Dr=0.0026（In△E-2）。结果表明坡面土壤侵蚀的发生具有一定的临界能量条件，当径流能耗大于7.38（J）时坡面开始有细沟侵蚀发生。     

美国TH─100A型车载式钻机基本功能初评

通过在黑龙江省两个勘探区９个钻孔的试生产，对ＴＨ—１００Ａ型车载式钻机的基本功能、平台布置、钻进方法作了扼要评述，并提出了改进建议。     

Longitudinal velocity patterns and bed morphology interaction in a rill

Present‐day understanding of rill dynamics is hampered by a lack of detailed data on velocity distributions in rills. The latter are dif?cult to collect with traditional techniques due to the very low water depths and the relatively high ?ow velocities in rills. The objectives of this paper were to investigate the feasibility of miniaturized acoustic Doppler velocimeter (mADV) measurements in rill ?ow and to explore longitudinal variations in ?ow velocities and their relationship with rill bed morphology. Detailed data on longitudinal ?ow velocity were required to achieve these objectives. A 1·8 m long rill was formed freely in a ?ume at 5° slope and 0·001 m³ s^?1 discharge. Rill topography was characterized by an alternation of steps and pools. The ?ume surface was then ?xed to preserve rill roughness. A topographical scanning of the entire ?ume surface was made. Velocity was measured with a mADV along the rill, and at different depths. Flow depth in a longitudinal direction was also measured using an elevation gauge. A strong relationship exists between rill topography and ?ow hydraulics. Over steps, ?ow was unidirectional and rapidly accelerating until a threshold Froude number (Fn) value between 1·3 and 1·7 was reached and a hydraulic jump occurred leading to the formation of a pool. In the pool, the ?ow pattern was multidirectional and complex. The ?ow was subcritical when leaving the pool and accelerated over the next step until the threshold Froude number value was again reached. Energy loss in the rill was concentrated in the pools, mainly due to the action of a hydraulic jump. This mechanism of energy dissipation appeared to be an essential factor in rill formation and bedform evolution. Copyright © 2004 John Wiley & Sons, Ltd.