In highway construction across the southeastern coastal regions of China, granite residual soil is widely used as subgrade fill material in pavement engineering. Its mechanical behaviour under dynamic loads warrants in-depth investigation. Dynamic events such as vehicular traffic and earthquakes are complex, involving multidirectional loads. The dynamic behaviour of soil under bidirectional cyclic loading differs significantly from that under cyclic loading in one direction. A large-scale bidirectional cyclic direct shear apparatus was utilised to carry on a series of horizontal cyclic direct shear tests on granite residual soil with water contents of 14% and 24% at different normal stress amplitudes (sigma a) (0, 100, 200 kPa). Based on these tests, discrete element method (DEM) models were developed to simulate the laboratory tests. The test results revealed that cyclic normal stress increases dynamic shear strength during forward shear but reduces it during reverse shear. The energy dissipation capacity increases with rising sigma a. The dynamic behaviour of granite residual soil is more significantly affected by cyclic normal stress when the water content is higher. The DEM simulation results indicated that as cyclic shearing progresses, the location of the maximum principal stress (sigma 1) shifts from the top of the specimen toward the shear interface. The distribution of the angle between sigma 1 and the x-axis, as well as sigma 1 and the z-axis, transitions from 'M' distribution to 'Arch' distribution. With increasing sigma a, during forward shear, the magnitude of the maximum principal stress increases, and the orientation of sigma 1 rotates toward the normal direction. Conversely, during reverse shear, the magnitude of the maximum principal stress decreases, and its orientation shifts toward the horizontal shear direction. The material fabric anisotropy coefficient decreases with increasing sigma a, while the anisotropy orientation increases with increasing normal stress.

In order to investigate the frost-heaving characteristics of wintering foundation pits in the seasonal frozen ground area, an outdoor in-situ test of wintering foundation pits was carried out to study the changing rules of horizontal frost heave forces, vertical frost heave forces, vertical displacement, and horizontal displacement of the tops of the supporting piles under the effect of groundwater and natural winterization. Based on the monitoring condition data of the in-situ test and the data, a coupled numerical model integrating hydrothermal and mechanical interactions of the foundation pit, considering the groundwater level and phase change, was established and verified by numerical simulation. The research results show that in the silty clay-sandy soil strata with water replenishment conditions and the all-silty clay strata without water replenishment conditions, the horizontal frost heave force presents a distribution feature of being larger in the middle and smaller on both sides in the early stage of overwintering. With the extension of freezing time, the horizontal frost heave force distribution of silty clay-sand strata gradually changes from the initial form to the Z shape, while the all-silty clay strata maintain the original distribution characteristics unchanged. Meanwhile, the peak point of the horizontal frost heave force in the all-silty clay stratum will gradually shift downward during the overwintering process. This phenomenon corresponds to the stage when the horizontal displacement of the pile top enters a stable and fluctuating phase. Based on the monitoring conditions of the in-situ test, a numerical model of the hydro-thermo-mechanical coupling in the overwintering foundation pit was established, considering the effects of the groundwater level and ice-water phase change. The accuracy and reliability of the model were verified by comparison with the monitoring data of the in-situ test using FLAC3D finite element analysis software. The evolution of the horizontal frost heaving force of the overwintering foundation pit and the change rule of its distribution pattern under different groundwater level conditions are revealed. This research can provide a reference for the prevention of frost heave damage and safety design of foundation pit engineering in seasonal frozen soil areas.

Tree destruction induced by heavy rainfall, an overlooked type of forest degradation, has been exacerbated along with global climate change. On the Chinese Loess Plateau, especially in afforested gully catchments dominated by Robinia pseudoacacia, destructive rainfall events have increasingly led to widespread forest damage. Previous study has manifested the severity of heavy rainfall-induced tree destruction and its association with topographic change, yet the contributions of tree structure and forest structure remain poorly understood. In this study, we quantified the destroyed trees induced by heavy rainfall using light detection and ranging (LiDAR) techniques. We assessed the influence of tree structure (tree height, crown diameter, and crown area), forest structure (tree density, gap fraction, leaf area index, and canopy cover), and terrain parameters (elevation, slope, and terrain relief) using machine learning models (random forest and logistic regression). Based on these, we aimed to clarify the respective and combined contributions of structural and topographic factors to rainfall-induced tree destruction. Key findings revealed that when considered in isolation, greater tree height, crown diameter, crown area, leaf area index (LAI), and canopy cover suppressed tree destruction, whereas higher gap fractions increased the probability of tree destruction. However, the synergistic increases of tree structural factors (tree height, crown diameter, and crown area) and forest structural factors (LAI and canopy cover) significantly promoted tree destruction, which can counteract the inhibitory effect of terrain on destruction. In addition, increases in tree structure or canopy density (LAI and canopy cover) also increased the probability of tree destruction at the same elevation. Our findings challenge conventional assumptions in forest management by demonstrating the interaction of tree structure and canopy density can significantly promote tree destruction during heavy rainfall. This highlights the need to avoid overly dense afforestation in vulnerable landscapes and supports more adaptive, climate-resilient restoration strategies.

In cold regions, the strength and deformation characteristics of frozen soil change over time, displaying different mechanical properties than those of conventional soils. This often results in issues such as ground settlement and deformation. To analyze the rheological characteristics of frozen soil in cold regions, this study conducted triaxial creep tests under various creep deviatoric stresses and established a corresponding Discrete Element Method (DEM) model to examine the micromechanical properties during the creep process of frozen clay. Additionally, the Burgers creep constitutive model was used to theoretically validate the creep deformation test curves. The research findings indicated that frozen clay primarily exhibited attenuated creep behavior. Under low confining pressure and relatively high creep deviatoric stress, non-attenuated creep was more likely to occur. The theoretical model demonstrated good fitting performance, indicating that the Burgers model could effectively describe and predict the creep deformation characteristics of frozen clay. Through discrete element numerical simulations, it was observed that with the increase in axial displacement, particle displacement mainly occurs at both ends of the specimen. Additionally, with the increase in creep deviatoric stress, the specimen exhibits different deformation characteristics, transitioning from volumetric contraction to expansion. At the same time, the vertical contact force chains gradually increase, the trend of particle sliding becomes more pronounced, and internal damage in the specimen progresses from the ends toward the middle.

Soybean urease-induced calcium carbonate precipitation (SICP) is an innovative and eco-friendly approach with demonstrated potential for mitigating soil liquefaction. However, the specific impacts of the concentrations of soybean urease and salt solutions require further elucidation. The research examines how the two compositions influence calcium carbonate formation. Dynamic characteristics of one-cycle SICP-treated clean and silty sand were analyzed based on cyclic triaxial tests. It was revealed that SICP-treated specimens of both liquefied sand and silty sand exhibit reduced accumulation of excess pore pressure and diminished strain growth under cyclic loading, thereby delaying liquefaction failure. Although higher concentrations of both soybean urease and salt solution can enhance liquefaction resistance, salt solution concentration has a more pronounced effect on improving liquefaction resistance due to the more production of calcium carbonate. Scanning electron microscopy observations confirmed the presence of calcium carbonate crystals at the interfaces between sand particles and between sand and fine particles. These crystals effectively bond the loose sand and fine particles into a cohesive matrix, reinforcing soil structure. A direct linear correlation was established between the liquefaction resistance improvement and precipitated calcium carbonate content. Notably, the one-cycle SICP treatment method adopted in this study demonstrates a better biocementation effect compared to cement mortar or multi-cycle MICP-treated sand under the same content of cementitious materials. These findings provide valuable insights for optimizing SICP treatments, aiming to reduce the risk of soil liquefaction in potential field applications.

In this paper, through extensive on-site research of the plain concrete composite foundation for the Jiuma Expressway, the study conducted proportional scaling tests. This study focused on the temperature, moisture, pile-soil stress, and deformation of this foundation under freeze-thaw conditions. The findings indicate that the temperature of the plain concrete pile composite foundation fluctuates sinusoidally with atmospheric temperature changes. As the depth increases, both temperature and lag time increase, while the fluctuation range decreases. Furthermore, the effect of atmospheric temperature on the shoulder and slope foot is more significant than on the interior of the road. During the freeze-thaw cycle, the water content and pore-water pressure in the foundation fluctuate periodically. The pile-soil stress fluctuates periodically with the freeze-thaw cycle, with the shoulder position exhibiting the most significant changes. Finally, the road displays pronounced freeze-thaw deformations at the side ditch and slope toe. This study provides a valuable basis for the construction of highway projects in cold regions.

In the loess tableland, gully slope instability induces severe soil erosion and land degradation, yet the synergistic effects of dominant vegetation under varying restoration modes combined with dynamic rainfall regimes and topographic variations on gully slope stabilization mechanisms remain inadequately quantified. Therefore, the dominant vegetation species under natural (NR) and artificial restoration (AR) was chosen as the object. Through field sampling, root-soil complex mechanical experiments, and numerical simulations, the protection effect of dominant vegetation under different restoration modes combination with rainfall and topographic variations was investigated. The result revealed significant differences in basic soil physical properties, root morphological characteristics, root and root-soil complex mechanical properties among five dominant vegetated plots under the different restoration modes (P < 0.05). The soil properties in the Scop plot under AR were slightly better than those in the other plots. The roots in the Spp plot developed better under NR. The shear strength of Lespedeza bicolor Turcz. was the highest under NR. The tensile strength of Digitaria sanguinalis (L.) Scop. was greatest under AR. The tensile force and tensile strength of single roots exhibited a significant positive linear correlation and a significant negative exponential correlation, with root diameter, respectively (P < 0.01). For the unstable gully slopes (F-s < 1.0), maximum displacement occurred at the slope foot, where tensile shear failure dominated, while the interior experienced compressive yielding. The grey relational analysis identified rainfall intensity as the primary destabilizing factor, followed by dominant vegetation species, slope height, and slope gradient. Notably, when rainfall intensity reaches or exceeds 0.06 m/h, or when slope height exceeds 20 m combined with long-duration rainfall, the regulatory impacts of dominant vegetation under different restoration modes on the gully slope stability are substantially diminished and become negligible. This study provides a theoretical basis for gully slope protection and ecological environmental construction in loess tableland.

Landslides, which are a type of process-based geological hazard, exhibit stagewise characteristics that serve as important guidance for the prevention and mitigation of slope engineering disasters. The cross-correlation and randomness of soil parameters can influence the evolution of landslide characteristics. This paper, based on the spatial variability of slope soil parameters, combines copula theory and the material point method (MPM) to establish a Monte Carlo-random material point method considering the cross-correlation of soil parameters. This resulting method is called copula-RMPM. It investigates the probability distributions of slope instability and landslide large deformation characteristics, such as sliding distance, landslide thickness, collapse range, and volume of sliding mass. The results indicated that in the study of soil parameter characteristics, failure probability increases with increased correlation coefficient. Also, failure probability showed a positive correlation with the variability coefficient of cohesion and internal friction angle, with failure probability being more sensitive to the variability coefficient of the internal friction angle. The landslide large deformation characteristics generally follow the normal distribution; they exhibit significant fluctuations in sliding distance and sliding mass area despite the relatively small variability coefficient. Compared with the results of random field simulation of soil parameters, the probability of landslide large deformation characteristics obtained by deterministic soil parameters is often lower. Therefore, the probability distribution of landslide large deformation characteristics obtained by the Monte Carlo-random material point method considering the cross-correlation of soil parameters is more meaningful for engineering guidance.

Soil-rock mixtures are composed of a complex heterogeneous medium, and its mechanical properties and mechanism of failure are intermediate between those of soil and rock, which are difficult to determine. To consider the influence of different particle groups on soil-rock mixture's shear strengths, based on the mesomotion properties of the particles of different particle groups when the soil-rock mixture is deformed, it is classified into two-phase composites, matrix and rock mass. In this paper, based on the representative volume element model of soil-rock mixtures and the Eshelby-Mori-Tanaka equivalent contained mean stress principle, a model of shear constitutive of the accumulation considering the mesoscopic characteristics of the rock is established, the influence of different factors on the shear strength of the accumulation is investigated, and the mesoscopic strengthening mechanism of the rock on the shear strength of the accumulation is discussed. The results show that there is a positive correlation between the rock content, the surface roughness of the rock, the stress concentration coefficient, coefficient of average shear displacement, and the accumulation's shear strength. When the accumulation is deformed, it stores or releases additional energy than the pure soil material, so it shows an increase in deformation resistance and shear strength on a macroscopic scale.

Volume changes in soil caused by freeze-thaw cycles can affect the shear performance of the saline soil-geotextile interface. To investigate this issue, the study examined changes in shear strength, deformation characteristics, and failure modes of the saline soil-geotextile interface under different numbers of freeze-thaw cycles. The experimental results indicate that with the increase in freeze-thaw cycles, the shear stiffness of the interface initially increases and then decreases, demonstrating the reduction in elasticity and resistance to deformation caused by freeze-thaw cycles. And the enhancement of normal stress can effectively increase the density of the soil and the adhesion at the interface, thereby improving shear stiffness. Meanwhile, the salt content in the soil also significantly impacts the mechanical properties, with notable changes in the dynamic characteristics of the interface as the salt content varies. Furthermore, after freeze-thaw actions, the soil becomes loose, reduces in integrity, features uneven surfaces, and sees increased internal porosity leading to slip surfaces. Trend analysis from this study provides new insights into the failure mechanisms at the saline soil-geotextile interface.