1. Introduction
Landslide mobility is often connected to the water-table fluctuations owing to weather, and therefore, it is distinguished by an alternation of stages of movement and inactivity. Particularly, a reactivation of the landslide or an acceleration of the motion (if it is moving) can be caused by a raising in the water table throughout rainy periods. On the other hand, a decreasing in the groundwater table occurring in arid periods causes a lowering of the displacement rate of the unstable soil mass, leading eventually to the stop of the motion. According to the classification system of Cruden and Varnes [
1], these landslides are categorized as very slow-moving landslides and are characterized by a displacement rate of approximately a few centimeters per year. A further common characteristic of such landslides is that the soil strains are usually located inside a narrow shear band situated at the base of the unstable soil mass in which the shear strength of the soil is under residual conditions owing to the important gained deformations [
2]. Contrarily, the landslide body is subjected to low level of deformation and is generally outlined by horizontal displacements that are substantially steady along the depth. In other words, the unstable soil mass practically slides rigidly.
In engineering practice, the problem is commonly dealt with by separately analyzing the changes in groundwater regime and slope stability by means of an uncoupled approach. In particular, groundwater pressures at the level of the slip surface due to rain infiltration are firstly evaluated, and the latter are afterwards used in a limit equilibrium approach to calculate a slope-safety factor. However, a rational forecasting of the landslide mobility is more useful from a practical point of view than a safety factor calculation, especially when some structures are located on the slope. In these circumstances, the designer should assess whether the calculated displacements exceed the tolerable values established to safeguard the serviceability of the existing structures or prevent their damage [
3,
4]. Advanced numerical techniques can undeniably accomplish this target, providing a good comprehension of the complicated deformation processes that characterize the slope behaviour [
5,
6,
7,
8,
9,
10,
11]. Since the soil viscosity exerts a fundamental role in slow-moving landslides [
12,
13,
14,
15,
16,
17,
18], a constitutive model based on the elasto-viscoplastic theory is generally used when performing such a type of simulation [
19,
20,
21,
22]. This may sensibly increase the computational costs as well as need many input parameters that are often of difficult experimental evaluation. Therefore, the usage of numerical techniques might not be rational when dealing with routine applications.
Several methods of practical interest were also published in the literature [
23,
24,
25,
26,
27,
28,
29,
30,
31,
32,
33]. Such methods are generally able to calculate the landslide velocity when groundwater level measurements are available. Therefore, their use is unsuitable for predictive purposes. To overcome this limitation, Conte et al. [
34] developed a simple method based on the infinite slope model to relate landslide mobility to rainfall. In this paper, the approach originally introduced by Conte et al. [
34] is extended to more complex geometries. Specifically, it is assumed that the unstable soil mass slides on a bi-planar slip surface and can be schematized by means of a two-wedges model. As highlighted by Chowdhury [
35], such a schematization may reasonably represent several real situations. Using some analytical solutions derived in this study and rainfall recorded during a sufficiently long time period, the groundwater level fluctuations are evaluated and the landslide velocity predicted. The proposed solution is simple to use and necessitates few data as input, most of which are determined by means of traditional geotechnical tests, whereas the remaining ones are evaluated using a calibration procedure. After this calibration, the method may be adopted to perform a prediction of future displacements of a given landslide body based on weather forecast.
2. Method of Analysis
The proposed method uses a two-wedges model to schematize the landslide body. The wedges are separated by an ideal vertical wall and can slide on a bi-planar slip surface, making angles β
1 and β
2 with the horizontal plane (
Figure 1). It is assumed that seepage takes place with flow lines parallel to the slope and a stationary water table located at a depth
zw below the ground surface. The water level is subjected to changes
h(
t) due to rain infiltration.
By assuming that groundwater level oscillates simultaneously with rainfall,
h(
t) might be calculated as follows [
34]:
where
ho defines the position of water table at the beginning of the considered period of observation (detected from the stationary position),
N is the number of rain events,
is a model parameter depending on the soil permeability,
where
n and
Sr are the porosity and the saturation degree of the portion of soil located above the water table, respectively, and
hi is the volume of water that can infiltrate in a unity area of the ground surface during the
i-th event, whose beginning is defined by
toi. Following [
36], it is assumed
, where
p is the potential infiltration rate and represents the maximum water volume that can permeate through a unity area of soil. In the present study, the values of
n,
Sr, and
p are assumed to be constant for the sake of simplicity.
The forces acting on wedge
i (with
i = 1, 2) are weight
Wi along with the forces developed on the failure surface. The latter are the shear force
Ti, effective normal force
, the force due to water pressure
Ui, and the viscous force
Fi activated at the base of the unstable soil mass during motion (
Figure 2).
Additionally, in correspondence of the vertical interface between the blocks, the resultant of pore water pressure
U and the force indicated by
act. The latter is inclined of
, with respect to the horizontal plane (with
, where
is the friction angle corresponding to the material forming the landslide body). Forces
Wi,
U, and
Ui can be calculated using the following equations:
and
in which
γ is the unit weight of the soil (changes of
γ are neglected for the sake of simplicity),
Ai represents the volume per unit length of the considered block,
hw defines the length of the vertical interface between the two wedges included from the stationary water table to the failure surface,
h(
t) is provided by Equation (1), and
ϑi and
ωi are the angles shown in
Figure 2. The viscous forces activated at the base of the landslide body during motion are expressed as:
where
μs is a coefficient of viscosity of the soil in the shear zone, the expression of which is provided in
Appendix A, and
vi(
t) is the velocity of the i-th block at time
t in the movement direction. Forces
T1,
T2,
,
, and
are unknown. However, the following expressions relating these forces can be written using the Mohr–Coulomb failure criterion:
where
L1 and
L2 define the extension of the failure surface involving the considered wedges, and
and
are the shear strength parameters of the soil in the shear band (intercept cohesion and friction angle, respectively) that are generally under residual condition for the considered type of landslides.
Since the stability of slopes is fervently influenced by forces
U1 and
U2 that, in turn, depend on force
U, a critical value of the latter force,
Ucrit, can be calculated by solving the equilibrium equations of the wedges under a condition of incipient failure. The resulting expression for
Ucrit takes the following form:
where
Force
Ucrit is useful to establish whether the landslide body moves or is at rest. Specifically, no motion occurs if
. On the contrary, the landslide is activated when
. In the latter case, it is convenient to invoke the motion equations of the wedges. Considering the direction of the failure surface and the normal one, the equations of motion take the following form for the first wedge:
and for the second wedge:
where
g is the gravity acceleration, and
v1(
t) and
v2(
t) are the velocities of wedge 1 and wedge 2, respectively, referring to the direction of the failure surface. In this paper, it is assumed that
v1(
t) and
v2(
t) are characterized by the same component in the horizontal direction, i.e.,
As shown in
Figure 3, the blocks slide rigidly on the slip surface, avoiding separation and overlapping at the vertical wall. However, the soil mass indicated as AA’A” in
Figure 3 passes over the slip surface, breaching the congruence. Such a violation could be considered acceptable from an engineering point of view if the landslide displacements are small in comparison with the size of the unstable soil mass [
37]. Therefore, after the computation is completed, it necessary to check if the area AA’A” can be neglected in comparison with that of the unstable soil mass. In this connection, area AA’A” can be calculated as
, where
d is the horizontal displacement of the unstable soil mass (
Figure 3).
Lastly, substituting Equations (6), (7), and (13) in Equations (9)–(12), the following differential equation is attained:
where
v(
t) coincides with
v1(
t) (it is related to
v2(
t) by Equation (13)), and the expressions of the coefficients
λ and
χ are:
where
and
Equation (14) can be solved by employing the Duhamel’s theorem, imposing the initial condition that the unstable soil mass is initially quiescent, i.e.,
v = 0 at
t = 0. Using a procedure similar to that adopted by Conte and Troncone [
25], the following analytical solution is derived:
In which
τ is a fictitious time employed to solve the integral,
ts is the time when the motion starts, and
is the solution of Equation (14) with
at any time. This solution is
After replacing
t by (
t-
τ) in Equation (19) and making the derivate of
v(
t) with respect to
t, Equation (18) assumes the following expression:
Considering that force
U depends on the groundwater level, a solution for
v(
t) can be obtained when measurements of
h(
t) are available. Alternatively, Equation (1) may be used to this end. Following Conte and Troncone [
25],
U(
t) could be conveniently expanded in a finite number of
K harmonic functions using the Fourier series:
where
ωj is the frequency of the
j-th function:
in which
T is the period of
U(
t) that needs to be taken larger than the time interval considered in the simulation, and
Ao,
Aj, and
Bj are the series amplitudes calculated using the following equations:
Closed-form expressions for
Ao,
Aj, and
Bj may be obtained when
U(
t) is expressed as a series of linear functions (Conte and Troncone [
25]). Finally, the analytical expression of
v(
t) is
where
with
Velocity
v(
t) is computed if
. Positive values of
v(
t) corresponding to movements in the downhill direction are only considered. On the contrary, velocity is nil when
. To the end, displacements of the unstable soil mass may be calculated by numerically integrating Equation (24). A similar approach was also proposed by Troncone et al. [
38] to study a different problem concerning the effects of stabilizing piles on rainfall-induced landslide mobility.
Summarizing, the input data necessary to employ the presented approach are soil unit weight,
γ; angle of shearing resistance of the landslide body,
; and the residual cohesion and friction angle of the soil belonging to the failure surface,
and
. They can be determined by means of traditional geotechnical tests. The other parameters are the potential infiltration rate,
p; the model parameters,
A and
,; and the coefficient
μs. The latter are more difficult to be obtained experimentally, but they may be assessed by performing a calibration procedure based on the available recordings of water-table and soil displacements. In particular,
p,
A, and
are determined by fitting the available measurements of water-table changes with the theoretical results obtained using Equation (1). Likewise,
μs is estimated by matching the measured landslide displacements with those obtained from Equation (24). To this purpose, the calibration method described in Conte et al. [
34] is included in the proposed approach. After this calibration is performed, the proposed approach may be used to carry out an approximate assessment of future displacement of the unstable soil mass directly from weather forecast.
3. Application of the Proposed Approach
The approach presented in this study is employed in the present section to analyze the mobility of a landslide located near Palermo, in Sicily (Italy), which is named Cerda landslide (
Figure 4a). A detailed documentation of this landslide can be found in Rosone et al. [
39]. The Cerda landslide was triggered by an earthquake in 2002 and subsequently was seasonally reactivated owing to rainfall-induced fluctuations of water table. From the morphological viewpoint, the unstable soil mass was subdivided in three different bodies characterized by a different mobility (
Figure 4b), as documented by Rosone et al. [
39]. In the present study, attention is focused on the body named A (
Figure 4b), which is characterized by an inclination angle of the ground surface of about 8°. The subsoil consists of a layer of weathered silty-clay overlying a formation of varicoloured clay (
Figure 5a). The slip surface is located at the base of the upper layer and is characterized by a maximum depth of approximately 20 m below the ground surface. For the purposes of the present study, the unstable soil mass is schematized by two wedges bounded at their base by a bi-planar failure surface with β
1 = 5° and β
2 = 12° (
Figure 5b). In autumn 2008, some piezometers and inclinometers along with a pluviometer were set up in the site affected by movements [
39]. Rain recordings and groundwater level and horizontal displacements measurements are available from September 2009 to May 2011. These data are shown in
Figure 6a–c from which it is clear that a noticeable simultaneity exists among rainfall, water-table fluctuations, and horizontal displacements.
Figure 7 shows the available horizontal displacement profiles of inclinometer BH4 located at the toe of the landslide body (
Figure 4b), as documented by Rosone et al. [
39]. It is evident from
Figure 7 that the soil deformations are essentially located within a well-defined zone situated below the layer of weathered silty-clay. This observation is consistent with the hypotheses on which the proposed approach is based. Shear strength parameters of the weathered silty-clay were obtained by Rosone et al. [
39] by performing a large number of direct shear tests. Peak intercept cohesion
and shearing resistance angle
fall into the range 5–50 kPa and 22°–23°, respectively. An intercept cohesion
equal to zero, and a friction angle
were found at residual conditions. Considering that the Cerda landslide is periodically activated by rainfall, the residual strength parameters are employed in this work to characterize the soil shear strength in the shear zone. On the contrary, the peak shear strength parameters are employed for the soils forming the landslide body, considering that these soils are affected by moderate deformations. Therefore, it is assumed
for the soil in correspondence of the interface between the two blocks (
Figure 2). Soil unit weight is
γ = 20.5 kN/m
3 [
40].
The remaining parameters involved in the proposed approach are calibrated using the above-mentioned procedure by matching the available recordings of water-table and soil displacement with the outcomes obtained using the proposed approach. In order to perform a reliable prediction of the landslide mobility, this calibration procedure should refer to a sufficiently prolonged time interval. To this purpose, the period from September 2009 to August 2010 is considered in the present study and from which the following parameters are obtained using Equation (1) along with rain recordings and groundwater level measurements:
p = 161.9 mm/month,
A = 0.27, and
= 2.9 month
−1. In addition, a coefficient of viscosity
kN·month/m
2 is evaluated on the basis of the soil displacements measured in the same time interval. The calibration results are reported in
Figure 8a–c. It is clear from
Figure 8b that the groundwater level calculated using the above-mentioned parameters well resembles the measured one. The critical level
zw,crit (associated with force
Ucrit) is also indicated in this figure for completeness. A reasonably fine accord between theoretical and observed soil displacements is also gained (
Figure 8c). As a further requirement of the proposed method, area AA’A” (
Figure 3) is calculated in the order of 10
−3 m
2. This value can be neglected in comparison with the total area of the entire unstable soil mass (12,000 m
2).
Once the calibration procedure is completed, the method is used to predict the displacement of the unstable soil mass in the subsequent period from August 2010 to May 2011. The beginning of the prediction is defined by a dashed line in
Figure 9. A comparison between the calculated displacements and the recorded ones is reported in
Figure 9b. The results of this comparison corroborate the predictive capacity of the approach proposed in the present study.
For the sake of completeness, a comparison between the present method and that developed by Conte et al. [
34] using the simple infinite slope model is shown in
Figure 10a,b. In this last method, it is assumed that the slip surface is located at a depth of 20 m from the ground surface and is parallel to the slope. Specifically,
Figure 10a documents the results of the calibration stage, and
Figure 10b shows the horizontal displacements predicted by the two methods on the basis of rainfall recorded from August 2010 to May 2011. It is evident from these figures that, for the considered case study, the predictive capacity of the proposed method is better than that of the method based on an infinite slope model that more roughly approximates the real geometry of the landslide body.