Energy piles, which serve the dual functions of load-bearing and geothermal energy exchange, are often modeled with surrounding soil assumed to be either fully saturated or completely dry in existing design and computational methods. These simplifications neglect soil saturation variability, leading to reduced predictive accuracy of the thermomechanical response of energy piles. This study proposes a novel theoretical framework for predicting the thermo-hydro-mechanical (THM) behavior of energy piles in partially saturated soils. The framework incorporates the effects of temperature and hydraulic conditions on the mechanical properties of partially saturated soils and pile-soil interface. A modified cyclic generalized nonlinear softening model and a cyclic hyperbolic model were developed to describe the interface shear stress-displacement relationship at the pile shaft and base, respectively. Governing equations for the load-settlement behavior of energy piles in partially saturated soils were derived using the load transfer method (LTM) and solved numerically using the matrix displacement method. The proposed approach was validated against experimental data from both field and centrifuge tests, demonstrating strong predictive performance. Specifically, the average relative error (ARE) was less than 15% for saturated soils and below 23% for unsaturated soils when evaporation effects were considered. Finally, parametric analyses were conducted to assess the effects of flow rate, groundwater table position, and softening parameters on the THM behavior of energy piles. This framework can offer a valuable tool for predicting THM behavior of energy piles in partially saturated soils, supporting their broader application as a sustainable foundation solution in geotechnical engineering.

The parameters of the soil water characteristic curve (SWCC) play a pivotal role in the examination of unsaturated soil behavior. This study employs three machine learning models-random forest (RF), extreme gradient boosting (XGBoost), and multiexpression programming (MEP)-to predict the SWCC using key soil properties. Among them, the RF model demonstrated the most robust performance in SWCC prediction. The Shapley Additive Explanation (SHAP) analysis further reveals that suction is the most influential factor affecting SWCC predictions, with other input parameters also contributing significantly. Additionally, the MEP model offers a straightforward expression for SWCC estimation and, thus, proved practical for predicting embankment responses and exhibited superior accuracy over traditional methods, such as the Arya and Paris model (ACAP). For a precise assessment of the hydromechanical response of the embankment subjected to infiltration, an increase in pore pressure is observed when employing the MEP model compared to the ACAP model for fine-grained soils. The findings emphasize the potential of RF and MEP in enhancing SWCC prediction and their practical implications for soil engineering applications.

The permanent displacement of earth slopes during earthquake shaking is a key indicator for landslide hazard assessment. Previous studies mostly attempt to evaluate the earthquake-induced displacement of dry or saturated soil slopes, while it is less common to deal with partially saturated soils. In the present study, a simplified procedure is proposed to account for the seismic-induced excess pore pressure in slopes with partially saturated sandy soils. The effect of matric suction, suction stress, and excess pore pressure on the yield acceleration of partially saturated sandy slopes is investigated, and the coupled Newmark sliding block method, known as the flexible soil columns with dynamic shear modulus and damping ratio, is modified to estimate the seismic slope displacement. Detailed discussions are made about the effect of different degrees of saturation on the excess pore pressure ratio, yield acceleration, and slope displacement. The numerical results show that the excess pore pressure ratio tends to exponentially increase with saturation, and the change of yield acceleration and displacement with saturation can be divided into suction stress dominant and excess pore water pressure dominant stages.

The stability of the slopes is critical for ensuring the safety and longevity of soil structures including, embankments, stockpiles, and retaining walls. Many of these soil structures remain in partially saturated conditions throughout their design life. In such cases, the matric suction (i4.'m) becomes a key parameter influencing the stability of the soil. However, i4.'m in soil varies during the water infiltration, leading to potential instability and collapse. To better understand this behaviour, constant shear drained (CSD) tests, which replicate the stress path experienced by the soil during water infiltration were conducted. A series of CSD experiments on silty sand in the partially saturated state was performed to examine the effects of i4.'m on the onset of instability. The results show that the onset of instability in silty sand increases with i4.'m. The volume change behaviour for partially saturated silty sand exhibits dilative behaviour compared to the fully saturated condition for the density considered in this study. Moreover, the methodology for the onset of instability for fully saturated soils was extended to analyse the onset of instability in partially saturated soil. A unique onset of instability was obtained for partially saturated soil using various methods considered for this study.

The occurrence of earthquake-induced soil liquefaction poses a significant threat, leading to extensive damage to building foundations and other structures, resulting in substantial economic repercussions. The seismic performance of geotechnical systems is markedly influenced by the saturation level of the soil. This study examines the impact of dynamic response on Palar sand. Cyclic triaxial tests were conducted on partially saturated finegrained loose sand with a relative density of 35 % and a degree of saturation ranging from 65 % to 75 %. These tests were carried out at a strain rate of 0.1 % and confining pressures of 50 and 75 kPa. The study findings reveal that an increase in back pressure corresponds to a rise in the excess pore water pressure ratio of the sand. Additionally, the sand undergoes liquefaction as the number of cycles increases, and the degree of saturation decreases for different confining pressures at frequencies of 0.75 and 1 Hz. It was observed that soil liquefies more rapidly at lower strain rates with an increase in effective confining pressure. Conversely, at higher frequencies, soil liquefaction occurs in a smaller number of cycles. Comparing the effects of confining pressure and frequency, a damping ratio of 13 % and a shear modulus of 40 MPa were achieved at a frequency of 0.75 Hz and a confining pressure of 50 kPa. The shear modulus of partially saturated sand decreases with an increase in the initial degree of saturation due to specific characteristics of the Palar sand and the loading conditions.

Earthquakes and groundwater are pivotal factors affecting slope stability. However, the majority of previous studies have focused on these factors individually, neglecting their combined effects. Hence, this paper aims to develop a framework using the kinematic approach of limit analysis to investigate the stability of slopes in partially saturated soils under the combined effects of seismic force and pore-water pressure. The pseudodynamic method (PDM) was employed to capture the temporal-spatial effect of horizontal and vertical seismic waves. Variations in suction and effective unit weight profiles with moisture content under steady-state unsaturated flow were considered. External rates arising from both static pore-water pressure and earthquake-induced excess pore-water pressure were incorporated into the energy-balance equation. With the aid of gravity increase method (GIM), an explicit expression of safety factor (FS) was derived and optimized using a genetic algorithm (GA). The validity of this approach was verified through a comparison with existing solutions. Parametric analyses were conducted to explore the influence of varying groundwater level, seismic coefficients, suction, threedimensional effects, excess pore water pressure, unsaturated flow types, and pseudo-dynamic parameters, on the FS and critical sliding surface of slopes in partially saturated slopes. This framework can provide a good reference for the safety design of reservoir slope under the combined effects of earthquakes and groundwater.

Earthquake-induced liquefaction is a relevant natural hazard due to the damages caused in numerous buildings, facilities and infrastructures worldwide. The damages caused to the infrastructure by this phenomenon are caused by the loss of stiffness and strength in granular soils, which leads to settlements and lateral spreading. Earthquake-induced liquefaction typically occurs in saturated deposits composed of non-plastic soils. Hence, the degree of saturation reduction is considered one of the most favourable and optimistic methods for liquefaction resistance mitigation. This paper explores the earthquake-induced liquefaction in saturated and gassy sands, varying their degree of saturation and state parameters. The state parameter was used to analyse the mechanical behaviour by combining the effects of relative density (or initial void ratio) with confinement pressure. Results show that liquefaction resistance improvement caused by the reduction in the degree of saturation is higher as the state parameter increases. This improvement can be described and quantified by multivariate models integrating the effects of degree of saturation and state parameter on liquefaction resistance. This provides a potential solution for improving the resilience of infrastructures susceptible to earthquake-induced liquefaction.

A fully coupled micro-hydromechanical (micro-HM) model is developed for partially saturated soils in this study by integrating two-dimensional pore morphology (PM) approach and discrete element method (DEM). In the proposed model, the PM approach is employed to predict the tentative water distribution. The porous media marching cubes (PMMC) algorithm is adopted to evaluate the interphase interfaces and to further calculate the capillary forces. The combined effects of interparticle contact forces and the capillary forces on the motion of particles are handled by DEM. The developed model was then employed to conduct a series of numerical biaxial shear tests on a partially saturated soil with real particle shapes. The typical macroscopic responses such as stress-strain relationship, volume change, and saturation change can be well simulated by the micro-HM model. Based on the micro-HM model, a novel equation is proposed to directly evaluate the effective stress from the pore water distribution. The effective stress parameter and the suction contribution to effective stress calculated by the new equation well matches the experimental data, thus confirming the validity of the micro-HM model and the new equation of effective stress. The microscopic responses are then revealed and discussed through the proposed model.

In this study, a generalized hypoplastic constitutive model for unsaturated soils is presented. The constitutive model is formulated in terms of effective stresses, degree of saturation, and suction. A feature of the model is the introduction of a limiting surface (LS) in hypoplasticity, allowing for the description of the maximum achievable void ratio as a function of mean effective stress and degree of saturation. The LS allows to capture the wetting-induced collapse of initially unsaturated soils. Contrary to other models, a concept of curved normal compression lines without limiting the range of applicable stresses is proposed. The performance of the proposed model is demonstrated by back-calculation of a well-documented experimental study on over 30 samples of compacted Pearl clay under isotropic as well as triaxial loading conditions over a wide range of stresses and void ratios. For this purpose, the proposed model is coupled with a hypoplastic model for the soil-water retention curve, which interrelates the effective degree of saturation with the suction and the void ratio.

Earthquake-induced soil liquefaction is a catastrophic phenomenon that can damage existing building foundations and other structures, resulting in significant economic losses. Traditional mitigation techniques against liquefaction present critical aspects, such as high construction costs, impact on surrounding infrastructure and effects on the surrounding environment. Therefore, research is ongoing in order to develop new approaches and technologies suitable to mitigate liquefaction risk. Among the innovative countermeasures against liquefaction, Induced Partial Saturation (IPS) is considered one of the most promising technologies. It consists of introducing gas/air bubbles into the pore water of sandy soils in order to increase the compressibility of the fluid phase and then enhance liquefaction resistance. IPS is economical, eco-friendly and suitable for urbanised areas, where the need to reduce the risk of liquefaction must be addressed, taking into account the integrity of existing buildings. However, IPS is still far from being a routine technology since more aspects should be better understood. The main aim of this review is to raise some important questions and encourage further research and discussions on this topic. The review first analyses and discusses the effects of air/gas bubbles on the cyclic behaviour of sandy soils, focusing on the soil volume element scale and then extending the considerations to the real scale. The use of useful design charts is also described. Moreover, a will be devoted to the effect of IPS under shallow foundations. The readers will fully understand the research trend of IPS liquefaction mitigation and will be encouraged to further explore new practical aspects to overcome the application difficulties and contribute to spreading the use of this technology.