Stochastic modelling of relative water permeability in vegetative soils with implications on stability of bioengineered slope

Vegetation is known to influence the hydrological state variables, suction (ψ) and volumetric water content (θ_w) of soil. In addition, vegetation induces heterogeneity in the soil porous structure and consequently the relative permeability (k_r) of water under unsaturated conditions. The indirect method of utilising the soil water characteristic curve (SWCC) is commonly adopted for the determination of k_r. In such cases, it is essential to address the stochastic behaviour of SWCC, in order to conduct a robust analysis on the k_r of vegetative cover. The main aim of this study is to address the uncertainties associated with k_r, using probabilistic analysis, for vegetative covers (i.e., grass and tree species) with bare cover as control treatment. We propose two approaches to accomplish the aforesaid objective. The univariate suction approach predicts the probability distribution functions of k _r, on the basis of identified best probability distribution of suction. The bivariate suction and water content approach deals with the bivariate modelling of the water content and suction (SWCC), in order to capture the randomness in the permeability curves, due to presence of vegetation. For this purpose, the dependence structure of ψ and θ_w is established via copula theory, and the k_r curves are predicted with respect to varying levels of ψ- θ_w correlation. The results showed that the k_r of vegetative covers is substantially lower than that in bare covers. The reduction in k_r with drying is more in tree cover than grassed cover, since tree roots induce higher levels of suction. Moreover, the air entry value of the soil depends on the magnitude of ψ- θ_w correlation, which in turn, is influenced by the type of vegetation in the soil. k_r is found to be highly uncertain in the desaturation zone of the relative permeability curve. The stochastic behaviour of k_r is found to be most significant in tree covers. Finally, a simplified case study is also presented in order to demonstrate the impact of the uncertainty in k_r, on the stability of vegetates slopes. With an increment in the parameter α, factor of safety (FS) is found to decrease. The trend of FS is reverse of this with parameter n. Overall FS is found to vary around 4–5%, for both bare and vegetative slopes.