Qingci Qin, Chaoquan Luo, Kegang Li, Fa Dong, Naeem Abbas
Fractures are one of the most critical factors influencing the mechanical properties of rocks. To explore the quantitative description method of fracture networks in natural rock masses and achieve precise mathematical reconstruction of the stochastic distribution characteristics of natural fractures, the apparent fracture distribution characteristics of 499 fractures in 12 standard specimens were investigated. The results show that there is a significant correlation among the density, length, and dip angle of fractures in different directions. Based on this, a quantitative description method for fracture networks that considers the relationship between fracture density and orientation is proposed. Combined with the Fisher model, the probability density distribution function of natural fracture orientation density was established. Using this method, the equivalent numerical analysis model of rock specimens containing natural fractures was reconstructed on the RFPA platform, and numerical experiments of triaxial loading and unloading were performed on the reconstructed equivalent fractured rock numerical model. This method can quantitatively describe the geometric distribution characteristics of fractures in natural rocks. The modified Fisher model enables the visual reconstruction of natural fracture networks, offering an effective technical approach for building equivalent numerical analysis models of rocks with natural fractures. It is highly valuable for studying the coupled mechanical behavior of multi‐physical fields in natural fractured rocks and provides an equivalent analysis method for visualizing and analyzing the damage process of natural fractured rocks.
{"title":"Mathematical Reconstruction Analysis of a Natural Fracture Network in Rock Based on a Modified Fisher Model and Its Numerical Realization","authors":"Qingci Qin, Chaoquan Luo, Kegang Li, Fa Dong, Naeem Abbas","doi":"10.1002/nag.70205","DOIUrl":"https://doi.org/10.1002/nag.70205","url":null,"abstract":"Fractures are one of the most critical factors influencing the mechanical properties of rocks. To explore the quantitative description method of fracture networks in natural rock masses and achieve precise mathematical reconstruction of the stochastic distribution characteristics of natural fractures, the apparent fracture distribution characteristics of 499 fractures in 12 standard specimens were investigated. The results show that there is a significant correlation among the density, length, and dip angle of fractures in different directions. Based on this, a quantitative description method for fracture networks that considers the relationship between fracture density and orientation is proposed. Combined with the Fisher model, the probability density distribution function of natural fracture orientation density was established. Using this method, the equivalent numerical analysis model of rock specimens containing natural fractures was reconstructed on the RFPA platform, and numerical experiments of triaxial loading and unloading were performed on the reconstructed equivalent fractured rock numerical model. This method can quantitatively describe the geometric distribution characteristics of fractures in natural rocks. The modified Fisher model enables the visual reconstruction of natural fracture networks, offering an effective technical approach for building equivalent numerical analysis models of rocks with natural fractures. It is highly valuable for studying the coupled mechanical behavior of multi‐physical fields in natural fractured rocks and provides an equivalent analysis method for visualizing and analyzing the damage process of natural fractured rocks.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"42 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-12-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145731499","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Giulia Capati, Diana Salciarini, Alessandro F. Rotta Loria
The subsurface undergoes temperature variations in many situations due to anthropogenic and natural causes, which strongly influence the long‐term behavior of soils. These phenomena can involve temperature anomalies in the ground across distances of a few meters, as in the case of nuclear waste repositories, to distances encompassing entire cities, as in the case of subsurface urban heat islands. To date, a number of constitutive models have been proposed to capture the mechanics of soils under non‐isothermal conditions, with particular attention to fine‐grained soils due to their renowned sensitivity to temperature variations. However, most of the available models suffer from many constitutive parameters that hinder their applicability to the analysis of large and complex problems involving thermally induced deformations of fine‐grained soils. This study extends the classical Kelvin–Voigt model with a temperature‐dependent formulation for capturing the reversible or irreversible thermally induced deformations of fine‐grained soils, interpreted through the theory of thermally accelerated creep. Implemented in a finite element software and validated against experimental data, the model shows that the thermally induced deformations of fine‐grained soils are highly sensitive to the temperature variation rate, while they are little influenced by the magnitude of the applied mechanical loads. The proposed model effectively captures the complex, time‐dependent deformations of fine‐grained soils with only a few easily calibrated parameters, making it a practical tool for the long‐term analysis of thermally induced creep in such materials.
{"title":"Extended Kelvin–Voigt Model for Simulating Thermally Accelerated Creep in Fine‐Grained Soils","authors":"Giulia Capati, Diana Salciarini, Alessandro F. Rotta Loria","doi":"10.1002/nag.70190","DOIUrl":"https://doi.org/10.1002/nag.70190","url":null,"abstract":"The subsurface undergoes temperature variations in many situations due to anthropogenic and natural causes, which strongly influence the long‐term behavior of soils. These phenomena can involve temperature anomalies in the ground across distances of a few meters, as in the case of nuclear waste repositories, to distances encompassing entire cities, as in the case of subsurface urban heat islands. To date, a number of constitutive models have been proposed to capture the mechanics of soils under non‐isothermal conditions, with particular attention to fine‐grained soils due to their renowned sensitivity to temperature variations. However, most of the available models suffer from many constitutive parameters that hinder their applicability to the analysis of large and complex problems involving thermally induced deformations of fine‐grained soils. This study extends the classical Kelvin–Voigt model with a temperature‐dependent formulation for capturing the reversible or irreversible thermally induced deformations of fine‐grained soils, interpreted through the theory of thermally accelerated creep. Implemented in a finite element software and validated against experimental data, the model shows that the thermally induced deformations of fine‐grained soils are highly sensitive to the temperature variation rate, while they are little influenced by the magnitude of the applied mechanical loads. The proposed model effectively captures the complex, time‐dependent deformations of fine‐grained soils with only a few easily calibrated parameters, making it a practical tool for the long‐term analysis of thermally induced creep in such materials.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"15 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-12-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145731502","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
This study proposes an efficient modeling approach for shield tunnel joints based on user‐defined elements (UELs), which is applied to the analysis of transverse deformation in shield tunnels. The element stiffness matrix of longitudinal joints with irregular geometric configurations was derived under arbitrary loading conditions. The solution procedure of the stiffness matrix was then incorporated into the UEL algorithm framework. The developed element simultaneously accounts for the varying load state at longitudinal joints, geometric configurations of the joint, and material nonlinearity, as well as their coupling effect on the load‐bearing capacity of the joint. Subsequently, the feasibility and accuracy of the UEL‐based modeling approach were evaluated across multiple modeling scales. The results demonstrate that this element not only significantly simplifies the modeling process but also effectively captures the nonlinear characteristics of tunnel convergence deformation observed in practical engineering scenarios. In addition, the UEL‐based modeling approach was employed to investigate the transverse deformation of shield tunnels subjected to an extreme surcharge–unloading–grouting rehabilitation sequence. Lining displacement–subgrade reaction and in situ stress state were appropriately incorporated into the numerical model. The effects of initial soil stiffness and the coefficient of earth pressure at rest on tunnel transverse deformation were quantitatively examined, along with the relationship between tunnel convergence and joint rotation. The results reveal a flag‐shaped hysteresis in the moment‐rotation response of tunnel joints under unloading and grouting, underscoring the importance of precisely capturing the multi‐degree‐of‐freedom coupling in joint element for practical engineering applications.
{"title":"An Efficient User‐Defined Element Framework for Joint Behavior Simulation and Transverse Deformation Analysis in Shield Tunnels","authors":"Hanwen Ji, Ximin Hu, Yu Miao, Hongjun He","doi":"10.1002/nag.70206","DOIUrl":"https://doi.org/10.1002/nag.70206","url":null,"abstract":"This study proposes an efficient modeling approach for shield tunnel joints based on user‐defined elements (UELs), which is applied to the analysis of transverse deformation in shield tunnels. The element stiffness matrix of longitudinal joints with irregular geometric configurations was derived under arbitrary loading conditions. The solution procedure of the stiffness matrix was then incorporated into the UEL algorithm framework. The developed element simultaneously accounts for the varying load state at longitudinal joints, geometric configurations of the joint, and material nonlinearity, as well as their coupling effect on the load‐bearing capacity of the joint. Subsequently, the feasibility and accuracy of the UEL‐based modeling approach were evaluated across multiple modeling scales. The results demonstrate that this element not only significantly simplifies the modeling process but also effectively captures the nonlinear characteristics of tunnel convergence deformation observed in practical engineering scenarios. In addition, the UEL‐based modeling approach was employed to investigate the transverse deformation of shield tunnels subjected to an extreme surcharge–unloading–grouting rehabilitation sequence. Lining displacement–subgrade reaction and in situ stress state were appropriately incorporated into the numerical model. The effects of initial soil stiffness and the coefficient of earth pressure at rest on tunnel transverse deformation were quantitatively examined, along with the relationship between tunnel convergence and joint rotation. The results reveal a flag‐shaped hysteresis in the moment‐rotation response of tunnel joints under unloading and grouting, underscoring the importance of precisely capturing the multi‐degree‐of‐freedom coupling in joint element for practical engineering applications.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"144 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-12-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145731501","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Frequent occurrences of shield clogging during slurry shield tunneling emphasize the necessity of optimizing cutterhead scouring systems. Conventional studies largely rely on empirical design and simplified analyses, which fail to capture complex flow–solid interactions or enable quantitative optimization. To overcome these limitations, a computational fluid dynamics (CFD) model was developed to simulate jet flow formation in shield nozzles and the flow field in the cutterhead excavation zone. It allows to systematically investigate the effects of nozzle inlet velocity, slurry chamber pressure, and cutterhead rotation speed, promoting the optimization of nozzle geometry and arrangement. The optimized configuration has been validated in the Qingdao Second Submarine Tunnel. Results show that increasing inlet velocity enhances jet strength without altering the overall velocity distribution pattern. Cutterhead rotation generates a rotating flow field that intensifies scouring but causes jet deflection, while interference between central and adjacent nozzles limits the effective scouring area. By modifying the outlet geometry and eliminating ineffective main arm nozzles, both the effective scouring area ratio and scouring efficiency index were improved. The optimized configuration ensures stable and efficient advancement when the nozzle inlet velocity exceeds 4 m/s, effectively preventing shield clogging. This work provides a validated modeling framework and practical optimization strategy for enhancing slurry shield performance in complex geological conditions.
{"title":"Optimizing Cutterhead Scouring Systems for Large‐Diameter Slurry Shields: A Case Study of the Qingdao Second Submarine Tunnel","authors":"Hanbiao Zhu, Shuying Wang, Zihao Jin, Jiazheng Zhong, Xiangcou Zheng, Pengfei Liu","doi":"10.1002/nag.70201","DOIUrl":"https://doi.org/10.1002/nag.70201","url":null,"abstract":"Frequent occurrences of shield clogging during slurry shield tunneling emphasize the necessity of optimizing cutterhead scouring systems. Conventional studies largely rely on empirical design and simplified analyses, which fail to capture complex flow–solid interactions or enable quantitative optimization. To overcome these limitations, a computational fluid dynamics (CFD) model was developed to simulate jet flow formation in shield nozzles and the flow field in the cutterhead excavation zone. It allows to systematically investigate the effects of nozzle inlet velocity, slurry chamber pressure, and cutterhead rotation speed, promoting the optimization of nozzle geometry and arrangement. The optimized configuration has been validated in the Qingdao Second Submarine Tunnel. Results show that increasing inlet velocity enhances jet strength without altering the overall velocity distribution pattern. Cutterhead rotation generates a rotating flow field that intensifies scouring but causes jet deflection, while interference between central and adjacent nozzles limits the effective scouring area. By modifying the outlet geometry and eliminating ineffective main arm nozzles, both the effective scouring area ratio and scouring efficiency index were improved. The optimized configuration ensures stable and efficient advancement when the nozzle inlet velocity exceeds 4 m/s, effectively preventing shield clogging. This work provides a validated modeling framework and practical optimization strategy for enhancing slurry shield performance in complex geological conditions.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"367 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-12-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145731500","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
This study presents a semi‐analytical solution for modeling suspended sediment distribution in turbulent flows within ice‐covered channels under unsteady, non‐equilibrium conditions. The solution is derived using the generalized integral transform technique (GITT). Validation was performed against the cell‐centered finite volume method and existing experimental data. The results confirm high accuracy, supported by error analysis. Optimized parameter values were obtained through a hybrid genetic and interior point algorithm. Several underlying phenomena of particle‐turbulence interactions in ice‐covered channels are explored. The focus is on the influence of key sediment transport parameters on the time‐dependent evolution of vertical concentration profiles of suspended sediment particles. Key findings indicate that increasing the settling‐velocity correction coefficient raises sediment concentration profiles over time. In contrast, greater ice‐cover roughness reduces sediment suspension. Sensitivity analysis highlights the inverse of the Schmidt number as a critical factor. This novel application of GITT and variance‐based sensitivity analysis (VBSA) provides a detailed solution library, and serves as a benchmark for numerical models.
{"title":"Suspended Sediment Transport in Ice‐Covered Turbulent Flow: Semi‐Analytical Solution and Parametric Sensitivity","authors":"Sourav Hossain, Ashok Das, Sanjib Naskar, Sweta Narayan Sahu, Koeli Ghoshal","doi":"10.1002/nag.70192","DOIUrl":"https://doi.org/10.1002/nag.70192","url":null,"abstract":"This study presents a semi‐analytical solution for modeling suspended sediment distribution in turbulent flows within ice‐covered channels under unsteady, non‐equilibrium conditions. The solution is derived using the generalized integral transform technique (GITT). Validation was performed against the cell‐centered finite volume method and existing experimental data. The results confirm high accuracy, supported by error analysis. Optimized parameter values were obtained through a hybrid genetic and interior point algorithm. Several underlying phenomena of particle‐turbulence interactions in ice‐covered channels are explored. The focus is on the influence of key sediment transport parameters on the time‐dependent evolution of vertical concentration profiles of suspended sediment particles. Key findings indicate that increasing the settling‐velocity correction coefficient raises sediment concentration profiles over time. In contrast, greater ice‐cover roughness reduces sediment suspension. Sensitivity analysis highlights the inverse of the Schmidt number as a critical factor. This novel application of GITT and variance‐based sensitivity analysis (VBSA) provides a detailed solution library, and serves as a benchmark for numerical models.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"9 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-12-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145731503","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Yazhen Sun, Lifan Yuan, Jinchang Wang, Longyan Wang, Youlin Ye
A novel analytical approach based on the state‐space method (SSM) is proposed to facilitate the rapid and accurate analysis of the mechanical behavior of pipe umbrellas in shallow‐buried tunnels. First, a mechanical analysis model and its governing equations are developed based on the Winkler elastic foundation beam theory, incorporating key construction‐related factors such as the lag effect of primary support, differential stress release in the surrounding rock, and the elasto‐plastic behavior of the stratum. The governing equations are then transformed into a concise matrix form using the SSM and solved based on matrix theory and the continuity conditions between adjacent beam segments. Analytical solutions for deformation and internal forces at any cross‐section of the pipe umbrella are derived under the nonlinear interactions among the pipe umbrella, surrounding rock, and primary support. The accuracy and applicability of the proposed method are verified through comparisons with existing field monitoring data, analytical solutions, and numerical simulation results from other researchers. On this basis, taking the Aketepu Tunnel as a case study, the influences of the surrounding rock stiffness ahead of the tunnel face, excavation footage, steel pipe diameter, and other factors on the deformation and internal forces of the pipe umbrella were investigated, and a recommended design scheme was proposed accordingly. The scheme was subsequently applied in the field, and monitoring results showed that the maximum deformations at two sections were 26.1 and 22.3 mm, respectively, both within the acceptable limits for surrounding rock deformation control.
{"title":"A Novel Approach to Analyzing the Mechanical Response of Pipe Umbrella Support in Shallow Buried Tunnels Considering Construction Characteristics","authors":"Yazhen Sun, Lifan Yuan, Jinchang Wang, Longyan Wang, Youlin Ye","doi":"10.1002/nag.70196","DOIUrl":"https://doi.org/10.1002/nag.70196","url":null,"abstract":"A novel analytical approach based on the state‐space method (SSM) is proposed to facilitate the rapid and accurate analysis of the mechanical behavior of pipe umbrellas in shallow‐buried tunnels. First, a mechanical analysis model and its governing equations are developed based on the Winkler elastic foundation beam theory, incorporating key construction‐related factors such as the lag effect of primary support, differential stress release in the surrounding rock, and the elasto‐plastic behavior of the stratum. The governing equations are then transformed into a concise matrix form using the SSM and solved based on matrix theory and the continuity conditions between adjacent beam segments. Analytical solutions for deformation and internal forces at any cross‐section of the pipe umbrella are derived under the nonlinear interactions among the pipe umbrella, surrounding rock, and primary support. The accuracy and applicability of the proposed method are verified through comparisons with existing field monitoring data, analytical solutions, and numerical simulation results from other researchers. On this basis, taking the Aketepu Tunnel as a case study, the influences of the surrounding rock stiffness ahead of the tunnel face, excavation footage, steel pipe diameter, and other factors on the deformation and internal forces of the pipe umbrella were investigated, and a recommended design scheme was proposed accordingly. The scheme was subsequently applied in the field, and monitoring results showed that the maximum deformations at two sections were 26.1 and 22.3 mm, respectively, both within the acceptable limits for surrounding rock deformation control.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"20 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-12-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145711206","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Bosong Ding, Ping Lou, You Wang, Can Huang, Rui Wang
Shield tunnels in soft clay are highly susceptible to long‐term uplift and differential settlement. To precisely control these deformations, this study develops a novel anisotropic cavity expansion model integrated with Mindlin's solution, forming a coupled 2D/3D segment–soil interaction framework. Validated against field data from Fuzhou Metro Line 4, the model demonstrates superior accuracy over traditional isotropic approaches. The results quantitatively reveal that: (1) in single‐ring grouting, increasing the grouting burial depth from 5 to 15 m reduces invert uplift by 28%; (2) in planar multi‐hole grouting, symmetric‐spacing sequences effectively control segment displacements within 5 mm; and (3) in longitudinal multi‐ring grouting, symmetric injection reduces maximum stress fluctuations by 12.5% compared to staggered grouting under a 5% grouting volume. The study concludes that the anisotropic soil–grout interaction is the primary mechanism governing deformation, and the proposed model provides a foundation for optimizing grouting parameters to ensure tunnel stability.
{"title":"Grouting Parameter Effects on Structural Response: Precision Modeling of Soft Soil Tunnel‐Base Grouting via Cavity Expansion Theory and Mindlin Solution","authors":"Bosong Ding, Ping Lou, You Wang, Can Huang, Rui Wang","doi":"10.1002/nag.70203","DOIUrl":"https://doi.org/10.1002/nag.70203","url":null,"abstract":"Shield tunnels in soft clay are highly susceptible to long‐term uplift and differential settlement. To precisely control these deformations, this study develops a novel anisotropic cavity expansion model integrated with Mindlin's solution, forming a coupled 2D/3D segment–soil interaction framework. Validated against field data from Fuzhou Metro Line 4, the model demonstrates superior accuracy over traditional isotropic approaches. The results quantitatively reveal that: (1) in single‐ring grouting, increasing the grouting burial depth from 5 to 15 m reduces invert uplift by 28%; (2) in planar multi‐hole grouting, symmetric‐spacing sequences effectively control segment displacements within 5 mm; and (3) in longitudinal multi‐ring grouting, symmetric injection reduces maximum stress fluctuations by 12.5% compared to staggered grouting under a 5% grouting volume. The study concludes that the anisotropic soil–grout interaction is the primary mechanism governing deformation, and the proposed model provides a foundation for optimizing grouting parameters to ensure tunnel stability.","PeriodicalId":13786,"journal":{"name":"International Journal for Numerical and Analytical Methods in Geomechanics","volume":"1 1","pages":""},"PeriodicalIF":4.0,"publicationDate":"2025-12-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145711057","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Piao Miao, Fusheng Zha, Qiao Wang, Shijin Feng, Hongxin Chen, Lin Mu
The cover image is based on the article Innovative Sealing Approaches for Karst Flowing Groundwater: Numerical Analysis of Secondary Grouting Mechanism After Grout-Filled Bag Sealed by Piao Miao et al., https://doi.org/10.1002/nag.70092