STABILITY ANALYSIS OF RIVERBANK SUBJECT TO SEEPAGE

The stability of riverbanks subject to seepage is studied experimentally and theoretically in this paper. By including seepage in a 3-dimensional theoretical analysis, the study first shows how the critical slope or angle of repose of a cohesionless material is related to the ratio of the hydraulic gradient of seepage to its critical value under the fluidization condition. The critical stable slope is shown to be related to not only the hydraulic gradient but also the seepage direction. Measured laboratory data reasonably fit well with the theoretical relationship for the case of injection and suction. The data reveal that the slope is reduced with injection and increased with suction, respectively. Additionally, the study identifies the seepage direction which results in a minimum critical stable slope for a certain hydraulic gradient of seepage.     

Effect of seepage on initiation of cohesionless sediment transport

This paper presents theoretical analyses and experimental results of seepage effects, especially downward seepage, on the initiation of cohesionless sediment particles. The theoretical analysis examines how the additional seepage force acts to modify the critical shear stress for sediment entrainment. Laboratory experiments were conducted using medium sand with diameter of 0.9 mm with downward seepage to quantitatively show suction effects on sediment entrainment. The critical shear stresses with different suction rates were calculated using the experimental results. The measured data together with published results provide an overall view on seepage effects on the initiation of cohesionless sediment transport. Depending on whether seepage is in the form of injection or suction, it will either increase or decrease the critical shear stress. The result reveals that the ratio of drag force at the threshold condition with seepage to that without seepage is dependent on the ratio of the hydraulic gradient with seepage to its value at the quick condition.     

Hydrodynamic characteristics of rill flow on steep slopes

Rill erosion is a dominant sediment source on sloping lands. However, the amount of soil loss from rills on steep slopes is vastly more than that on gentle slopes because of differences in rill shape and hydraulic patterns. The aims of this paper are to determine the hydrodynamic characteristics of rills and the friction coefficients in steep slope conditions and to propose modifications of some hydraulic parameters used in soil loss prediction models. A series of inflow experiments was conducted on loess slopes. The results show that the geometric and hydraulic properties of rill on the steep loess slopes, which are characterized by the mean width of cross sections, mean velocity and mean depth of flow, are related to discharge and slope gradient in power functions. However, the related exponents to discharge are 0.26, 0.48 and 0.26, respectively, which are different from the exponents derived in previous studies, which were conducted on gentle slopes. The Manning roughness coefficient ranged from 0.035 to 0.071, with an average of 0.0536, and the Darcy–Weisbach friction coefficients varied from 0.4 to 1.9. The roughness coefficients are closely related to the Reynolds numbers and flow volumes; however, the correlations vary with slope gradient. The roughness coefficients are directly proportional to the Reynolds number and the flow volume on steep slopes, in contrast with the roughness coefficients found on gentle slopes, which decrease as the Reynolds number and flow volume increase. This difference is caused by the interactions among the hydraulics of the flow, the shape of the rills and the sediment concentrations on steep slopes. The results indicate that parameters used in models to predict rill erosion have to be modified according to slope gradient. Copyright © 2015 John Wiley & Sons, Ltd.     

Response of soil detachment rate to the hydraulic parameters of concentrated flow on steep loessial slopes on the Loess Plateau of China

Using hydraulic parameters is essential for describing soil detachment and developing physically based erosion prediction models. Many hydraulic parameters have been used, but the one that performs the best for describing soil detachment on steep slopes when the lateral expansion (widening) of rills is not limited has not been identified. An indoor concentrated flow scouring experiment was performed on steep loessial slopes to investigate soil detachment rates for different flow rates and slope gradients. The experiments were conducted on a slope‐adjustable plot (5 m length, 1 m width, 0.5 m depth). Sixteen combinations of 4 flow rates (10, 15, 20, and 25 L/min) and 4 slope gradients (17.6%, 26.8%, 36.4%, and 46.6%) were investigated. The individual and combined effects of slope gradient and flow hydraulic parameters on soil detachment rate were analysed. The results indicated that soil detachment rate increased with flow rate and slope gradient. Soil detachment rate varied linearly and exponentially with flow rate and slope gradient, respectively. Multivariate, nonlinear regression analysis indicated that flow depth exerted the greatest influence on the soil detachment rate, followed by unit discharge per unit width, slope gradient, and flow rate in this study. Shear stress and stream power could efficiently describe the soil detachment rate using a power equation. However, the unit stream power and unit energy of the water‐carrying section changed linearly with soil detachment rate. Stream power was an optimal hydraulic parameter for describing soil detachment. These findings improve our understanding of concentrated flow erosion on steep loessial slopes.     

Development of an electronic seepage chamber for extended use in a river

Seepage chambers have been used to characterize the flux of water across the water-sediment interface in a variety of settings. In this work, an electronic seepage chamber was developed specifically for long-term use in a large river where hydraulic gradient reversals occur frequently with river-stage variations. A bidirectional electronic flowmeter coupled with a seepage chamber was used to measure temporal changes in the magnitude and direction of water flux across the water-sediment interface over an 8-week period. The specific discharge measured from the seepage chamber compared favorably with measurements of vertical hydraulic gradient and previous specific discharge calculations. This, as well as other supporting data, demonstrates the effectiveness of the electronic seepage chamber to accurately quantify water flux in two directions over a multimonth period in this setting. The ability to conduct multimonth measurements of water flux at a subhourly frequency in a river system is a critical capability for a seepage chamber in a system where hydraulic gradients change on a daily and seasonal basis.     

Multiple failure modes-based practical calculation model on comprehensive risk for levee structure

Risk calculation of levee in the complex environment has significant theoretical values and practical meanings. As there still exists problems in present risk calculation models, a simple and efficient calculation model of the comprehensive levee risk is needed to be established. This paper studies the comprehensive risk calculation model of levee with multiple failure modes based on the analysis of levee instability and seepage failure. Firstly, coupling calculation of the seepage field and stress field is made using the finite element method, and safety factor of levee slope and critical failure hydraulic gradient of levee foundation are determined according to the strength reduction and piping theory. Then, particle swarm optimization is applied to conduct a comparative study on potential impact factors of levee stability and seepage, thus to determine explicit expressions for instability of levee body and seepage failure of levee foundation. Finally, Monte Carlo method is introduced to simulate the levee structure stochastically and calculated the comprehensive levee risk. Calculation results of the example show that the risk calculation method proposed in this paper has higher computational efficiency and can provide references for the decision of levee reinforcement.     

Landscape position,surface hydraulic gradients and erosion processes

Different hydraulic gradients, especially due to seepage or drainage, at different locations on a hillslope profile may have a profound effect on the dominant erosion processes. A laboratory study was designed to simulate hillslope processes and quantify effects of surface hydraulic gradients on erosion for a Glynwood clay loam soil (fine, illitic, mesic Aquic Hapludalf). A 5 m long, 1·2 m wide soil pan was used at 5 and 10 per cent slopes with an external watering tube to vary the soil bed's hydrological conditions. Different combinations of slope steepness with seepage or drainage gradients were used to simulate the hydrologic conditions on a 5 m segment of a hillslope profile. Runoff samples were taken during rainfall-only and rainfall with added inflow. Results showed that, under drainage conditions, interrill processes dominated and rilling was limited. The surface contained scattered crescent-shaped pits after the run. Under seepage conditions, rilling processes dominated and the inflow introduced at the top of the soil pan further accelerated the headward erosion of the rills. Erosion rates increased by as much as 60 times under seepage conditions representative of the lower backslope when compared to drainage conditions that generally occur at the upper backslope. This indicated that rills and gullies on backslopes and footslopes may be catalysed or enhanced by seepage conditions rather than form from flow hydraulic shear stress alone. An understanding of spatial and temporal changes that affect both hillslope hydrology and erosional processes is needed to develop accurate process-based erosion prediction models. This knowledge may lead to different management practices on landscape positions where seepage occurs. © 1998 John Wiley & Sons, Ltd.     

Morphology,structure and flow phases in soil pipes developing in forested hillslopes underlain by a Quaternary sand–gravel formation,Hokkaido, northern main island in Japan

Preferential flow pathways, such as soil pipes, are usually present in the soil of slopes. Subsurface flow through the soil pipes affects the subsurface drainage system and is responsible for sediment removal from slopes. However, a record of the inner structure of soil pipes has rarely been reported for slopes. A fibrescope examination of the morphology and flow phases in soil pipes in hillslopes underlain by a Quaternary sand–gravel formation provided the following information: the main pores of the soil pipes ran mostly parallel with the slope gradient; the cross‐sections of the soil pipes were approximately circular; and occurred on a few occasions; with some triple junctions being present. In addition, both full flow and partly full‐depth conditions occurred simultaneously in the soil pipe. The full flow condition has long been used in hydrological studies to model the pipe flow mechanism. Both the full flow condition and the partly full‐depth condition, however, must be examined closely in order to evaluate the subsurface hydrology in heterogeneous soil and the hydrogeomorphological processes of subsurface hydraulic erosion. Copyright © 2000 John Wiley & Sons, Ltd.     

Velocities of shallow overland flow on semiarid hillslopes

A series of 188 rainfall plot simulations was conducted on grass, shrub, oak savanna, and juniper sites in Arizona and Nevada. A total of 897 flow velocity measurements were obtained on 3.6% to 39.6% slopes with values ranging from 0.007 m s^‐1 to 0.115 m s^‐1. The experimental data showed that shallow flow velocity on rangelands was related to discharge and ground litter cover and was largely independent of slope gradient or soil characteristics. A power model was proposed to express this relationship. These findings support the slope–velocity equilibrium hypothesis. Namely, eroding soil surfaces evolve such that steeper areas develop greater hydraulic roughness. As a result overland flow velocity becomes independent of the slope gradient over time. Our findings have implications for soil erosion modeling suggesting that hydraulic friction is a dynamic, slope and discharge dependent property. Copyright © 2018 John Wiley & Sons, Ltd.     

Seepage erosion and its implication to the formation of amphitheatre valley heads: A case study at Obara,Japan

Seepage erosion was investigated in an amphitheatre with a semicircular valley head, steep slopes, and a flat bottom developed in granodiorite hills at Obara, Aichi prefecture, Japan. A high sediment yield occurred where the measuring sites were located at the base of the landslide debris in the base of the convex slopes, whereas sediment outflows were small where the measuring sites were located at the base of the strong convex slopes. This implies that the seepage erosion was an effective agent for removal of debris deposited at the base of the slope. Small landslides can be found at the lower slopes within the area of the observed amphitheatre. The slope stability analysis and subsurface water observation of the lower slope suggest that the small landslides in this amphitheatre are due to over-steepened slopes, and relatively insensitive to subsurface water status. Colluvium in the flat valley bottom thinly covers the bedrock surface. Therefore the topography of the amphitheatre was found to be formed by parallel retreat of slopes by the repetition of basal seepage erosion and subsequent small landslides.     

How hydrological factors initiate instability in a model sandy slope

Knowledge of the mechanisms of rain‐induced shallow landslides can improve the prediction of their occurrence and mitigate subsequent sediment disasters. Here, we examine an artificial slope's subsurface hydrology and propose a new slope stability analysis that includes seepage force and the down‐slope transfer of excess shear forces. We measured pore water pressure and volumetric water content immediately prior to a shallow landslide on an artificial sandy slope of 32°: The direction of the subsurface flow shifted from downward to parallel to the slope in the deepest part of the landslide mass, and this shift coincided with the start of soil displacement. A slope stability analysis that was restricted to individual segments of the landslide mass could not explain the initiation of the landslide; however, inclusion of the transfer of excess shear forces from up‐slope to down‐slope segments improved drastically the predictability. The improved stability analysis revealed that an unstable zone expanded down‐slope with an increase in soil water content, showing that the down‐slope soil initially supported the unstable up‐slope soil; destabilization of this down‐slope soil was the eventual trigger of total slope collapse. Initially, the effect of apparent soil cohesion was the most important factor promoting slope stability, but seepage force became the most important factor promoting slope instability closer to the landslide occurrence. These findings indicate that seepage forces, controlled by changes in direction and magnitude of saturated and unsaturated subsurface flows, may be the main cause of shallow landslides in sandy slopes. Copyright © 2013 John Wiley & Sons, Ltd.     

The effects of rock fragments on hydrologic and hydraulic responses along a slope

A model describing the three‐dimensional matrix flow along a slope with rock fragments or impermeable blocks was developed. The model was combined with modified Picard's iteration to ensure mass conservation in the unsaturated flow. We found that rock fragments obstruct water flow along the slope. The groundwater table must be raised to provide a sufficient pore water pressure gradient to facilitate water flow, but higher pore water pressure may induce slope failure. We also conducted a bench‐scale laboratory flume experiment to examine the effects of impermeable blocks on downstream seepage flow. In addition, a numerical experiment was conducted to examine how different arrangements of impermeable blocks affect downstream seepage flow and pore water pressure. This research demonstrated that the hydraulic phenomena were affected when impermeable blocks were present, and pore water pressure increased as the position of impermeable blocks was lowered. Copyright © 2006 John Wiley & Sons, Ltd.     

A Tube Seepage Meter for In Situ Measurement of Seepage Rate and Groundwater Sampling

We designed and evaluated a “tube seepage meter” for point measurements of vertical seepage rates (q), collecting groundwater samples, and estimating vertical hydraulic conductivity (K) in streambeds. Laboratory testing in artificial streambeds show that seepage rates from the tube seepage meter agreed well with expected values. Results of field testing of the tube seepage meter in a sandy‐bottom stream with a mean seepage rate of about 0.5 m/day agreed well with Darcian estimates (vertical hydraulic conductivity times head gradient) when averaged over multiple measurements. The uncertainties in q and K were evaluated with a Monte Carlo method and are typically 20% and 60%, respectively, for field data, and depend on the magnitude of the hydraulic gradient and the uncertainty in head measurements. The primary advantages of the tube seepage meter are its small footprint, concurrent and colocated assessments of q and K, and that it can also be configured as a self‐purging groundwater‐sampling device.     

Hillslope response to tectonic forcing in threshold landscapes

Hillslopes are thought to poorly record tectonic signals in threshold landscapes. Numerous previous studies of steep landscapes suggest that large changes in long‐term erosion rate lead to little change in mean hillslope angle, measured at coarse resolution. New LiDAR‐derived topography data enables a finer examination of threshold hillslopes. Here we quantify hillslope response to tectonic forcing in a threshold landscape. To do so, we use an extensive cosmogenic beryllium‐10 (¹⁰Be)‐based dataset of catchment‐averaged erosion rates combined with a 500 km² LiDAR‐derived 1 m digital elevation model to exploit a gradient of tectonic forcing and topographic relief in the San Gabriel Mountains, California. We also calibrate a new method of quantifying rock exposure from LiDAR‐derived slope measurements using high‐resolution panoramic photographs. Two distinct trends in hillslope behavior emerge: below catchment‐mean slopes of 30°, modal slopes increase with mean slopes, slope distribution skewness decreases with increasing mean slope, and bedrock exposure is limited; above mean slopes of 30°, our rock exposure index increases strongly with mean slope, and the prevalence of angle‐of‐repose debris wedges keeps modal slopes near 37°, resulting in a positive relationship between slope distribution skewness and mean slope. We find that both mean slopes and rock exposure increase with erosion rate up to 1 mm/a, in contrast to previous work based on coarser topographic data. We also find that as erosion rates increase, the extent of the fluvial network decreases, while colluvial channels extend downstream, keeping the total drainage density similar across the range. Our results reveal important textural details lost in 10 or 30 m resolution digital elevation models of steep landscapes, and highlight the need for process‐based studies of threshold hillslopes and colluvial channels. Copyright © 2012 John Wiley & Sons, Ltd.     

Local‐scale variability of seepage and hydraulic conductivity in a shallow gravel‐bed river

Seepage rate and direction measured with a seepage metre modified for use in flowing water were greatly variable along a 300‐m reach of a shallow, gravel‐bed river and depended primarily on the local‐scale bed topography. The median value of seepage measured at 24 locations was 24 cm/day, but seepage measured at specific sites ranged from ?340 to +237 cm/day. Seepage also varied substantially over periods of hours to days and occasionally reversed direction in response to evolution of the sediment bed. Vertical hydraulic conductivity was related to seepage direction and was larger during upward seepage than during downward seepage; with differences ranging from 4 to 40% in areas of active sediment transport to more than an order of magnitude in areas where current was too slow to mobilize bed sediment. Seepage was poorly related to hydraulic gradient measured over vertical distances of 0·3 m and appeared to be opposite the hydraulic gradient at 18% of the locations where both parameters were measured. Results demonstrate the scale dependence of these measurements in coarse‐grained hyporheic settings and indicate that hydraulic gradients should be determined over a much shorter vertical increment if used to indicate exchange across the sediment–water interface. Published in 2009 by John Wiley & Sons, Ltd.     

温度作用下煤层瓦斯解吸渗流规律数值模拟

如何提高煤层气渗透率是目前煤层气开采研究中的重要课题。基于煤层瓦斯渗流规律数学模型,利用COMSOL Multiphysics软件,对流-固-热耦合条件下的非等温煤层气解吸、渗流变化规律进行了数值模拟。结果表明,在注热条件下,煤层气渗流压力随着温度的增加而下降,且下降速度加剧,压力差越大,气体从高压区域流向低压区域的渗流速度越快。气体在煤层中径向流向井口,井口附近压力的梯度增大,气体渗流速度较快;在未受到加热影响的区域,煤层气不受外加热量影响,煤层气解吸速率保持不变;注热后煤层温度升高,可以加快煤层气渗流速度、提高渗透率、增加煤层气产量。研究成果可为煤层中注热开采煤层气的工程实践提供相应的理论依据。     

Flow and deformation failure of sandy slopes

The effects of earthquake induced pore pressure on seismic and post seismic stability conditions of cohesionless slopes are investigated with reference to the infinite slope scheme. In cohesionless slopes the shear strength reduction caused by pore pressure build-up may lead the slope to a deformation failure or to a flow failure if liquefaction conditions are approached. Two critical values of the seismic induced pore pressure ratio are introduced to evaluate the effect of shear strength reduction on the slope failure mechanism. The results are given in the form of stability charts and a procedure for the evaluation of the seismic stability condition is described. The procedure gives useful information about the failure mechanism that slopes may exhibit and the displacement analysis which should be carried out.     

Controls on the stability and inclinations of hillslopes formed on hard rock

In identifying controls on rock slope form a distinction is made between: (1) rock slopes with joints which dip steeply out of a cliff and hence are subject to mass failure of the rock mass above a critical joint; and (2) rock slopes with inclinations which are either in equilibrium with the mass strength of their rocks, or have profiles which will develop towards strength equilibrium as cross joints open. In the first class of slope, stability results not just from the basic frictional resistance of the rock but also from the frictional roughness along the critical joint and from the normal stress acting across that joint. Stability may be reduced by weathering and loss of strength of the joint wall rock. As a result of normal stress variations with depth, induced by overburdens, high cliffs which are not undercut have a concave profile. The second group of slopes includes those with inclinations controlled at the scale of individual joint blocks, buttressed slopes and those on unjointed rock masses. Buttressed and unjointed rock masses develop towards a condition of mass strength equilibrium as cross joints open. Strength equilibrium slopes may be recognized by application of a rock mass strength classification proposed for geomorphic purposes. Eleven propositions are formulated which identify controls on rock slope development and some consequences of these controls.     

Bank undercutting and tension failure by groundwater seepage: predicting failure mechanisms

Groundwater seepage can lead to the erosion and failure of streambanks and hillslopes. Two groundwater instability mechanisms include (i) tension failure due to the seepage force exceeding the soil shear strength or (ii) undercutting by seepage erosion and eventual mass failure. Previous research on these mechanisms has been limited to non‐cohesive and low cohesion soils. This study utilized a constant‐head, seepage soil box packed with more cohesive (6% and 15% clay) sandy loam soils at prescribed bulk densities (1.30 to 1.70 Mg m^?3) and with a bank angle of 90° to investigate the controls on failure mechanisms due to seepage forces. A dimensionless seepage mechanism (SM) number was derived and evaluated based on the ratio of resistive cohesion forces to the driving forces leading to instability including seepage gradients with an assumed steady‐state seepage angle. Tension failures and undercutting were both observed dependent primarily on the saturated hydraulic conductivity, effective cohesion, and seepage gradient. Also, shapes of seepage undercuts for these more cohesive soils were wider and less deep compared to undercuts in sand and loamy sand soils. Direct shear tests were used to quantify the geotechnical properties of the soils packed at the various bulk densities. The SM number reasonably predicted the seepage failure mechanism (tension failure versus undercutting) based on the geotechnical properties and assumed steady‐state seepage gradients of the physical‐scale laboratory experiments, with some uncertainty due to measurement of geotechnical parameters, assumed seepage gradient direction, and the expected width of the failure block. It is hypothesized that the SM number can be used to evaluate seepage failure mechanisms when a streambank or hillslope experiences steady‐state seepage forces. When prevalent, seepage gradient forces should be considered when analyzing bank stability, and therefore should be incorporated into commonly used stability models. Copyright © 2013 John Wiley & Sons, Ltd.     

Influence of a thin veneer of low‐hydraulic‐conductivity sediment on modelled exchange between river water and groundwater in response to induced infiltration

A thin layer of fine‐grained sediment commonly is deposited at the sediment–water interface of streams and rivers during low‐flow conditions, and may hinder exchange at the sediment–water interface similar to that observed at many riverbank‐filtration (RBF) sites. Results from a numerical groundwater‐flow model indicate that a low‐permeability veneer reduces the contribution of river water to a pumping well in a riparian aquifer to various degrees, depending on simulated hydraulic gradients, hydrogeological properties, and pumping conditions. Seepage of river water is reduced by 5–10% when a 2‐cm thick, low‐permeability veneer is present on the bed surface. Increasing thickness of the low‐permeability layer to 0·1 m has little effect on distribution of seepage or percentage contribution from the river to the pumping well. A three‐orders‐of‐magnitude reduction in hydraulic conductivity of the veneer is required to reduce seepage from the river to the extent typically associated with clogging at RBF sites. This degree of reduction is much larger than field‐measured values that were on the order of a factor of 20–25. Over 90% of seepage occurs within 12 m of the shoreline closest to the pumping well for most simulations. Virtually no seepage occurs through the thalweg near the shoreline opposite the pumping well, although no low‐permeability sediment was simulated for the thalweg. These results are relevant to natural settings that favour formation of a substantial, low‐permeability sediment veneer, as well as central‐pivot irrigation systems, and municipal water supplies where river seepage is induced via pumping wells. Published in 2011 by John Wiley & Sons, Ltd.