the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Technical note: The beneficial role of a natural permeable layer in slope stabilization by drainage trenches
Gianfranco Urciuoli
Marianna Pirone
Luciano Picarelli
Slope stabilization through drainage trenches is a classic approach in geotechnical engineering. Considering the low hydraulic conductivity of the soils in which this measure is usually adopted, a major constraint to the use of trenches is the time required to obtain a significant pore pressure decrease, here called “time lag”. In fact, especially when the slope safety factor is small, the use of drainage trenches may be a risky approach due to the probability that slope deformations will damage the system well before it will become fully operative.
However, this paper shows that the presence of persistent permeable natural soil layers can provide a significant benefit by increasing drainage efficiency and reducing time lag. As a matter of fact, any permeable layer that is intercepted by trenches may operate as part of the global hydraulic system, reducing the drainage paths.
A simplified approach to designing a drainage system that accounts for the presence of a persistent permeable layer is proposed. This approach, which can exploit solutions available in the literature for parallel drainage trenches, has been validated by numerical analyses.
The stabilization of deep landslides in clay is one of the greatest challenges to engineers due to the high cost and the unreliability of many structural solutions. Often, the only available approach is deep drainage, which can lead to some shear strength increase through a generalized pore pressure decrease. Available solutions (Hutchinson, 1977; Bromhead, 1984; Stanic, 1984; Desideri et al., 1997; Pun and Urciuoli, 2008; Urciuoli and Pirone, 2013) include the use of deep parallel trenches (and of deeper drainage panels as well), which are also dealt with in this paper, and tubular drains in a homogeneous soil.
Considering the finegrained nature of the soil, a major constraint to slope stabilization by draining trenches is the long time required to obtain a significant pore pressure decrease (time lag). In particular, when the slope is characterized by a small safety factor or is subjected to slow movements (Urciuoli, 1998), the use of draining trenches is in fact problematic due to the probability that slope deformations will damage the system well before it will become fully operative, thus cancelling out its potential effectiveness. However, the higher the depth of the trenches (or of drainage panels), the higher the probability that these intercept even thin soil layers of higher hydraulic conductivity at an intermediate depth between the ground surface and the slip surface. This would be a favourable condition since the incorporation of such layers in the drainage system may play a highly beneficial role in both the time to attain the final steadystate condition and the system efficiency.
The scope of this paper is restricted to an examination of the influence on the drainage system of a preexisting permeable soil layer parallel to the ground surface.
The solutions presented below are based on the following assumptions:

The groundwater flow is twodimensional.

Each soil layer is homogeneous, isotropic and characterized by a linear elastic constitutive law.

Total stresses are constant during the consolidation process (this allows the analysis of the hydraulic and of the mechanical soil responses to be uncoupled).
The governing equation of the problem (i.e. soil consolidation induced by the draining elements) is the following:
where $h=\mathit{\zeta}+\frac{u}{{\mathit{\gamma}}_{\mathrm{w}}}$ and ${c}_{v}^{\mathrm{2}\mathrm{D}}=\frac{KE}{\mathrm{2}\left(\mathrm{1}+\mathit{\nu}\right)\left(\mathrm{1}\mathrm{2}\mathit{\nu}\right){\mathit{\gamma}}_{\mathrm{w}}}$.
The technical literature reports solutions concerning the case of parallel draining trenches and of tubular drains in homogeneous soils, which are generally presented in the form of dimensionless design charts, providing the average efficiency, $\stackrel{\mathrm{\u203e}}{E}\left(t,\mathrm{\Gamma}\right)$, along the slip surface Γ:
In Eq. (2), u(0,Γ) is the initial pore pressure on the slip surface, Γ, and $\stackrel{\mathrm{\u203e}}{u}\left(t,\mathrm{\Gamma}\right)$ is the average pore pressure at time t modified by the draining elements; u(0,Γ)is generally assumed to be hydrostatic. During the consolidation phase, pore pressures decrease towards the minimum steadystate value $\stackrel{\mathrm{\u203e}}{u}(\mathrm{\infty},\mathrm{\Gamma})$, which is attained at time t→∞ when the efficiency $\stackrel{\mathrm{\u203e}}{E}\left(\mathrm{\infty},\mathrm{\Gamma}\right)$ reaches the highest value.
The available solutions for parallel trenches, featured by a thickness H_{0} and a width b, consider the soil volume between the two axes of symmetry, which respectively coincide with the middle of a trench and the centre line between two adjacent trenches (Fig. 1a). This volume is delimited by the ground surface and by an impermeable bottom surface located at the distance H from the ground surface. The ground and bottom surfaces are both horizontal: the slope angle is indeed assumed to play a negligible role in the hydraulic process (Aloi et al., 2019). The slip surface Γ is a horizontal plane as well, located at depth D. In this paper it is assumed to be coincident with the base of the trenches (D=H_{0}).
A key hypothesis, which strongly affects the solution, is the presence of a permanent film of water at the ground surface (Burghignoli and Desideri, 1987; D'Acunto and Urciuoli, 2006, 2010; D'Acunto et al., 2007). However, due to local formation of water ponding and saturation of vertical cracks in the ground, often this is not far from the truth, at least during the wet season. Based on this assumption, the pore pressure decrease is uniquely due to rotation of the flow lines towards the drainage trenches. Pore pressures in the zone between parallel trenches are then at any time less than hydrostatic (Fig. 1b). In contrast, beyond the bottom of the trenches, the upward direction of the flow lines leads to a pore pressure distribution higher than hydrostatic. It is just for this reason that the drains should always reach a depth close to the slip surface.
3.1 Time of consolidation
As outlined above, the presence of one or more persistent permeable layers in the soil body to be stabilized (a not unlikely situation in deep clay deposits to be stabilized with draining panels) may play a highly beneficial role in time lag and the effectiveness of the drainage system.
The influence of a layer parallel to the ground surface, here indicated as the “draining layer d”, featured by a thickness H_{d} as in Fig. 1c, has been investigated by FEM analyses using the code SEEP^{®} (GEOSLOPE Int. Ltd., 2012). The cases examined in this paper are indicated in Table 1; the results of the analyses are shown in a dimensionless form.
The data show that the presence of the permeable layer allows a considerable shortening of the time lag, here represented by the time factor, T_{90}, corresponding to an efficiency $\stackrel{\mathrm{\u203e}}{E}\left({t}_{\mathrm{90}},\mathrm{\Gamma}\right)=\mathrm{90}$ % of the average steadystate efficiency $\stackrel{\mathrm{\u203e}}{E}\left(\mathrm{\infty},\mathrm{\Gamma}\right)$:
The value of t_{90} in Eq. (3) has been obtained by a numerical integration of Eq. (1), being the time at which $\stackrel{\mathrm{\u203e}}{E}(t,\mathrm{\Gamma})=\mathrm{0.90}\cdot \stackrel{\mathrm{\u203e}}{E}(\mathrm{\infty},\mathrm{\Gamma})$ (see Eq. 2).
Figure 2a and b, which report some results concerning the horizontal plane located at depth D=H_{0}, suggest quite a rapid attainment of $\stackrel{\mathrm{\u203e}}{E}({t}_{\mathrm{90}},\mathrm{\Gamma})$, which is a crucial issue of the design. For significant values of trench spacing in the practice (i.e. $s/{H}_{\mathrm{0}}<\mathrm{3}$), the following considerations may be drawn: (i) for ${H}_{\mathrm{1}}/{H}_{\mathrm{0}}=\mathrm{0.75}$ and ${K}_{d}/K=\mathrm{100}$ (Fig. 2a), the dimensionless time T_{90} ranges between onehalf and onethird of the value that would be obtained in the absence of the draining layer; (ii) for ${K}_{d}/K=\mathrm{1000}$ (Fig. 2b), T_{90} significantly decreases with depth of the layer d (for ${H}_{\mathrm{1}}/{H}_{\mathrm{0}}=\mathrm{0.75}$, it drops to about 20 % of the value obtainable in homogenous soils).
3.2 Steadystate condition
The presence of a permeable layer allows higher values of $\stackrel{\mathrm{\u203e}}{E}\left(\mathrm{\infty},\mathrm{\Gamma}\right)$ to be obtained, and sooner than in homogeneous soils. Some significant data are provided:
 i.
Figure 2c, shows the steadystate efficiency for ${H}_{\mathrm{1}}/{H}_{\mathrm{0}}=\mathrm{0.75}$ reported as a function of the ratio K_{d}∕K and of trench spacing. As shown, the higher the hydraulic conductivity of the draining layer, the higher the efficiency (as an example, for ${K}_{d}/K=\mathrm{1000}$ and $s/{H}_{\mathrm{0}}=\mathrm{3}$ it practically doubles). A major effect of layer d is in fact diversion of a significant part of water coming from the ground surface towards the trench, thus strongly reducing water flow towards the slip surface.
 ii.
Figure 2d shows the efficiency for ${K}_{d}/K=\mathrm{1000}$ reported as a function of depth of layer d and trench spacing. The figure shows that it increases as the dimensionless distance, H_{1}∕H_{0}, increases. The effect of layer d is a strong pore pressure reduction at depth H_{1}. As a consequence, pore pressure decrease, due to the action of layer d, increases with its depth.
 iii.
Figure 2e shows the nondimensional pore pressure distribution, u_{d}, along the lower boundary of the draining layer d plotted as a function of trench spacing for ${H}_{\mathrm{1}}/{H}_{\mathrm{0}}=\mathrm{0.75}$ and ${K}_{d}/K=\mathrm{1000}$. Near the trench boundary, the pressure head is less than H_{d}; hence, a free water surface forms in the layer d (here water can move towards the trench only below this surface, where pore pressures are positive).
A simplified approach to predict the steadystate condition
In the following, a simplified model for the optimization of the design is briefly described. A very efficient working condition is achieved if, at the centre line between two adjacent trenches, the atmospheric pressure is attained at the uppermost point of layer d.
The first step in the design of the drainage system is just creating the conditions for full layer activation. This is obtained when the spacing, s, of the trenches is equal to the value s_{d}, according to the following expression:
The values of s_{d}∕H_{0} in Eq. (4) have been obtained from the results of the numerical integration of Eq. (1). These have been reported in Fig. 2f, which shows the dependency of s_{d}∕H_{0} on K_{d}∕K and H_{1}∕H_{0}, having fixed H_{d}∕H_{0}.
In case of full activation of layer d, the response of the entire draining system may be analysed by a simplified approach. Since the fluid pressure at the uppermost boundary of the layer d is equal to the atmospheric pressure (or to a small amount of suction, especially near the trench boundary), a water film may be fictitiously assumed at the same depth (Fig. 1d). This obviously leads to a generalized pore pressure decrease in the lowermost soil. In the following, any parameter referring to this fictitious condition will be indicated by the apex ^{∗}.
The values of ${\stackrel{\mathrm{\u203e}}{u}}^{\ast}\left(\mathrm{\infty},\mathrm{\Gamma}\right)$ and ${\stackrel{\mathrm{\u203e}}{E}}^{\ast}\left(\mathrm{\infty},\mathrm{\Gamma}\right)$ may be obtained from the well known dimensionless solutions for the case of parallel trenches in homogeneous soil, as a function of spacing (see the simplified scheme in Fig. 1d). The steadystate efficiency is
thus
It is worth noting that
and
The values obtained from Eq. (8), through the value of ${\stackrel{\mathrm{\u203e}}{u}}^{\ast}\left(\mathrm{\infty},\mathrm{\Gamma}\right)$ provided by mentioned solutions, have been compared to those obtained by FEM (Table 2). The good agreement allows the proposed method to be validated. It is worth mentioning that the solid lines in Fig. 2c and d for H_{d}=0 are just those that are reported in the design charts.
The aim of this paper is to demonstrate that the presence of soil layers of higher permeability, a not unlikely condition in some deep landslides in clay, may be exploited to improve the efficiency of systems of drainage trenches for slope stabilization. Once the depth of the trenches, which should reach the slip surface, is established, the selection of proper spacing may create a hydraulic system in which such layers can work as additional drains. The problem has been examined for the case that a unique permeable layer is present at an elevation higher than the bottom of the trenches.
The results of numerical analyses show that it significantly speeds up the consolidation process triggered by drainages, also leading to higher steady efficiency of the system. However, as mentioned in the Introduction, in many practical cases the design must take into account the time requested to achieve an adequate effective stress and safety factor increase. In these cases, trench spacing should be governed by the T_{90} value.
If pore pressures in the draining layer do not exceed the atmospheric pressure, a hydraulic disconnection forms between the two parts of the landslide body located above and below the layer, respectively. In such a way, the water film which is normally assumed at the ground surface ideally moves to the depth of the draining layer. This simple consideration allows the employment of the design charts available for the design of drainage trenches in homogeneous soils in the equivalent scheme characterized by groundwater level located at the depth of the draining layer, in order to calculate the final system efficiency.
It is worth mentioning that the hydraulic continuity of layer d is a fundamental condition for the design. Considering the variability and the unpredictability of many natural situations, proper investigations to check the validity of such an assumption are warmly recommended. In particular, the adoption of such a stabilization measure should always be managed through the “observational method”, i.e. by monitoring the system response in order to (i) check the validity of the design and (ii) adopt proper modifications to it due to unexpected or neglected factors. The installation of piezometers is an obvious measure to check in real time the efficiency of the drainage system (especially during the critical rainy season). The piezometers should be installed both in proximity to the slip surface (near and far from the trenches) and, if possible depending on thickness, in the permeable layer. This will allow the full activation of the permeable layer to be verified.
b  width of the trench 
${c}_{v}^{\mathrm{2}\mathrm{D}}$  2D coefficient of consolidation 
d  draining layer 
D  depth of the slip surface 
E  Young modulus of the soil 
$\stackrel{\mathrm{\u203e}}{E}(t,\mathrm{\Gamma})$  average efficiency of the draining trenches at time t along the sliding surface Γ 
$\stackrel{\mathrm{\u203e}}{E}(\mathrm{\infty},\mathrm{\Gamma})$  average steadystate efficiency of the draining trenches along the sliding surface Γ 
${\stackrel{\mathrm{\u203e}}{E}}^{\ast}\left(\mathrm{\infty},\mathrm{\Gamma}\right)$  average steadystate efficiency of the draining trenches along the sliding surface Γ according to the simplified approach (full activation of layer d) 
γ_{w}  unit weight of water 
Γ  slip surface 
h  total head 
h_{t}  first derivative of total head h with respect to time t 
h_{xx}  second derivative of total head h with respect to abscissa x 
h_{yy}  second derivative of total head h with respect to ordinate y 
H  depth of the impermeable bottom surface 
H_{0}  depth of the base of trench 
H_{1}  depth of the draining layer d 
H_{d}  thickness of the draining layer d 
K  coefficient of hydraulic conductivity 
K_{d}  coefficient of hydraulic conductivity of the draining layer d 
θ  soil moisture 
s  spacing between trenches 
s_{d}  spacing between trenches creating the conditions for full activation of the draining layer d 
s^{′}  distance between the boundaries of the trenches 
ν  Poisson ratio of the soil 
t  time 
t_{90}  dimensional time corresponding to $\stackrel{\mathrm{\u203e}}{E}\left(t,\mathrm{\Gamma}\right)=\mathrm{90}$ % of the average steadystate efficiency $\stackrel{\mathrm{\u203e}}{E}\left(\mathrm{\infty},\mathrm{\Gamma}\right)$ 
T_{90}  dimensionless time factor for $\stackrel{\mathrm{\u203e}}{E}\left(t,\mathrm{\Gamma}\right)=\mathrm{90}$ % of the average steadystate efficiency $\stackrel{\mathrm{\u203e}}{E}\left(\mathrm{\infty},\mathrm{\Gamma}\right)$ 
u  pore pressure 
u_{d}  pore pressure at the base of the draining layer d 
u(0,Γ)  initial pore pressure (time t=0) on the slip surface Γ 
${u}^{\ast}\left(\mathrm{0},\mathrm{\Gamma}\right)$  initial pore pressure (time t=0) on the slip surface Γ according to the simplified approach (full activation of the draining layer d) 
$\stackrel{\mathrm{\u203e}}{u}\left(t,\mathrm{\Gamma}\right)$  average pore pressure at time t on the slip surface Γ, modified by drainage trenches 
$\stackrel{\mathrm{\u203e}}{u}(\mathrm{\infty},\mathrm{\Gamma})$  average steadystate pore pressure on the slip surface Γ, modified by drainage trenches 
${\stackrel{\mathrm{\u203e}}{u}}^{\ast}\left(\mathrm{\infty},\mathrm{\Gamma}\right)$  average steadystate pore pressure on the slip surface Γ, modified by drainage trenches according to the simplified approach (full activation of the draining layer d) 
ζ  elevation head 
The datasets are freely downloadable from https://doi.org/10.5281/zenodo.3518143 (Urciuoli et al., 2019).
GU and LP jointly conceived and set up the research, discussing the issues with the other two authors, and wrote the general parts of the paper in Sects. 1 and 2. LC and MP carried out numerical analyses and interpreted the results, discussing them with the other two authors, and wrote the parts of the paper devoted to commenting on the results, in Sect. 3. The contributions of the authors are equal.
The authors declare that they have no conflict of interest.
This research has been developed within the framework of the PRIN 2015 project titled “Innovative Monitoring and Design Strategies for Sustainable Landslide Risk Mitigation”, supported by the Italian Ministry of Education, University and Research (MIUR) (grant no. 201572YTLA).
This paper was edited by Nadia Ursino and reviewed by two anonymous referees.
Aloi, F., Pirone, M., and Urciuoli, G.: Numerical investigation of small and mediumdiameter drain wells to stabilise deep landslides, Acta Geotech., 14, 1065–1080, https://doi.org/10.1007/s1144001806888, 2019.
Bromhead, E. N.: An analytical solution to the problem of seepage into countfort drains, Can. Geotech. J., 21, 657–662, 1984.
Burghignoli, A. and Desideri, A.: On the effectiveness of tubular drains, in: Proceedings of the IX ECSMFE, 31 August–3 September 1987, Dublin, Ireland, vol. 1, 121–124, 1987.
D'Acunto, B. and Urciuoli, G.: Groundwater regime in a slope stabilised by drain trenches, Math. Comput. Model., 43, 754–765, 2006.
D'Acunto, B. and Urciuoli, G.: Heavy rainfalls long term effects and mitigation with drains, in: Advances in environmental research, edited by: Aider, M. and Muller, T., Nova Science Publishers, New York, USA, vol. 2., 323–347, 2010.
D'Acunto, B., Parente, F., and Urciuoli, G.: Numerical models for 2D free boundary analysis of groundwater in slopes stabilized by drain trenches, Comput. Math. Appl., 53, 1615–1626, 2007.
Desideri, A., Miliziano, S., and Rampello, S.: Drenaggi a Gravità per la Stabilizzazione dei Pendii, Hevelius Edizioni, Benevento, Italy, 1997.
GEOSLOPE International Ltd.: Seepage modelling with SEEP/W, Calgary, Alberta, Canada, 2012.
Hutchinson, J. N.: Assessment of the effectiveness of corrective measures in relation to geological conditions and types of slope movement, Bulletin of the International Association of Engineering Geology, 16, 131–155, 1977.
Pun, W. K. and Urciuoli, G.: Soil nailing and subsurface drainage for slope stabilization, in: Proceedings of 10th International Symposium on Landslides and Engineered Slopes, 30 June–4 July 2008, Xi'an, China, Keynote paper, 85–126, 2008.
Stanic, B.: Influence of drainage trenches on slope stability, J. Geotech. Eng.ASCE, 110, 1624–1635, 1984.
Urciuoli, G.: Pore pressures in unstable slopes constituted by fissured clay shales, Proceedings of the Second International Symposium on Hard SoilsSoft Rocks, 12–14 October 1998, Naples, Italy, 1177–1185, 1998.
Urciuoli, G. and Pirone, M.: Subsurface drainage for slope stabilization, in: Landslide Science and Practice, edited by: Margottini, C., Canuti, P., and Sassa, K., Springer, Berlin, Heidelberg, Germany, vol. 6, 577–585, 2013.
Urciuoli, G., Comegna, L., Pirone, M., and Picarelli, L.: The beneficial role of stratigraphy on slope stabilization by drainage trenches: results of FEM numerical analyses, Zenodo, https://doi.org/10.5281/zenodo.3518143, 2019.
 Abstract
 Introduction
 The basic model
 Influence of a permeable layer located at an intermediate depth between ground and slip surface
 Conclusions and final considerations
 Appendix A: List of symbols
 Data availability
 Author contributions
 Competing interests
 Financial support
 Review statement
 References
 Abstract
 Introduction
 The basic model
 Influence of a permeable layer located at an intermediate depth between ground and slip surface
 Conclusions and final considerations
 Appendix A: List of symbols
 Data availability
 Author contributions
 Competing interests
 Financial support
 Review statement
 References