Enhanced application of root‐reinforcement algorithms for bank‐stability modeling

Riparian vegetation is known to exert a number of mechanical and hydrologic controls on bank stability. In particular, plant roots provide mechanical reinforcement to a soil matrix due to the different responses of soils and roots to stress. Root reinforcement is largely a function of the strength of the roots crossing potential shear planes, and the number and diameter of such roots. However, previous bank stability models have been constrained by limited field data pertaining to the spatial and temporal variability of root networks within stream banks. In this paper, a method is developed to use root‐architecture data to derive parameters required for modeling temporal and spatial changes in root reinforcement. Changes in root numbers over time were assumed to follow a sigmoidal curve, which commonly represents the growth rates of organisms. Regressions for numbers of roots crossing potential shear planes over time showed small variations between species during the juvenile growth phase, but extrapolation led to large variations in root numbers by the time the senescent phase of the sigmoidal growth curve had been reached. In light of potential variability in the field data, the mean number of roots crossing a potential shear plane at each year of tree growth was also calculated using data from all species and an additional sigmoidal regression was run. After 30 years the mean number of roots predicted to cross a 1 m shear plane was 484, compared with species‐specific curves whose values ranged from 240 roots for black willow trees to 890 roots for western cottonwood trees. In addition, the effect of spatial variations in rooting density with depth on stream‐bank stability was modeled using the bank stability and toe erosion model (BSTEM). Three root distributions, all approximating the same average root reinforcement (5 kPa) over the top 1 m of the bank profile, were modeled, but with differing vertical distributions (concentrated near surface, non‐linear decline with depth, uniform over top meter). It was found that stream‐bank F_S varied the most when the proportion of the failure plane length to the depth of the rooting zone was greatest. Copyright © 2008 John Wiley & Sons, Ltd.     

A new model to analyse the impact of woody riparian vegetation on the geotechnical stability of riverbanks

We present a geotechnical stability analysis for the planar failure of riverbanks, which incorporates the effects of root reinforcement and surcharge for mature stands of woody riparian vegetation. The analysis relies on a new method of representing the root distribution in the soil, which evaluates the effects of the vegetation's position on the bank. The model is used in a series of sensitivity analyses performed for a wide range of bank morphological (bank slope and height) and sedimentological (bank cohesion and friction angle) conditions, enabling discrimination of the types of bank environment for which vegetation has an effect on bank stability. The results indicate that woody vegetation elements have a maximal impact on bank stability when they are located at the ends of the incipient failure plane (i.e. at the bank toe or at the intersection of the failure plane with the floodplain) and that vegetation has a greater effect on net bank stability when it is growing on low, shallow, banks comprised of weakly cohesive sediments. However, the magnitude of these effects is limited, with vegetation typically inducing changes (relative to non‐vegetated banks) in simulated factors of safety of less than 5%. If correct, this suggests that the well documented effects of vegetation on channel morphology must be related to alternative process mechanisms (such as the interaction of vegetation with river flows) rather than the mechanical effects of vegetation on bank failure, except in special cases where the equivalent non‐vegetated bank has a highly marginal stability status. Copyright © 2007 John Wiley & Sons, Ltd.     

Measuring streambank erosion due to ground water seepage: correlation to bank pore water pressure,precipitation and stream stage

Limited information exists on one of the mechanisms governing sediment input to streams: streambank erosion by ground water seepage. The objective of this research was to demonstrate the importance of streambank composition and stratigraphy in controlling seepage flow and to quantify correlation of seepage flow/erosion with precipitation, stream stage and soil pore water pressure. The streambank site was located in Northern Mississippi in the Goodwin Creek watershed. Soil samples from layers on the streambank face suggested less than an order of magnitude difference in vertical hydraulic conductivity (K_s) with depth, but differences between lateral K_s of a concretion layer and the vertical K_s of the underlying layers contributed to the propensity for lateral flow. Goodwin Creek seeps were not similar to other seeps reported in the literature, in that eroded sediment originated from layers underneath the primary seepage layer. Subsurface flow and sediment load, quantified using 50 cm wide collection pans, were dependent on the type of seep: intermittent low‐flow (LF) seeps (flow rates typically less than 0·05 L min^?1), persistent high‐flow (HF) seeps (average flow rate of 0·39 L min^?1) and buried seeps, which eroded unconsolidated bank material from previous bank failures. The timing of LF seeps correlated to river stage and precipitation. The HF seeps at Goodwin Creek began after rainfall events resulted in the adjacent streambank reaching near saturation (i.e. soil pore water pressures greater than ?5 kPa). Seep discharge from HF seeps reached a maximum of 1·0 L min^?1 and sediment concentrations commonly approached 100 g L^?1. Buried seeps were intermittent but exhibited the most significant erosion rates (738 g min^?1) and sediment concentrations (989 g L^?1). In cases where perched water table conditions exist and persistent HF seeps occur, seepage erosion and bank collapse of streambank sediment may be significant. Copyright © 2007 John Wiley & Sons, Ltd.     

Physical environmental controls on riparian root profiles associated with black poplar (Populus nigra L.) along the Tagliamento River,Italy

Tree roots contribute to the resistance of riparian sediments to physical deformation and disintegration. Understanding reinforcement by roots requires information on root distributions within riparian soils and sediments. Continuous‐depth models or curves have been proposed to describe vertical root density variations, providing useful indicators of the types of function that may be appropriate to riparian trees, but have generally been estimated for terrestrial species or broad vegetation types rather than riparian species or environments. We investigated vertical distributions of roots >0.1 mm diameter of a single riparian tree species (Populus nigra L.) along the middle reaches of a single river (Tagliamento River, Italy), where Populus nigra dominates the riparian woodland. Root density (hundreds m^?2) and root area ratio (RAR in cm² m^?2) were measured within 10 cm depth increments of 24 excavated bank profiles across nine sites. Sediment samples, extracted from distinct strata within the profiles, were analysed for moisture content, organic matter content and particle size. Statistical analyses identified two groups of wetter and drier profiles and five sediment types. Following log_e‐transformation of root density and RAR, linear regression analysis explored their variation with depth and, using dummy variables, any additional influence of moisture and sediment type. Significant linear regression relationships were estimated between both root density and RAR and depth which explained only 15% and 8% of the variance in the data. Incorporating moisture and then sediment characteristics into the analysis increased the variance explained in root density to 29% and 36% and in RAR to 14% and 26%. We conclude that riparian tree root density and RAR are highly spatially variable and are poorly explained by depth alone. Complex riparian sedimentary structures and moisture conditions are important influences on root distributions and so need to be incorporated into assessments of the contribution of roots to river bank reinforcement. Copyright © 2016 John Wiley & Sons, Ltd.     

Computer program for stability analysis of steep,cohesive riverbanks

The ability to predict the stability of eroding riverbanks is a prerequisite for modelling alluvial channel width adjustment and a requirement for predicting bank erosion rates and sediment yield associated with bank erosion. Mass‐wasting of bank materials under gravity occurs through a variety of specific mechanisms, with a separate analysis required for each type of failure. This paper presents a computer program for the analysis of the stability of steep, cohesive riverbanks with respect to planar‐type failures. Planar‐type failures are common along stream channels destabilized by severe bed degradation. Existing stability analyses for planar‐type failures have a number of limitations that affect their physical basis and predictive ability. The computer program presented here is based on an analysis developed by Darby and Thorne. The software takes account of the geotechnical characteristics of the bank materials, the shape of the bank profile, and the relative elevations of the groundwater and surface water to estimate stability with respect to mass failure along a planar‐type failure surface. Results can be displayed either in terms of a factor of safety (ratio of resisting to driving forces), or probability of failure. The computer analysis is able to determine the relative amounts of bed degradation and bank‐toe erosion required to destabilize an initially stable bank. Data for the analysis are supplied in the form of either HEC‐2 hydrographic survey data files or user‐supplied bank profile data, in conjunction with user‐supplied geotechnical parameter values. Some examples, using data from the Upper Missouri River in Montana, are used to demonstrate potential applications of the software. Copyright © 2000 John Wiley & Sons, Ltd.     

Sediment transferring function of the lower reaches of the Yellow River influenced by drainage basin factors and human activities

Since 1986, with a sharp decrease in water dis-charges, the Yellow River has entered a period charac-terized by low discharges and seasonally occurring dry-ups[1,2]. Since 1999, more strict management of water diversion has been imposed, and therefore the dry-ups have been well under control. However, the lower reaches of the Yellow River is still predominated by low-discharges, and has become a man-induced shrinking river. In the past 40 years, significant effect of soil and water conservat…     

Effects of wet meadow riparian vegetation on streambank erosion. 2. Measurements of vegetated bank strength and consequences for failure mechanics

We measured the effect of wet meadow vegetation on the bank strength and failure mechanics of a meandering montane meadow stream, the South Fork of the Kern River at Monache Meadow, in California's Sierra Nevada. Streambanks colonized by ‘wet’ graminoid meadow vegetation were on average five times stronger than those colonized by ‘dry’ xeric meadow and scrub vegetation. Our measurements show that strength is correlated with vegetation density indicators, including stem counts, standing biomass per unit area, and the ratio of root mass to soil mass. Rushes appear better than sedges at stabilizing coarse bar surfaces, while sedges are far more effective at stabilizing actively eroding cut banks. Wet meadow floodplain vegetation creates a composite cut bank configuration (a cohesive layer overlying cohesionless materials) that erodes via cantilever failure. Field measurements and a geotechnical model of cantilever stability show that by increasing bank strength, wet meadow vegetation increases the thickness, width, and cohesiveness of a bank cantilever, which, in turn, increases the amount of time required to undermine, detach, and remove bank failure blocks. At Monache Meadow, it takes approximately four years to produce and remove a 1 m wide wet meadow bank block. Wet meadow vegetation limits bank migration rates by increasing bank strength, altering bank failure modes, and reducing bank failure frequency. Copyright © 2002 John Wiley & Sons, Ltd.     

Decreasing trend of sediment transfer function of the Upper Yellow River,China, in response to human activity and climate change

Abstract

An index (F_s) for sediment transfer function is introduced, based on the sediment budget at the channel scale. The purpose of this study is two-fold: to gain a deeper insight into how F_s is influenced by natural and human factors, and to provide some new knowledge for decision making in the management of the Upper Yellow River, China. Since 1960, the F_s of the Lanzhou to Toudaoguai reach of the Upper Yellow River shows a decreasing trend. At the drainage basin level, the decreased F_s can be explained by changes in precipitation and air temperature, as well as by a number of variables describing human activity, such as reservoir regulation, water diversion, and soil and water conservation. The higher temperature reduces the transfer function, while the larger runoff coefficient increases it. At the channel level, the decreased F_s can be explained by a number of variables of flow and sediment input. Three countermeasures for restoration of the F_s are suggested.

Editor Z.W. Kundzewicz     

Factors controlling the chemical evolution of travertine‐depositing rivers of the Barkly karst,northern Australia

Groundwaters feeding travertine‐depositing rivers of the northeastern segment of the Barkly karst (NW Queensland, Australia) are of comparable chemical composition, allowing a detailed investigation of how the rate of downstream chemical evolution varies from river to river. The discharge, pH, temperature, conductivity and major‐ion concentrations of five rivers were determined by standard field and laboratory techniques. The results show that each river experiences similar patterns of downstream chemical evolution, with CO₂ outgassing driving the waters to high levels of calcite supersaturation, which in turn leads to widespread calcium carbonate deposition. However, the rate at which the waters evolve, measured as the loss of CaCO₃ per kilometre, varies from river to river, and depends primarily upon discharge at the time of sampling and stream gradient. For example, Louie Creek (Q = 0·11 m³ s^?1) and Carl Creek (Q = 0·50 m³ s^?1) have identical stream gradients, but the loss of CaCO₃ per kilometre for Louie Creek is twice that of Carl Creek. The Gregory River (Q = 3·07 m³ s^?1), O'Shanassy River (Q = 0·57 m³ s^?1) and Lawn Hill Creek (Q = 0·72 m³ s^?1) have very similar gradients, but the rate of hydrochemical evolution of the Gregory River is significantly less than either of the other two systems. The results have major implications for travertine deposition: the stream reach required for waters to evolve to critical levels of calcite supersaturation will, all others things being equal, increase with increasing discharge, and the length of reach over which travertine is deposited will also increase with increasing discharge. This implies that fossil travertine deposits preserved well downstream of modern deposition limits are likely to have been formed under higher discharge regimes. Copyright © 2002 John Wiley & Sons, Ltd.     

Plant root and rhizome strength: Are there differences between and within species and rivers?

The strength and architecture of roots and other below-ground organs of riparian and aquatic plants affect plant resistance to uprooting and contribute to reinforcing river bank, bar and bed materials. Therefore, root properties are an important element in models for estimating river bank stability and such models may focus on the role of plants by using root strength–diameter relationships for the particular plant species that are present. Here we explore the degree to which there appear to be significant differences in strength–diameter relationships between and within species-specific data sets obtained for two riparian tree/shrub (Populus nigra, Salix alba) and two emergent aquatic macrophyte (Sparganium erectum, Phalaris arundinacea) species in different European river environments. While the analysed data sets were not specifically collected to answer these research questions, the results are sufficiently compelling to make the case for the collection of a more comprehensive data set and its rigorous analysis. This would allow recommendations to be made on the degree to which (i) species-specific or more general relationships between root/rhizome strength and diameter are appropriate, (ii) such relationships are applicable within and between rivers in different geographical regions and subject to different local environmental conditions, and (iii) further (minimalist) field observations are needed to calibrate such relationships for investigations of new locales or species. © 2018 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Controls on effective settling velocities of mud flocs in the Changjiang Estuary,China

Theoretical and experimental studies were undertaken to gain insight into physical parameters controlling the flocculation and settling properties of mud flocs in the Changjiang Estuary, China. The Rouse equation is applied to vertical profiles of suspended sediment concentration to determine the bulk mean settling velocity (w_s) of sediment suspended in the Changjiang Estuary. Both in situ point‐sampled and acoustically measured profiles of suspended mud concentrations were fit selectively. The calculated settling velocities w_s mainly ranged from 0·4 to 4·1 mm s^?1 for the point‐sampled data set, and from 1·0 to 3·0 mm s^?1 for the acoustically measured data set. Furthermore, the settling velocities of mud flocs increased with mean concentration (C?) of mud flocs in suspension and were proportional to increasing bottom shear stress (τ_b) of tidal flow. The best equation for the field settling velocity of mud flocs in the Changjiang Estuary can be expressed by the power law: w_s = mC?ⁿ (m, 1·14–2·37; n, 0·84–1·03). It is suggested that C? and τ_b were the dominant physical parameters controlling the flocculation and w_s of mud flocs in suspension. Copyright © 2004 John Wiley & Sons, Ltd.     

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.     

Snow accumulation in forests from ground and remote‐sensing data

Winter‐forest processes affect global and local climates. The interception‐sublimation fraction (F) of snowfall in forests is a substantial part of the winter water budget (up to 40%). Climate, weather‐forecast and hydrological modellers incorporate increasingly realistic surface schemes into their models, and algorithms describing snow accumulation and snow‐interception sublimation are now finding their way into these schemes. Spatially variable data for calibration and verification of wintertime dynamics therefore are needed for such modelling schemes. The value of F was determined from snow courses in open and forested areas in Hokkaido, Japan. The value of F was related to species and canopy‐structure measures such as closure, sky‐view fraction (SVF) and leaf‐area index (LAI). Forest structure was deduced from fish‐eye photographs. The value of F showed a strong linear correlation to structure: F = 0·44 ? 0·6 × SVF for SVF < 0·72 and F = 0 for SVF > 0·72, and F = 0·11 LAI. These relationships seemed valid for evergreen conifers, larch trees, alder, birch and mixed deciduous stands. Forest snow accumulation (S_F) could be estimated from snowfall in open fields (S_o) and to LAI according to S_F = S_o (1 ? 0·11 LAI) as well as from SVF according to S_F = S_o (0·56 + 0·6 SVF) for SVF < 0·72. The value of S_F was equal to S_o for SVF values above 0·72. The value of sky‐view fraction was correlated to the normalized difference snow index (NDSI) using a Landsat‐TM image for observation plots exceeding 1 ha. Variables F and S_F were related to NDSI for these plots according to: F = ?0·37NDSI + 0·29 and S_F = S_o (0·81 + 0·37NDSI). These relationships are somewhat hypothetical because plot‐size limitation only allowed one sparse‐forest observation of NDSI to be used. There is, therefore, a need to confirm these relationships with further studies. Copyright © 2004 John Wiley & Sons, Ltd.     

Soil moisture causes dynamic adjustments to root reinforcement that reduce slope stability

In steep soil‐mantled landscapes, the initiation of shallow landslides is strongly controlled by the distribution of vegetation, whose roots reinforce the soil. The magnitude of root reinforcement depends on the number, diameter distribution, orientation and the mechanical properties of roots that cross potential failure planes. Understanding how these properties vary in space and time in forests remains a significant challenge. Here we test the hypothesis that spatio‐temporal variations in root reinforcement along a hillslope occur as a function of topographic soil moisture gradients. To test this hypothesis we compared root reinforcement measurements from relatively dry, divergent noses to relatively wet, convergent hollows in the southern Appalachian Mountains, North Carolina, USA. Our initial results showed that root reinforcement decreased in areas of higher soil moisture because the tensile strength of roots decreased. A post hoc laboratory experiment further demonstrated that root tensile strength decreased as root moisture content increased. This effect is consistent with other experiments on stem woods showing that increased water content in the cell wall decreases tensile strength. Our experimental data demonstrated that roots can adjust to changes in the external root moisture conditions within hours, suggesting that root moisture content will change over the timescale of large storm events (hours–days). We assessed the effects of the dynamic changes in root tensile strength to the magnitude of apparent cohesion within the infinite slope stability model. Slopes can be considerably less stable when precipitation‐driven increases in saturated soil depth both increase pore pressures and decrease root reinforcement. Copyright © 2016 John Wiley & Sons, Ltd.     

Snowpack disrupts relationship between young water fraction and isotope amplitude ratio; approximately one fifth of mountain streamflow less than one year old

Previous “fraction of young water” (F_yw) estimates based on relative annual isotopic amplitudes in precipitation (A_p) and streamflow (A_s) produced low F_yw values in mountain catchments, which is contrary to extensive research that reports rapid water transmission in mountains. This study investigated this discrepancy by testing the effect of snow accumulation on the model that underpins the F_yw method. A Monte-Carlo analysis of simulations for 20,000 randomly-generated catchment model configurations used 10 years of precipitation inputs for the Upper Elbow River catchment in the Rocky Mountains (Alberta, Canada) to model discharge with and without snowpack storage of winter precipitation. Neither direct nor modified precipitation input produced a 1:1 relationship between A_s/A_p and F_yw, undermining the applicability of the original F_yw method in mountain watersheds with large seasonal snow accumulation. With snowpack-modified input a given A_s/A_p ratio corresponds to a range of F_yw values, which can still provide semi-quantitative information. In the small (435 km²) Elbow River catchment a F_yw range of 7–23% supports previous findings of rapid transmission in mountain catchments. Further analysis showed that the improved discharge prediction (Nash–Sutcliffe efficiency > 0.9) correlates with higher F_yw values and demonstrated that the interannual shifts in δ¹⁸O can be used to estimate of new water (<1 year) fraction in winter streamflow, and the estimate of 20% for the Elbow River further supports rapid transmission in mountain catchments.     

Factors controlling specific sediment yield in the upper Indus River basin,northern Pakistan

Estimates of sediment yield are essential in water resources analysis, modelling and engineering, in investigations of continental denudation rates, and in studies of drainage basin response to changes in climate and land use. The availability of high resolution, global environmental datasets offers an opportunity to examine the relationships between specific sediment yield (SY_sp) and drainage basin attributes in a geographical information system (GIS) environment. This study examines SY_sp at 14 long‐term gauging stations within the upper Indus River basin. Twenty‐nine environmental variables were derived from global datasets, the majority with a 1 × 1 km resolution. The SY_sp ranges from 194 to 1302 t km^?2 yr^?1 for sub‐basins ranging from 567 to 212 447 km². The high degree of scatter in SY_sp is greatly reduced when the stations are divided into three groups: upper, glacierized sub‐basins; lower, monsoon sub‐basins; and the main Indus River. Percentage snow/ice cover (LC_s) emerges as the single major land cover control for SY_sp in the high mountainous upper Indus River basin. A regression model with percentage snow/ice cover (LC_s) as the single independent variable explains 73·4% of the variance in SY_sp for the whole Indus basin. A combination of percentage snow/ice cover (LC_s), relief and climate variables explains 98·5% of the variance for the upper, glacierized sub‐basins. For the lower monsoon region, a regression model with only mean annual precipitation (P) explains 99·4% of the variance. Along the main Indus River, a regression model including just basin relief (R) explains 92·4% of the variance in SY_sp. Based on the R²_adj and P‐value statistics, the variables used are capable of explaining the majority of variance in the upper Indus River basin. Copyright © 2007 John Wiley & Sons, Ltd.     

On the frequency distribution of land slope

Slope frequency distributions are computed from Digital Elevation Models (DEMs) of 18 areas from five different geographic and physiographic regions. Over 160,000 slope estimates are used in the definition of the slope frequency distribution for each area. Analysis of several transformations of the slope data indicates that no single transformation is capable of normalizing all slope distributions. This is, in part, due to the arbitrary placement of the boundaries of the map sheets used in the data collection. The square root of sine, however, provides better results than sine, log-tangent, or several other transformations. The degree of (positive) skewness of gradients in degrees correlates with the proportion of slopes below 2° or 5° (r_s = + 0·93), as does the kurtosis (r_s = + 0·91).     

Magnitude and processes of bank erosion at a small stream in Denmark

River banks are important sources of sediment and phosphorus to fluvial systems, and the erosion processes operating on the banks are complex and change over time. This study explores the magnitude of bank erosion on a cohesive streambank within a small channelized stream and studies the various types of erosion processes taking place. Repeat field surveys of erosion pin plots were carried out during a 4‐year period and observations were supplemented by continuous monitoring of volumetric soil water content, soil temperature, ground water level and exposure of a PEEP sensor. Bank erosion rates (17·6–30·1 mm year^?1) and total P content on the banks were relatively high, which makes the bank an important source of sediment and phosphorus to the stream, and it was estimated that 0·27 kg P_tot year^?1 ha^?1 may potentially be supplied to the stream from the banks. Yearly pin erosion rates exceeding 5 cm year^?1 were mainly found at the lower parts of the bank and were associated with fluvial erosion. Negative erosion pin readings were widespread with a net advance of the bank during the monitoring period mainly attributed to subaerial processes and bank failure. It was found that dry periods characterized by low soil water content and freeze–thaw cycles during winter triggered bank failures. The great spatial variability, in combination with the temporal interaction of processes operating at different scales, requires new tools such as 3‐D topographical surveying to better capture bank erosion rates. An understanding of the processes governing bank erosion is required for riparian management using vegetational measures as root size and structure play different roles when it comes to controlling bank erosion processes. Copyright © 2010 John Wiley & Sons, Ltd.     

Influence of seepage undercutting on the stability of root‐reinforced streambanks

Several mechanisms contribute to streambank failure including fluvial toe undercutting, reduced soil shear strength by increased soil pore‐water pressure, and seepage erosion. Recent research has suggested that seepage erosion of noncohesive soil layers undercutting the banks may play an equivalent role in streambank failure to increased soil pore‐water pressure. However, this past research has primarily been limited to laboratory studies of non‐vegetated banks. The objective of this research was to utilize the Bank Stability and Toe Erosion Model (BSTEM) in order to determine the importance of seepage undercutting relative to bank shear strength, bank angle, soil pore‐water pressure, and root reinforcement. The BSTEM simulated two streambanks: Little Topashaw Creek and Goodwin Creek in northern Mississippi. Simulations included three bank angles (70° to 90°), four pore‐water pressure distributions (unsaturated, two partially saturated cases, and fully saturated), six distances of undercutting (0 to 40 cm), and 13 different vegetation conditions (root cohesions from 0·0 to 15·0 kPa). A relative sensitivity analysis suggested that BSTEM was approximately three to four times more sensitive to water table position than root cohesion or depth of seepage undercutting. Seepage undercutting becomes a prominent bank failure mechanism on unsaturated to partially saturated streambanks with root reinforcement, even with undercutting distances as small as 20 cm. Consideration of seepage undercutting is less important under conditions of partially to fully saturated soil pore‐water conditions. The distance at which instability by undercutting became equivalent to instability by increased soil pore‐water pressure decreased as root reinforcement increased, with values typically ranging between 20 and 40 cm at Little Topashaw Creek and between 20 and 55 cm at Goodwin Creek. This research depicts the baseline conditions at which seepage undercutting of vegetated streambanks needs to be considered for bank stability analyses. Copyright © 2008 John Wiley & Sons, Ltd.