# Study on the Stability of Slopes Reinforced by Composite Vegetation Combined with a Geogrid under Rainfall Conditions

Authors: Qizhi Hu; Yong Zhou; Gaoliang Tao
Journal: Advances in Civil Engineering
 (2021)

---

## Abstract

The planting of shrubs and trees on geogrid-reinforced slopes is an important ecological slope protection method that is frequently implemented in the rainy areas of southern China. First, this paper analyzes the soil-fixing principle of the geogrid and root system and demonstrates the feasibility of using composite vegetation of shrubs and trees to reinforce the slope with a geogrid. Using the Yushi Expressway project in Guizhou, we conducted a stability analysis of slopes under different working conditions and different reinforcement modes. We determined that the ecological protection method of combining composite vegetation with a geogrid can effectively increase the stability of slopes. The maximum displacement of the ecological slope under rainfall conditions was reduced by 82% compared with the original slope, and the overall stability was improved by 35%. Four factors affect the slope stability: the depth of shrub reinforcement, depth of anchorage of trees, distribution of trees, and spacing of the geogrids. An orthogonal analysis considering these 4 factors with 3 levels was implemented. The following optimal combination was obtained to ensure ecological protection under rainfall conditions: a shrub reinforcement depth of 0.6 m, a tree anchorage depth of 3 m, a grid spacing of 0.4 m, and a top-sparse and bottom-dense tree distribution. The combined slope protection schematic was applied to the Yushi Expressway project in Guizhou, and a strong reinforced slope protection effect was observed.


---

## Body

## 1. Introduction

Under rainfall conditions, the erosion resistance of exposed slopes decreases, and soil erosion increases, which can easily induce major natural disasters such as landslides [1, 2]. Ecological slope protection efforts consider both engineering protection and environmental beautification and are widely used in China. Under the action of plant roots, the surface soil of an exposed slope is consolidated, and the shear strength is increased. The transpiration of plant leaves and suction of roots also improve the suction of the surface soil, stop ground runoff, and prevent scouring, which helps soil and water conservation. However, the contribution of root consolidation to slope stability is limited and not applicable to unstable or understable slopes. With the development of geosynthetics, geogrids have been widely used in slope engineering due to their obvious reinforcement effect and relatively economical nature [3]. Therefore, the planting of shrubs and trees on geogrid-reinforced slopes can improve the stability of slopes and highlight the ecological benefits of such conservation efforts since it is an important method for slope management in rainy areas.At present, many scholars have conducted various studies on geogrids and root reinforced slopes. Among them, numerical simulation technology can intuitively analyze the failure mechanism inside the rock and soil and is widely used in slope stability analysis. For example, Sun et al. used the finite element software Midas-GTS to simulate the performance of reinforced soil slopes, and they compared the numerical simulation results with monitoring data to verify that numerical simulation techniques can reliably analyze geogrid-reinforced slopes [4, 5]. Zhang et al. took the actual engineering of geogrid reinforcement of the slope as an example, analyzed the stability of the slope, and optimized the geogrid reinforcement scheme [6, 7]. Yang et al. discussed the mechanical mechanism of plant roots for soil consolidation and compared the stability of grass, shrub, arbor, and composite vegetation reinforcement under the slope. They noted that composite vegetation had significantly higher reinforcement effect than single vegetation [8]. Waldron et al. conducted a shear test on reinforced soil and noted that reinforced soil mainly exerted the reinforcement effect by improving the original cohesion of the soil [9–11]. Su et al. conducted a pull test of the root system of the arbor. The influence of the root diameter on the anchorage of arbor roots was quantitatively analyzed [12]. Abdi et al. calculated the spatial distribution of plant roots in the soil; they found that the roots of shrubs were mainly distributed 0.6 m below the slope, and the roots of trees could reach 3 m in length [13, 14]. However, these studies did not consider the influence of geogrids and root systems on the hydraulic characteristics of slopes.Rainfall infiltration is an important factor that induces slope instability [15, 16]. Based on the theory of unsaturated seepage, Zhang et al. established a hydraulic-mechanical coupling numerical calculation model for reinforced soil slopes and studied the effects of rainfall infiltration on the stress, displacement, and pore pressure of geogrid-reinforced soil slopes [17, 18]. Zhang et al. analyzed the seepage field distribution of root-fixed soil under rainfall conditions and found that when the root system continued to absorb water, the shear strength of unsaturated soil also increased [19, 20]. Ng et al. studied the size and distribution of suction generated by plant roots in soil and noted that plants could significantly increase soil suction under both dry and rainy conditions [21]. Wang et al. conducted a slope model test under rainfall conditions and believed that plant roots could improve the hydraulic properties of the soil and reduce the impact of rainfall on the slope runoff [22–24]. Song et al. proposed a composite soil treatment and slope protection method using the geocell structures and the wheat straw reinforcement, which can effectively reduce soil erosion on the slope [25].These studies mainly focused on the influence of a single reinforcement mode and a single factor on the slope deformation and stability. There are few studies on slope reinforcement by composite vegetation with geogrids. With the improvement of the dual requirements of economic and ecological benefits, joint ecological slope protection technology has widely become the optimal choice for slope treatment. However, its work performance is more complicated, and there are many influencing factors. There is no clear understanding of the design plan for planting shrubs and trees on slopes reinforced by geogrids. Especially under rainfall conditions, it is of great significance to quantify the influence of various factors on the mechanical and hydraulic characteristics of ecological slopes. Therefore, to popularize this technology in engineering, it is necessary to study the stability of composite vegetation with geogrid reinforcement under rainfall conditions.


## 2. Soil Fixation Mechanism of Geogrid and Root Systems

As common reinforcement materials, geogrid and root systems can affect the mechanical properties of the soil and root system as an active medium. They can absorb water from different soil depths through the transpiration of plant leaves, increase the soil suction, and change the hydraulic characteristics of unsaturated soil. We analyzed the soil fixation mechanism of a geogrid and a root system based on mechanical and hydraulic soil properties as follows.


### 2.1. Mechanical Activation Behavior

Geogrids and root systems have high tensile strengths. After being implanted in the soil, roots can interact with the soil to improve its tensile and shear strength and form a reinforced soil. Therefore, its mechanical properties can be quantitatively analyzed by reinforced soil theory, which mainly includes two types: the friction reinforcement principle and quasi-cohesive force mechanism.The principle of friction reinforcement is that the reinforcing material is firmly fixed inside the soil body under the action of soil pressure; when the relative displacement of the soil body and reinforcing material occurs, the slip force will be transferred to the reinforcing material. Due to the friction between soil and tendons and the high tensile strength of the tendons, the soil deformation can be limited, which improves the soil stability. Sectiondl of the reinforced soil is selected for analysis, as shown in Figure 1:(1)F=2σnfbdl,dT=T1−T2,where F is the frictional resistance, σn is the effective stress that acts on the reinforced material, f is the friction coefficient, b and dl are the width and length of the reinforced material, and dT is the tensile force on the reinforced material. When F > dT, the reinforced soil is in a stable state, and vice versa for destabilization.Figure 1
Principle of friction reinforcement.The mechanism behind quasi-cohesion is to regard the reinforced soil as a composite material. Reinforcements play a similar role in the confining pressure, which limits the soil deformation. Reinforced soil mainly exerts the reinforcement effect by improving the original cohesion of the soil. As shown in Figure2, the strength envelope of reinforced soil is parallel to that of unreinforced soil, and only the point of intersection with the vertical axis is significantly different. When the confining pressure σ3 is constant and the soil reaches the ultimate equilibrium state, the maximum principal stress σ1′ of the reinforced soil is obviously greater than the maximum principal stress σ1 of unreinforced soil, and the strength of the reinforced soil is substantially improved. When the maximum principal force σ1′ is constant, the minimum principal stress σ3′ of unreinforced soil can be calculated using its strength envelope. Figure 2 shows that σ3′ is significantly larger than σ3, which indicates that the tendons can limit the soil deformation and increase the enclosing pressure. The increase in confining pressure is defined as follows:(2)Δσ3=σ3′−σ3=2Δc⋅tan45°−φ2,(3)Δσ3=σ1′−σ1tan245°−φ2=Δσ1⋅tan245°−φ2.Figure 2
Diagram of the strength analysis of reinforced soil.From formulas (2) and (3), we obtain(4)Δc=12Δσ1⋅tan45°−φ2,where Δc is the quasi-cohesive force, φ is the friction angle in the soil body, and Δσ1 is the maximum principal stress difference.


### 2.2. Hydraulic Performance

Unlike other reinforced materials, the root system, as a type of active reinforcement, can play the mechanical role of reinforcement and anchoring, absorb water from different soil depths through the transpiration of plant leaves, reduce the pore water pressure of the soil, improve soil suction, and change the permeability characteristics of unsaturated soil. Therefore, the influence of the root system on soil hydraulic characteristics can be quantitatively analyzed using the unsaturated permeability theory, which is mainly reflected in two aspects.On one hand, based on unsaturated theory, a shear strength formula considering matrix suction is established [26]:(5)τ=c+σ−μatanφ+μa−μwtantanφθ−θrθs−θr,where τ, σ, c, and φ are the shear stress, principal stress, cohesion, and internal friction angle of the soil in sequence, respectively, μa−μw is the matrix suction, μa is the pore-air pressure, μw is the pore water pressure, θ is the volumetric water content, θs is the saturated water content, and θr is the residual water content. This formula shows that plant roots can increase soil suction, reduce the pore water pressure, and consequently increase the shear strength of the soil.On the other hand, the relationship between suction and permeability coefficient of unsaturated soil can be described as follows:(6)Ks=Ks∫Insbθey−es/eyθ′eydy∫Insavebθey−θs/eyθ′eydy,where Ks is the saturated permeability coefficient, s = μa − μw, b = ln (106), y is the dummy variable of the suction integral, and θ′ is the derivative of θ. From this formula, the soil suction is negatively related to the permeability coefficient. The change curve of various soil permeability coefficients with suction in Figure 3 shows that the permeability coefficient of soil decreases with increasing suction. Therefore, the ecological slope can increase the soil suction and reduce the permeability coefficient under the action of root suction. During rainfall, it can reduce the infiltration of rainwater and reduce the risk of soil erosion.Figure 3
Diagram of the change in soil permeability coefficient with suction.After the above analysis, as reinforced materials, the geogrid and root system play the role of reinforced anchors in terms of their mechanical properties, and the principle of soil fixation is identical. In actual engineering practices, due to the diversity of the plant root distribution, the root system can be intertwined with the grid through the grid gap to form a joint entity, which can limit the relative displacement of the geogrid and soil and prevent the root system from being pulled out to improve the shear strength of the soil. Simultaneously, the suction effect of the root system can improve the hydraulic properties of the soil, increase the shear strength of the soil, reduce the infiltration coefficient, and reduce the harm of rainwater to the soil. In summary, the ecological slope protection method of composite vegetation with geogrids can fully utilize the reinforcing and anchoring effects of root systems and geogrids, effectively reduce the risk of soil erosion, and provide a theoretical basis for promotion.


## 2.1. Mechanical Activation Behavior

Geogrids and root systems have high tensile strengths. After being implanted in the soil, roots can interact with the soil to improve its tensile and shear strength and form a reinforced soil. Therefore, its mechanical properties can be quantitatively analyzed by reinforced soil theory, which mainly includes two types: the friction reinforcement principle and quasi-cohesive force mechanism.The principle of friction reinforcement is that the reinforcing material is firmly fixed inside the soil body under the action of soil pressure; when the relative displacement of the soil body and reinforcing material occurs, the slip force will be transferred to the reinforcing material. Due to the friction between soil and tendons and the high tensile strength of the tendons, the soil deformation can be limited, which improves the soil stability. Sectiondl of the reinforced soil is selected for analysis, as shown in Figure 1:(1)F=2σnfbdl,dT=T1−T2,where F is the frictional resistance, σn is the effective stress that acts on the reinforced material, f is the friction coefficient, b and dl are the width and length of the reinforced material, and dT is the tensile force on the reinforced material. When F > dT, the reinforced soil is in a stable state, and vice versa for destabilization.Figure 1
Principle of friction reinforcement.The mechanism behind quasi-cohesion is to regard the reinforced soil as a composite material. Reinforcements play a similar role in the confining pressure, which limits the soil deformation. Reinforced soil mainly exerts the reinforcement effect by improving the original cohesion of the soil. As shown in Figure2, the strength envelope of reinforced soil is parallel to that of unreinforced soil, and only the point of intersection with the vertical axis is significantly different. When the confining pressure σ3 is constant and the soil reaches the ultimate equilibrium state, the maximum principal stress σ1′ of the reinforced soil is obviously greater than the maximum principal stress σ1 of unreinforced soil, and the strength of the reinforced soil is substantially improved. When the maximum principal force σ1′ is constant, the minimum principal stress σ3′ of unreinforced soil can be calculated using its strength envelope. Figure 2 shows that σ3′ is significantly larger than σ3, which indicates that the tendons can limit the soil deformation and increase the enclosing pressure. The increase in confining pressure is defined as follows:(2)Δσ3=σ3′−σ3=2Δc⋅tan45°−φ2,(3)Δσ3=σ1′−σ1tan245°−φ2=Δσ1⋅tan245°−φ2.Figure 2
Diagram of the strength analysis of reinforced soil.From formulas (2) and (3), we obtain(4)Δc=12Δσ1⋅tan45°−φ2,where Δc is the quasi-cohesive force, φ is the friction angle in the soil body, and Δσ1 is the maximum principal stress difference.


## 2.2. Hydraulic Performance

Unlike other reinforced materials, the root system, as a type of active reinforcement, can play the mechanical role of reinforcement and anchoring, absorb water from different soil depths through the transpiration of plant leaves, reduce the pore water pressure of the soil, improve soil suction, and change the permeability characteristics of unsaturated soil. Therefore, the influence of the root system on soil hydraulic characteristics can be quantitatively analyzed using the unsaturated permeability theory, which is mainly reflected in two aspects.On one hand, based on unsaturated theory, a shear strength formula considering matrix suction is established [26]:(5)τ=c+σ−μatanφ+μa−μwtantanφθ−θrθs−θr,where τ, σ, c, and φ are the shear stress, principal stress, cohesion, and internal friction angle of the soil in sequence, respectively, μa−μw is the matrix suction, μa is the pore-air pressure, μw is the pore water pressure, θ is the volumetric water content, θs is the saturated water content, and θr is the residual water content. This formula shows that plant roots can increase soil suction, reduce the pore water pressure, and consequently increase the shear strength of the soil.On the other hand, the relationship between suction and permeability coefficient of unsaturated soil can be described as follows:(6)Ks=Ks∫Insbθey−es/eyθ′eydy∫Insavebθey−θs/eyθ′eydy,where Ks is the saturated permeability coefficient, s = μa − μw, b = ln (106), y is the dummy variable of the suction integral, and θ′ is the derivative of θ. From this formula, the soil suction is negatively related to the permeability coefficient. The change curve of various soil permeability coefficients with suction in Figure 3 shows that the permeability coefficient of soil decreases with increasing suction. Therefore, the ecological slope can increase the soil suction and reduce the permeability coefficient under the action of root suction. During rainfall, it can reduce the infiltration of rainwater and reduce the risk of soil erosion.Figure 3
Diagram of the change in soil permeability coefficient with suction.After the above analysis, as reinforced materials, the geogrid and root system play the role of reinforced anchors in terms of their mechanical properties, and the principle of soil fixation is identical. In actual engineering practices, due to the diversity of the plant root distribution, the root system can be intertwined with the grid through the grid gap to form a joint entity, which can limit the relative displacement of the geogrid and soil and prevent the root system from being pulled out to improve the shear strength of the soil. Simultaneously, the suction effect of the root system can improve the hydraulic properties of the soil, increase the shear strength of the soil, reduce the infiltration coefficient, and reduce the harm of rainwater to the soil. In summary, the ecological slope protection method of composite vegetation with geogrids can fully utilize the reinforcing and anchoring effects of root systems and geogrids, effectively reduce the risk of soil erosion, and provide a theoretical basis for promotion.


## 3. Numerical Simulation of Combined Ecological Slope Protection

### 3.1. Engineering Background

The 10th section of the Guizhou Yushi Expressway starts from Majiao, Tianxing Township, Cengong County and ends at the Yangtang River, Tianxing Township, Cengong County, with a total length of 8.155 km. There are many deep excavation and high filling subgrades along the entire line of the tender section, with 3.47 million square meters of excavation and 3.06 million square meters of filling and rainfall conditions throughout the year. Therefore, the main technical difficulties in the construction process are abundant rainfall, difficulties associated with soil and water conservation, and more high-side and larger slopes. In addition, an urgent blockade of this project is how to adapt to local conditions in the context of green transportation, ensure the safety, stability, and green premise of the integration, optimize various protection methods, and achieve maximum economic and ecological benefits.


### 3.2. Establishment of the Numerical Model and Validation

MIDAS/GTS provides a slope stability analysis under various working conditions, which is widely used in slope stability calculations, and the calculation accuracy has been verified. Therefore, the finite element software Midas-GTS is used to simulate the joint ecological slope protection. A section of the slope of the Yushi Expressway in Guizhou was selected for analysis, which has a slope fill height of 6 m and a slope of 45°. To achieve maximum economic and ecological benefits, ecological protection of the composite vegetation with geogrids is performed on the slopes. The composite vegetation consists of a combination of shrubs and trees. Shrub roots are mainly distributed within 0.6 m below the slope surface. To facilitate the simulation, the arbor root system is simplified. The main roots of the arbor are vertically distributed with a length of 3 m and a diameter of 0.1 m. The lateral roots are distributed symmetrically about the main root, parallel to the slope, with a length of 0.6 m, a diameter of 0.02 m, and a horizontal planting interval of 2 m. The geogrid is 6 m long. The layout spacing is 0.6 m. To reduce the influence of the size effect, the grid of the slope model is extended and converted according to the research of Zheng et al. on the calculation accuracy of numerical simulation [27]. The specific slope form and model dimensions are shown in Figure 4.Figure 4
Slope model.Due to the well-developed root system of shrubs with diverse and dense distribution patterns, they can interlock with the soil and form a joint mass. Therefore, the soil in the shrub reinforcement area is simulated in the form of a root soil complex. The Mohr–Coulomb model is used for plain soil and reinforced soil. Both geogrid and arbor root system were simulated using a separated elastic model. A normal constraint was applied around the model with a fixed constraint at the bottom and a rainfall boundary with an intensity of 5 mm/h at the top for 24 h, which was recorded every 6 h. The initial water level was 4.5 m. According to the engineering geological survey report and corresponding references, the mechanical parameters and hydraulic characteristic parameters of plain soil and reinforced soil are provided. Among them, the mechanical parameters of both geogrid and arboreal root system were obtained from indoor tensile tests, which are shown in Tables1 and 2.Table 1
Mechanical parameters of the materials.
MaterialVolumetric weightr (kN/m3)Cohesionc (kPa)Friction angleφ (°)Elastic modulusE (MPa)Poisson’s ratioμPlain soil1810.322.8540.3Shrub reinforced soil1824.628.394.80.3Arbor lateral roots12.85——5000.3Geogrid———15000.23Table 2
Hydraulic characteristic parameters of the soil.
MaterialFitting parametera (kPa−1)Fitting parameternθr (m3/m3)θs (m3/m3)ks (10−6 m/s)Plain soil0.1141.8060.0170.453.29Shrub reinforced soil0.1761.7980.0180.4563.18To verify the accuracy of the built numerical model, further research was performed on the joint ecological slope. The strength reduction method and stress limit equilibrium method are used to calculate the stability of the original slope. The safety factor of the slope obtained by the strength reduction method is 1.890, which is closer to the safety factor of 1.906 calculated by the stress limit equilibrium method. Thus, the numerical model can be used for the stability calculation of combined ecological slope protection.


### 3.3. Result Analysis

Figure5 shows the displacement field of the original slope and ecological slope under the combined reinforcement after 24 h of rainfall. As seen from the figure, the slope displacement field considerably changes under combined ecological protection. First, in terms of displacement values, the maximum displacements of the original slope and ecological slope are 186.5 mm and 33.4 mm, respectively, i.e., the slope displacement is reduced by 82%. Second, the maximum displacement of the original slope is located at the foot of the slope, the stress is relatively concentrated, and the slope exhibits an overall slippage effect. Meanwhile, the sensitive area of the maximum displacement of the ecological slope has changed and is located in the middle and upper parts of the slope. The reason is that, under the action of rainfall, the weight of the soil sharply rises; the soil at the foot of the slope is most severely eroded by rainwater, and the original slope has an integral sliding from top to bottom. However, under the combined action of the geogrid and root system of the ecological slope, the displacement at the toe of the slope is fully restricted. Simultaneously, the suction effect of the root system reduces the pore water pressure at the foot of the slope and improves the shear strength of the soil at the foot of the slope. However, the effective stress on the tendons in the middle and upper parts of the slope is relatively small, the tendon reinforcement effect is limited, and the connection force between shrub reinforcement area and plain soil is insufficient, which can easily produce relative slip and cause shallow damage to the slope.Figure 5
Displacement field of the original slope and ecological slope after 24 h of rainfall. (a) Original slope. (b) Ecological slope protection by composite vegetation with geogrid reinforcement.
(a)(b)According to research on the mechanism of root soil fixation, considering the influence of the root system on the hydraulic characteristics of the slope, the unsaturated strength reduction method is used to calculate the slope stability under different reinforcement modes. As shown in Figure6, under natural working conditions, the shrub reinforcement has limited effect on the slope, but the tree and geogrid reinforcement has an obvious effect. Compared with the original slope, the slope safety factor increases by 6% and 9%, respectively. In the combined form, the ecological protection method of composite vegetation with the geogrid has the highest safety factor, which is represented by an increase of 12% relative to the reinforcement effect of a single vegetation with geogrid. Under rainfall conditions, as rainfall continued, the effect of the shrub reinforcement considerably improved, and the slope safety factor increased by 6% after 24 h of rainfall. Although the shrub reinforcement has limited contribution to the slope stability, it can effectively reduce the hazard of rainfall on the slope. Among them, the ecological protection method of composite vegetation with a geogrid still obtains the highest safety factor, which is 35% higher than that of the original slope.Figure 6
Slope stability under various protection methods.Regardless of whether rainfall occurs, composite vegetation with the geogrid reinforcement mode has the best performance. Under the anchoring action of the geogrid and arbor root system, the reinforced soil in the shrub reinforcement area can be connected with the stabilizing layer, which improves the stability of the shallow soil. Simultaneously, the root system passes through the grid gap and interweaves with the grid to form a united entity. This unity effectively plays the role of geogrid, root reinforcement, and anchor, limits slope deformation, and changes the sensitive area of maximum ecological slope displacement. It can improve the stability of the slope under natural working conditions and greatly reduce the hazard of rainfall to the slope.


## 3.1. Engineering Background

The 10th section of the Guizhou Yushi Expressway starts from Majiao, Tianxing Township, Cengong County and ends at the Yangtang River, Tianxing Township, Cengong County, with a total length of 8.155 km. There are many deep excavation and high filling subgrades along the entire line of the tender section, with 3.47 million square meters of excavation and 3.06 million square meters of filling and rainfall conditions throughout the year. Therefore, the main technical difficulties in the construction process are abundant rainfall, difficulties associated with soil and water conservation, and more high-side and larger slopes. In addition, an urgent blockade of this project is how to adapt to local conditions in the context of green transportation, ensure the safety, stability, and green premise of the integration, optimize various protection methods, and achieve maximum economic and ecological benefits.


## 3.2. Establishment of the Numerical Model and Validation

MIDAS/GTS provides a slope stability analysis under various working conditions, which is widely used in slope stability calculations, and the calculation accuracy has been verified. Therefore, the finite element software Midas-GTS is used to simulate the joint ecological slope protection. A section of the slope of the Yushi Expressway in Guizhou was selected for analysis, which has a slope fill height of 6 m and a slope of 45°. To achieve maximum economic and ecological benefits, ecological protection of the composite vegetation with geogrids is performed on the slopes. The composite vegetation consists of a combination of shrubs and trees. Shrub roots are mainly distributed within 0.6 m below the slope surface. To facilitate the simulation, the arbor root system is simplified. The main roots of the arbor are vertically distributed with a length of 3 m and a diameter of 0.1 m. The lateral roots are distributed symmetrically about the main root, parallel to the slope, with a length of 0.6 m, a diameter of 0.02 m, and a horizontal planting interval of 2 m. The geogrid is 6 m long. The layout spacing is 0.6 m. To reduce the influence of the size effect, the grid of the slope model is extended and converted according to the research of Zheng et al. on the calculation accuracy of numerical simulation [27]. The specific slope form and model dimensions are shown in Figure 4.Figure 4
Slope model.Due to the well-developed root system of shrubs with diverse and dense distribution patterns, they can interlock with the soil and form a joint mass. Therefore, the soil in the shrub reinforcement area is simulated in the form of a root soil complex. The Mohr–Coulomb model is used for plain soil and reinforced soil. Both geogrid and arbor root system were simulated using a separated elastic model. A normal constraint was applied around the model with a fixed constraint at the bottom and a rainfall boundary with an intensity of 5 mm/h at the top for 24 h, which was recorded every 6 h. The initial water level was 4.5 m. According to the engineering geological survey report and corresponding references, the mechanical parameters and hydraulic characteristic parameters of plain soil and reinforced soil are provided. Among them, the mechanical parameters of both geogrid and arboreal root system were obtained from indoor tensile tests, which are shown in Tables1 and 2.Table 1
Mechanical parameters of the materials.
MaterialVolumetric weightr (kN/m3)Cohesionc (kPa)Friction angleφ (°)Elastic modulusE (MPa)Poisson’s ratioμPlain soil1810.322.8540.3Shrub reinforced soil1824.628.394.80.3Arbor lateral roots12.85——5000.3Geogrid———15000.23Table 2
Hydraulic characteristic parameters of the soil.
MaterialFitting parametera (kPa−1)Fitting parameternθr (m3/m3)θs (m3/m3)ks (10−6 m/s)Plain soil0.1141.8060.0170.453.29Shrub reinforced soil0.1761.7980.0180.4563.18To verify the accuracy of the built numerical model, further research was performed on the joint ecological slope. The strength reduction method and stress limit equilibrium method are used to calculate the stability of the original slope. The safety factor of the slope obtained by the strength reduction method is 1.890, which is closer to the safety factor of 1.906 calculated by the stress limit equilibrium method. Thus, the numerical model can be used for the stability calculation of combined ecological slope protection.


## 3.3. Result Analysis

Figure5 shows the displacement field of the original slope and ecological slope under the combined reinforcement after 24 h of rainfall. As seen from the figure, the slope displacement field considerably changes under combined ecological protection. First, in terms of displacement values, the maximum displacements of the original slope and ecological slope are 186.5 mm and 33.4 mm, respectively, i.e., the slope displacement is reduced by 82%. Second, the maximum displacement of the original slope is located at the foot of the slope, the stress is relatively concentrated, and the slope exhibits an overall slippage effect. Meanwhile, the sensitive area of the maximum displacement of the ecological slope has changed and is located in the middle and upper parts of the slope. The reason is that, under the action of rainfall, the weight of the soil sharply rises; the soil at the foot of the slope is most severely eroded by rainwater, and the original slope has an integral sliding from top to bottom. However, under the combined action of the geogrid and root system of the ecological slope, the displacement at the toe of the slope is fully restricted. Simultaneously, the suction effect of the root system reduces the pore water pressure at the foot of the slope and improves the shear strength of the soil at the foot of the slope. However, the effective stress on the tendons in the middle and upper parts of the slope is relatively small, the tendon reinforcement effect is limited, and the connection force between shrub reinforcement area and plain soil is insufficient, which can easily produce relative slip and cause shallow damage to the slope.Figure 5
Displacement field of the original slope and ecological slope after 24 h of rainfall. (a) Original slope. (b) Ecological slope protection by composite vegetation with geogrid reinforcement.
(a)(b)According to research on the mechanism of root soil fixation, considering the influence of the root system on the hydraulic characteristics of the slope, the unsaturated strength reduction method is used to calculate the slope stability under different reinforcement modes. As shown in Figure6, under natural working conditions, the shrub reinforcement has limited effect on the slope, but the tree and geogrid reinforcement has an obvious effect. Compared with the original slope, the slope safety factor increases by 6% and 9%, respectively. In the combined form, the ecological protection method of composite vegetation with the geogrid has the highest safety factor, which is represented by an increase of 12% relative to the reinforcement effect of a single vegetation with geogrid. Under rainfall conditions, as rainfall continued, the effect of the shrub reinforcement considerably improved, and the slope safety factor increased by 6% after 24 h of rainfall. Although the shrub reinforcement has limited contribution to the slope stability, it can effectively reduce the hazard of rainfall on the slope. Among them, the ecological protection method of composite vegetation with a geogrid still obtains the highest safety factor, which is 35% higher than that of the original slope.Figure 6
Slope stability under various protection methods.Regardless of whether rainfall occurs, composite vegetation with the geogrid reinforcement mode has the best performance. Under the anchoring action of the geogrid and arbor root system, the reinforced soil in the shrub reinforcement area can be connected with the stabilizing layer, which improves the stability of the shallow soil. Simultaneously, the root system passes through the grid gap and interweaves with the grid to form a united entity. This unity effectively plays the role of geogrid, root reinforcement, and anchor, limits slope deformation, and changes the sensitive area of maximum ecological slope displacement. It can improve the stability of the slope under natural working conditions and greatly reduce the hazard of rainfall to the slope.


## 4. Optimization of the Combined Ecological Slope Protection

At present, most studies on the effect of ecological slope reinforcement are single-factor analyzes, while the ecological protection method of composite vegetation with geogrids has multiple factors and multiple levels, which affect one another. Therefore, compared with the single-factor analysis, the orthogonal model test can quantitatively calculate the degree of influence of each factor on the stability of ecological slopes under different working conditions. This process can also obtain the optimal combination form to guide the application of joint ecological slope protection technology in practical engineering.


### 4.1. Establishment of the Orthogonal Model

Based on the existing research, taking the shrub reinforcement depth, arbor anchoring depth, arbor distribution form, and geogrid spacing as influencing factors, we designed an orthogonal model test with 4 factors and 3 levels, as shown in Table3. The distribution forms of trees include uniform (a), upper dense and lower sparse (b), and upper sparse and lower dense (c) distributions. Specific planting methods are shown in Figure 7.Table 3
Orthogonal design level of optimized parameters (unit: m).
LevelOptimization parametersShrub reinforcement depthArbor anchoring depthArbor distribution formGeogrid spacing10.21a0.420.42b0.630.63c0.8Figure 7
Schematic diagram of the tree planting pattern (unit: m).
(a)(b)(c)


### 4.2. Orthogonal Test Results

Nine combinations of 4 factors and 3 levels under different working conditions were numerically calculated, and the results are shown in Table4. The calculation results are analyzed in the range in Table 5. A larger extreme difference indicates that this factor has a greater influence on the slope stability, which can be used to obtain the degree of influence of different factors on the ecological slope stability under different working conditions and derive the optimal combination form, which is reported in Table 5.Table 4
Orthogonal test results of the slope safety factor.
Combining formsOptimization parametersSafety factor of slopeShrub reinforcement depthArbor anchoring depthArbor distribution formGeogrid spacingRainfall 0 hRainfall 24 h10.21a0.42.0751.52220.22b0.62.0941.52930.23c0.82.1001.56340.41b0.82.0631.53850.42c0.42.1011.55660.43a0.62.1091.58170.61c0.62.0781.56380.62a0.82.0881.56690.63b0.42.1051.613Table 5
Range analysis of the slope safety factor.
Rainfall duration (h)Average waterShrub reinforcement depthArbor anchoring depthArbor distribution formGeogrid spacing0K12.0902.0722.0912.094K22.0912.0942.0872.094K32.0902.1052.0932.084Range0.0010.0330.0060.01024K11.5381.5411.5561.564K21.5581.5501.5601.558K31.5811.5861.5611.556Range0.0430.0450.0050.008Under natural working conditions, trees have the largest extreme difference in anchorage depth, and shrubs have the smallest extreme difference in reinforcement depth. Therefore, the degree of influence of each factor on the stability of ecological slopes is as follows: depth of anchorage of trees > spacing of geogrid > form of distribution of trees > depth of reinforcement of shrubs. Meanwhile, it can be seen from Figure8(a), the slope safety coefficient first increases and subsequently decreases with increasing reinforcement depths; the slope safety increases with increasing arbor anchorage depths; the slope safety first increases and then stabilizes with the decrease of the geogrid spacing. The size of the slope safety coefficient under different forms of arbor distribution follows the following order: upper sparse and lower dense > uniform distribution > upper dense and lower sparse. The above result shows that the stability of shallow slopes can be improved with increasing reinforcement depth. However, when the influence depth of reinforced soil is too large, the intermediate connection force is insufficient, which increases the risk of sliding along shallow slopes and induces shallow slope instability. Increasing the anchoring depth and reducing the spacing of the geogrid can improve the reinforced anchoring effect of the reinforced material. However, when the spacing of the grids is reduced, the reinforcement effect will not improve. Considering the economic benefits of the project, the spacing of grating should be 0.6 m. In addition, the sensitive area of the slope safety factor is at the foot of the slope, and the arbor adopts the upper sparse and lower dense distribution, which can effectively play the role of anchoring. In summary, under natural working conditions, the optimal combination form is as follows: the reinforcement depth is 0.4 m, the anchorage depth is 3 m, the arbor is sparse on the top and dense at the bottom, and the grid spacing is 0.6 m.Figure 8
Curves of slope safety factor varying with the average level of various factors under different working conditions. (a) Rainfall 0 h. (b) Rainfall 24 h.
(a)(b)Under rainfall conditions, the degree of influence of each factor on the stability of ecological slopes is as follows: depth of anchorage of trees > depth of reinforcement of shrubs > spacing of geogrid > form of distribution of trees. The extreme difference in shrub reinforcement depth was greatly enhanced. Thus, the effect of shrub reinforcement on the hydraulic properties of the slope was considerably greater than its mechanical properties. Under joint ecological protection, shrub reinforcement can effectively reduce the hazard of rainfall on the slope stability and improve the hydraulic characteristics of the slope. It can be seen from Figure8(b) that slope safety factor, shrub reinforcement depth, arbor anchorage depth, and geogrid spacing are linearly related, which obviously differs from the trend of natural working conditions. A greater depth of shrub reinforcement corresponds to a greater degree of improvement in the hydraulic properties of the slope and better protection against rainfall hazards. In other words, when planting shrubs in rainfall areas, vegetation with more developed root systems should be used. Simultaneously, rainfall increases the tendency of relative displacement between the soil and the grid. This reinforcement effect can be effectively improved by reducing the grid spacing, so the grid spacing should be 0.4 m. Under different arbor distribution forms, the slope safety factor follows the following order: upper sparse and lower dense > upper dense and lower sparse > evenly distributed. The optimal reinforcement effect can be obtained with the upper sparse and lower dense planting method. The root system in the lower part more easily passes through the plastic zone of the slope and gives full play to the anchoring effect of the root system. However, under the action of rainfall, the displacement-sensitive area of the ecological slope is changed, and the maximum displacement is located in the middle and upper parts of the slope. Therefore, unlike the natural working condition, the planting method of upper density and lower sparsity can improve the connection force between the shallow soil and the stable layer in the middle and upper parts of the slope, prevents shallow damage to the slope, and is better than the uniform distribution. In summary, under rainfall conditions, the optimal combination form is as follows: reinforcement depth of 0.6 m, anchorage depth of 3 m, sparse arbor at the top, dense arbor at the bottom, and grid spacing of 0.4 m.The above analysis shows that, under natural and rainfall conditions, the degree of influence of various factors on the stability of ecological slopes is very different, and the obtained optimal combination form is also different. Therefore, for the combination forms of ecological slopes in different regions, the corresponding optimal combination forms should be formulated according to the local hydrogeological conditions.


### 4.3. Practical Engineering Application

We applied the optimal combination form obtained from the orthogonal model test to the corresponding actual project to verify the superiority of the combination.According to research on the reinforcement mode of composite vegetation with geogrids, the combined ecological protection method can be applied to the side slope project of the Yushi Expressway in Guizhou. Therefore, a geogrid was laid on the filling section of a subgrade side slope, and a vegetation composite consisting of shrubs and trees was planted on the slope. The grid spacing was 0.4 m, and trees were planted with sparse tops and dense bottoms. Figure9 shows a site picture of ecological protection in the bid section. By monitoring the finished slope project, we found that, in the early stage of planting, the reinforcement of the geogrid improved the slope stability. As the vegetation growth period progressed, the root system passed through the mesh of the geogrid. The vegetation and geogrid were interweaved and jointly reinforced and anchored. Hence, the reinforcement of composite vegetation and ecological benefits considerably improved. Under rainfall conditions, shrubs could effectively weaken splash erosion and reduce the scouring of slopes by rainwater. No obvious soil erosion occurred, and the local ecological environment was quickly restored, which achieved good reinforcement and demonstrates the ecological benefits of the project.Figure 9
Practical application of joint ecological protection.


## 4.1. Establishment of the Orthogonal Model

Based on the existing research, taking the shrub reinforcement depth, arbor anchoring depth, arbor distribution form, and geogrid spacing as influencing factors, we designed an orthogonal model test with 4 factors and 3 levels, as shown in Table3. The distribution forms of trees include uniform (a), upper dense and lower sparse (b), and upper sparse and lower dense (c) distributions. Specific planting methods are shown in Figure 7.Table 3
Orthogonal design level of optimized parameters (unit: m).
LevelOptimization parametersShrub reinforcement depthArbor anchoring depthArbor distribution formGeogrid spacing10.21a0.420.42b0.630.63c0.8Figure 7
Schematic diagram of the tree planting pattern (unit: m).
(a)(b)(c)


## 4.2. Orthogonal Test Results

Nine combinations of 4 factors and 3 levels under different working conditions were numerically calculated, and the results are shown in Table4. The calculation results are analyzed in the range in Table 5. A larger extreme difference indicates that this factor has a greater influence on the slope stability, which can be used to obtain the degree of influence of different factors on the ecological slope stability under different working conditions and derive the optimal combination form, which is reported in Table 5.Table 4
Orthogonal test results of the slope safety factor.
Combining formsOptimization parametersSafety factor of slopeShrub reinforcement depthArbor anchoring depthArbor distribution formGeogrid spacingRainfall 0 hRainfall 24 h10.21a0.42.0751.52220.22b0.62.0941.52930.23c0.82.1001.56340.41b0.82.0631.53850.42c0.42.1011.55660.43a0.62.1091.58170.61c0.62.0781.56380.62a0.82.0881.56690.63b0.42.1051.613Table 5
Range analysis of the slope safety factor.
Rainfall duration (h)Average waterShrub reinforcement depthArbor anchoring depthArbor distribution formGeogrid spacing0K12.0902.0722.0912.094K22.0912.0942.0872.094K32.0902.1052.0932.084Range0.0010.0330.0060.01024K11.5381.5411.5561.564K21.5581.5501.5601.558K31.5811.5861.5611.556Range0.0430.0450.0050.008Under natural working conditions, trees have the largest extreme difference in anchorage depth, and shrubs have the smallest extreme difference in reinforcement depth. Therefore, the degree of influence of each factor on the stability of ecological slopes is as follows: depth of anchorage of trees > spacing of geogrid > form of distribution of trees > depth of reinforcement of shrubs. Meanwhile, it can be seen from Figure8(a), the slope safety coefficient first increases and subsequently decreases with increasing reinforcement depths; the slope safety increases with increasing arbor anchorage depths; the slope safety first increases and then stabilizes with the decrease of the geogrid spacing. The size of the slope safety coefficient under different forms of arbor distribution follows the following order: upper sparse and lower dense > uniform distribution > upper dense and lower sparse. The above result shows that the stability of shallow slopes can be improved with increasing reinforcement depth. However, when the influence depth of reinforced soil is too large, the intermediate connection force is insufficient, which increases the risk of sliding along shallow slopes and induces shallow slope instability. Increasing the anchoring depth and reducing the spacing of the geogrid can improve the reinforced anchoring effect of the reinforced material. However, when the spacing of the grids is reduced, the reinforcement effect will not improve. Considering the economic benefits of the project, the spacing of grating should be 0.6 m. In addition, the sensitive area of the slope safety factor is at the foot of the slope, and the arbor adopts the upper sparse and lower dense distribution, which can effectively play the role of anchoring. In summary, under natural working conditions, the optimal combination form is as follows: the reinforcement depth is 0.4 m, the anchorage depth is 3 m, the arbor is sparse on the top and dense at the bottom, and the grid spacing is 0.6 m.Figure 8
Curves of slope safety factor varying with the average level of various factors under different working conditions. (a) Rainfall 0 h. (b) Rainfall 24 h.
(a)(b)Under rainfall conditions, the degree of influence of each factor on the stability of ecological slopes is as follows: depth of anchorage of trees > depth of reinforcement of shrubs > spacing of geogrid > form of distribution of trees. The extreme difference in shrub reinforcement depth was greatly enhanced. Thus, the effect of shrub reinforcement on the hydraulic properties of the slope was considerably greater than its mechanical properties. Under joint ecological protection, shrub reinforcement can effectively reduce the hazard of rainfall on the slope stability and improve the hydraulic characteristics of the slope. It can be seen from Figure8(b) that slope safety factor, shrub reinforcement depth, arbor anchorage depth, and geogrid spacing are linearly related, which obviously differs from the trend of natural working conditions. A greater depth of shrub reinforcement corresponds to a greater degree of improvement in the hydraulic properties of the slope and better protection against rainfall hazards. In other words, when planting shrubs in rainfall areas, vegetation with more developed root systems should be used. Simultaneously, rainfall increases the tendency of relative displacement between the soil and the grid. This reinforcement effect can be effectively improved by reducing the grid spacing, so the grid spacing should be 0.4 m. Under different arbor distribution forms, the slope safety factor follows the following order: upper sparse and lower dense > upper dense and lower sparse > evenly distributed. The optimal reinforcement effect can be obtained with the upper sparse and lower dense planting method. The root system in the lower part more easily passes through the plastic zone of the slope and gives full play to the anchoring effect of the root system. However, under the action of rainfall, the displacement-sensitive area of the ecological slope is changed, and the maximum displacement is located in the middle and upper parts of the slope. Therefore, unlike the natural working condition, the planting method of upper density and lower sparsity can improve the connection force between the shallow soil and the stable layer in the middle and upper parts of the slope, prevents shallow damage to the slope, and is better than the uniform distribution. In summary, under rainfall conditions, the optimal combination form is as follows: reinforcement depth of 0.6 m, anchorage depth of 3 m, sparse arbor at the top, dense arbor at the bottom, and grid spacing of 0.4 m.The above analysis shows that, under natural and rainfall conditions, the degree of influence of various factors on the stability of ecological slopes is very different, and the obtained optimal combination form is also different. Therefore, for the combination forms of ecological slopes in different regions, the corresponding optimal combination forms should be formulated according to the local hydrogeological conditions.


## 4.3. Practical Engineering Application

We applied the optimal combination form obtained from the orthogonal model test to the corresponding actual project to verify the superiority of the combination.According to research on the reinforcement mode of composite vegetation with geogrids, the combined ecological protection method can be applied to the side slope project of the Yushi Expressway in Guizhou. Therefore, a geogrid was laid on the filling section of a subgrade side slope, and a vegetation composite consisting of shrubs and trees was planted on the slope. The grid spacing was 0.4 m, and trees were planted with sparse tops and dense bottoms. Figure9 shows a site picture of ecological protection in the bid section. By monitoring the finished slope project, we found that, in the early stage of planting, the reinforcement of the geogrid improved the slope stability. As the vegetation growth period progressed, the root system passed through the mesh of the geogrid. The vegetation and geogrid were interweaved and jointly reinforced and anchored. Hence, the reinforcement of composite vegetation and ecological benefits considerably improved. Under rainfall conditions, shrubs could effectively weaken splash erosion and reduce the scouring of slopes by rainwater. No obvious soil erosion occurred, and the local ecological environment was quickly restored, which achieved good reinforcement and demonstrates the ecological benefits of the project.Figure 9
Practical application of joint ecological protection.


## 5. Conclusions

For the ecological slope protection technology of composite vegetation with a geogrid, theoretical analysis, numerical simulation, optimization design, and practical application were performed. The mechanical and hydraulic characteristics of the combined ecological slope were comprehensively evaluated. The impacts of the shrub shallow root reinforcement, arbor deep root anchoring, and combined geogrid reinforcement on the stability of the slope were quantified. The design of the combined ecological slope was optimized. All of this work will contribute to the design and application of ecological slope engineering. The study indicates that(1)
Based on the reinforced soil theory and unsaturated seepage theory, we demonstrated the feasibility of using composite vegetation of shrubs and trees to reinforce the slope with a geogrid.(2)
The ecological protection method of composite vegetation with geogrids is obviously better than the single reinforcement method. The maximum displacement of the ecological slope under rainfall conditions decreased by 82% compared to the original slope, and the overall stability increased by 35%.(3)
The shrub reinforcement has a much greater effect on the hydraulic characteristics of the slope than its mechanical characteristics. The anchoring depth of the arbor root system plays a key role in ensuring slope stability. The encryption of geogrid spacing can reduce the hazard of rainfall to the slope. According to the degree of influence of various factors on the stability of ecological slopes, the optimal combination form under different working conditions is obtained. The optimal combination form is applied to the actual project, and good reinforcement and ecological benefits have been achieved.(4)
Due to the diversity of root distribution morphology, this study simplified the arbor root system, which can only reflect the spatial relationship between the root system and the geogrid. Therefore, future research will examine how to more truly reflect the actual effect of the joint reinforcement of the root system and grid.



---

*Source: 8058009-2021-08-19.xml*