Pub Date : 2025-12-17DOI: 10.1016/j.enggeo.2025.108516
Zhengbin Liu , Shuai Wang , Shuwei Wu , Jianbo Guo , Yiwei Mao , Zeren Chen , Qingxue Huang
Accurately setting the friction coefficient between rock particles is a critical prerequisite for ensuring the validity of dynamic mechanical behavior simulations of rocks. The geometric and physical parameters of rock particles have complex effects on the friction coefficient. However, existing calibration methods often have limitations in terms of precision, efficiency, and applicability. To address these issues, this study proposes a novel calibration method for the friction coefficient of rock particles, which integrates sphero-polyhedron modeling techniques with a data-driven strategy. The method uses the angle of repose (AOR) as a reference for quantitative analysis, considering the influence of the particle geometric parameters and material physical properties on the friction coefficient. By constructing a discrete element simulation database and generating a sample dataset, a mapping relationship is established with AOR and vertical aspect ratios as inputs, and the static friction coefficient, dynamic friction coefficient, and rolling resistance coefficient as outputs. This enables rapid calibration of the friction coefficient through a data-driven approach. The experimental results show that the proposed method not only achieves excellent accuracy but also demonstrates strong generalizability, providing a new approach for determining the friction coefficient in rock particle simulation analysis and offering valuable support for geotechnical engineering analysis.
{"title":"A data-driven calibration method for the friction coefficients between rock particles","authors":"Zhengbin Liu , Shuai Wang , Shuwei Wu , Jianbo Guo , Yiwei Mao , Zeren Chen , Qingxue Huang","doi":"10.1016/j.enggeo.2025.108516","DOIUrl":"10.1016/j.enggeo.2025.108516","url":null,"abstract":"<div><div>Accurately setting the friction coefficient between rock particles is a critical prerequisite for ensuring the validity of dynamic mechanical behavior simulations of rocks. The geometric and physical parameters of rock particles have complex effects on the friction coefficient. However, existing calibration methods often have limitations in terms of precision, efficiency, and applicability. To address these issues, this study proposes a novel calibration method for the friction coefficient of rock particles, which integrates sphero-polyhedron modeling techniques with a data-driven strategy. The method uses the angle of repose (AOR) as a reference for quantitative analysis, considering the influence of the particle geometric parameters and material physical properties on the friction coefficient. By constructing a discrete element simulation database and generating a sample dataset, a mapping relationship is established with AOR and vertical aspect ratios as inputs, and the static friction coefficient, dynamic friction coefficient, and rolling resistance coefficient as outputs. This enables rapid calibration of the friction coefficient through a data-driven approach. The experimental results show that the proposed method not only achieves excellent accuracy but also demonstrates strong generalizability, providing a new approach for determining the friction coefficient in rock particle simulation analysis and offering valuable support for geotechnical engineering analysis.</div></div>","PeriodicalId":11567,"journal":{"name":"Engineering Geology","volume":"361 ","pages":"Article 108516"},"PeriodicalIF":8.4,"publicationDate":"2025-12-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145785691","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-17DOI: 10.1016/j.enggeo.2025.108515
Ruian Wu , Yongshuang Zhang , Chang Qi , Wenbo Zhao , Xiang Li , Deguang Song , Haishan Ma , Qijun Zou
Large-scale ancient landslides in the Himalayan region are increasingly susceptible to reactivation due to climate change and intensifying engineering activities, posing catastrophic geohazard risks. This study deciphers the complete failure chain of the Pangcun ancient landslide (∼18.9 × 106 m3) in Tibet, employing a multi-methodological approach that integrates remote sensing, field investigation, geotechnical testing, and numerical modeling. Our findings reveal a composite failure mechanism characterized by initial retrogressive deformation followed by thrust-style propagation. The reactivation manifests as a creep-slip process within the accumulation mass at depths of 6–25 m, where toe excavation induced early-stage retrogressive cracking, while subsequent rainfall infiltration triggered a thrust-style failure pushing from the rear. Stability analysis quantitatively confirms this vulnerability, showing the Factor of Safety (FoS) decreasing from a marginally stable 1.043 under natural conditions to an unstable 0.951 during heavy rainfall. Furthermore, post-failure simulations predict that a shallow failure could evolve into a high-speed event, reaching peak velocities of up to 17.8 m/s and a runout distance of 840 m, thereby directly endangering the G219 National Highway and downstream communities. Ultimately, this study provides a robust mechanistic framework for assessing similar ancient landslides, facilitating a critical shift in hazard management from reactive response to proactive, mechanism-based prevention.
{"title":"Ancient landslide on the Tibet Plateau(China): Reactivation mechanism and post-failure behavior prediction","authors":"Ruian Wu , Yongshuang Zhang , Chang Qi , Wenbo Zhao , Xiang Li , Deguang Song , Haishan Ma , Qijun Zou","doi":"10.1016/j.enggeo.2025.108515","DOIUrl":"10.1016/j.enggeo.2025.108515","url":null,"abstract":"<div><div>Large-scale ancient landslides in the Himalayan region are increasingly susceptible to reactivation due to climate change and intensifying engineering activities, posing catastrophic geohazard risks. This study deciphers the complete failure chain of the Pangcun ancient landslide (∼18.9 × 10<sup>6</sup> m<sup>3</sup>) in Tibet, employing a multi-methodological approach that integrates remote sensing, field investigation, geotechnical testing, and numerical modeling. Our findings reveal a composite failure mechanism characterized by initial retrogressive deformation followed by thrust-style propagation. The reactivation manifests as a creep-slip process within the accumulation mass at depths of 6–25 m, where toe excavation induced early-stage retrogressive cracking, while subsequent rainfall infiltration triggered a thrust-style failure pushing from the rear. Stability analysis quantitatively confirms this vulnerability, showing the Factor of Safety (FoS) decreasing from a marginally stable 1.043 under natural conditions to an unstable 0.951 during heavy rainfall. Furthermore, post-failure simulations predict that a shallow failure could evolve into a high-speed event, reaching peak velocities of up to 17.8 m/s and a runout distance of 840 m, thereby directly endangering the G219 National Highway and downstream communities. Ultimately, this study provides a robust mechanistic framework for assessing similar ancient landslides, facilitating a critical shift in hazard management from reactive response to proactive, mechanism-based prevention.</div></div>","PeriodicalId":11567,"journal":{"name":"Engineering Geology","volume":"361 ","pages":"Article 108515"},"PeriodicalIF":8.4,"publicationDate":"2025-12-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145785695","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-16DOI: 10.1016/j.enggeo.2025.108509
Furong Liu , Wei Ma , Yanhu Mu , Zhi Wen , Mingde Shen , Pengfei He
Under the background of global climate changing, the warming of permafrost has led to numerous engineering infrastructures being operated on warm permafrost foundations with diminishing bearing capacity. Meanwhile, infrastructure construction not only increases the overburden load on permafrost foundations but also induces directional deviation of the principal stress axis relative to the vertical direction. Therefore, conducting study on the stress-strain behavior and strength characteristics along different principal stress directions in warm frozen soils is imperative for accurately assessing deformation evolution patterns and bearing capacity of warm permafrost foundations. Thus, the stress-strain relationships respond, non-coaxiality evolution and strength distribution characteristics during directional loading along different principal stress directions were systematically investigated. The results indicated that the influence of principal stress direction on the strength intensifies with decreasing initial mean principal stress (when p₀ = 500 kPa, the strength at α = 45° exhibits a 27.3 % reduction compared to the α = 0°). Concurrently, increasing initial mean principal stress diminishes both the stress-strain non-coaxiality angle and the directional dependence of strength. Furthermore, a novel strength model incorporating principal stress direction is proposed for warm frozen silt. These findings elucidate the correlation mechanisms between non-coaxiality evolution and strength anisotropy in warm frozen silt under fixed principal stress direction, providing theoretical foundations for optimizing engineering designs in permafrost regions under warming scenarios.
{"title":"Strength and non-coaxiality behavior of warm frozen silt under inclined principal stress axes","authors":"Furong Liu , Wei Ma , Yanhu Mu , Zhi Wen , Mingde Shen , Pengfei He","doi":"10.1016/j.enggeo.2025.108509","DOIUrl":"10.1016/j.enggeo.2025.108509","url":null,"abstract":"<div><div>Under the background of global climate changing, the warming of permafrost has led to numerous engineering infrastructures being operated on warm permafrost foundations with diminishing bearing capacity. Meanwhile, infrastructure construction not only increases the overburden load on permafrost foundations but also induces directional deviation of the principal stress axis relative to the vertical direction. Therefore, conducting study on the stress-strain behavior and strength characteristics along different principal stress directions in warm frozen soils is imperative for accurately assessing deformation evolution patterns and bearing capacity of warm permafrost foundations. Thus, the stress-strain relationships respond, non-coaxiality evolution and strength distribution characteristics during directional loading along different principal stress directions were systematically investigated. The results indicated that the influence of principal stress direction on the strength intensifies with decreasing initial mean principal stress (when <em>p</em><sub><em>₀</em></sub> = 500 kPa, the strength at <em>α</em> = 45° exhibits a 27.3 % reduction compared to the <em>α</em> = 0°). Concurrently, increasing initial mean principal stress diminishes both the stress-strain non-coaxiality angle and the directional dependence of strength. Furthermore, a novel strength model incorporating principal stress direction is proposed for warm frozen silt. These findings elucidate the correlation mechanisms between non-coaxiality evolution and strength anisotropy in warm frozen silt under fixed principal stress direction, providing theoretical foundations for optimizing engineering designs in permafrost regions under warming scenarios.</div></div>","PeriodicalId":11567,"journal":{"name":"Engineering Geology","volume":"361 ","pages":"Article 108509"},"PeriodicalIF":8.4,"publicationDate":"2025-12-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145785697","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-16DOI: 10.1016/j.enggeo.2025.108512
Ming Wei , Jinlai Zhu , Zhen Guo , Wen Zhang , Linpeng Qin , Zongzheng Li , Xiaoyan Wang , Qi Sun
Traditional microseismic location methods face severe limitations in complex mountainous terrain due to oversimplified velocity assumptions and neglect of topographic effects, often yielding location errors exceeding 20–30 m. This case study demonstrates how high-precision 3D seismic event (SE) location can be achieved in such challenging environments through two key methodological innovations: (1) incorporation of complex stratigraphic structures using high-resolution 3D velocity models derived from dense array surface wave tomography (SWT), capturing velocity variations from ∼200–2500 m/s characteristic of weathered and fractured slope masses; and (2) integration of topographic effects through fast marching ray tracing within DEM-constrained domains, computing physically realistic wave paths that honor both velocity structure and terrain geometry. Application to actively deforming slopes in the Hengduan Mountains of eastern Tibet—where extreme topographic relief (>700 m) and complex geological structures exemplify the challenges confronting conventional methods—demonstrates location accuracies of 3 m overall and 1.5 m within dense array coverage areas. The excellent agreement between SWT-derived velocity structures and independent geological observations from boreholes and field mapping confirms the physical validity of the wave propagation models. Furthermore, analysis of 1470 SEs located over one year reveals shallow microseismic activity (0–43 m depth) concentrated within zones of maximum surface deformation identified by interferometric synthetic aperture radar (InSAR), with characteristic frequencies of 4–9 Hz and balanced energy distributions indicative of continuous creeping behavior. The strong spatial correlation between located SE clusters and independently measured surface deformation validates that our dual consideration of complex strata and topographic effects successfully captures the true subsurface source distribution. This methodology provides the spatial resolution essential for reliable slope stability assessment in complex geological settings.
{"title":"High-precision 3D seismic event (SE) location method for slopes incorporating complex strata and topographic effects: A case study of creeping slopes in the Hengduan Mountains, Eastern Tibet","authors":"Ming Wei , Jinlai Zhu , Zhen Guo , Wen Zhang , Linpeng Qin , Zongzheng Li , Xiaoyan Wang , Qi Sun","doi":"10.1016/j.enggeo.2025.108512","DOIUrl":"10.1016/j.enggeo.2025.108512","url":null,"abstract":"<div><div>Traditional microseismic location methods face severe limitations in complex mountainous terrain due to oversimplified velocity assumptions and neglect of topographic effects, often yielding location errors exceeding 20–30 m. This case study demonstrates how high-precision 3D seismic event (SE) location can be achieved in such challenging environments through two key methodological innovations: (1) incorporation of complex stratigraphic structures using high-resolution 3D velocity models derived from dense array surface wave tomography (SWT), capturing velocity variations from ∼200–2500 m/s characteristic of weathered and fractured slope masses; and (2) integration of topographic effects through fast marching ray tracing within DEM-constrained domains, computing physically realistic wave paths that honor both velocity structure and terrain geometry. Application to actively deforming slopes in the Hengduan Mountains of eastern Tibet—where extreme topographic relief (>700 m) and complex geological structures exemplify the challenges confronting conventional methods—demonstrates location accuracies of 3 m overall and 1.5 m within dense array coverage areas. The excellent agreement between SWT-derived velocity structures and independent geological observations from boreholes and field mapping confirms the physical validity of the wave propagation models. Furthermore, analysis of 1470 SEs located over one year reveals shallow microseismic activity (0–43 m depth) concentrated within zones of maximum surface deformation identified by interferometric synthetic aperture radar (InSAR), with characteristic frequencies of 4–9 Hz and balanced energy distributions indicative of continuous creeping behavior. The strong spatial correlation between located SE clusters and independently measured surface deformation validates that our dual consideration of complex strata and topographic effects successfully captures the true subsurface source distribution. This methodology provides the spatial resolution essential for reliable slope stability assessment in complex geological settings.</div></div>","PeriodicalId":11567,"journal":{"name":"Engineering Geology","volume":"361 ","pages":"Article 108512"},"PeriodicalIF":8.4,"publicationDate":"2025-12-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145785696","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-15DOI: 10.1016/j.enggeo.2025.108510
Jinxi Liang, Wanghua Sui, Ming Ye, Sara Kasmaeeyazdi, Francesco Tinti
Water-sand mixture inrush (WSMI) events pose severe threats to mining safety, infrastructure stability, and subsurface operations. This study first develops a pathway loss model to integrating frictional and expansion-induced hydraulic head losses, and then applies the Sobol-based global sensitivity analysis (GSA) to the model to evaluate WSMI risk for the following two scenarios (1) direct pathway-induced WSMI (with short, gravity-driven pathways) and (2) indirect or combined pathway-induced WSMI (with long, complex, pressure-driven pathways). For the two scenarios, GSA identifies fluid velocity as the dominant parameter, with pathway expansion loss governing direct inrush and friction loss dominating indirect inrush. Hydraulic head loss is markedly higher in the indirect inrush scenario than in the direct inrush scenario. Accordingly, tailored mitigation strategies are developed. For the direct inrush scenario (simple pathways), the priority is to cut off the energy conversion chain; for indirect inrush scenario (complex pathways), the focus is on dissipating excess energy. These findings advance the mechanistic understanding of WSMI and offer scenario-specific guidance for hazard control.
{"title":"Water-sand mixture inrush in underground pathways: Risk factors and mitigation strategies","authors":"Jinxi Liang, Wanghua Sui, Ming Ye, Sara Kasmaeeyazdi, Francesco Tinti","doi":"10.1016/j.enggeo.2025.108510","DOIUrl":"https://doi.org/10.1016/j.enggeo.2025.108510","url":null,"abstract":"Water-sand mixture inrush (WSMI) events pose severe threats to mining safety, infrastructure stability, and subsurface operations. This study first develops a pathway loss model to integrating frictional and expansion-induced hydraulic head losses, and then applies the Sobol-based global sensitivity analysis (GSA) to the model to evaluate WSMI risk for the following two scenarios (1) direct pathway-induced WSMI (with short, gravity-driven pathways) and (2) indirect or combined pathway-induced WSMI (with long, complex, pressure-driven pathways). For the two scenarios, GSA identifies fluid velocity as the dominant parameter, with pathway expansion loss governing direct inrush and friction loss dominating indirect inrush. Hydraulic head loss is markedly higher in the indirect inrush scenario than in the direct inrush scenario. Accordingly, tailored mitigation strategies are developed. For the direct inrush scenario (simple pathways), the priority is to cut off the energy conversion chain; for indirect inrush scenario (complex pathways), the focus is on dissipating excess energy. These findings advance the mechanistic understanding of WSMI and offer scenario-specific guidance for hazard control.","PeriodicalId":11567,"journal":{"name":"Engineering Geology","volume":"23 1","pages":""},"PeriodicalIF":7.4,"publicationDate":"2025-12-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145753363","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-15DOI: 10.1016/j.enggeo.2025.108511
Chuanxiang Qu , Yutong Liu , Haowen Guo , Leilei Liu
Probabilistic stability analysis of unsaturated soil slope with spatial variability under rainfall infiltration is computationally intensive due to highly non-linear behaviour and numerous repeated computations. In the field, unsaturated soil typically experiences specific stress states, and these stress levels can influence soil water capacity, thereby affecting slope stability. However, such stress effects have rarely been considered in previous probabilistic analyses of unsaturated soil slope stability. The relative importance of stress effects and spatial variability on slope stability remains unclear. To tackle these issues, a convolutional neural network with Bayesian optimisation (CNNB) is proposed as a surrogate algorithm. A completely decomposed tuff (CDT) slope, which is commonly observed in Hong Kong, serves as an example. Stress effects are characterised by a stress-dependent water retention model that effectively captures the influence of stress on water capacity at any given stress level. The spatially varying soil hydraulic and mechanical parameters of the slope are simulated by multivariate cross-correlated random fields. It is found that the proposed CNNB considerably enhances computational efficiency by at least 7.7 times compared to the random finite element method combined with the random limit equilibrium method (RFEM-RLEM). Meanwhile, it maintains a reliable probability of failure (pf) assessment with a prediction error as low as 2.9 %. Ignoring stress effects underestimates pf of the slope by up to 90 % under rainfall in Hong Kong with a 100-year return period. Stress effects have a more significant influence than spatial variability when computing the factor of safety (FOS) of the slope. Utilising deterministic analysis without stress effects as a benchmark, the difference in FOS due to stress effects is about 3.5 times that of spatial variability. Additionally, without considering spatial variability can also lead to unsafe assessments, as evidenced by a mean FOS value of 1.04 corresponding to a 22.6 % pf, indicating a hazardous performance level.
{"title":"Probabilistic analysis of stress effects on an unsaturated soil slope stability using convolutional neural networks and Bayesian optimisation","authors":"Chuanxiang Qu , Yutong Liu , Haowen Guo , Leilei Liu","doi":"10.1016/j.enggeo.2025.108511","DOIUrl":"10.1016/j.enggeo.2025.108511","url":null,"abstract":"<div><div>Probabilistic stability analysis of unsaturated soil slope with spatial variability under rainfall infiltration is computationally intensive due to highly non-linear behaviour and numerous repeated computations. In the field, unsaturated soil typically experiences specific stress states, and these stress levels can influence soil water capacity, thereby affecting slope stability. However, such stress effects have rarely been considered in previous probabilistic analyses of unsaturated soil slope stability. The relative importance of stress effects and spatial variability on slope stability remains unclear. To tackle these issues, a convolutional neural network with Bayesian optimisation (CNNB) is proposed as a surrogate algorithm. A completely decomposed tuff (CDT) slope, which is commonly observed in Hong Kong, serves as an example. Stress effects are characterised by a stress-dependent water retention model that effectively captures the influence of stress on water capacity at any given stress level. The spatially varying soil hydraulic and mechanical parameters of the slope are simulated by multivariate cross-correlated random fields. It is found that the proposed CNNB considerably enhances computational efficiency by at least 7.7 times compared to the random finite element method combined with the random limit equilibrium method (RFEM-RLEM). Meanwhile, it maintains a reliable probability of failure (<em>p</em><sub>f</sub>) assessment with a prediction error as low as 2.9 %. Ignoring stress effects underestimates <em>p</em><sub>f</sub> of the slope by up to 90 % under rainfall in Hong Kong with a 100-year return period. Stress effects have a more significant influence than spatial variability when computing the factor of safety (FOS) of the slope. Utilising deterministic analysis without stress effects as a benchmark, the difference in FOS due to stress effects is about 3.5 times that of spatial variability. Additionally, without considering spatial variability can also lead to unsafe assessments, as evidenced by a mean FOS value of 1.04 corresponding to a 22.6 % <em>p</em><sub>f</sub>, indicating a hazardous performance level.</div></div>","PeriodicalId":11567,"journal":{"name":"Engineering Geology","volume":"361 ","pages":"Article 108511"},"PeriodicalIF":8.4,"publicationDate":"2025-12-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145753480","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-12DOI: 10.1016/j.enggeo.2025.108508
Fanfan Yang , Renguang Zuo , Oliver P. Kreuzer
Data-driven deep learning approaches have exhibited promising performance in engineering geological mapping. However, existing methods face challenges in geological mapping based on multimodal data fusion due to their limited ability to exploit the complementary features among geoscience data. Moreover, the poor interpretability of deep learning methods limits their applicability for downstream engineering decision-making. To address these issues, this study designed a novel interpretable framework combining a contrastive multimodal graph attention network (CMGAT) with GNNExplainer (generating explanations for graph neural networks) for geological mapping. CMGAT was developed to extract discriminative features from multimodal graphs and align cross-modal representations via contrastive learning, while GNNExplainer was applied to quantify the influence of graph structure and geological features on the identification of geological units. The proposed CMGAT outperformed other unimodal models, achieving overall accuracies of 91 % and 82.9 % in lithological and fault mapping, respectively, in southwestern Fujian Province of China. Moreover, the GNNExplainer analysis identified key graph structure and geological indicators for geological unit delineation, strengthening the credibility of the predictive results. The framework can be further extended to diverse engineering geological mapping tasks.
{"title":"Interpretable regional-scale geological mapping using a contrastive graph attention network for multimodal data fusion and recognition of controlling factors","authors":"Fanfan Yang , Renguang Zuo , Oliver P. Kreuzer","doi":"10.1016/j.enggeo.2025.108508","DOIUrl":"10.1016/j.enggeo.2025.108508","url":null,"abstract":"<div><div>Data-driven deep learning approaches have exhibited promising performance in engineering geological mapping. However, existing methods face challenges in geological mapping based on multimodal data fusion due to their limited ability to exploit the complementary features among geoscience data. Moreover, the poor interpretability of deep learning methods limits their applicability for downstream engineering decision-making. To address these issues, this study designed a novel interpretable framework combining a contrastive multimodal graph attention network (CMGAT) with GNNExplainer (generating explanations for graph neural networks) for geological mapping. CMGAT was developed to extract discriminative features from multimodal graphs and align cross-modal representations via contrastive learning, while GNNExplainer was applied to quantify the influence of graph structure and geological features on the identification of geological units. The proposed CMGAT outperformed other unimodal models, achieving overall accuracies of 91 % and 82.9 % in lithological and fault mapping, respectively, in southwestern Fujian Province of China. Moreover, the GNNExplainer analysis identified key graph structure and geological indicators for geological unit delineation, strengthening the credibility of the predictive results. The framework can be further extended to diverse engineering geological mapping tasks.</div></div>","PeriodicalId":11567,"journal":{"name":"Engineering Geology","volume":"361 ","pages":"Article 108508"},"PeriodicalIF":8.4,"publicationDate":"2025-12-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145731549","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-12DOI: 10.1016/j.enggeo.2025.108482
Philipp Frieß , Hervé Vicari , Brian McArdell , Amanda Åberg , Johan Gaume
When debris and mud flows pass through curved channels, centrifugal forces lead to a height difference – known as superelevation – between the inner and outer banks. Analytical models describe this phenomenon by relating the superelevation angle to flow speed. However, these models assume simplified flow dynamics, a linear flow free surface, and do not explicitly account for solid–fluid interactions, requiring an empirical correction factor. In this study, we perform fully depth-resolved SPH-DEM numerical experiments to investigate the influence of water content on superelevation in curved channels. DEM represents the coarse solid particles, while SPH models the fluid phase, including both fines and water. The model is first validated against laboratory-scale experiments of debris flow superelevation. A parametric study is then conducted by varying the water content in debris and mud flows. The results show that increased water content leads to higher flow velocity and thus greater superelevation. The transverse flow surface depends strongly on material composition: mud flows typically exhibit convex-downward profiles, whereas granular flows display concave-downward profiles. By balancing centrifugal forces with basal normal stresses, we establish a correlation between the empirical correction factor, water content, and flow-surface curvature. However, the numerical experiments also reveal significant spatial variability in the correction factor along the bend, indicating additional mechanisms – specifically, a run-up impact that promotes superelevation, and subsequent alternating transverse motions – that limit the applicability of this analytical approach. Finally, SPH-DEM simulations of a real debris flow event at Illgraben successfully reproduce the observed field data, demonstrating the ability of the model for large-scale applications.
{"title":"Two-phase SPH-DEM modeling of the superelevation phenomenon of debris and mud flows","authors":"Philipp Frieß , Hervé Vicari , Brian McArdell , Amanda Åberg , Johan Gaume","doi":"10.1016/j.enggeo.2025.108482","DOIUrl":"10.1016/j.enggeo.2025.108482","url":null,"abstract":"<div><div>When debris and mud flows pass through curved channels, centrifugal forces lead to a height difference – known as superelevation – between the inner and outer banks. Analytical models describe this phenomenon by relating the superelevation angle to flow speed. However, these models assume simplified flow dynamics, a linear flow free surface, and do not explicitly account for solid–fluid interactions, requiring an empirical correction factor. In this study, we perform fully depth-resolved SPH-DEM numerical experiments to investigate the influence of water content on superelevation in curved channels. DEM represents the coarse solid particles, while SPH models the fluid phase, including both fines and water. The model is first validated against laboratory-scale experiments of debris flow superelevation. A parametric study is then conducted by varying the water content in debris and mud flows. The results show that increased water content leads to higher flow velocity and thus greater superelevation. The transverse flow surface depends strongly on material composition: mud flows typically exhibit convex-downward profiles, whereas granular flows display concave-downward profiles. By balancing centrifugal forces with basal normal stresses, we establish a correlation between the empirical correction factor, water content, and flow-surface curvature. However, the numerical experiments also reveal significant spatial variability in the correction factor along the bend, indicating additional mechanisms – specifically, a run-up impact that promotes superelevation, and subsequent alternating transverse motions – that limit the applicability of this analytical approach. Finally, SPH-DEM simulations of a real debris flow event at Illgraben successfully reproduce the observed field data, demonstrating the ability of the model for large-scale applications.</div></div>","PeriodicalId":11567,"journal":{"name":"Engineering Geology","volume":"361 ","pages":"Article 108482"},"PeriodicalIF":8.4,"publicationDate":"2025-12-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145731842","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-11DOI: 10.1016/j.enggeo.2025.108504
Ming-Jen Lo , Tso-Ren Wu , Kenji Satake
On January 1, 2024, a powerful earthquake (M 7.6) struck the Noto Peninsula, Japan, triggering a tsunami in the Sea of Japan. In Toyama Bay, the tsunami arrived earlier than expected. This study investigates the 2024 Noto tsunami event by separately modeling three potential tsunami generation mechanisms: vertical displacement from fault motion, horizontal displacement, and submarine landslides. To enhance the accuracy of submarine landslide-induced tsunami modeling, a computational fluid dynamics model, SPLASH3D, is utilized to simulate the landslide dynamics and determine its duration. Subsequently, a temporally variable seabed motion is used as the initial condition for a tsunami simulation code, COMCOT, to generate a dynamic tsunami source. The simulation results indicate that the sliding process has a significant influence on the observed tsunami in Toyama Bay, producing waveforms that better match observations than those derived from the equivalent instantaneous initial free surface displacement method. The combined simulation of dynamic submarine landslides, vertical displacements from fault motion, and horizontal displacements of the Noto Peninsula closely matches the observed data, enabling a detailed analysis of each source's contribution to the anomalous tsunami. Simulation results indicate that the submarine landslide was responsible for the early arrival of the tsunami. The contributions of the vertical fault displacement and submarine landslide each account for approximately 45 % of the maximum wave height, elucidating the unexpectedly high tsunami wave height. Therefore, the risks posed by landslide-generated tsunamis constitute a critical issue that must be addressed in tsunami early warning and coastal engineering risk assessment.
{"title":"Contribution of time-evolving landslide sources to the anomalous tsunami observed in the 2024 Noto earthquake","authors":"Ming-Jen Lo , Tso-Ren Wu , Kenji Satake","doi":"10.1016/j.enggeo.2025.108504","DOIUrl":"10.1016/j.enggeo.2025.108504","url":null,"abstract":"<div><div>On January 1, 2024, a powerful earthquake (M 7.6) struck the Noto Peninsula, Japan, triggering a tsunami in the Sea of Japan. In Toyama Bay, the tsunami arrived earlier than expected. This study investigates the 2024 Noto tsunami event by separately modeling three potential tsunami generation mechanisms: vertical displacement from fault motion, horizontal displacement, and submarine landslides. To enhance the accuracy of submarine landslide-induced tsunami modeling, a computational fluid dynamics model, SPLASH3D, is utilized to simulate the landslide dynamics and determine its duration. Subsequently, a temporally variable seabed motion is used as the initial condition for a tsunami simulation code, COMCOT, to generate a dynamic tsunami source. The simulation results indicate that the sliding process has a significant influence on the observed tsunami in Toyama Bay, producing waveforms that better match observations than those derived from the equivalent instantaneous initial free surface displacement method. The combined simulation of dynamic submarine landslides, vertical displacements from fault motion, and horizontal displacements of the Noto Peninsula closely matches the observed data, enabling a detailed analysis of each source's contribution to the anomalous tsunami. Simulation results indicate that the submarine landslide was responsible for the early arrival of the tsunami. The contributions of the vertical fault displacement and submarine landslide each account for approximately 45 % of the maximum wave height, elucidating the unexpectedly high tsunami wave height. Therefore, the risks posed by landslide-generated tsunamis constitute a critical issue that must be addressed in tsunami early warning and coastal engineering risk assessment.</div></div>","PeriodicalId":11567,"journal":{"name":"Engineering Geology","volume":"361 ","pages":"Article 108504"},"PeriodicalIF":8.4,"publicationDate":"2025-12-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145731845","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-11DOI: 10.1016/j.enggeo.2025.108506
Li Fei, Michel Jaboyedoff, Tiggi Choanji, Marc-Henri Derron
Over the past two decades, accelerated rock wall retreat has become a growing concern due to its link to global warming. While most research has focused on high-altitude cryosphere and deglacial regions, rock wall retreat in low-elevation areas remains understudied, despite posing higher risks to infrastructure and public safety. To address this gap, we investigated a molasse rock wall at La Cornalle located in the subalpine region (Vaud, Switzerland), composed of interbedded marl and sandstone layers. Using monthly Structure from Motion (SfM) photogrammetry and terrestrial laser scanning (TLS), we established a detailed four-year rockfall inventory and examined it with meteorological factors, including precipitation (including the snow melting), air temperature, and evapotranspiration (ET), collected from two nearby weather stations. A total of 4051 rockfall events, with a cumulative volume of 285 m3, were recorded. The annual retreat rates for sandstones and marls were 35.6 mm/yr and 26.0 mm/yr, respectively, with newly exposed rock faces showing a higher retreat rate (43.8 mm/yr) for marls. Spatially, rockfalls were concentrated in steep, thinly bedded, and highly fractured zones, as well as around large sandstone overhangs. Temporally, rockfall frequency peaked during winter and wet spring-summer periods, with duration of rainfall emerging as the primary driver, as prolonged rain facilitates deep water infiltration and weakens the water-sensitive marl layers. Following an extreme heatwave in August 2022, a notable spike in small rockfall events was observed at the early autumn (from Mid-September to Mid-October), indicating that local climatic shifts, such as extreme heatwave (coupled drying and heating) followed by effective water input (wetting), can significantly destabilize rock walls. This study highlights the importance of understanding temporal variations in rockfall activity and rock wall retreat by incorporating geological and climatic factors to improve rockfall hazard assessments in low-elevation regions.
{"title":"Analysis of rockfall-induced retreat and influencing factors in a sandstone-marl interbedded rock wall in a low-elevation environment","authors":"Li Fei, Michel Jaboyedoff, Tiggi Choanji, Marc-Henri Derron","doi":"10.1016/j.enggeo.2025.108506","DOIUrl":"10.1016/j.enggeo.2025.108506","url":null,"abstract":"<div><div>Over the past two decades, accelerated rock wall retreat has become a growing concern due to its link to global warming. While most research has focused on high-altitude cryosphere and deglacial regions, rock wall retreat in low-elevation areas remains understudied, despite posing higher risks to infrastructure and public safety. To address this gap, we investigated a molasse rock wall at La Cornalle located in the subalpine region (Vaud, Switzerland), composed of interbedded marl and sandstone layers. Using monthly Structure from Motion (SfM) photogrammetry and terrestrial laser scanning (TLS), we established a detailed four-year rockfall inventory and examined it with meteorological factors, including precipitation (including the snow melting), air temperature, and evapotranspiration (ET), collected from two nearby weather stations. A total of 4051 rockfall events, with a cumulative volume of 285 m<sup>3</sup>, were recorded. The annual retreat rates for sandstones and marls were 35.6 mm/yr and 26.0 mm/yr, respectively, with newly exposed rock faces showing a higher retreat rate (43.8 mm/yr) for marls. Spatially, rockfalls were concentrated in steep, thinly bedded, and highly fractured zones, as well as around large sandstone overhangs. Temporally, rockfall frequency peaked during winter and wet spring-summer periods, with duration of rainfall emerging as the primary driver, as prolonged rain facilitates deep water infiltration and weakens the water-sensitive marl layers. Following an extreme heatwave in August 2022, a notable spike in small rockfall events was observed at the early autumn (from Mid-September to Mid-October), indicating that local climatic shifts, such as extreme heatwave (coupled drying and heating) followed by effective water input (wetting), can significantly destabilize rock walls. This study highlights the importance of understanding temporal variations in rockfall activity and rock wall retreat by incorporating geological and climatic factors to improve rockfall hazard assessments in low-elevation regions.</div></div>","PeriodicalId":11567,"journal":{"name":"Engineering Geology","volume":"361 ","pages":"Article 108506"},"PeriodicalIF":8.4,"publicationDate":"2025-12-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145732132","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}