The reasonable value of good gradation characteristic parameters is key in designing and optimising soil-rock mixed high fill embankment materials. Firstly, the DJSZ-150 dynamic-static large-scale triaxial testing instrument was used for triaxial compression shear tests on compacted skeleton structure soil-rock mixture standard specimens. The changes in strength and deformation indicators under different gradation parameters and confining pressure were analysed. Then, based on the Janbu empirical formula, relationships between parameters K, n, and (sigma 1-sigma 3)ult and the coefficient of uniformity Cu and coefficient of curvature Cc were explored. Empirical fitting formulas for Duncan-Chang model constants a and b were proposed, establishing an improved Duncan-Chang model for soil-rock mixtures considering gradation characteristics and stress states. Finally, based on significant differences in particle spatial distribution caused by gradation changes, three generalised models of matrix-block stone motion from different particle aggregation forms were proposed. Results indicate the standard specimen's strength and deformation indicators exhibit significant gradation effects and stress-state correlations. The improved Duncan-Chang model effectively simulates the stress-strain relationship curve under different gradations and confining pressure, with its characteristics explainable based on the matrix block stone motion generalised model.

Flash floods are often responsible for deaths and damage to infrastructure. The objective of this work is to create a data-driven model to understand how predisposing factors influence the spatial variation of the triggering factor (rainfall intensity) in the case of flash floods in the continental area of Portugal. Flash floods occurrences were extracted from the DISASTER database. We extracted the accumulated precipitation from the Copernicus database by considering two days of duration. The analysed predisposing factors for flooding were extracted considering the whole basin where each occurrence is located. These factors include the basin area, the predominant lithology, drainage density, and the mean or median values of elevation, slope, stream power index (SPI), topographic wetness index (TWI), roughness, and four soil properties. The Random Forest algorithm was used to build the models and obtained mean absolute percentage error (MAPE) around 19%, an acceptable value for the objectives of the work. The median of SPI, mean elevation and the area of the basin are the top three most relevant predisposing factors interpreted by the model for defining the rainfall input for flash flooding in mainland Portugal.

Canopy reflectance (CR) models describe the transfer and interaction of radiation from the soil background to the canopy layer and play a vital role in the retrieval of biophysical variables. However, few efforts have focused on estimating soil background scattering operators, resulting in uncertainties in CR modelling, especially over sloping terrain. This study developed a canopy reflectance model for simulating CR over sloping terrain, which combines the general spectral vector (GSV) model, the PROSPECT model, and 4SAIL model coupled with topography (GSV-PROSAILT). The canopy reflectance simulated by GSV-PROSAILT was validated against two datasets: discrete anisotropic radiative transfer (DART) simulations and remote sensing observations. A comparison with DART simulations under various conditions revealed that the GSV-PROSAILT model captures terrain-induced CR distortion with high accuracy (red band: coefficient of determination $\lpar {\rm R 2} \rpar = 0.731$(R2)=0.731, root-mean-square error (RMSE) = 0.007; near infrared (NIR) band: $\rm R2 = 0.8319$R2=0.8319, RMSE = 0.0098). The results of remote sensing observation verification revealed that the GSV-PROSAILT model can be successfully used in CR modelling. These validations confirmed the performance of GSV-PROSAILT in soil and canopy reflectance modelling over sloping terrain, indicating that it can provide a potential tool for biophysical variable retrieval over mountainous areas.

The root-knot nematode, Meloidogyne javanica, is one of the most damaging plant-parasitic nematodes, affecting chickpea and causing substantial yield losses worldwide. The damage potential and population dynamics of this nematode in chickpea in Ethiopia have yet to be investigated. In this study, six chickpea cultivars were tested using 12 ranges of initial population densities (Pi) of M. javanica second-stage juveniles (J2): 0, 0.125, 0.25, 0.5, 1, 2, 4, 8, 16, 32, 64 and 128 J2 (g dry soil)-1 in a controlled glasshouse pot experiment. The Seinhorst yield loss and population dynamics models were fitted to describe population development and the effect on different measured growth variables. The tolerance limit (TTFW) for total fresh weight ranged from 0.05 to 1.22 J2 (g dry soil)-1, with corresponding yield losses ranging from 31 to 64%. The minimum yield for seed weight (mSW) ranged from 0.29 to 0.61, with estimated yield losses of 71 and 39%. The 'Haberu' and 'Geletu' cultivars were considered good hosts, with maximum population densities (M) of 16.27 and 5.64 J2 (g dry soil)-1 and maximum multiplication rate (a) values of 6.25 and 9.23, respectively. All other cultivars are moderate hosts for M. javanica; therefore, it is crucial to initiate chickpea-breeding strategies to manage the tropical root-knot nematode M. javanica in Ethiopia.

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.

The existence of rock weathering products has an important effect on the infiltration of water in the soil. Understanding the mechanism of water infiltration in a mixed soil and weathered rock debris medium is highly important for soil science and hydrology. The purpose of this study is to explore the effects of mudstone hydrolysis on water infiltration in the soil under different mixing ratios (0-70 %) of weathered mudstone contents. Soil column experiments and numerical modelling were used to study the processes of hydrolysis of weathered mudstone and water infiltration in the mixed medium. The results revealed that water immersion can cause the dense mudstone surface to fall off, thus forming pores, and that the amount of these pores first increase but then decrease over time. The disintegration of post-hydrolysis mudstone debris occurs mainly among particles ranging from 2-2000 mu m, predominantly transforming sand particles into finer fractions. Increasing the mudstone content in the soil from 0 % to 50 % enhances the infiltration rate and cumulative infiltration volume. However, when the mudstone content exceeds 50 %, these parameters decrease. The mudstone weathering products promote water infiltration in the soil within a certain range of mudstone contents, but as the ratio of weathered products increases, excessive amounts of mudstone hinder the movement of water in the soil. The identified transformation phenomenon suggests that the infiltration capacity of mixed soil will not scale linearly with mudstone content. The findings enable some mitigation strategies of geologic hazards based on the hydrological stability in heterogeneous environments.

The cyclic behavior of clay significantly influences the dynamic response of offshore wind turbines (OWTs). This study presents a practical bounding surface model capable of describing both cyclic shakedown and cyclic degradation. The model is characterized by a simple theoretical framework and a limited number of parameters, and it has been numerically implemented in ABAQUS through a user-defined material (UMAT) subroutine. The yield surface remains fixed at the origin with isotropic hardening, while a movable projection center is introduced to capture cyclic hysteresis behavior. Cumulative plastic deviatoric strain is integrated into the plastic modulus to represent cyclic accumulation. Validation against undrained cyclic tests on three types of clay demonstrates its capability in reproducing stress-strain hysteresis, cyclic shakedown, and cyclic degradation. Additionally, its effectiveness in solving finite element boundary value problems is verified through centrifuge tests on large-diameter monopiles. Furthermore, the model is adopted to analyze the dynamic response of monopile OWTs under seismic loading. The results indicate that, compared to cyclic shakedown, cyclic degradation leads to a progressive reduction in soil stiffness, which diminishes acceleration amplification, increases settlement accumulation, and results in higher residual excess pore pressure with greater fluctuation. Despite its advantages, this model requires a priori specification of the sign of the plastic modulus parameter cd to capture either cyclic degradation or shakedown behavior. Furthermore, under undrained conditions, the model leads pstabilization of the effective stress path, which subsequently results in underestimation of the excess pore pressure.

Accurately modeling soil-fluid coupling under large deformations is critical for understanding and predicting phenomena such as slope failures, embankment collapses, and other geotechnical hazards. This topic has been studied for decades and remains challenging due to the nonlinear responses of geotechnical structures, which typically result from plastic yielding and finite deformation of the soil skeleton. In this work, we comprehensively summarize the theory involved in the soil-fluid coupling problem. Within a finite strain framework, we employ an elasto-plastic constitutive model with linear hardening to represent the solid skeleton and a nearly incompressible model for water. The water content influences the behavior of the solid skeleton by affecting its cohesion. The governing equations are discretized by material point method and two sets of material points are employed to independently represent solid skeleton and fluid, respectively. The proposed method is validated by comparing simulation results with experimental results for the impact of water on dry soil and wet soil. The capability of the method is further demonstrated through two cases: (1) the impact of a rigid body on saturated soil, causing water seepage, and (2) the filling of a ditch, which considers the erosion of the foundation. This work may provide a versatile tool for analyzing the dynamic responses of fluid and solid interactions, considering both mixing and separation phenomena.

This paper presents a rigorous, semi-analytical solution for the drained cylindrical cavity expansion in transversely isotropic sand. The constitutive model used for the sand is the SANISAND-F model, which is developed within the anisotropic critical state theory framework that can account for the essential fabric anisotropy of soils. By introducing an auxiliary variable, the governing equations of the cylindrical expansion problem are transformed into a system of ten first-order ordinary differential equations. Three of these correspond to the stress components, three are associated with the kinematic hardening tensor, three describe the fabric tensor, and the last one represents the specific volume. The solution is validated through comparison with finite element analysis, using Toyoura sand as the reference material. Parametric analyses and discussion on the impact of initial void ratio, initial mean stress level, at-rest earth pressure coefficient and initial fabric anisotropy intensity are presented. The results demonstrate that the fabric anisotropy of sand significantly influences the distribution of stress components and void ratio around the cavity. When fabric anisotropy is considered, the solution predicts lower values of radial, circumferential and vertical stresses near the cavity wall compared to those obtained without considering fabric anisotropy. The proposed solution is expected to enhance the accuracy of cavity expansion predictions in sand, which will have significant practical applications, including interpreting pressuremeter tests, predicting effects of driven pile installation, and improving the understanding of sand mechanics under complex loading scenarios.

In geotechnical engineering, the development of efficient and accurate constitutive models for granular soils is crucial. The micromechanical models have gained much attention for their capacity to account for particle-scale interactions and fabric anisotropy, while requiring far less computational resources compared to discrete element method. Various micromechanical models have been proposed in the literature, but none of them have been conclusively shown to agree with the critical state theory given theoretical proof, despite the authors described that their models approximately reach the critical state. This paper modifies the previous CHY micromechanical model that is compatible with the critical state theory based on the assumption that the microscopic force-dilatancy relationship should align with the macroscopic stress-dilatancy relationship. Moreover, under the framework of the CHY model, the fabric anisotropy can be easily considered and the anisotropic critical state can be achieved with the introduction of the fabric evolution law. The model is calibrated using drained and undrained triaxial experiments and the results show that the model reliably replicates the mechanical behaviors of granular materials under both drained and undrained conditions. The compatibility of the model with the critical state theory is verified at both macroscopic and microscopic scales.