Significant climate change has increased rainfall intensity in some parts of the world. Recently, a catastrophic landslide occurred at Cimanggung village (Figure 1- Catastrophic landslide_Courtesy to PT. Geotechnical Engineering Consultant) due to significant rainfall intensity. Rainfall infiltration affects the unsaturated condition of the soil near the ground surface and it increases the pore-water pressure (uw) . As a result, it reduces the soil shear strength.
Based on the geological map of Cimanggung (Figure 2), it is known that most of the soil at the site consists of residual soil from volcanic rock (and hence referred as volcanic soil). Soil consists mostly of tuff, breccia, lava, and agglomerate. Groundwater table varies from 8.5m to 20m depth. Hence, it is important to consider the unsaturated soil mechanics to understand the mechanism of rainfall-induced slope failure in the analyses of slope stability.
(Source: Badan Geologi, Pusat vulkanologi dan Mitigasi Bencana Geologi)
Pore‑water pressure (uw) in the unsaturated region is usually negative as described in Figure 3. Under steady-state condition and during the periods where there are no evapotranspiration or infiltration, uw is hydrostatic.
However, the presence of evapotranspiration results in the subsequent decrease in pore-water pressure. Hence, no longer hydrostatic pore-water pressure is observed above the groundwater table. On the other hand, uw will decrease and it might become positive when there is an infiltration occurred at the ground surface. The depth at which the pore pressure is affected by infiltration is referred to as the depth of the wetting front. If the infiltration is significant enough, it is possible that a perched water table will occur in the unsaturated zone.
Negative pore-water pressure is usually referred to as matric suction (s) which is commonly described as:
"It is important to note that soil may also have osmotic suction and hence the total suction is equal to the summation between matric suction and osmotic suction. However, it is quite common to use soil suction as generic term to describe any of the matric suction, osmotic suction, or total suction. In modelling unsaturated slope, it is important to consider the change in permeability and shear strength due to the change in soil suction (Wijaya and Leong 2018) which is shown in Figure 4."
Permeability is commonly estimated from the so-called soil-water characteristic curve (SWCC) which is shown in Figure 5. The Y-axis of the SWCC can be either gravimetric water content (SWCC‑w), volumetric water content (SWCC‑
In MIDAS GTS NX, change in shear strength is considered in the numerical analyses using Bishop effective stress which is given as follow:
Where
1. The value of
2. The collapse phenomenon which was observed in the double oedometer test by Jennings and Burland (1962) invalidates the concept of effective stress as the inundation which decreases the value of effective stress causes the soil to decrease in volume while it is expected to increase in volume.
3. The
Thus, while the concept of effective stress can be useful in stability problems, it might not be suitable especially when dealing with collapsible soil. Thus, care must be taken when it is used.
👇 Watch the popular case study webinar
Unsaturated Slope Stability Analysis using GTS NX
In order to assess the validity of the model, and to ensure the stability of the unsaturated slope, it is important to conduct suction and moisture measurements. Suction measurement is usually carried by using a tensiometer while moisture measurement is usually carried by using moisture sensors such as frequency domain reflectometry (FDR) or time domain reflectometry (TDR). Figure 5 shows the typical TDR sensor and Figure 6 shows the typical arrangement for the moisture sensor and suction probe for slope monitoring.
Figure 7. Field Monitoring (Hamdany et al. 2021)
The problem in conducting a 3D slope model is in constructing the shape of the slope surface especially on the natural slope which has a complex geometry. MIDAS GTS NX has the so-called terrain geometry maker (TGM) which helps to convert contour into a surface as shown in Figure 8.
Figure 9 shows the 3D slope model where rainfall only occurred in the blue area. Instrumentation is placed on the observation point in order to validate the 3D model with field observation.
Figure 10 shows the verification of the 3D slope model. The change in pore pressure between MIDAS 3D matches well with the field measurement. The difference between the measured and modeled pore‑pressure is due to the fact that the input for SWCC usually ignores the hysteresis of SWCC. It is interesting to note that the phreatic surface is developed on 11/02/2020 due to the significant rainfall.
Figure 11. Phreatic surface
Modeling unsaturated slope stability requires a number of additional input parameters such as soil-water characteristic curve, permeability function of unsaturated soils, and rainfall/evaporation rate. On the other hand, the finite element model is the simplification of the actual conditions which are rarely homogeneous. Moreover, when the field condition is no longer can be represented by a 2D model, 3D analysis is required and hence further complicate the issues as 3D geometry is always challenging.
Thus, validation of the finite element model is important and can be carried by comparing the change in pore-water pressure or moisture in the model and from the field monitoring.
MIDAS GTS NX is a powerful tool that can be used to model 3D slope with ease due to its powerful terrain geometry maker (TGM) and it allows coupled analysis between stress and seepage in 3D. The pore-pressure responses from MIDAS GTS NX also appear to agree with the field measurement.