In this study, an analytical model for the three-dimensional (3D) dynamic stability analysis of vegetation-rooted slopes is first developed under steady-state unsaturated flow conditions. Root reinforcement, defined as the increase in the soil shear strength produced by the mechanical and hydrological effects of vegetation roots, is included in the proposed analytical model. By combining the modified pseudo-dynamic approach (MPDA) and the kinematic theory of limit analysis to the 3D discretized failure model, the most critical failure surface and the corresponding factor of safety (FS) are derived to examine the stability of vegetation-rooted slopes with the aid of the optimization algorithm of particle swarm. The proposed approach is verified by comparing with published analytical solutions and numerical results. A series of parametric analysis are then conducted to examine the influence of seismic-related parameters, vegetation properties, possible surcharge and slope geometry parameters on the slope stability. Finally, a comparison between the slope stability under different root architectures is provided and discussed. The results show that, for these selected cases, the stability of vegetation-rooted slopes is significantly improved by approximately 45% compared to bare soil slopes, and the divergences of reinforcement effects between different root architectures can be negligible.

Important unsaturated soil mechanics topics for all geotechnical engineers and geotechnical engineering students are reviewed. These key topics include: (1) Soil is an elastoplastic material for which the macro-level response, in general, is controlled by two separate stress variables: total stress (net stress) and negative pore water pressure (suction). (2) Pore water pressures are always negative above the groundwater table-and should not be conservatively assumed zero; (3) shear strength and volume change of unsaturated soils are dependent on soil suction, as well as confining stress, and therefore geotechnical site investigations and testing must account for both stress variables; (4) water flow follows Darcy's law, but hydraulic conductivity is a strong function of water content such that fine-grained soil can have a higher conductivity than course-grained soil, leading to unexpected results when using saturated flow thinking processes; (5) unsaturated soil response is complex and difficult to intuit in the absence of laboratory testing and simulation. Features of unsaturated soil behavior most frequently encountered in geotechnical practice are highlighted, with discussion and demonstration from existing literature. Suggestions are given for relatively simple approaches for first steps in taking unsaturated soil mechanics principles into consideration in site investigation, laboratory testing, and design-related decisions.

Alpine permafrost environments are highly vulnerable and sensitive to changes in regional and global climate trends. Thawing and degradation of permafrost has numerous adverse environmental, economic, and societal impacts. Mathematical modeling and numerical simulations provide powerful tools for predicting the degree of degradation and evolution of subsurface permafrost as a result of global warming. A particularly significant characteristic of alpine environments is the high variability in their surface geometry which drives large lateral thermal and fluid fluxes along topographic gradients. The combination of these topography-driven fluxes and unsaturated ground makes alpine systems markedly different from Arctic permafrost environments and general geotechnical ground freezing applications, and therefore, alpine permafrost demands its own specialized modeling approaches. In this work, we present a multi-physics permafrost model tailored to subsurface processes of alpine regions. In particular, we resolve the ice-water phase transitions, unsaturated conditions, and capillary actions, and account for the impact of the evolving pore space through freezing and thawing processes. Moreover, the approach is multi-dimensional, and therefore, inherently resolves the topography-driven horizontal fluxes. Through numerical case studies based on the elevation profiles of the Zugspitze (DE) and the Matterhorn (CH), we show the strong influence of lateral fluxes in 2D on active layer dynamics and the distribution of permafrost.