Skip to main content Skip to main navigation menu Skip to site footer

Modeling of Root-reinforced Soil Slope under Rainfall Condition


Download data is not yet available.


This paper presents the results of geotechnical centrifuge and numerical modeling study of root-reinforced systems on soil slopes. The centrifuge models were designed to simulate a soil slope reinforced by vetiver grass root system commonly grows in Southeast Asian countries. Unreinforced and root-reinforced soil slope models were subjected to heavy rainfall using a rainfall simulator designed for a geotechnical centrifuge. Results of the study show that the root system causes a reduction of rainwater infiltration rate, delay in the response of groundwater table, and increasing in soil shear strength. The results were validated by comparing centrifuge modeling test results with numerical modeling analyses based on limit equilibrium and finite element methods. The transient seepage analysis results were employed in the slope stability analysis. The numerical analysis shows a good agreement in failure mechanisms with the observations from centrifuge modeling tests.

1. Introduction

Landslides induced by rainfall are one of most serious and common disasters worldwide. In recent years, natural slope instability has increased especially in the tropical monsoon zone such as Southeast Asian countries including Thailand. There are several factors that could cause the natural slope failure, such as geological activity, hydrological influence and human interference, but rainfall is usually the main factor. This is because the change of pore water pressures can be induced by ground water rising or infiltration from the rain water. The increase in pore water pressure shall reduce the matric suction and hence decrease the shear strength of soil in unsaturated zone. The physical process of rainfall infiltration into the ground and its seepage through the soil layers have been studied by hydrogeologists, soil scientist and geotechnical engineers [ 1 ].

Vegetation along a slope plays a major role in preventing failure as it can provide both hydrological and mechanical benefits to the stability of natural soil slopes. The role of vegetation on slope stability has been defined by Greenway (19 78) [ 2 ], Coppin and Richard (1990) [ 3 ], Wu (1995) [ 4 ] and Gray and Sotir (1996) [ 5 ]. For the hydrological process, changing the soil moisture regime and drawing the water from the soil via evapotranspiration [ 6 ] could increase the soil suction. For mechanical process, the roots of vegetation could enhance slope stability by increasing the shear strength of the soil [ 5 , 7 ]. The benefits of vegetation along natural soil slopes are to prevent shallow failure and soil surface erosion. However, research to understand fundamental mechanisms of root-reinforced soil slopes is limited, particularly for the study of root reinforcement effects on slope stability.

Vetiver grass ( Vetiveria nemoralis A. Camus ), has been implemented to prevent soil erosion, water runoff, and infiltration by the World Bank in the 1980s [ 8 ]. In Thailand, the Chaipattana Foundation and the Office of the Royal Development Projects Board has promoted the use of vetiver grass for soil and water conservation for many royal projects. Vetiver grass is fast growing, cost effective, and requires low maintenance. Chinapan et al. (1 997) [ 9 ] observed the length of vetiver roots to be up to 2-3.5 m. The vetiver roots can penetrate deep into the ground to form a net-like barrier capable of filtering silt and preserving top soil. In Thailand, slope failure due to prolonged and heavy rainfall typically occurs at shallow depth of approximately 1-1.5 m from the ground surface. As a result, the use of the vetiver grass root system as slope reinforcement is suitable for preventing slope failure. Root area ratio defined as the ratio between the cross-sectional area of root in the shear plane area and the total area of shear plane, plays an important role contributing to the increase of soil shear strength [5]. Eab et al. (2015) [ 10 ] correlated the root area ratio of the vetiver roots with the increase of soil shear strength. The root area ratios of the 4-6 month growth group vetiver were in the range of 2.44 – 4.56%. The shear strength of soil with vetiver root system was measured by direct shear tests. Results showed that the cohesion of the reinforced soil was increased approximately 11 kPa compared to unreinforced soil.

Centrifuge modeling was implemented in this study because the centrifuge allows modeling of soil slopes with full-scale stress conditions. Few researchers have investigated the mechanism of vegetation for slope stability through the use of a centrifuge. These models simulated the slope failure by continuously raising the groundwater table (e.g. Sonnenberg et al., 20 10, 2012) [ 11 , 12 ]. Some previous researches have investigated the mechanism of rainfall-induced slope failure based on laboratory tests and filed monitoring [ 13 ]. However, it is difficult to model the rainfall in a centrifuge. Recently, centrifuge modeling has been used to model rainfall induced landslides to provide a better understanding of the triggering mechanisms of landslides due to rainfall [ 14 , 15 , 16 .] Eab et al. (2015) [ 10 ] focused on the application of root-reinforced system to slope stabilization using a rainfall simulator in the centrifuge. In addition, Ng et al., (2014, 2015) [ 17 , 18 ] and Kamchoom et al. (2014) [ 19 ] studied the effects of root on transpiration-induced soil suction using centrifuge model.

This paper presents the results of centrifuge and numerical modeling of root-reinforced systems on soil slopes. The results of centrifuge tests carried out by Eab et al. (2014, 2015) [ 16 , 10 ] were validated by comparing with numerical modeling simulations. The pore water pressure generation due to transient seepage was used for the slope stability analysis. Input parameters used for the numerical analysis, including volumetric water content, hydraulic conductivity, and shear strength parameters, were obtained from laboratory testing results. Discussions of results from centrifuge and numerical modeling are provided in this paper.

2. Centrifuge Model Tests

Centrifuge modelling is a physical model test and is now widely used in many geotechnical applications. Since slope stability is a gravity-dependent problem, the major advantage of using centrifuge modelling is to enable researchers to test reduced-scale physical models at the correct stress level by increasing g-level [ 20 ]. The centrifuge model tests were conducted at Tokyo Institute of Technology in Japan. The tests were performed at a centrifugal acceleration of 50g. To simulate a heavy rainfall in the centrifuge, a pneumatic nozzle system was designed to spray water over the model slope. Figure 1 shows the schematic of the rainfall simulation system for this study. A brief summary of the experimental details is explained below. More detailed information of centrifuge model tests can be found in Eab et al. ( 2014, 2015) [ 16 , 10 ].

Fig 1 Schematic centrifuge model with a rainfall simulator system (After Eab et al, 2014) [ 16 ]

2.1 Material Properties

The Japan Geotecnical Society (JGS) Standard fine sand called Edosaki Sand was used in the model slopes. The engineering properties of the Edosaki sand are summarized in Table 1. The polyester fibers with approximately 39μm in diameter and 10mm in length were used to model the root fibers. In the model tests, the 2% by mass of polyester fibers were randomly mixed with the sandy soil to model the vegetated slope surface. The 2% fiber mixing by mass is approximately 7% by volume, which is slightly higher than the root area ratio of vetiver root at 6-months growth reported by Eab et al. (2015) [ 10 ]. It is noted that the typical values of the root area ratio of 3 – 5% for small vegetation are reported by Gray and Sotir (1996) [5]. The engineering properties of soil mixed with 2% of polyester fiber are summarized in Table 1.

Table 1 Properties of compacted soils (After Eab et al., 2014 [ 16 ]

2.2 Slope Model Tests

The slope model was constructed inside a 450mm×150mm×270mm steel box. The model design is schematically illustrated in Fig. 1. Three model cases, i.e., Case 1: without reinforcement, Case 2: with 20-mm fiber reinforcement and Case 3: with 40-mm fiber reinforcement, were conducted and summarized in Table 2.

Table 2 Summary of centrifuge tests (After Eab et al., 2014)[ 16 ]

2.3 Results of Centrifuge Tests

Results showed that the unreinforced slope failed due to the rising groundwater table. It was observed that the failure started to deform from the toe and progressively propagate upward to the top of the slope. Both root-reinforced slopes deformed but did not collapse. The tests were terminated when the failure of the slope was observed or the pressure head at the toe slope reached the critical value. The centrifuge test results are summarized in Table 2 . Results indicated that the existence of root fibers helps increase the soil strength and reduce the infiltration rate of rainfall into the ground, resulting in a delayed response in the groundwater table. In other words, roots can delay slope failure and slow movement of a soil slope could provide a first early warning of failure. More detailed discussions of the centrifuge results can be found in Eab et al. (2014) [ 16 ].

3. Slope Stability Analysis

Numerical analysis is commonly carried out to validate the results of physical modeling tests and to better understand the effect of suction loss due to rainfall infiltration on slope stability. In this study, the three cases similar to the centrifuge tests were modeled. Each model was performed by using transient analysis to define the pore water pressure change in the soil slope. Based on the pore water pressure information, the factor of safety of the soil slope was calculated using the limit equilibrium methods.

3.1 Seepage and Slope Stability Modeling

Two-dimensional seepage and slope stability analyses were performed using GEO-SLOPE software [ 22 ]. Simplified homogenous slope profiles for each soil layer were used in the study. For a transient analysis in seepage, the soil-water characteristic curve (SWCC) and the hydraulic conductivity function are required input parameters. The estimation method to determine the SWCC proposed by Aubertin et al. (2003) [ 23 ] was chosen in this analysis. The proposed model assumes that soil-water characteristic results from the combined effect of capillary and adhesion forces. The SWCC can be evaluated from the negative pore water pressure and material property information such as particle-size distribution. In addition, the hydraulic conductivity function can also be estimated using several predictive methods that utilize either a grain-size distribution or a measured volumetric water content function and the saturated hydraulic conductivity (Ksat). In SEEP/W, a built-in predictive method can be used to estimate the hydraulic conductivity function once the volumetric water content function and a Ksat value have been specified. More details of the assumptions of the predictive methods to estimate the SWCC and the hydraulic conductivity function can be seen in the SEEP/W manual [ 22 ]. Figures 2(a) and 2(b) present the SWCC and the hydraulic conductivity function used in this study.

The initial total head was defined by specifying the initial water table position at the base of the toe slope. The negative pore water pressure in the unsaturated zone above the water table varies in a linear manner from zero at the base of the toe slope to the maximum value at the surface. For the rainfall simulation, the flux boundary, q equal to the desired rainfall intensity of 20 mm/hr was applied to the surface of the slope. The potential seepage face was also applied at the toe of the slope model. Figure 3 presents the geometry boundary and the initial water pressure conditions for the seepage and slope stability analyses. The pore water pressures were calculated for every time step at each node of the finite-element mesh. The pore water pressure results from SEEP/W were then incorporated into the SLOPE/W model. Slope stability analyses were performed using Bishop’s simplified and Morgenstern-Price methods. The input parameters from the direct shear tests shown in Table 1 were used for the slope stability analyses.

Fig 2 Input soil water characteristic curve and hydraulic conductivity for seepage analysis

Fig 3 Boundary and initial conditions for seepage and stability analyses

4. Results and Discussion

4.1 Pore Water Pressure Distribution

Figure 4 shows the pore water pressure contour distribution in the soil slope simulated by the seepage analysis for the three different cases. Results showed that during rainfall, pore water pressure in the soil slope initially started from negative and then moved gradually towards positive values. Figure 4(a) shows the result of the unreinforced case, where the pore water pressure increased by 10kPa at the toe of the slope while at the top slope pressure remained negative after an 18hr duration of rainfall. For the case of 1-m depth of vegetation, shown in Fig. 4(b) , the large amount of rainwater infiltrated into the slope and changed the pore water pressure from negative to positive value of 20kPa after 35hr. For the case of 2-m depth of vegetation, shown in Fig. 4(c) , pore water pressure increased up to 25kPa after 40hr. Based on these results, it is observed that pore water pressure increased to the critical values when the slope started to fail or deform significantly.

4.2 Slip Surface for Factor of Safety

Figure 5 shows the potential slip surface and ground water table simulated by SLOPE/W. Due to the increase of cohesion by the root reinforcement, the slip surface for the root-reinforced slope was mostly developed near the bedrock, especially for the case of 2-m depth of vegetation ( Fig. 5(c) ), in comparison with the slip surface of the unreinforced slope that was developed at shallower depth ( Fig. 5(a) ). The reason of this occurs is because the shear strength of the soil and root fiber at the top layer is higher than the soil only at the bottom layer. In addition, it was found that the factor of safety (FS) is slightly increased for the reinforced soil slopes as summarized in Table 3 . According to the results, the shear strength of unsaturated zone decreased and consequently the FS of the slope reduced during the rainfall. Comparing the FS at the same duration of rainfall, the FS increases with the depth of roots for the root-reinforced cases.

Table 3 Factors of safety from numerical analysis

Fig 4 Results pore water pressure change at the end of test.

Fig 5 Results slip surface at the end of test.

Fig 6 Comparison of pore water pressure change between centrifuge test and numerical analysis

4.3 Comparison of Centrifuge and Numerical Modeling Results

Fig ure 6 shows a comparison of pore water pressure between the centrifuge test and numerical analysis for the unreinforced case, 1-m depth of root case, and 2-m depth of root case. Pore water pressure measured by PWP 13 and PWP 14, which are located near the toe of the slope were plotted with time during the continuous rainfall. The increase of the pore water pressure was used to indicate the triggering time for slope failure as summarized in Table 2. Buildup of pore water pressure for the unreinforced slope is faster than for the root-reinforced slopes. For the unreinforced case, the pore water pressures measured from the centrifuge are much higher than the numerical values. As a result, the factor of safety calculated from the numerical model is higher than expected. This might be due to the fact that the numerical analysis results rely on the input parameters and boundary conditions. In this analysis, assumptions of the SWCC and hydraulic conductivity were made. Despite this, the comparisons show a reasonable agreement between the centrifuge and numerical results with the exception of the negative value measured at PWP 14. This is due to the limitation of pore water pressure transducer installed in the centrifuge test which cannot accurately measure negative pore water pressure. Based on this study, results from both centrifuge and numerical modeling provide good insight into the mechanisms of vegetated root reinforced slope during rainfall. The root-reinforced system delays the rising of the groundwater table and increases the factor of safety of the slope.

5. Conclusions

This research evaluated the effects of vegetation root on slope stability. A series of centrifuge model tests on slopes reinforced by the root model were conducted to illustrate the advantage of the vegetation on the rainfall-induced shallow failure. The results of centrifuge tests were verified by the numerical simulations. Conclusions from this study are as follows.

  1. The centrifuge model with rainfall simulator was employed to illustrate the slope instability due to heavy rainfall. The slope failure is triggered by rising of the water table and starts around the toe of the slope. The increase in pore water pressures from water table causes the decrease in the effective stress and shear strength of the soil.

  2. The polyester fiber mixed soil can be reasonably used to model the vegetated soil slope in the centrifuge. The results showed that roots in the surface of a slope tie the soil particles together and prevent the formation of cracks on the slope surface. The reinforcement efficiency increases with the root depth and the thicker reinforced zone reduces the slope deformation. In addition, the roots can delay the infiltration of rainfall into the ground to delay the rising of the groundwater table. The delay of slope failure or movement could provide an early warning for evacuation.

  3. Results from numerical analyses show an agreement with the results observed in centrifuge model tests. The factor of safety of the slope increases as the depth of root reinforcement increases. The mixture of soil and fiber roots increase the shear strength of soil, hence it provides overall improvement of the slope stability.


This research was supported by the Thailand Research Fund under the TRF Research Scholar Award Grant No. RSA-5880023. The work was also carried out under the research and development project on landslide prevention and protection according to Royal Initiatives of the Chaipattana Foundation. The authors would like to thank Dr. Inthuorn Sasanakul for her comments to improve the manuscript.


  1. A numerical investigation of the stability of unsaturated soil slopes subjected to transient seepage Ng CWW, Shi Q. Computers and Geotechnics.1998-jan;:1-28.
  2. Stress Analysis and Slope Stability Developments in Geotechnical Engineering.1978;:173-212.
  3. Bio-engineering: the use of vegetation in civil engineering Freer R. Construction and Building Materials.1991-mar;:23-26.
  5. Biotechnical and Soil Bioengineering Slope Stabilization: A Practical Guide for Erosion Control. Lal R. Soil Science.1998-jan;:83-85.
  7. Strength of tree roots and landslides on Prince of Wales Island, Alaska Wu TienH, McKinnell III WilliamP, Swanston DouglasN. Canadian Geotechnical Journal.1979-feb;:19-33.
  8. Soil Moisture Conservation Shaxson TF. Conservation Agriculture.2003;:317-326.
  9. Application of airborne gamma-ray imagery to assist soil survey: A case study from Thailand Moonjun Ruamporn, Shrestha DhrubaPikha, Jetten VictorG, Ruitenbeek FrankJAvan. Geoderma.2017-mar;:196-212.
  10. Laboratory and modelling investigation of root-reinforced system for slope stabilisation Eab KrengHav, Likitlersuang Suched, Takahashi Akihiro. Soils and Foundations.2015-oct;:1270-1281.
  11. Centrifuge modelling of soil slopes reinforced with vegetation Sonnenberg R, Bransby MF, Hallett PD, Bengough AG, Mickovski SB, Davies MCR. Canadian Geotechnical Journal.2010-dec;:1415-1430.
  12. Centrifuge modelling of soil slopes containing model plant roots Sonnenberg R, Bransby MF, Bengough AG, Hallett PD, Davies MCR. Canadian Geotechnical Journal.2012-jan;:1-17.
  13. Causative Mechanisms of Rainfall-Induced Fill Slope Failures Chen H, Lee CF, Law KT. Journal of Geotechnical and Geoenvironmental Engineering.2004-jun;:593-602.
  14. Centrifuge Model Simulations of Rainfall-Induced Slope Instability Ling Henry, Ling HoeI. Journal of Geotechnical and Geoenvironmental Engineering.2012-sep;:1151-1157.
  15. On the seepage-induced failure of vegetation-stabilised slopes Takahashi A, Nakamura K, Liktilersuang S. ICPMG2014 – Physical Modelling in Geotechnics.2013;:1233-1239.
  16. Centrifuge modelling of root-reinforced soil slope subjected to rainfall infiltration Eab KH, Takahashi A, Likitlersuang S. Géotechnique Letters.2014-jul;:211-216.
  17. A Novel Root System for Simulating Transpiration-Induced Soil Suction in Centrifuge Ng CharlesWangWai, Leung AnthonyKwan, Kamchoom Viroon, Garg Ankit. Geotechnical Testing Journal.2014-jun.
  18. Centrifuge modelling of the effects of root geometry on transpiration-induced suction and stability of vegetated slopes Ng CWW, Kamchoom V, Leung AK. Landslides.2015-oct;:925-938.
  19. Effects of root geometry and transpiration on pull-out resistance Kamchoom V, Leung AK, Ng CWW. Géotechnique Letters.2014-sep;:330-336.
  20. Geotechnical Centrifuge Technology TAYLOR RN. .1988.
  21. An objective classification of rainfall events on the basis of their convective features: application to rainfall intensity in the northeast of spain Llasat Mar�a-Carmen. International Journal of Climatology.2001;:1385-1400.
  22. Monitoring and evaluating of slope stability for setting out of critical limit at slope stability radar Kumar Ajay, Rathee Ritika. International Journal of Geo-Engineering.2017-sep.
  23. A model to predict the water retention curve from basic geotechnical properties Aubertin M, Mbonimpa M, Bussière B, Chapuis RP. Canadian Geotechnical Journal.2003-dec;:1104-1122.
View article in other formats
Author Biographies
Suched Likitlersuang

Geotechnical Research Unit, Department of Civil Engineering, Faculty of Engineering, Chulalongkorn University, 254 Phayathai Road, Wang-Mai, Pathumwan, Bangkok 10330, Thailand

Akihiro Takahashi

Department of Civil Engineering, Tokyo Institute of Technology, M1-3, 2-12-1, Ookayama Meguro-ku Tokyo 152-8552, Japan

Kreng Hav Eab

Geotechnical Research Unit, Department of Civil Engineering, Faculty of Engineering, Chulalongkorn University, 254 Phayathai Road, Wang-Mai, Pathumwan, Bangkok 10330, Thailand

Vol 21 No 3, Jun 15, 2017
How to Cite
S. Likitlersuang, A. Takahashi, and K. Eab, “Modeling of Root-reinforced Soil Slope under Rainfall Condition”, Eng. J., vol. 21, no. 3, pp. 123-132, Jun. 2017.

Authors who publish with Engineering Journal agree to transfer all copyright rights in and to the above work to the Engineering Journal (EJ)'s Editorial Board so that EJ's Editorial Board shall have the right to publish the work for nonprofit use in any media or form. In return, authors retain: (1) all proprietary rights other than copyright; (2) re-use of all or part of the above paper in their other work; (3) right to reproduce or authorize others to reproduce the above paper for authors' personal use or for company use if the source and EJ's copyright notice is indicated, and if the reproduction is not made for the purpose of sale.

Article Statistics
Total PDF downloads: 614
Total XML downloads: 40