Impacts of the Degraded Alpine Swamp Meadow on Tensile Strength of Riverbank: A Case Study of the Upper Yellow River

: In the meandering riverbank of the Upper Yellow River (UYR), the native alpine swamp meadow (AS) has continuously degenerated into an alpine meadow (AM) due to climate change and intensiﬁed grazing. Its implication on river morphology is still not well known. This study examined this e ﬀ ect by in situ measurings of (1) physical properties of roots and their distribution in the soil-root mixture of the upper bank layer, and (2) the tensile strength in terms of excavating tests for triggering cantilever collapses of AS and AM riverbanks. The results showed that the root number in AS was signiﬁcantly greater than that in AM, though the root distribution in both was similar. Also, the average tensile strength of individual roots in AS was 31,310 kPa, while that in AM was only 16,155 kPa. For the soil-root mixture, it decreased from 67.39 to 21.96 kPa. The weakened mechanical property was mainly ascribed to the lessened root number and the simpler root structure in the soil-root mixture of AM that reduces its ability to resist the external force. These ﬁndings conﬁrmed that healthy AS can enhance bank stability and delay the development of tensile cracks in the riverbank of the meandering rivers in the UYR. Our study provides ﬁrsthand evidence of weakened bank strength in the soil-root mixture due to degradation of vegetation and will serve as a benchmark for future investigation of riverbank strength, not only in the UYR, but also in other alpine regions in the world.


Introduction
Banks of meandering rivers are often composed of silts and sand that have significantly higher compressive strength than tensile strength or cohesion [1][2][3][4]. There has been a consensus that riparian vegetation can reinforce bank strength [5][6][7][8][9] and stability through the interaction between soil and roots [10][11][12][13][14][15][16]. Mechanisms of riverbank failure are quite different from those found on hillslopes because steeper and shorter riverbanks tend to have a more variable profile and relatively small size of the failed block [17]. Accordingly, vegetation types and their root distributions throughout the bank profile play a critical role in resisting riverbank failure. In addition, the lower soil layer of the composite riverbank is subject to fluvial erosion, often resulting in a cantilever upper layer that includes the mixture of cohesive soil and vegetation. Once the gravity moment generated by the upper cantilever Grassland is the major type of land cover in the region. It is dominated by alpine swamp meadow (AS), which provides important ecosystem services to the regional environment. Affected by global warming and intensified grazing, AS in the region has gradually degraded to alpine meadows (AM) and alpine steppe meadows (ASM) [40]. Spatially, AS and AM are often seen along riverbanks, the Lanmucuo River, while ASM is mostly distributed on the upper sloping parts of the Grassland is the major type of land cover in the region. It is dominated by alpine swamp meadow (AS), which provides important ecosystem services to the regional environment. Affected by global warming and intensified grazing, AS in the region has gradually degraded to alpine meadows (AM) and alpine steppe meadows (ASM) [40]. Spatially, AS and AM are often seen along riverbanks, the Lanmucuo River, while ASM is mostly distributed on the upper sloping parts of the piedmonts (Figure 1d). Because the area covered by ASM is far away from river banks, it is not a concern in this study. Based on a set of qualitative and semi-qualitative indicators used for pasture degradation classification [33,35], a survey for riparian vegetation along a 22 km reach of the upper Lanmucuo River was conducted in the 2017-2019 period [16,41]. During the survey, the coverage and number of species, dominant species, and underground biomass of AS and AM were determined within each sampling plot with the size of 1 × 1 m. The underground biomass was represented by the root mass within the soil layer that is 0.3 m in depth. The results (Table 1) showed that the values of all measured metrics were different with statistical significance between AS and AM. In general, distributions of surficial plant species, coverage, and rooting depth and lateral root spread in underground root systems are determined by the prevailing physical habitats [42]. In particular, AS tends to appear in areas with topographic lows and occupied by seasonally saturated water (i.e., in swales), whereas AM is typically developed around the apex of river bends (Figure 1d). Vegetation communities of AS are mainly composed of cold-tolerant hygrophytes and hydromesophytes, which have a simple community structure. Blysmus sinocompressus is the dominant and healthy specie ( Figure 1e and Table 1) and its coverage reaches nearly 98%. The remaining 2% is contributed from other species, such as Ranunculus nephelogenes and Pedicularis longiflora. The dominant species of AM are Kobresia pygmaea, Elymus nutans, and Potentilla saundersiana, accounting for 30% of the coverage (Table 1). Other herbaceous species, such as Poa annua, Nardostachys jatamansi, Saxifraga montana, Aconitum tanguticum, and degraded species, such as Leontopodium pusillum, Oxytropis ochrocephala, take up to 30% of the coverage. This means that AM only takes 60% of the habitat (Figure 1f). Compared with AS, the number of species and the mesophytes in AM are significantly high (Table 1). These surficial ecological differences between AS and AM must affect their underground properties and the associated mechanical characteristics, which were investigated in this study.

Field Experiments and Measurements
Field experiments were conducted at the two selected sites that are about 200 m apart from each other in the study area ( Figure 1c). These sites have natural vegetation communities and are not disturbed by livestock and human activities. Therefore, they are representative of the general vegetation distribution in the study area. Site 1 is covered by AS, while site 2 is topped by AM. The height of the riverbank to the water level of the channel flow at site 1 is about 0.2-0.4 m lower than that at site 2, suggesting that the groundwater level at site 1 is higher than that at site 2. This ·difference is the main factor that determines the spatial distribution of AS and AM in the study area. The riverbank at both sites had been undercut by river flows, forming a cantilever arm with a width of about 0.2 m and thickness of about 0.3 m at site 1, while 0.3 and 0.35 m, respectively at site 2 ( Figure 2). Both cantilever arms were stable with no tensile cracks developed on the top. They were underlined by a layer of silt and sand with a thickness of approximately 0.3 to 0.6 m, whose bottoms were gradually intergraded into fine gravels deposited at the toes of both banks ( Figure 2). Field experiments were conducted during an October storm event in 2018, such that the entire bank profiles were exposed, facilitating subsequent bank excavation and measurements.  At each site, experiments were performed at the face of a bank profile stretched longitudinally for 8 m (Figure 3). An area of 1 m long on each side (i.e., E1 and E2 in Figure 3) was selected for extending the existing cantilever arm by excavating the lower part of the bank profile. The two areas were selected apart from each other by 3 m to ensure that the excavation of the first one would not affect the stability of the second. During each experiment, excavation was executed gradually, such that the development of tensile cracks on top of the cantilever arm could be observed. Their emergence indicated the status of arm stability. Excavation ended when the cantilever arm reached the threshold that triggered the failure of the cantilever arm. Three types of measurements were subsequently performed after bank collapse, determining tensile strengths of individual roots, measuring root diameter, number, and distribution, and sampling soil for both in situ and laboratory analyses.
The tensile strength of a single root for a given plant species (Tr, MPa) may be determined by [2], where F is the maximum pullout force of measured individual roots (N) and d is the diameter of the corresponding single root (10 −3 m). The value of F was measured using a HP-500 digital push-pull meter with a maximum load of 500 N and an indication error of ±0.5% (Leqing Aidebao Instruments Co. Ltd., Leqing, China). This measurement was taken at three plots within the experimental area; two (i.e., R1 and R3) were on the new face of the upper soil-root layer after cantilever failure and the other (i.e., R2) was at the face of the original upper bank ( Figure 3). This design assured that the measured tensile strength accounted for spatial variability. At each plot, vegetation roots were separated from their surrounding soils by brushing soil particles away. Then, every single root was connected to the tension meter by a clamp. The value of F was recorded after applying a horizontal tensile force at a uniform speed until the root is broken or pulled out ( Figure 4). The diameter of the same root was also measured using a Vernier caliper with an accuracy of 2 × 10 −5 m. If the broken position of the root was close to the clamp, then the recorded F value was biased and not counted. There were more than 30 roots measured in each plot. At each site, experiments were performed at the face of a bank profile stretched longitudinally for 8 m (Figure 3). An area of 1 m long on each side (i.e., E1 and E2 in Figure 3) was selected for extending the existing cantilever arm by excavating the lower part of the bank profile. The two areas were selected apart from each other by 3 m to ensure that the excavation of the first one would not affect the stability of the second. During each experiment, excavation was executed gradually, such that the development of tensile cracks on top of the cantilever arm could be observed. Their emergence indicated the status of arm stability. Excavation ended when the cantilever arm reached the threshold that triggered the failure of the cantilever arm. Three types of measurements were subsequently performed after bank collapse, determining tensile strengths of individual roots, measuring root diameter, number, and distribution, and sampling soil for both in situ and laboratory analyses.
The tensile strength of a single root for a given plant species (T r , MPa) may be determined by [2], where F is the maximum pullout force of measured individual roots (N) and d is the diameter of the corresponding single root (10 −3 m). The value of F was measured using a HP-500 digital push-pull meter with a maximum load of 500 N and an indication error of ±0.5% (Leqing Aidebao Instruments Co. Ltd., Leqing, China). This measurement was taken at three plots within the experimental area; two (i.e., R1 and R3) were on the new face of the upper soil-root layer after cantilever failure and the other (i.e., R2) was at the face of the original upper bank ( Figure 3). This design assured that the measured tensile strength accounted for spatial variability. At each plot, vegetation roots were separated from their surrounding soils by brushing soil particles away. Then, every single root was connected to the tension meter by a clamp. The value of F was recorded after applying a horizontal tensile force at a uniform speed until the root is broken or pulled out ( Figure 4). The diameter of the same root was also measured using a Vernier caliper with an accuracy of 2 × 10 −5 m. If the broken position of the root was close to the clamp, then the recorded F value was biased and not counted. There were more than 30 roots measured in each plot.
Water 2020, 12, x FOR PEER REVIEW 6 of 19  Root diameters and numbers were measured within the front faces of the soil-root layer (0.3 m in thickness) in the collapsed bank block at R1 and R3, and the face of the same layer in the original bank at R2 (Figure 3). The surface of the face was washed by water to expose the root system and then divided into 30 small zones by laying a grid whose cells had the size of 0.10 × 0.10 m on the top. Within each of the 30 cells, root diameters were measured using the same Vernier caliper and the number of roots was recorded. All the roots that passed through the soil-root layer were measured. These values were classified into three groups based on the root diameter (i.e., <0.5, 0.5-1.0, and >1.0 × 10 −3 m), and the depths with the soil-root layer (i.e., 0-0.10, 0.11-0.20, and 0.21-0.30 m). The root area ratio (RAR), defined as the ratio of the total area of all roots passing through the collapsed soilroot layer to the area of the plot (i.e., 1 × 0.3 m), was then calculated using the measured values. The measurements were repeated three times at each plot (i.e., R1, R2, and R3).
Soil samples were taken from two small plots in the soil-root layer around the depth of 0.10-0.20 m beneath the ground surface (i.e., S1 and S3 in Figure 3) and from two plots in the lower layer between 0.3 and 0.5 m from the ground surface (i.e., S2 and S4 in Figure 3) at both sites. These samples were carefully stored and transported back to the soil mechanics laboratory of Geological Engineering Department, Qinghai University for further analysis. The oven-drying method [43] was used to determine the soil moisture content, and the sieving method [43] was applied to determine the grain size distribution. Each measurement was repeated three times for accuracy. Additionally, a cutting ring with a capacity volume of 6.0 × 10 −5 m 3 was used in situ to determine the soil bulk density at each plot ( Figure 3).

Calculation
In the process of cantilever bank collapse, the weight of the cantilever body is originally balanced   Root diameters and numbers were measured within the front faces of the soil-root layer (0.3 m in thickness) in the collapsed bank block at R1 and R3, and the face of the same layer in the original bank at R2 (Figure 3). The surface of the face was washed by water to expose the root system and then divided into 30 small zones by laying a grid whose cells had the size of 0.10 × 0.10 m on the top. Within each of the 30 cells, root diameters were measured using the same Vernier caliper and the number of roots was recorded. All the roots that passed through the soil-root layer were measured. These values were classified into three groups based on the root diameter (i.e., <0.5, 0.5-1.0, and >1.0 × 10 −3 m), and the depths with the soil-root layer (i.e., 0-0.10, 0.11-0.20, and 0.21-0.30 m). The root area ratio (RAR), defined as the ratio of the total area of all roots passing through the collapsed soilroot layer to the area of the plot (i.e., 1 × 0.3 m), was then calculated using the measured values. The measurements were repeated three times at each plot (i.e., R1, R2, and R3).
Soil samples were taken from two small plots in the soil-root layer around the depth of 0.10-0.20 m beneath the ground surface (i.e., S1 and S3 in Figure 3) and from two plots in the lower layer between 0.3 and 0.5 m from the ground surface (i.e., S2 and S4 in Figure 3) at both sites. These samples were carefully stored and transported back to the soil mechanics laboratory of Geological Engineering Department, Qinghai University for further analysis. The oven-drying method [43] was used to determine the soil moisture content, and the sieving method [43] was applied to determine the grain size distribution. Each measurement was repeated three times for accuracy. Additionally, a cutting ring with a capacity volume of 6.0 × 10 −5 m 3 was used in situ to determine the soil bulk density at each plot ( Figure 3).

Calculation
In the process of cantilever bank collapse, the weight of the cantilever body is originally balanced Root diameters and numbers were measured within the front faces of the soil-root layer (0.3 m in thickness) in the collapsed bank block at R1 and R3, and the face of the same layer in the original bank at R2 ( Figure 3). The surface of the face was washed by water to expose the root system and then divided into 30 small zones by laying a grid whose cells had the size of 0.10 × 0.10 m on the top. Within each of the 30 cells, root diameters were measured using the same Vernier caliper and the number of roots was recorded. All the roots that passed through the soil-root layer were measured. These values were classified into three groups based on the root diameter (i.e., <0.5, 0.5-1.0, and >1.0 × 10 −3 m), and the depths with the soil-root layer (i.e., 0-0.10, 0.11-0.20, and 0.21-0.30 m). The root area ratio (RAR), defined as the ratio of the total area of all roots passing through the collapsed soil-root layer to the area of the plot (i.e., 1 × 0.3 m), was then calculated using the measured values. The measurements were repeated three times at each plot (i.e., R1, R2, and R3).
Soil samples were taken from two small plots in the soil-root layer around the depth of 0.10-0.20 m beneath the ground surface (i.e., S1 and S3 in Figure 3) and from two plots in the lower layer between 0.3 and 0.5 m from the ground surface (i.e., S2 and S4 in Figure 3) at both sites. These samples were carefully stored and transported back to the soil mechanics laboratory of Geological Engineering Department, Qinghai University for further analysis. The oven-drying method [43] was used to determine the soil moisture content, and the sieving method [43] was applied to determine the grain size distribution. Each measurement was repeated three times for accuracy. Additionally, a cutting ring with a capacity volume of 6.0 × 10 −5 m 3 was used in situ to determine the soil bulk density at each plot ( Figure 3).

Calculation
In the process of cantilever bank collapse, the weight of the cantilever body is originally balanced by the tensile strength from the soil-root mixture of the body. Extension of the cantilever arm (by artificial excavation in this case or by fluvial erosion under the natural condition) will increase its weight. Once the tensile strength is insufficient to balance the weight, tensile cracks develop from the top of the cantilever body. At this time, both tensile stress and compressive stress exhibit a triangle distribution, and the center axis of the failing cantilever block is located in the stress center of the cantilever body below the crack [44] ( Figure 5). Further extension of the cantilever arm catalyzes the development of the tensile crack until the failure of the cantilever body occurs. Under the critical condition of the failure, the external moment of the cantilever body is balanced by the resistance moment of the soil-root mixture, which may lead to [44], where W = ρgb c d 1 l (N) is the weight of the cantilever body; ρ is the bulk density of the soil-root mixture (kg m −3 ); b c is the critical width of the cantilever arm (m); d 1 is the thickness of the cantilever layer (m), d t is the depth of the crack from the top of the bank (m); l is the unit length of the cantilever layer (m), that is 1 m; a = σ t / σ c , σ t and σ c are tensile and compressive stress of the soil-root mixture (N m −2 ), following Ajaz's results [45], a = 0.1. Substituting W = ρgb c d 1 l into Equation (2), the tensile strength σ t of the soil-root mixture of the cantilever body may be expressed as Water 2020, 12, x FOR PEER REVIEW 7 of 19 artificial excavation in this case or by fluvial erosion under the natural condition) will increase its weight. Once the tensile strength is insufficient to balance the weight, tensile cracks develop from the top of the cantilever body. At this time, both tensile stress and compressive stress exhibit a triangle distribution, and the center axis of the failing cantilever block is located in the stress center of the cantilever body below the crack [44] ( Figure 5). Further extension of the cantilever arm catalyzes the development of the tensile crack until the failure of the cantilever body occurs. Under the critical condition of the failure, the external moment of the cantilever body is balanced by the resistance moment of the soil-root mixture, which may lead to [44], is the weight of the cantilever body; ρ is the bulk density of the soil-root mixture (kg m −3 ); bc is the critical width of the cantilever arm (m); d1 is the thickness of the cantilever layer (m), dt is the depth of the crack from the top of the bank (m); l is the unit length of the cantilever layer (m), that is 1 m; a = σt / σc, σt and σc are tensile and compressive stress of the soil-root mixture (N m −2 ), following Ajaz's results [45], a = 0.1. Substituting (2), the tensile strength σt of the soil-root mixture of the cantilever body may be expressed as Using the measured d1 and dt in our tested AS and AM bank blocks, we calculated their σt values. The product of Tr and RAR also reflects the tensile strength of a soil-root mixture. We compared these two types of calculations for the tensile strength.

Characteristics of Soils and Root Distribution
Our results demonstrated that the composition and physical properties of the two layers in the vertical profile of the riverbank are significantly different between the two sites (Table 2). Generally, the upper soil-root layer was composed of silt, while the lower layer consisted of silty sand with some poorly graded fine gravel. In the upper layer, the average bulk density and moisture content at the two sites were 1260 kg·m −3 and 1560 kg·m −3 ; 39.36% and 41.72%, respectively. They were about 180 kg·m −3 Using the measured d 1 and d t in our tested AS and AM bank blocks, we calculated their σ t values. The product of T r and RAR also reflects the tensile strength of a soil-root mixture. We compared these two types of calculations for the tensile strength.

Characteristics of Soils and Root Distribution
Our results demonstrated that the composition and physical properties of the two layers in the vertical profile of the riverbank are significantly different between the two sites (Table 2). Generally, the upper soil-root layer was composed of silt, while the lower layer consisted of silty sand with some poorly graded fine gravel. In the upper layer, the average bulk density and moisture content at the two sites were 1260 kg·m −3 and 1560 kg·m −3 ; 39.36% and 41.72%, respectively. They were about 180 kg·m −3 and 405 kg·m −3 ; 17.12% and 31.42% higher than those in the lower layer. These two soil properties in the soil-root layer of AS were 21.96% and 3.57% higher than those in that of AM. For the lower layer, the average bulk density and moisture content at site 1 were greater than those at site 2 by 118.67% and 3.27%, respectively. It is obvious that the physical parameters (i.e., bulk density and moisture content) of the soil-root mixture of AS and AM are disparate. These differences could have different influences on the tensile strength of the root system in the soil [27], which lend the support for our analyses of tensile strength for both single roots and the soil-root mixture in this study. Although the distribution of roots over different classes of root diameters is highly variable for different species [46], it may still be characterized by vertical patterns of root number and diameters along the depth of a riverbank. In this study, AS is mainly composed of Blysmus sinocompressus, belonging to the Cyperaceous family. This species has a typical dense and fibrous root system that may be up to 0.8 m long, and its rhizomes are typically about 0.25-0.60 m long. These roots mix with the surrounding soil, forming a soil-root layer that extends from the bank surface to the depth of 0.30 m. Our results showed that the total root number within the experimental block of this layer, which was 1.0 m long and 0.3 m deep (Figure 3), was 4345, with 2505, 1160, and 680 in the depth ranges of 0-0.10, 0.11-0.20, and 0.21-0.30 m, respectively. Among these roots, those with the root diameter <0.5, 0.5-1.0, and >1.0 × 10 −3 m took about 66%, 25%, and 9%, respectively. AS is dominated by finer roots, which is evidenced by the fact that within each of the three depth ranges, they took 66%, 67%, and 76%, respectively (Figure 6a). Roots with medium diameters (0.5-1.0 × 10 −3 m) in all three depth ranges took about 25%, while those with the diameters >1.0 × 10 −3 m only existed in the depths of 0-0.10 and 0.11-0.20 m, taking merely 9% and 8%, respectively. More fine roots (<0.5 × 10 −3 m) developed in the shallow depth range rather than in the deep depth range, supported by their distributions of 56%, 26%, and 18% in the depth ranges of 0-0.10, 0.11-0.20, and 0.21-0.30 m, respectively (Figure 6a). The roots with the mean diameter (0.5-1.0 × 10 −3 m) followed a similar vertical distribution, featured by 58%, 27%, and 15% in the three depth ranges, respectively. Roots with the diameters >1.0 × 10 −3 m only existed in the first two depth ranges with 71% and 29%, respectively. 0.11-0.20 m, taking merely 9% and 8%, respectively. More fine roots (<0.5 × 10 −3 m) developed in the shallow depth range rather than in the deep depth range, supported by their distributions of 56%, 26%, and 18% in the depth ranges of 0-0.10, 0.11-0.20, and 0.21-0.30 m, respectively (Figure 6a). The roots with the mean diameter (0.5-1.0 × 10 −3 m) followed a similar vertical distribution, featured by 58%, 27%, and 15% in the three depth ranges, respectively. Roots with the diameters >1.0 × 10 −3 m only existed in the first two depth ranges with 71% and 29%, respectively. AM. The number within each color of the three columns represents the percentage of roots in each sub-layer of the soil-root layer (i.e., the sum of the numbers in each color is 100%). The number outside of the columns represents the percentages of three root diameters in each sub-layer (i.e., the sum of these numbers along each column is 100%). The scaling of the y-axis is different between (a,b). The number within each color of the three columns represents the percentage of roots in each sub-layer of the soil-root layer (i.e., the sum of the numbers in each color is 100%). The number outside of the columns represents the percentages of three root diameters in each sub-layer (i.e., the sum of these numbers along each column is 100%). The scaling of the y-axis is different between (a,b).
The AM is dominated by Kobresia pygmaea, Elymus nutans, and Potentilla saundersiana. Their root number took about 93% of the total in AM. The first two plants have dense and fibrous roots that are about 0.15-0.55 m long. Besides the three dominant plants, AM contains degraded species such as Nardostachys chinensis, Oxytropis ochrocephala, Saxifraga montana, and Leontopodium pusillum, whose roots are generally sparse and short. They have a typical tap root system with the length ranging between 0.04 and 0.13 m. The total number of roots in AS was only 2355 with 1420, 580, and 355 in the three depth ranges downward, respectively. Similarly to those in the AS, roots with the diameters <0.5 × 10 −3 m dominated in each of the three depth ranges, taking 71%, 71%, and 76%, respectively (Figure 6b). The remaining roots at each depth were those with the mean diameters (0.5-1.0 × 10 −3 m). They took roughly the same percentage of the total roots within each depth range, which was 20%, 22%, and 24%, respectively. Again, the coarse roots (i.e., diameter >1.0 × 10 −3 m) only occupied the first two depth ranges, taking a small portion of the total roots in each (i.e., 9% and 7%, respectively). Along the vertical direction, roots with each diameter class demonstrated a similar distribution to their AS peers. For example, the fine roots (diameter <0.5 × 10 −3 m) took 60%, 24%, and 16% from the top to the bottom depth ranges (Figure 6b). A similar vertical distribution appeared for the roots with the mean diameter (i.e., diameters were between 0.1 and 1.0 × 10 −3 m) with 57%, 26%, and 17% in the three depth ranges, respectively. Also, the coarse roots only stayed in the first two depth ranges. Most of them (75%) occupy in the 0-0.10 m depth range (Figure 6b).
Although roots in AS and AM showed similar spatial and diameter distributions, which led to their similar average root diameters (i.e., 0.46 ± 0.31 ×10 −3 m and 0.41 ± 0.39 ×10 −3 m, respectively) in the entire mixed layer of 0.30 m, their root numbers were greatly different, which can be proved by a two-sample difference test (p < 0.05). The root number of AS was 46% higher than that of AM. Consequently, the RAR of the AS experimental block was 0.225% on average, which was about twice that of the AM block (Table 2). These results showed that the root number, which plays an important role in resisting tensile force imposed to the experimental block, is reduced greatly when AS is degraded to AM.

Tensile Trength of Individual Roots and Soil-Root Mixture
The effect of roots on soil strength of the riverbank does not only depend on the root number, but also the tensile strength of a single root (T r ) (i.e., Equation (1)). The average T r of the dominant plant Blysmus sinocompressus, whose root diameters ranged between 0.20 and 1.76 × 10 −3 m at site 1, was 31,310 kPa, and the maximum value can reach up to 128,000 kPa. The T r for roots with the diameters ranging between 0.24 and 1.88 × 10 −3 m at site 2 (i.e., AM) was 16,160 kPa. The average T r of the AS block was 48.4% higher than that of the AM block. There existed a strong relationship between the root diameter and T r (Figure 7), which indicated that T r decreases as the root diameter increases for both AS and AM. This finding is consistent with those reported in earlier studies [11,[46][47][48][49]. This relationship may be described by a power function, whose exponent was different between the two sites (Figure 7). For the AS block, when the root diameter was less than 0.70 × 10 −3 m, T r decreased sharply as the root diameter increased, nearly following a linear relationship. Specifically, the value of T r decreased drastically from 127,930 to 30,810 kPa, as the root diameter only increased from 0.20 to 0.70 × 10 −3 m. As the root diameter continuously increased from 0.71 to 1.76 × 10 −3 m, T r only decreased from 26,480 to 7050 kPa. The differences of T r were 97,120 kPa and 19,430 kPa before and after the threshold root diameter (i.e., 0.7 × 10 −3 m). For the AM block, when the diameter was less than 0.4 × 10 −3 m, T r decreased greatly as the root diameter increased, and its value decreased from 73,380 to 23,590 kPa with the average of 39,050 kPa. However, the difference of T r was only 25,860 kPa when the root diameter increased from 0.4 to 1.88 × 10 −3 m. This change was much less than that for root diameters less than 0.4 × 10 −3 m (i.e., 49,790 kPa). Therefore, the root diameter of 0.7 and 0.4 × 10 −3 m can be taken as a threshold of the AS and AM blocks, respectively. Values of T r were more sensitive to the changes of root diameters less than the threshold while remaining less variable for coarse roots whose diameters are greater than the threshold. Because there were about 93% and 90% fine roots with diameters less than 1 × 10 −3 m in the mixed layer of AS and AM blocks, respectively, the difference in the root number should not be the main cause for the different Tr values between AS and AM blocks. Rather, their difference represented the true mechanical discrepancy between the two. In other words, the tensile strength of individual roots in healthy meadow plants (i.e., AS) is generally higher than that in the degraded meadow plants (i.e., AM).
The tensile strength of the soil-root mixture should be relevant to the physical properties of the mixture. At site 1, the moisture content and volume of the two tested bank blocks were different. They were 40.92% and 0.175 m 3 for block 1, and 41.71% and 0.187 m 3 fo block 2, respectively (Tables 2 and 3). The tensile strength of the soil-root mixture (σt) between the two tested sites was also slightly different. It was 66.86 kPa for Site 1 and 67.93 kPa for Site 2 (Table 3). At site 2, the moisture content and volume of the two tested experimental blocks were different by 1.06% and 0.004 m 3 (Tables 2 and 3), respectively. However, the variation of σt between the two was still minor, merely 1.34 kPa (Table 3). This shows that for either AS or AM, repetitive measurements of σt for different experimental blocks were consistent with each other, though the moisture content and volume of the experimental blocks might be slightly different.  Because there were about 93% and 90% fine roots with diameters less than 1 × 10 −3 m in the mixed layer of AS and AM blocks, respectively, the difference in the root number should not be the main cause for the different T r values between AS and AM blocks. Rather, their difference represented the true mechanical discrepancy between the two. In other words, the tensile strength of individual roots in healthy meadow plants (i.e., AS) is generally higher than that in the degraded meadow plants (i.e., AM).
The tensile strength of the soil-root mixture should be relevant to the physical properties of the mixture. At site 1, the moisture content and volume of the two tested bank blocks were different. They were 40.92% and 0.175 m 3 for block 1, and 41.71% and 0.187 m 3 fo block 2, respectively (Tables 2 and 3). The tensile strength of the soil-root mixture (σ t ) between the two tested sites was also slightly different. It was 66.86 kPa for Site 1 and 67.93 kPa for Site 2 (Table 3). At site 2, the moisture content and volume of the two tested experimental blocks were different by 1.06% and 0.004 m 3 (Tables 2 and 3), respectively. However, the variation of σ t between the two was still minor, merely 1.34 kPa (Table 3). This shows that for either AS or AM, repetitive measurements of σ t for different experimental blocks were consistent with each other, though the moisture content and volume of the experimental blocks might be slightly different. These results suggest that variation of the moisture content and volume of the experimental blocks for the same type of plants has a negligible impact on σ t . Nonetheless, between the AS and AM, the mean tensile strength of the rooted soil was significantly different, which was 67.39 and 21.96 kPa, respectively ( Table 3). The value of σ t for the AS at site 1 was about 3.07 times higher than that of the AM at site 2 ( Table 3). The dramatic difference clearly reflected the differences in the mechanical characteristics of the root systems between the AS and AM blocks. It follows that the characteristics of root distribution and tensile strength of individual roots are the important factors that influence the tensile strength of the soil-root mixture.
The product of RAR and T r is often used to evaluate the contribution of the root system to the soil tensile strength [9,27,[49][50][51][52][53][54]. For a soil-vegetation mixture, the tensile strength is often viewed as that from the soil and root system, and the latter is much greater than the former [9,51]. In this study, the product of RAR and T r for the AS and AM was 70.45 and 18.58 kPa, respectively. Their ratio was 3.79, which was greater than that of σ t for the AS and AM blocks (i.e., 3.07) ( Table 3). This indicates that conditions of the root system in a soil-root mixture are critical for determining the mechanical characteristic of the mixture.
In the two excavation tests for AS, the depths of the developed cracks in the collapsed block were 0.035 and 0.037 m, respectively, with the average of 0.036 m (Table 3). However, for AM, these depths were 0.048 and 0.050 m, respectively, giving rise to the average of 0.049 m ( Table 3). The crack developed in the AM collapsed block was deeper by about 27% than that in the AS block.

Effect of Degraded Riparian Vegetation on Tensile Strength of Individual Roots and the Soil-Root Mixture
Continuous vegetation degradation has forced alpine meadow, dominated by members of the Cyperaceae, to transform into bare land and subsequently Heitutan on the Qinghai-Tibet Plateau in western China. Areas of severely degraded alpine meadow on the Qinghai-Tibet Plateau are referred to as Heitutan and are characterized by increased proportions of bare land, reduced edible herbage, and commensurate increases in the dominance of less palatable species [55]. Specifically, when AS degenerates into AM, the distribution of dominant plants reduced significantly, from an original coverage of 98% to only about 30%. This change of the surface vegetation communities has further led to changes in the plant's underground biomass. Our results showed that the root number of AM block was 46% less than that of the AS block within the depth of 0-0.30 m (Figure 6). For plants of AS, many roots were much longer than this depth, some of which may have a length of up to 0.80 m (Figure 2a). For the plants of AM, however, the lengths of most of their roots are less than 0.30 m (Figure 2b). Thus, transforming from AS to AM means that deep-rooted plants with dense root systems are gradually replaced by plants with short and sparse root systems [35,[55][56][57][58]. This change indicates that the shallow-rooted plants of AM are more vulnerable to high evaporation, livestock trampling, and human activities in the study area, resulting in further degradation [35].
In addition to their different densities and lengths, the two types of plants also have different root structures. The dominant plant of AS (i.e., Blysmus sinocompressus) has developed rhizomes. Roots derived from them often grow laterally (Figure 8a). Moreover, some roots have a wave shape (Figure 8a). In AM, however, most roots of the dominant plant (i.e., Kobresia pygmaea) grow vertically with no rhizome (Figure 8b). The other two dominant plants (i.e., Elymus nutans and Potentilla saundersiana) share a similar root structure. Compared with that of plant roots in AM, the structure of the dominant plant in AS effectively enhances the tensile strength of the soil-root mixture because (1) the laterally distributed roots may increase root density and (2) the wave-shaped roots have greater contact area with the surrounding soils and may resist higher external force by straightening its shape [59].
( Figure 2b). Thus, transforming from AS to AM means that deep-rooted plants with dense root systems are gradually replaced by plants with short and sparse root systems [35,[55][56][57][58]. This change indicates that the shallow-rooted plants of AM are more vulnerable to high evaporation, livestock trampling, and human activities in the study area, resulting in further degradation [35].
In addition to their different densities and lengths, the two types of plants also have different root structures. The dominant plant of AS (i.e., Blysmus sinocompressus) has developed rhizomes. Roots derived from them often grow laterally (Figure 8a). Moreover, some roots have a wave shape (Figure 8a). In AM, however, most roots of the dominant plant (i.e., Kobresia pygmaea) grow vertically with no rhizome (Figure 8b). The other two dominant plants (i.e., Elymus nutans and Potentilla saundersiana) share a similar root structure. Compared with that of plant roots in AM, the structure of the dominant plant in AS effectively enhances the tensile strength of the soil-root mixture because (1) the laterally distributed roots may increase root density and (2) the wave-shaped roots have greater contact area with the surrounding soils and may resist higher external force by straightening its shape [59]. The above-mentioned characteristics of root diameters and branching, and the tortuosity of plant roots, mainly affected the mechanical behavior of individual roots [28]. Our study showed that the The above-mentioned characteristics of root diameters and branching, and the tortuosity of plant roots, mainly affected the mechanical behavior of individual roots [28]. Our study showed that the influence of vegetation degradation is on not only the root number and structure but also on the mechanical characteristics of roots. In the study area, effective roots (i.e., roots with diameters less than 1 × 10 −3 m) account for about 93% in the layer of 0-0.30 m (Figure 6) below the surface, indicating that most root systems mainly play the role of reinforcement [60]. When the cantilever arm of the meandering riverbank is formed, the soil and root system is subjected to the external load from the weight of the arm, giving rise to the deformation of the arm. The root system can convert part of the external load into the tensile stress and dissipate it to the surrounding soil through the soil-root interface. In this way, the root system and the surrounding soil particles can work together to balance the load and enhance the soil-root tensile resistance [61].
The tensile strength of the root system is a critical factor that directly reflects the effect of rooted soil consolidation. In our study, the average tensile strength of the root for the dominant plants in the AS reached 31,310 kPa, which was 48% higher than that of the AM (16,155 kPa). Mattia et al. [47] and Li et al. [62] measured the root tensile strength of Gramineae plants of Lygeum spartum, Stipa purpurea and sedge plant of Kobresia pygmaea. They found that their tensile strength ranged from 36,260 to 45,670 kPa, close to that of the sedge plant, Blysmus sinocompressus (i.e., 31,310 kPa) in our study. The tensile strength of the root for degraded plants of Potentilla bifurca, Ajania tenuifolia and Saussurea salsa decreases significantly, from 5110 to 25,610 kPa [62]. In our study, this value for the degraded meadow plants (i.e., the AM) was within this range. Evidently, the mechanical characteristics of roots changed because of degeneration. The root tensile strength of healthy meadow plants (i.e., the AS) is much higher than that of the degraded plants (i.e., the AM).
The contribution of the root system to the tensile strength of the rooted soil may be appropriately evaluated using the product of RAR and T r [9,27,49,51,52]. To further analyze the reduction of the tensile strength in the soil-root mixture due to degradation, we compared our calculated values with those from earlier studies [12,47,62]. All of the products of RAR and root tensile strength of the four degraded herbaceous plants (i.e., 1, 2, 3, 5 in Figure 9) and the AM (i.e., 4) fell in the range of 2.04-38.4 kPa, which are lower than those of four healthy herbaceous plants (i.e., 7, 8, 9, 10 in Figure 9) and the AS (i.e., 6), which are within the range of 72.01-118.74 kPa. Clearly, there is a discrepancy of this product between healthy and degenerated plants, separated by the threshold of 50 kPa (Figure 9). strength in the soil-root mixture due to degradation, we compared our calculated values with those from earlier studies [12,47,62]. All of the products of RAR and root tensile strength of the four degraded herbaceous plants (i.e., 1, 2, 3, 5 in Figure 9) and the AM (i.e., 4) fell in the range of 2.04-38.4 kPa, which are lower than those of four healthy herbaceous plants (i.e., 7, 8, 9, 10 in Figure 9) and the AS (i.e., 6), which are within the range of 72.01-118.74 kPa. Clearly, there is a discrepancy of this product between healthy and degenerated plants, separated by the threshold of 50 kPa (Figure 9). Because not all of the tensile strength of the roots is mobilized instantaneously at the moment of bank failure [9,50,51,63], the product of RAR and Tr overestimates the true tensile strength ( Table 3). As more and more fiber materials have been used in engineering practices, their ability to improve the tensile strength of the composite has been tested for natural fibers [24,[64][65][66][67] and synthetic fibers [18,24,58,68]. It is generally believed that fiber materials can reinforce tensile strength, and the degree of reinforcement varies for different fiber materials. Tang et al. [3] proposed that the tensile strength of Figure 9. The product of T r and RAR for healthy and degraded grassland plants. 1 Because not all of the tensile strength of the roots is mobilized instantaneously at the moment of bank failure [9,50,51,63], the product of RAR and T r overestimates the true tensile strength ( Table 3). As more and more fiber materials have been used in engineering practices, their ability to improve the tensile strength of the composite has been tested for natural fibers [24,[64][65][66][67] and synthetic fibers [18,24,58,68]. It is generally believed that fiber materials can reinforce tensile strength, and the degree of reinforcement varies for different fiber materials. Tang et al. [3] proposed that the tensile strength of soil-fiber composite (σ composite ) includes two parts, which are the tensile strength of natural soils (σ soil ) and the increase of the tensile strength due to the fiber (∆σ fiber ): According to Zhu et al. [69], the tensile strength of unsaturated clays is 0.7~0.8 times that of the shear strength. In the current study, the ratio of tensile strength to the shear strength of soil without root is taken as 0.75, and the shear strength is 9.29 kPa [70]. Using these values and Equation (4), the reinforced tensile strengths of plant roots in AS and AM may be calculated as 60.43 and 15.0 kPa, respectively. The reinforcement of the healthy plant root system of AS is about 4 times higher than that of the AM. Also, they are less than the products of RAR and T r by 11.58 and 4.36 kPa, respectively.
The relationship between the value of ∆σ fiber for short (0.03 m) and long (0.05 m) fibers and different fiber contents [24] (Figure 10) was established by setting the root content of AS and AM as 0.48% and 0.18%, respectively [70]. It demonstrated a positive correlation between ∆σ fiber and the fiber content [1,20,21,67], and showed that the longer the fiber length, the greater the strengthening effect [3,24]. In Figure 10, the points representing AS and AM are in line with the established relationship for the fiber lengths of 0.05 and 0.03 m, respectively. This consistency indicates that both the root length and root content contribute to ∆σ fiber , and both variables are higher for plants of AS than those for plants of AM.
Furthermore, the tensile strength may also be affected by the initial water content and dry bulk density. In the fiber-reinforced soils, it is negatively related to water content and positively related to dry bulk density [3,24]. The average density of the natural soil-root mixture for AS is 18% higher than that of AM (Table 2). This may be attributed to the decreased root number and the dramatically increased pore volume of surface soil (i.e, soil in the 0-0.05 m layer) [58]. Higher bulk density should increase not only bonding forces between particles, which enhances the tensile strength, but also the interfacial contact area of the fiber-matrix structure, which improves the interfacial shear strength and associated friction [24]. In theory, the bonding forces between particles and friction of the soil-root mixture for AS is stronger than that of AM. Because the difference in the moisture content of the upper layer of soils between AS and AM is merely 1.43% (Table 2), its influence on the tensile strength of the soil-root mixture is relatively small. In this study, the influence of plant degradation on the soil tensile strength is mainly analyzed from a mechanical point of view, and the influences of soil properties (e.g., bulk density, moisture content, porosity, and particle size) on tensile strength and the impacts of degradation on hydraulic properties need to be further studied. It has been well known that the root system can change the hydraulic characteristics of soil [71][72][73][74], because it occupies the pore spaces of soil, thus reducing the porosity and increasing the water-holding capacity of soil [75]. In this study, the root number of AM decreased by about 46% compared with that of AS in the depth of 0-30 cm on the upper part of the riverbank, suggesting that the proportion of pores in the soil mass occupied by roots was reduced, which is evidenced by the fact that the density of soil mass was relatively lower than that of AS. Therefore, the water holding capacity of soil mass should be lower in AM than that in AS, implying that its suction of soil mass should also be lower [76].
According to Equation (4), the contributions of the root systems of AS and AM to the tensile strength of the soil-root mixture are 89.67% and 68.3%, respectively. This means that the influence of degradation on the tensile strength of the riverbank is mainly derived from the roots. In our study, when the AS degenerated to the AM, the rooted soil tensile strength decreased by 67.42% (Table 3). Li et al. [58] showed that the shear strength of the root-soil mixture decreases by 36.0% and 52.3% from severely degraded alpine meadows to moderately and slightly degraded alpine meadows, respectively. Thus, degradation of alpine swamp meadows obviously reduces the mechanical properties of the soils and weakens their ability to resist bank failure, wind and water erosion, and other external forces.
According to Zhu et al. [69], the tensile strength of unsaturated clays is 0.7~0.8 times that of the shear strength. In the current study, the ratio of tensile strength to the shear strength of soil without root is taken as 0.75, and the shear strength is 9.29 kPa [70]. Using these values and Equation (4), the reinforced tensile strengths of plant roots in AS and AM may be calculated as 60.43 and 15.0 kPa, respectively. The reinforcement of the healthy plant root system of AS is about 4 times higher than that of the AM. Also, they are less than the products of RAR and Tr by 11.58 and 4.36 kPa, respectively.
The relationship between the value of Δσfiber for short (0.03 m) and long (0.05 m) fibers and different fiber contents [24] (Figure 10) was established by setting the root content of AS and AM as 0.48% and 0.18%, respectively [70]. It demonstrated a positive correlation between Δσfiber and the fiber content [1,20,21,67], and showed that the longer the fiber length, the greater the strengthening effect [3,24]. In Figure 10, the points representing AS and AM are in line with the established relationship for the fiber lengths of 0.05 and 0.03 m, respectively. This consistency indicates that both the root length and root content contribute to Δσfiber, and both variables are higher for plants of AS than those for plants of AM. Furthermore, the tensile strength may also be affected by the initial water content and dry bulk density. In the fiber-reinforced soils, it is negatively related to water content and positively related to dry bulk density [3,24]. The average density of the natural soil-root mixture for AS is 18% higher than that of AM (Table 2). This may be attributed to the decreased root number and the dramatically increased pore volume of surface soil (i.e, soil in the 0-0.05 m layer) [58]. Higher bulk density should increase not only bonding forces between particles, which enhances the tensile strength, but also the interfacial contact area of the fiber-matrix structure, which improves the interfacial shear strength and associated friction [24]. In theory, the bonding forces between particles and friction of the soilroot mixture for AS is stronger than that of AM. Because the difference in the moisture content of the upper layer of soils between AS and AM is merely 1.43% (Table 2), its influence on the tensile strength of the soil-root mixture is relatively small. In this study, the influence of plant degradation on the soil tensile strength is mainly analyzed from a mechanical point of view, and the influences of soil properties (e.g., bulk density, moisture content, porosity, and particle size) on tensile strength and the impacts of degradation on hydraulic properties need to be further studied. It has been well known Figure 10. The relationship between fiber (root) content and reinforced tensile strength. The letter l in the diagram stands for the fiber length. The hollow circles and triangles are the data from Meriem et al. [24], while the solid circle and triangle refer to the reinforced tensile strength of AS and AM in this study, respectively.

The Role of Root System in Preventing Development of Riverbank Cracks
The development of cracks greatly destroys the integrity of the soil structure, weakens the mechanical properties of soils, reduces stability, increases permeability, intensifies evaporation, and increases soil erosion, resulting in a series of subsequent adverse effects on geotechnical engineering and the environment [77][78][79][80]. The tensile strength of the soil is an important mechanical parameter that controls the initiation and propagation of tensile cracks [3]. Therefore, enhancing the tensile strength of the soil and preventing the occurrence or slowing the expansion of cracks are critical for riverbank protection. Because plants in the AS have dense roots and strong single-root pullout resistance, they may enhance the tensile ductility of the rooted soils and inhibit the initial formation of tensile cracks. In our study, the expansion rate of tension cracks in cantilever arms due to weights of the AS blocks and the duration of block failure are generally slower than those for the AM blocks. This indicates that types of riparian vegetation, fiber root number within the soils, and the root tensile strength may considerably constrain the number of cracks and their propagation rates [81].
When a tested block of cantilever arms is about to fail, a penetrating crack forms from the surface of AS or AM by following the path of least resistance (Figure 11). For the tested block of AS, the crack had a zigzag shape, and its length was about 1.24 m, which is longer than the crack length of AM by 0.14 m. The crack propagated with a preference angle until it was interrupted by roots. This angle seemed greater in the block with a higher root content. This property is consistent with the result of Meriem et al. [24]. At the depth of 0-0.3 m, the root number of AS that passing through the collapsed profile was 1.84 times that of the root number of AM ( Figure 6). As a result, a crack in AS propagated in a way that avoided most roots. Given that the length of the crack was longer than that of AM, the time needed to penetrate through the soil-root mixture is longer. Nonetheless, the crack pattern is smoother on the surface of the AM block because of its lower number of roots. mechanical properties of soils, reduces stability, increases permeability, intensifies evaporation, and increases soil erosion, resulting in a series of subsequent adverse effects on geotechnical engineering and the environment [77][78][79][80]. The tensile strength of the soil is an important mechanical parameter that controls the initiation and propagation of tensile cracks [3]. Therefore, enhancing the tensile strength of the soil and preventing the occurrence or slowing the expansion of cracks are critical for riverbank protection. Because plants in the AS have dense roots and strong single-root pullout resistance, they may enhance the tensile ductility of the rooted soils and inhibit the initial formation of tensile cracks. In our study, the expansion rate of tension cracks in cantilever arms due to weights of the AS blocks and the duration of block failure are generally slower than those for the AM blocks. This indicates that types of riparian vegetation, fiber root number within the soils, and the root tensile strength may considerably constrain the number of cracks and their propagation rates [81].
When a tested block of cantilever arms is about to fail, a penetrating crack forms from the surface of AS or AM by following the path of least resistance (Figure 11). For the tested block of AS, the crack had a zigzag shape, and its length was about 1.24 m, which is longer than the crack length of AM by 0.14 m. The crack propagated with a preference angle until it was interrupted by roots. This angle seemed greater in the block with a higher root content. This property is consistent with the result of Meriem et al. [24]. At the depth of 0-0.3 m, the root number of AS that passing through the collapsed profile was 1.84 times that of the root number of AM ( Figure 6). As a result, a crack in AS propagated in a way that avoided most roots. Given that the length of the crack was longer than that of AM, the time needed to penetrate through the soil-root mixture is longer. Nonetheless, the crack pattern is smoother on the surface of the AM block because of its lower number of roots. Considering the impact of the degraded alpine swamp meadow on tensile strength of riverbank in the UYR region, the grassland management should practice regional rotation grazing, rational distribution of herds, and balanced use of the riparian alpine swamp meadow. In addition, anthropogenic disturbances, in particular engineering construction, should be reduced as much as possible. Considering the impact of the degraded alpine swamp meadow on tensile strength of riverbank in the UYR region, the grassland management should practice regional rotation grazing, rational distribution of herds, and balanced use of the riparian alpine swamp meadow. In addition, anthropogenic disturbances, in particular engineering construction, should be reduced as much as possible.

Conclusions
Under the influence of climate change and increased anthropogenic disturbances, the alpine swamp meadow (AS), which is the main type of vegetation cover in the meandering riverbank of the Upper Yellow River (UYR), is subject to severe degradation and typically has transformed into the alpine meadow (AM). However, little is known about the tensile strength of the roots in soils and the influence of riparian vegetation degradation on the tensile strength of riverbanks. To reflect the actual interaction between soils and the roots with the natural root branching structure for AS and AM and its effect on bank strength, we measured properties of root vertical distribution and number and performed in situ root pullout and artificial excavation tests. Our results led to the following conclusions. First, though spatial and size distributions of roots in the soil-root mixture of AS and AM were similar, AS was characterized by a higher number of roots than AM. Second, the tensile strength of individual roots decreased with the diameter of roots in both AS and AM. Yet, the former always had higher tensile strength than the latter for any given root diameter. Similarly, the tensile strength of the soil-root mixture in AS was about three times higher than that in AM. This difference is mainly caused by the fact that the lateral extended and wave-shaped root structure in AS is more effective in enhancing the resistance of the soil-root mixture to the external force than the simple vertically distributed root structure in AM. Third, the tensile crack developed in the collapsed block for AM was deeper than that for AS, indicating the reduced resistance to the external force as AS gradually degraded to AM. The impact of the roots for the degraded vegetation (i.e., AM) was also reflected by the relatively smooth and shorter crack route on the surface of the collapsed block for AM. These findings call for better ecological management for preventing AS from being degraded to AM.
Although grassland degradation is a worldwide problem [82,83], how such degradation affects the mechanical properties of riparian riverbanks in the alpine environment has not been fully understood.
Our study provides firsthand evidence of weakened bank strength in the soil-root mixture due to degradation of vegetation and will serve as a benchmark for future investigation of riverbank strength, not only in the UYR, but also in other alpine regions in the world.