Next Article in Journal
The Influence of Precision Dripping Irrigation System on the Phenology and Yield Indices of Sweet Maize Hybrids
Previous Article in Journal
The Variations and Influences of the Channel Centerline of the Zhenjiang-Yangzhou Reach of the Yangtze River Based on Archival and Contemporary Data Sets
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Freeze-Thaw Effect on Riverbank Stability

1
College of Water Conservancy and Civil Engineering, Inner Mongolia Agricultural University, Hohhot 010018, China
2
Department of Civil and Environmental Engineering, Clarkson University, Potsdam, NY 13699, USA
*
Author to whom correspondence should be addressed.
Water 2022, 14(16), 2479; https://doi.org/10.3390/w14162479
Submission received: 2 July 2022 / Revised: 26 July 2022 / Accepted: 9 August 2022 / Published: 12 August 2022
(This article belongs to the Section Water Erosion and Sediment Transport)

Abstract

:
The stability of riverbanks in cold regions is affected by the freeze-thaw action. The freeze-thaw process causes changes in the moisture content, friction angle, and cohesiveness of the bank material. Together with the freeze-thaw effect, seepage pressure influenced by the changing water levels, and the bank slope are the key factors affecting bank stability. A limit equilibrium bank stability model considering the infiltration water pressure effect is developed and applied to the Shisifenzi section of Inner Mongolia reach of the Yellow River. Laboratory tests of field samples with moisture contents of 15%, 18%, 21%, and saturated showed that the freeze-thaw action reduced the degree of saturation by 34.37 %, 30.71%, 32.48%, and 46.23%, respectively, accompanied by reductions in the internal friction angles by 1.78%, 2.74%, 6.33%, and 5.32%. These changes resulted in a 24.35% to 29.13% reduction in the safety factor of bank stability. Together with seasonal variations in the water levels the field data showed that the bank stability safety factor in the study site increases gradually through the melting period, dry period, wet period, flooding period, and low flow period. The slope stability safety factor increases with the stage of the river but decreases with the groundwater level.

1. Introduction

Bank stability is an important phenomenon affecting river channel morphology. For the Inner Mongolia reach of the Yellow River, the stability of riverbanks is essential for flood control and bank protection projects. Bank erosion also contributes to sediment supply to the river. The bank stability is mainly influenced by channel flow and the composition and properties of the riverbank material [1,2,3,4,5]. Erosion at the bank toe could lead to bank instability. Under the action of gravity, cracks will develop along the banks, especially when there is a significant difference between the groundwater level and water level in the channel due to the difference in the infiltration water pressures [6]. These cracks could be further enhanced by heavy rainfall and freeze-thaw action [7].
Bank stability is influenced by numerous factors [8,9]. It is mainly affected by the channel flow condition, embankment groundwater level, and sediment transport condition [3,10,11,12]. They are mostly manifested in the scouring and seepage effects, resulting in gravity failure of the bank [13,14]. The composition and properties of riparian soils are internal factors affecting the stability of riverbanks [4,15,16,17,18,19]. The type of bank instability may vary with the channel and cross-sectional profiles and bank material properties [20,21,22]. Regarding the slope and material properties, the more significant the proportion of cohesive components in riparian soils, the stronger their ability to withstand scouring [18]. However, soil cohesion and shear strength are affected by water content and freeze-thaw cycling [4,15,17,23]. The former affects the pore water pressure and matrix suction within the soil, thereby changing the scour resistance of the bank [24]. The latter alters the soil’s mechanical properties, i.e., cohesion and internal friction angle, by weakening its internal structure [23,25].
Bank instability is often the result of multiple factors, and the pattern of bank instability can differ depending on the bank and hydraulic conditions [4,11,26,27,28]. Osman and Thorne [11] proposed the classical steep slope plane instability model for homogeneous clayey soil banks, which considers the influence of scouring on the bank boundary. Darby and Thorne [26] improved the Osman and Thorne [11] bank stability model by considering the effects of hydrostatic and pore water pressure in the force analysis of the damaged surface. In addition, other influencing factors, e.g., matrix suction, groundwater level, soil properties, etc., and dynamic water pressure have been analyzed [15,17,29,30,31]. For complex non-homogeneous multilayered soil banks, Amiri–Tokaldany [32,33] proposed a bank stability model for non-single soil layers based on Dandy and Thorne [26]. Lai [34] coupled bank stability analysis with a two-dimensional streambed dynamics model to analyze and predict the collapse of multilayered clayey soil bank slopes. This study will examine the stability of seasonal-frost banks with coupled hydraulic effects of channel water level and the groundwater level in the Shisifenzi Bend section of the Inner Mongolia reach of the Yellow River, as shown in Figure 1.
The weather of the study area is controlled by the Mongolian high pressure and has a typical continental monsoon climate. The annual average temperature is 7.3 °C, the highest temperature is 38.4 °C and the lowest temperature reaches −36.3 °C. The average annual rainfall is 359.8 mm and the evaporation is 1849 mm. The freezing and thawing period is 4 months long, with the temperature turning negative in mid-November and positive in mid-March of the following year. The Inner Mongolia reach of the Yellow River has very gentle slopes. The bed slope is nearly 1/10,000, with many bends and broad, shallow cross-sections [35]. The riverbanks are mainly formed from the deposition of sediment from upstream and are prone to instability and collapse. Moreover, since it is situated in the seasonal frost zone, bank failures are more frequent due to the combined effect of hydraulic erosion and the freeze-thaw impact. This paper extends the Darby and Thorne bank stability model [26] by introducing the effect of infiltration water pressure. The model is applied to the Shisifenzi Bend for different river flow and groundwater conditions through the freeze-thaw period and other seasons of the year. Figure 2 shows an example of the bank failures in the Shisifenzi Bend section.

2. Slope Stability Analysis and Modeling

2.1. Bank Stability Mechanics

The critical height and limit equilibrium methods are commonly used in bank stability analysis [36]. The former applies to steep and vertical bank slopes, while the latter accounts for the actual geometric form of bank slopes in the stability analysis. The limit equilibrium method considers the stability safety factor, F S , of the bank slope. When F S is below the critical value, it is assumed that the bank will be destabilized and collapse. The safety factor can be calculated as:
F S = F R F D
In which, F D is the sum of the component of all forces acting on the potential sliding mass along the sliding surface; F R is the resistance force on the potential sliding surface, which is a function of the cohesive force and frictional force. Xia [5] suggested that the bank remains stable when F S > 1.3; and the bank is unstable when F S < 1.0. When 1.0 ≤ F S ≤ 1.3, the bank is at risk of possible failure. The threshold of instability safety factor is set to 1.3 in the present study.

Forces on the Potential Sliding Mass

The bank is subjected to the gravity force, the lateral water pressure from the channel water and the pore water pressure from the bank. The pore water pressure is generated by the groundwater level inside the bank, which lags behind the change in channel water level. Figure 3 shows a force diagram of riverbank, in which W is the gravity force, kN/m2; F c p is the lateral hydrostatic pressure, kN/m2; F x and F z are the horizontal and vertical components of F c p , respectively; ω is the angle between F c p and the horizontal direction, °; α is the angle between F c p and the potential sliding surface, °; U w is the pore water pressure, kN/m2; P d is the infiltration water pressure, kN/m2; e is the angle between the infiltration water pressure and the normal direction of the sliding surface, °; h w is the groundwater level, m; ξ is the river water level, m; β is the inclination of the potential sliding surface, °; i is the slope inclination, °; b is the top width of the potential sliding mass, m; H is the potential sliding depth, m; H 1 is the depth of the bank slope, m; H 2 is the depth of the bank tensile fracture, m; H 3 is the depth of the residual tensile fracture, m.
Gravity: Based on the channel water level and the groundwater level elevations, the sliding mass is divided into two parts by the groundwater table: above it is the unsaturated soil, area A 1 , m 2 ; below it is the saturated soil, area A 2 , m 2 . The gravity is calculated as:
W = γ A 1 + γ s a t A 2
In which, γ is the specific weight of the soil, kN/m3; γ s a t is the specific weight of saturated soil, kN/m3. and
A 1 = 1 2 H 2 2 + 2 H h w h w 2 tan β + H 3 2 + 2 H 1 h w + h w 2 tani A 2 = 1 2 [ H h w 2 tan β H 1 h w 2 tani ]
Lateral hydrostatic pressure: The hydrostatic pressure under the water surface is P = γ w h , where, P is the hydrostatic pressure, kN/m2; γ w is the specific weight water and is set to 10 kN/m3 consider the turbidity of the river water; h is the depth from the water surface, m. The hydrostatic pressure can be decomposed into horizontal hydrostatic pressure F x and vertical hydrostatic pressure F z , then the total hydrostatic pressure is:
F c p = F x 2 + F z 2
F x and F z are calculated by:
F x = 1 2 γ w h 2
F z = γ w S
ω = tan 1 F z F x
where S is the area of the water body above the sliding mass, m 2 .
Pore water pressure: The pore water pressure at any point on the sliding surface is U w = γ w h w + ξ . The total pore water pressure acting on the potential sliding body is:
U w = 0 x u w d x
Infiltration water pressure: The infiltration water pressure can be calculated by:
P d = γ w J A 2
where J is the infiltration gradient, J = Δ h l ; l is the seepage path length, m; Δ h is the head difference, m.
Sliding surface angle: The sliding surface angle, β , of the collapsed body can be estimated by Equation (10) [37]:
β = i + φ 2
where φ is the internal friction angle of the slope soil, °.

2.2. Modeling Bank Slope Instability

Depending on the interaction between the channel water and the bank, the bank failure mode varies, including plane sliding, arc sliding, and slumping. Specifically, arc sliding failure occurs when the bank slope is gentle, plane sliding failure occurs when the bank slope is steep, and slumping collapse occurs when the clay layer above a binary structure bank is suspended by erosion below [4,27,38,39]. Since the bank slope of the study reach is steep, and the soil is mainly low liquid limit powder, the bank can be considered monolithic. The main purpose of this paper is to analyze the stability of the riverbank and to facilitate the calculation according to the plane sliding analysis. Osman and Thorne [11] divided riverbank collapse into initial and secondary collapse based on plane sliding. The initial failure is the scouring of the bank toe, and the secondary collapse refers to the parallel retreating process of the riverbank and the sliding surface passing through the bank toe. Based on this assumption, they proposed a model for calculating the stability of the riverbank in the case of plane failure. However, the model failed to consider lateral water pressure, pore water pressure, infiltration water pressure, and the damaged surface must pass through the slope foot. Darby and Thorne [26] improved the Osman and Thorne [11] model so that the sliding surface is not necessary to pass through the bank toe, and the effects of lateral water pressure and pore water pressure in the force analysis of the sliding surface are considered. The established stability calculation model for riverbanks with clayey soils is:
F S = F R F D = C L + W U W cos β + F C P cos α tan φ W sin β F C P sin α
In which, C is the effective cohesive force of soil, kPa ; φ is the effective internal friction angle of soil, °.
Although the lateral and pore water pressures are considered in the model of Darby and Thorne [26], the model does not consider the role of infiltration water pressure. The infiltration water pressure acting on the sliding body is an important factor affecting the stability of the river embankment slope, especially during the freeze-thaw period when the channel water level declines and during the summer flood period when the river water level varies greatly.
Since the sliding force F D is a combination of the gravitational force, hydrostatic pressure, and infiltration water pressure acting on the potential sliding body along the direction of the sliding surface. Equation (12) modifies Darby and Thorne [26] formulation by including the infiltration water pressure in the sliding force calculation:
F D = W sin β F C P sin α ± P d sin e
The infiltration water pressure is negative when the river water level is higher than the groundwater level and positive when the groundwater level is higher.
By taking into account the effects of gravity, pore water pressure, lateral hydrostatic pressure, and infiltration water pressure, the slip resistance of the potential sliding surface is given by:
F R = C L + W U W cos β + F C P cos α ± P d cos e tan φ
Substituting Equations (12) and (13) into Equation (1) gives bank stability factor, F S :
F S = F R F D = C L + W U W cos β + F C P cos α ± P d cos e tan φ W sin β F C P sin α ± P d sin e ,

3. Model Application and Analysis

3.1. Study Site

The modified bank stability model is used to analyze the bank stability in the Shisifenzi Bend of the Yellow River in Inner Mongolia, located at 40°17′39″ N and 111°2′53″ E in Tuoketuo County, 4 km upstream of the Toudaoguai hydrological station. as shown in Figure 1. The river channel turns sharply from northwest to southeast. The channel slope is about 0.1%. The channel width varies between 200 m and 500 m.
Field samples indicate the bank material mainly consists of sandy loam and shallow clay. Table 1 lists the basic physical and mechanical parameters of the soils obtained by laboratory tests.

3.2. Bank Stability Analysis

The slope of the riverbank in the study area varies between 30° and 60°. An average value of 45° is adopted for calculation. Figure 4 shows the water level in the study reach changes rapidly over a wide range. The lowest water level was 987.14 m on May 7, while the highest was 990.65 m on September 6, with a variation of 3.51 m. The top elevation of the bank is about 992.0 m, and the height is 8 m. Hence the riverbank is modeled with H = 8   m , H 1 = 7   m , H 2 = H 3 = 2   m , and i = 45 o , as shown in Figure 5. Seven bank groundwater levels varying from 0 m to 3.0 m are used to calculate the safety coefficient of bank stability, with the river water level ranging from 1 to 8 m from the top of the bank under each working condition.
Figure 6 shows the variation of bank stability factors under different combinations of groundwater and river water levels. Table 2 shows the stability of the riverbank under different combinations of conditions and the critical stability safety factor of 1.3. From the calculated results shown in Figure 6 and Table 2, it can be seen that: the stability safety factor of the bank slope decreases rapidly with the lowering of the river water level, and finally tends to be stable. When the stability safety factor of bank slope is less than 1.3, the bank slope will be destabilized. The rise and fall of river water level will cause the change in groundwater level in the bank, which makes the water content, the cohesion, and the internal friction angle of the bank soil change.
The reach studied has a freezing period of nearly four months, during which the river water level changes significantly under the influence of ice conditions. In the meantime, the bank freeze and thaw affect the mechanical properties of the bank material. To determine the seasonal variation of the slope stability, according to the changes in temperature and river water level, the flow period is divided into the thawing period (early March–end of March), dry period (early April–end of June), wet period (early July–early September), flooding period (mid–September), and low flow period (late September–mid-November). The safety coefficient of the bank stability in each period is calculated. The calculated result is presented in Table 3.
As shown in Table 3, the safety factor increases from the thawing period to the high-water periods and decreases during the receding water period. The safety coefficient during the thawing period is only 0.75, far smaller than the threshold of 1.3, indicating a period of potential destabilization of the bank slope. The main reasons are: (1) After the freeze-thaw cycle, the cohesion and internal friction angle of the bank material is reduced, and the shear strength is weakened; (2) During this period, the river water level drops sharply after the melting of ice cover. The lateral water pressure decreases, while the infiltration water pressure towards the river rises. Under the combined action of these two, the bank slope is destabilized; (3) During the dry period, the soil is still impacted by the freeze-thaw effect. However, the difference between the groundwater and river water levels decreases, so the infiltration water pressure effect is weakened compared to the thawing period. Hence, the safety coefficient of the bank increased, but it is still lower than the 1.3 threshold, and the bank slope is still subjected to destabilization; (4) When the river water level rises during the wet water period, the lateral hydrostatic pressure increases. Although the change in groundwater level in the bank lags behind, and there is infiltration water pressure towards the bank. Hence, the safety factor is the highest; (5) During the flooding period, the groundwater level of the bank rises and is slightly higher than the river water level. The infiltration water pressure is small, and the safety factor is smaller than the wet water period. However, the bank slope is stable in this period; (6) The river water level and the lateral water pressure drop during the low water period. However, the groundwater level of the bank is still high, and there is a large infiltration water pressure toward the river. In this case, the safety factor of the bank slope is reduced to 1.29. These indicate that the bank slope may still be susceptible to instability during the low water period, which partly explains why the sand content of the river flow is still high in September and October after the flooding period.
Through the statistical analysis of the monthly average sediment content of the River at the Shizuishan and Toudaoguai hydrological stations from 2010 to 2019, as shown in Figure 7, it can be seen that the sand content of the Shizuishan section is higher than that of the Toudaoguai section during the freezing period (December–February). This means that the river channel is subjected to siltation during the freezing period. However, the sediment content of the Toudaoguai section is significantly higher during the open water period (March–November). In this period, the sediment supply is mainly from tributaries, i.e., the Ten Kungdui, Kundulun River, Wudangou, etc., and the collapsed bank material. During the summer flood period (June–September), the sand is mainly from the tributaries during flash floods. The spring flood period (March–April) is in the thawing and dry period, and the stability and safety coefficient of the bank slope is the smallest. The sediment mainly comes from the bank collapse under the effect of freeze-thaw erosion. From Figure 7, it can be seen that the average sediment content of the Toudaoguai section reaches 2.5 kg/m³ in March–April, the highest sediment content period second only to the summer flood season. Therefore, it can be concluded that the effect of freezing and thawing on the collapse and destabilization of the bank slope and on the sediment content in the river channel cannot be ignored.

3.3. Factors Affecting Bank Stability

The bank stability is governed by hydraulic and bank conditions and the bank condition. The hydraulic condition is mainly caused by the change in water level inside and outside the river channel. The bank condition includes the bank slope, soil composition, and physical and mechanical properties. The river water level changes a lot during the flood and summer flood periods, but the groundwater level of the bank lags behind, which produces the infiltration water pressure on the bank. Moreover, the freezing and thawing effect of the flood will change the soil properties. The effects of freeze-thaw action, infiltration water pressure, and bank slope on bank stability are discussed in this section.

3.3.1. Freeze-Thaw Effect

The Shisifenzi Bend has long, cold winters, with the air temperature remaining below freezing from mid-November to mid-March. The freeze-thaw period lasts for more than four months. The frozen zone surface will turn into a water-rich layer when it melts [41]. Hence, the groundwater level is set at 0 m at the end of the freeze-thaw period. The stability of the embankment soil before and after the freezing and thawing is calculated. Remodeled field soil samples with water contents of 15%, 18%, 21%, and fully saturated samples are used to investigate the effect of freeze-thaw action on their mechanical properties. To fully freeze-thaw the samples, they were frozen at −15 °C for 12 h then thawed at 15 °C for another 12 h. The number of freeze-thaw cycles was selected as 0, 1, 3, 5, 7, 9, 11, 13, and 15. The changes in cohesion and internal friction angle of the samples with varying moisture contents after different freeze-thaw cycles were obtained using the static triaxle test with and Mohr–Coulomb failure criterion, as shown in Figure 8. The cohesion c and internal friction angle φ of soil samples with each moisture content became stable after 15 freeze-thaw cycles. Figure 8a shows that the cohesion c of each moisture content soil sample tends to level off after 15 cycles and decreased by 34.37%, 30.71%, 32.48%, and 46.23%, respectively. Correspondingly, the internal friction angle φ declined by 1.78%, 2.74%, 6.33%, and 5.32%, respectively, as shown in Figure 8b.
The water content of the undisturbed soil sample measured on site is about 21%. The samples were taken on 15 July 2020, when the river water was rising. The effective cohesion and effective internal friction angle of soil samples with 21% moisture content after 0 and 15 freeze-thaw cycles are used to calculate the bank stability safety factor before and after freeze-thaw, as shown in Figure 9. It can be seen that the stability safety factor of the bank after freeze-thaw is significantly lower.

3.3.2. Effect of Infiltration Water Pressure

To study the effect of the infiltration water pressure, the safety factor of the bank stability with or without considering the infiltration water pressure when the groundwater level is 1 m, 2 m, and 3 m, and the river water level varies from 1 to 8 m from the top of the bank is calculated. The result is shown in Figure 10. When the groundwater level is lower than the river level, the infiltration water pressure is toward the bank, which inhibits the sliding of the slope. The bank stability safety factor is larger than when the infiltration water pressure is not considered. On the other hand, when the groundwater level is above the river water level, the infiltration water pressure is towards the river. As a result, the sliding force on the sliding mass increases, which promotes the destabilization of the slope. Therefore, the safety factor when infiltration water pressure is considered is smaller.

3.3.3. Effect of Bank Slope

To examine the effect of bank slope, the bank stability factor is calculated for the slope angles of 30°, 35°, 40°, 45°, 50°, 55°, and 60°, respectively, and the groundwater level is 1 m. As shown in Figure 11, the bank stability factor varies significantly with the bank slope. When the river water level is lower than the groundwater level, the infiltration water pressure is generated, which promotes the destabilization of the bank slope. When the river water level is below 4 m, there is a possibility of bank collapse under different bank slope angles. In addition, the safety factor of the bank decreases gradually as the bank slope increases.

4. Conclusions

This paper extended the Darby and Thorne [26] bank stability model by considering the effect of infiltration water pressure, which is important for rivers subjected to freeze-thaw actions in the banks. The model is used to analyze the bank slope stability of the Shisifenzi Bend of the Yellow River in Inner Mongolia. The effects of freeze-thaw action, infiltration water pressure, and bank slope are studied. The major conclusions are:
(1)
When the river water level is constant, the safety factor of bank stability decreases with the rising groundwater level. When the groundwater level is constant, the safety coefficient of bank stability declines with the decreasing river water level, with a trend of thawing period < dry period < low water period < flooding period < wet period.
(2)
The freeze-thaw action significantly changes the mechanical properties of the bank material, leading to a 24.35–29.13% reduction in the safety factor of bank stability, indicating the important effect of the freeze-thaw in reducing the bank stability.
(3)
When the groundwater level is lower than the river water level, the infiltration water pressure mainly shows an inhibitory effect on the bank stability. The safety factor of bank stability will increase. By contrast, when the groundwater level is higher, the infiltration water pressure plays a destabilizing role.
(4)
When the river water level is below 4 m, there is a possibility of bank collapse under different bank slope angles, and the safety factor decreases with larger bank slope angles.

Author Contributions

Conceptualization, C.L. and Z.Y.; consultation, H.T.S.; methodology, C.L. and Z.Y.; model development, Z.Y. and C.L.; field and laboratory experiments, C.L. and Z.Y.; writing, Z.Y., C.L. and H.T.S.; funding acquisition, C.L. and X.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by [National Natural Science Fund of China] grant number [No.51969025], [Natural Science Foundation of Inner Mongolia Autonomous Region] grant number [No.2019MS05006] And [Key project of Inner Mongolia Natural Science Foundation] grant number [No.2022ZD08].

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author or the first author.

Acknowledgments

Special thanks to Xinchuan Lu of the Yellow River Wanjiazhai Hydraulic Hub Co., Ltd. for providing the water level data of the Shisifenzi reach.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Hooke, J.M. An analysis of the processes of river bank erosion. J. Hydrol. 1979, 42, 39–62. [Google Scholar] [CrossRef]
  2. Chen, D.; Duan, J.G. Modeling width adjustment in meandering channels. J. Hydrol. 2006, 321, 59–76. [Google Scholar] [CrossRef]
  3. Henshaw, A.J.; Thorne, C.R.; Clifford, N.J. Identifying causes and controls of river bank erosion in a British upland catchment. Catena 2012, 100, 107–119. [Google Scholar] [CrossRef]
  4. Xia, J.Q.; Zong, Q.L.; Deng, S.S.; Xu, Q.X.; Lu, J.Y. Seasonal variations in composite riverbank stability in the Lower Jingjiang Reach, China. J. Hydrol. 2014, 519, 3664–3673. [Google Scholar] [CrossRef]
  5. Xia, J.Q.; Zong, Q.L. Mechanism of Bank Failure in Jingjiang Section of Yangtze River and Its Numerical Simulation; Science Press: Beijing, China, 2015; pp. 155–167. (In Chinese) [Google Scholar]
  6. Lu, Y.; Lu, Y.J.; Zhang, X.N. Influence of river water level landing on slope stability. Adv. Water Science. 2008, 19, 389–393. (In Chinese) [Google Scholar]
  7. Guo, Y.; Shen, W. Monitoring and Experiment on the Effect of Freeze-Thaw on Soil Cutting Slope Stability. Procedia Environ. Sci. 2011, 10, 1115–1121. [Google Scholar] [CrossRef]
  8. Xie, H.; Yu, K.Z. Analysis of factors affecting bank stability in construction of Yangtze River dry embankment protection project. People’s Yangtze River. 2006, 37, 93–94+101. (In Chinese) [Google Scholar]
  9. Hoomehr, S.; Akinola, A.; Wynn-Thompson, T.; Garnand, W.; Eick, M.J. Water Temperature, pH, and Road Salt Impacts on the Fluvial Erosion of Cohesive Streambanks. Water 2018, 10, 302. [Google Scholar] [CrossRef]
  10. Xu, W.; He, S.L. Analysis of influencing factors of riverbank collapse. Regul. Huai River. 2019, 4, 22–24. (In Chinese) [Google Scholar]
  11. Osman, A.M.; Thorne, C.R.; Member, A. Riverbank stability analysis. I: Theory. J. Hydraul. Eng. 1988, 114, 134–150. [Google Scholar] [CrossRef]
  12. Zhang, X.N.; Jiang, C.F.; Ying, Q.; Chen, C.Y. Review of research on bank collapse in natural rivers. Adv. Sci. Technol. Water Resour. 2008, 28, 80–83. (In Chinese) [Google Scholar]
  13. Zhang, X.N.; Jiang, C.F.; Chen, C.Y.; Ying, Q. Analysis of influencing factors of river bank failure. J. River Sea Univ. (Nat. Sci. Ed.) 2009, 37, 36–40. (In Chinese) [Google Scholar]
  14. Zhang, F.Z.; Chen, X.P. Study on the effect of river scour on seepage and deformation of embankment. Geotechnics 2011, 32, 441–447. (In Chinese) [Google Scholar]
  15. Simon, A.; Curini, A.; Darby, S.E.; Langendoen, E.J. Bank and near-bank processes in an incised channel. Geomorphology 2000, 35, 193–217. [Google Scholar] [CrossRef]
  16. Couper, P. Effects of silt–clay content on the susceptibility of river banks to subaerial erosion. Geomorphology 2003, 56, 95–108. [Google Scholar] [CrossRef]
  17. Chiang, S.W.; Tsai, T.L.; Yang, J.H. Conjunction effect of stream water level and groundwater flow for riverbank stability analysis. Environ. Earth Sci. 2011, 62, 707–715. [Google Scholar] [CrossRef]
  18. Karmaker, T.; Dutta, S. Erodibility of fine soil from the composite river bank of Brahmaputra in India. Hydrol. Processes. 2011, 25, 104–111. [Google Scholar] [CrossRef]
  19. Kimiaghalam, N.; Goharrokhi, M.; Clark, S.P.; Ahmari, H. A comprehensive fluvial geomorphology study of riverbank erosion on the Red River in Winnipeg, Manitoba, Canada. J. Hydrol. 2015, 529, 1488–1498. [Google Scholar] [CrossRef]
  20. Thoron, C.R.; Tovey, N.K. Stability of composite riverbanks. Earth Surf. Processes Landf. 1981, 6, 469–484. [Google Scholar] [CrossRef]
  21. Wang, Y.G.; Kuang, S.F. Study of critical collapse height of river banks. J. Water Resour. 2007, 10, 1158–1165. (In Chinese) [Google Scholar]
  22. Wang, J.; Zong, Q.L.; Yue, H.Y.; Liu, Z.X. Effects of alternating wet and dry on the mechanical properties of shoreland soils in a typical section of the Jingjiang River section of the Yangtze River. J. Agric. Eng. 2019, 35, 144–152. (In Chinese) [Google Scholar]
  23. Wang, X.D.; Li, S.Y.; Sun, Y.F.; Zhang, C.B.; Liu, G.W. Influence of freeze–thaw cycling on the soil mechanical properties of open-pit mine dump under different moisture contents. Environ. Earth Sci. 2021, 80, 279. [Google Scholar] [CrossRef]
  24. Simon, A.; Collison, A.J.C. Pore-water pressure effects on the detachment of cohesive streambeds: Seepage forces and matric suction. Earth Surf. Processes Landf. 2001, 26, 1421–1442. [Google Scholar] [CrossRef]
  25. Tang, L.Y.; Li, G.; Li, Z.; Jin, L.; Yang, G.S. Shear properties and pore structure characteristics of soil–rock mixture under freeze–thaw cycles. Bull. Eng. Geol. Environ. 2021, 80, 3233–3249. [Google Scholar] [CrossRef]
  26. Darby, S.E.; Thorne, C.R. Development and testing of riverbank-stability analysis. J. Hydraul. Eng. 1996, 122, 443–454. [Google Scholar] [CrossRef]
  27. Samadi, A.; Davoudi, M.H.; Amiri-Tokaldany, E. Experimental study of cantilever failure in the upper part of cohesive riverbanks. Res. J. Environ. Sci. 2011, 5, 444–460. [Google Scholar] [CrossRef]
  28. Wang, Y.G.; Kuang, S.F. Research on bank collapse types and collapse modes. Sediment Res. 2014, 1, 13–20. [Google Scholar]
  29. Rinaldi, M.; Casagli, N. Stability of streambanks formed in partially saturated soils and effect of negative pore water pressures: The Sieve River (Italy). Geomorphology 1999, 26, 253–277. [Google Scholar] [CrossRef]
  30. Deng, S.; Xia, J.Q.; Zhou, M.R.; Li, J.; Zhu, Y.H. Coupled modeling of bank retreat processes in the Upper Jingjiang Reach, China. Earth Surf. Processes Landf. 2018, 43, 2863–2875. [Google Scholar] [CrossRef]
  31. Zhang, Z.H.; Wang, W.D.; Zhang, Y.B. Improvement of transfer coefficient method under the consideration of hydrodynamic pressure. Eng. Fail. Anal. 2021, 124, 1167–1170. [Google Scholar] [CrossRef]
  32. Amiri-tokaldany, E.; Darby, S.E.; Tosswell, P. Bank stability analysis for predicting reach-scale land loss and sediment yield. J. Am. Water Resour. Assoc. 2003, 39, 897–909. [Google Scholar] [CrossRef]
  33. Amiri-tokaldany, E.; Darby, S.E.; Tosswell, P. Coupling bank stability and bed deformation models to predict equilibrium bed topography in river bends. J. Hydraul. Eng. 2007, 113, 1167–1170. [Google Scholar] [CrossRef]
  34. Lai, Y.G.; Thomas, R.E.; Ozeren, Y.; Simon, A.; Greimann, B.P.; Wu, K. Modeling of multilayer cohesive bank erosion with a coupled bank stability and mobile-bed model. Geomorphology 2015, 243, 116–129. [Google Scholar] [CrossRef]
  35. Li, C.; Quan, D.; Zhang, Y.; Shi, X.H.; Guo, Z.Y. Characteristics of water-sand movement process and evolution trend of the Yellow River (Inner Mongolia section). J. Soil Water Conserv. 2020, 34, 41–46+53. (In Chinese) [Google Scholar]
  36. Li, Y.; Su, H.F.; Wang, L. Advances in shore slope stability research. Constr. Technol. 2017, 46, 1225–1229. (In Chinese) [Google Scholar]
  37. Simon, A. Mass Wasting Algorithms in an Alluvial Channel Model. In Proceedings of the 5th Federal Inter-Agency Sedimentation Conference, Las Vegas, NV, USA, 22–29 August 1991. [Google Scholar]
  38. Zhang, X.N.; Jiang, C.F.; Chen, C.Y.; Ying, Q. Types and characteristics of river bank failure. Prog. Water Resour. Hydropower Sci. Technol. 2008, 28, 66–70. (In Chinese) [Google Scholar]
  39. Xia, J.Q.; Zong, Q.L.; Xu, Q.X.; Deng, C.Y. Characteristics of binary structure riverbank soils and bank failure mechanism in the lower Jingjiang River. Adv. Water Sci. 2013, 24, 810–820. (In Chinese) [Google Scholar]
  40. Ran, R.; Liu, Y.F. Analysis of the influence of reservoir bank slope morphology on its stability using BSTEM model. Groundwater 2011, 33, 162–165. (In Chinese) [Google Scholar]
  41. Yang, Z.; Li, C.; Ji, H.L.; Sun, C.; Shi, Z.F. Study on the characteristics of temperature and humidity changes in the embankment soil of the Yellow River during the freeze-thaw period and the influencing factors. Arid. Zone Resour. Environ. 2021, 35, 76–83. (In Chinese) [Google Scholar]
Figure 1. Inner Mongolia reach of the Yellow River, and the Shisifenzi Bend section.
Figure 1. Inner Mongolia reach of the Yellow River, and the Shisifenzi Bend section.
Water 14 02479 g001
Figure 2. Yellow Riverbank at Shisifenzi, Inner Mongolia.
Figure 2. Yellow Riverbank at Shisifenzi, Inner Mongolia.
Water 14 02479 g002
Figure 3. Schematic diagram of forces on the sliding surface of the riverbank (after Darby and Thorne with modification.).
Figure 3. Schematic diagram of forces on the sliding surface of the riverbank (after Darby and Thorne with modification.).
Water 14 02479 g003
Figure 4. Daily average water level at the Shisifenzi in 2020.
Figure 4. Daily average water level at the Shisifenzi in 2020.
Water 14 02479 g004
Figure 5. Calculated bank slope profile.
Figure 5. Calculated bank slope profile.
Water 14 02479 g005
Figure 6. Variation of bank stability factor with groundwater levels and river water levels.
Figure 6. Variation of bank stability factor with groundwater levels and river water levels.
Water 14 02479 g006
Figure 7. Monthly average sediment content from 2010–2019.
Figure 7. Monthly average sediment content from 2010–2019.
Water 14 02479 g007
Figure 8. Effect of freezing and thawing on soil cohesion and angle of internal friction. (a) Soil cohesion, (b) Angle of internal friction.
Figure 8. Effect of freezing and thawing on soil cohesion and angle of internal friction. (a) Soil cohesion, (b) Angle of internal friction.
Water 14 02479 g008
Figure 9. Effect of freezing and thawing on bank stability.
Figure 9. Effect of freezing and thawing on bank stability.
Water 14 02479 g009
Figure 10. Effect of infiltration water pressure on bank stability.
Figure 10. Effect of infiltration water pressure on bank stability.
Water 14 02479 g010
Figure 11. Effect of bank inclination on bank stability.
Figure 11. Effect of bank inclination on bank stability.
Water 14 02479 g011
Table 1. Basic physical and mechanical parameters of bank material.
Table 1. Basic physical and mechanical parameters of bank material.
Natural Specific Weight, (kN/m3)Saturated Specific Weight, (kN/m3) Median   Particle   Size ,   D 50 / ( μ m ) Dry Density, (g/cm3)Experimental ValuesValid Values
C ,   ( k P a ) φ C ,   ( k P a ) φ
18.119.333.551.5835.53024.8527
Note: The effective cohesive force C′ in the Table is 0.7 times the experimental value; the effective internal friction angle φ′ is 0.9 times the experimental value [40].
Table 2. Bank slope stability under different combinations of groundwater level and river water level.
Table 2. Bank slope stability under different combinations of groundwater level and river water level.
Groundwater   Level   h w / m 0−0.5−1−1.5−2.0−2.5−3.0
River   Water   Level   ξ / m
−1+++++++
−2+++++++
−3-++++++
−4--+++++
−5---++++
−6----+++
−7----+++
−8-----++
Note: + and - indicates stable and unstable conditions.
Table 3. Bank stability factors for different periods of the year.
Table 3. Bank stability factors for different periods of the year.
ParametersThawing
Period
Dry
Period
Wet
Period
Flooding PeriodLow Flow Period
Cohesion   C / k P a 16.816.824.8524.8524.85
Friction   angle   φ 25.225.2272727
Ground   water   level   h w /m0−3−3−0.5−0.5
River   water   level   ξ /m−5−5−1.5−1.5−4
Safety factor0.751.264.162.141.29
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, C.; Yang, Z.; Shen, H.T.; Mou, X. Freeze-Thaw Effect on Riverbank Stability. Water 2022, 14, 2479. https://doi.org/10.3390/w14162479

AMA Style

Li C, Yang Z, Shen HT, Mou X. Freeze-Thaw Effect on Riverbank Stability. Water. 2022; 14(16):2479. https://doi.org/10.3390/w14162479

Chicago/Turabian Style

Li, Chao, Zhen Yang, Hung Tao Shen, and Xianyou Mou. 2022. "Freeze-Thaw Effect on Riverbank Stability" Water 14, no. 16: 2479. https://doi.org/10.3390/w14162479

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop