In order to compensate for the shrinkage of geopolymer pastes uniformly, reactive MgO powders are evenly dispersed in the geopolymer. The deformation performance, mechanical properties, microstructure and components of geopolymer pastes with reactive MgO are characterized. The effects of the content and the activity of MgO are discussed. The results indicate that the chemical shrinkage, autogenous shrinkage and drying shrinkage decrease with the addition of reactive MgO. MgO reacted with water, and fine Mg(OH)₂ crystals forms as a geopolymer paste. Mg(OH)₂ produces uniform expansion, which refines the pore size of pastes and the compressive strength increases. The shrinkage of the geopolymer pastes is thus effectively compensated.

The change of spatial and temporal distribution of precipitation has an important impact on urban water security. The effect of land cover land use change (LCLUC) on the spatial and temporal distribution of precipitation needs to be further studied. In this study, transfer matrix, standard deviation ellipse and spatial autocorrelation analysis techniques were used. Based on the data of land cover land use and precipitation, this paper analyzed the land cover land use change and its influence on the spatial and temporal distribution pattern of precipitation in the Guangdong–Hong Kong–Macao Greater Bay Area (GBA). The results showed that from 2001 to 2019, the area of cropland, water, barren, forest/grassland in the GBA decreased by 44.03%, 8.05%, 50.22%, 0.43%, respectively, and the area of construction land increased by 20.05%. The precipitation in the GBA was mainly concentrated in spring and summer, and the precipitation in spring tended to increase gradually, while the precipitation in summer tended to decrease gradually, while the precipitation in autumn and winter has no obvious change. It was found that with the change of land cover land use, the spatial distribution of precipitation also changed. Especially in the areas where the change of construction land was concentrated, the spatial distribution of precipitation changed most obviously.

Planting vegetation is an environmentally friendly method for reducing landslides. Current vegetated slope analysis fails to consider the influence of different root architectures, and the accuracy and effectiveness of the numerical simulations need to be improved. In this study, an explicit smoothed particle finite element method (eSPFEM) was used to evaluate slope stability under the influence of vegetation roots. The Mohr–Coulomb constitutive model was extended by incorporating apparent root cohesion into the shear strength of the soil. The slope factors of safety (FOS) of four root architectures (uniform, triangular, parabolic, and exponential) for various planting distances, root depths, slope angles, and planting locations were calculated using the shear strength reduction technique with a kinetic energy-based criterion. The results indicated that the higher the planting density, the stronger the reinforcement effect of the roots on the slope. With increasing root depth, the FOS value first decreased and then increased. The FOS value decreased with an increase in slope angle. Planting on the entire ground surface had the best improvement effect on the slope stability, followed by planting vegetation with a uniform root architecture in the upper slope region or planting vegetation with triangular or exponential root architecture on the slope’s toe. Our findings are expected to deepen our understanding of the contributions of different root architectures to vegetated slope protection and guide the selection of vegetation species and planting locations.

A landslide is a common natural disaster that causes environmental damage, casualties and economic losses, which seriously affects the sustainable development of society. In geomechanics, it is one of the largest deformation problems. Herein, the GPU-accelerated explicit smoothed particle finite element method (eSPFEM) for large deformation analysis in geomechanics was developed on the CUDA platform based on high-performance computing using a self-designed eSPFEM program code. The eSPFEM combines the strain smoothing nodal integration techniques found in the particle finite element method (PFEM) framework, which allows for the use of low-order triangular elements without volume locking and avoids frequent information transfer and mapping errors between Gaussian points and particles in PFEM. A numerical simulation of slope instability using the eSPFEM and based on a strength reduction technique was conducted using various examples, including a cohesive homogeneous slope, a non-cohesive homogeneous slope, a non-homogeneous slope and a slope with a thin soft band. The calculation results show that the eSPFEM can be applied to slope stability analysis under different working conditions, simulating the entire process of slope instability initiation, sliding and reaccumulation, and obtaining reliable FOS values. A numerical simulation was conducted to analyse a landslide that occurred in the Zhangjiazhuang tunnel on the Lanzhou–Xinjiang high-speed railway line on 18 January 2016. A natural unsaturated soil slope, a soil slope with a high moisture content and a soil slope with a high moisture content subjected to an earthquake were analysed. The findings of this study are in good agreement with the actual slope failure conditions. The primary triggers identified for the landslide were heavy rainfall and earthquakes. The verification results indicate that the eSPFEM can effectively simulate an actual landslide case, showcasing high accuracy and applicability in simulating the large deformation behaviour of landslides.

Purpose: Understanding the fluid flow through rock masses, which commonly consist of rock matrix and fractures, is a fundamental issue in many application areas of rock engineering. As the equivalent porous medium approach is the dominant approach for engineering applications, it is of great significance to estimate the equivalent permeability tensor of rock masses. This study aims to develop a novel numerical approach to estimate the equivalent permeability tensor for fractured porous rock masses. Design/methodology/approach: The radial point interpolation method (RPIM) and finite element method (FEM) are coupled to simulate the seepage flow in fractured porous rock masses. The rock matrix is modeled by the RPIM, and the fractures are modeled explicitly by the FEM. A procedure for numerical experiments is then designed to determinate the equivalent permeability tensor directly on the basis of Darcy’s law. Findings: The coupled RPIM-FEM method is a reliable numerical method to analyze the seepage flow in fractured porous rock masses, which can consider simultaneously the influences of fractures and rock matrix. As the meshes of rock matrix and fracture network are generated separately without considering the topology relationship between them, the mesh generation process can be greatly facilitated. Using the proposed procedure for numerical experiments, which is designed directly on the basis of Darcy’s law, the representative elementary volume and equivalent permeability tensor of fractured porous rock masses can be identified conveniently. Originality/value: A novel numerical approach to estimate the equivalent permeability tensor for fractured porous rock masses is proposed. In the approach, the RPIM and FEM are coupled to simulate the seepage flow in fractured porous rock masses, and then a numerical experiment procedure directly based on Darcy’s law is introduced to estimate the equivalent permeability tensor.

During urbanisation, extensive production and construction activities encroach on ecological spaces, leading to changes in environmental structures and soil erosion. The issue of yellow muddy water caused by rainfall in cities with high construction intensity has garnered significant attention. Taking Guangzhou City as the research area, this study is the first to propose a risk assessment model for yellow muddy water in cities with high construction intensity, and the influence of construction sites on yellow muddy water was fully considered. Rainfall and construction sites were used as indicators to assess the hazards of yellow muddy water. Elevation, slope, normalised difference vegetation index (NDVI), soil erosion modulus, stream power index (SPI), surface permeability, and roads represent the exposure evaluation indicators. Population number and GDP (Gross Domestic Product) were used as vulnerability evaluation indicators. Based on the analytic hierarchy process (AHP) method, the weights of each evaluation indicator were determined, and a risk assessment system for yellow muddy water was established. By overlaying the weighted layers of different evaluation indicators on the geographic information system (GIS) platform, a risk degree distribution map of yellow muddy water disasters was generated. The evaluation results demonstrated that the disaster risk levels within the study area exhibited spatial differentiation, with areas of higher risk accounting for 14.76% of the total. The evaluation results were compared with historical yellow muddy water event information from Guangzhou, and the effectiveness of the model was verified by the receiver operating characteristic (ROC) curve. The validation results indicate that this model provides high accuracy in assessing the degree of risk of yellow muddy water in high-construction-intensity cities, offering effective technical support for precise disaster prevention and mitigation.