Black carbon (BC) is a major short-lived climate pollutant (SLCP) with significant climate and environmentalhealth impacts. This review synthesizes critical advancements in the identification of emerging anthropogenic BC sources, updates to global warming potential (GWP) and global temperature potential (GTP) metrics, technical progress in characterization techniques, improvements in global-regional monitoring networks, emission inventory, and impact assessment methods. Notably, gas flaring, shipping, and urban waste burning have slowly emerged as dominant emission sources, especially in Asia, Eastern Europe, and Arctic regions. The updated GWP over 100 years for BC is estimated at 342 CO2-eq, compared to 658 CO2-eq in IPCC AR5. Recent CMIP6-based Earth System Models (ESMs) have improved attribution of BC's microphysics, identifying a 22 % increase in radiative forcing (RF) over hotspots like East Asia and Sub-Saharan Africa. Despite progress, challenges persist in monitoring network inter-comparability, emission inventory uncertainty, and underrepresentation of BC processes in ESMs. Future efforts could benefit from the integration of satellite data, artificial intelligence (AI)assisted methods, and harmonized protocols to improve BC assessment. Targeted mitigation strategies could avert up to four million premature deaths globally by 2030, albeit at a 17 % additional cost. These findings highlight BC's pivotal roles in near-term climate and sustainability policy.

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.

Vegetation greening across the Tibetan Plateau, a critical ecological response to climate warming and land-cover change, affects soil hydrothermal regimes, altering soil moisture (SM) and soil temperature (ST) dynamics. However, its effects on SM-ST coupling remain poorly understood. Using integrated field measurements from a vegetation-soil (V-S) network, reanalysis, and physics-based simulations, we quantify responses of SM, ST, and their coupling to vegetation changes across the Upper Brahmaputra (UB) basin, southern Tibetan Plateau. Results show that strong positive SM-ST correlations occur throughout 0-289 cm soil layers across the basin, consistent with the monsoon-driven co-occurrence of rainy and warm seasons. Spatially, SM-ST coupling strength exhibits pronounced spatial heterogeneity, demonstrating strongest coupling in central basin areas with weaker intensities in eastern and western regions. Overall, vegetation greening consistently induces soil warming and drying: as leaf area index (LAI) increases from 20 % to 180 % of its natural levels, SM (0-160 cm) declines by 15 % to 29 % due to enhanced evapotranspiration and root water uptake. Mean ST simultaneously increases by 1.4 +/- 0.9 degrees C. Crucially, sparsely vegetated regions sustain warming (1.4-2.1 degrees C), while densely vegetated areas transition from initial warming to gradual cooling. These findings advance our understanding of soil hydrothermal dynamics and their broader environmental impacts, improving climate model parameterizations and informing sustainable land management strategies in high-altitude ecosystems.

Thawing-triggered slope failures and landslides are becoming an increasing concern in cold regions due to the ongoing climate change. Predicting and understanding the behaviour of frozen soils under these changing conditions is therefore critical and has led to a growing interest in the research community. To address this challenge, we present the first mesh-free smoothed particle hydrodynamics (SPH) computational framework designed to handle the multi-phase and multi-physic coupled thermo-hydro-mechanical (THM) process in frozen soils, namely the THM-SPH computational framework. The frozen soil is considered a tri-phase mixture (i.e., soil, water and ice), whose governing equations are then established based on u-p-T formulations. A critical-state elasto-plastic Clay and Sand Model for Frozen soils (CASM-F), formulated in terms of solid-phase stress, is then introduced to describe the transition response and large deformation behaviour of frozen soils due to thawing action for the first time. Several numerical verifications and demonstrations highlight the usefulness of this advanced THM-SPH computational framework in addressing challenging problems involving thawing-induced large deformation and failures of slopes. The results indicate that our proposed single-layer, fully coupled THM-SPH model can predict the entire failure process of thawing-induced landslides, from the initiation to post-failure responses, capturing the complex interaction among multiple coupled phases. This represents a significant advancement in the numerical modelling of frozen soils and their thawing-induced failure mechanisms in cold regions.

Thawing permafrost alters climate not only through carbon emissions but also via energy-water feedback and atmospheric teleconnections. This review focuses on the Tibetan Plateau, where strong freeze-thaw cycles, intense radiation, and complex snow-vegetation interactions constitute non-carbon climate responses. We synthesize recent evidence that links freeze-thaw cycles, ground heat flux dynamics, and soil moisture hysteresis to latent heat feedback, monsoon modulation, and planetary wave anomalies. Across these pathways, both observational and simulation studies reveal consistent signals of feedback amplification and nonlinear threshold behavior. However, most Earth system models underrepresent these processes due to simplifications in freezethaw processes, snow-soil-vegetation coupling, and cross-seasonal memory effects. We conclude by identifying priority processes to better simulate multi-scale cryosphere-climate feedback, especially under continued climate warming in high-altitude regions.

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.