Articles | Volume 16, issue 11
https://doi.org/10.5194/se-16-1335-2025
© Author(s) 2025. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
https://doi.org/10.5194/se-16-1335-2025
© Author(s) 2025. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Stress drop sequences in the simplest pressure-sensitive ideal elasto-plastic media: implications for earthquake cycles
Yury Alkhimenkov
CORRESPONDING AUTHOR
Department of Civil and Environmental Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
Lyudmila Khakimova
Skolkovo Institute of Science and Technology (Skoltech), Sikorskogo 11, Moscow 121205, Russia
Faculty of Mechanics and Mathematics, Lomonosov Moscow State University, Moscow 119991, Russia
Yury Y. Podladchikov
Institute of Earth Sciences, University of Lausanne, Lausanne, Switzerland
Faculty of Mechanics and Mathematics, Lomonosov Moscow State University, Moscow 119991, Russia
Related authors
Yury Alkhimenkov and Yury Y. Podladchikov
Solid Earth, 16, 1227–1247, https://doi.org/10.5194/se-16-1227-2025, https://doi.org/10.5194/se-16-1227-2025, 2025
Short summary
Short summary
We present a thermodynamically consistent derivation of extended Biot poroelasticity, showing that Gassmann, Brown–Korringa, and related models emerge as special cases. Our formulation clarifies the conditions under which Gassmann’s relation holds and extends it by incorporating off-diagonal Hessian terms. Symbolic Maple code with consistency checks ensures full transparency, reproducibility, and accessibility for further research.
Yury Alkhimenkov and Yury Y. Podladchikov
Geosci. Model Dev., 18, 563–583, https://doi.org/10.5194/gmd-18-563-2025, https://doi.org/10.5194/gmd-18-563-2025, 2025
Short summary
Short summary
The accelerated pseudo-transient (APT) method is an efficient way to solve partial differential equations, particularly well-suited for parallel computing. This paper explores the APT method's effectiveness in solving elastic, viscoelastic, and hydromechanical problems, focusing on quasi-static conditions in 1D, 2D, and 3D. The study examines the best numerical settings for fast and accurate solutions. The paper shows how the APT method can handle complex problems in high-resolution models.
Emmanuel Wyser, Yury Alkhimenkov, Michel Jaboyedoff, and Yury Y. Podladchikov
Geosci. Model Dev., 14, 7749–7774, https://doi.org/10.5194/gmd-14-7749-2021, https://doi.org/10.5194/gmd-14-7749-2021, 2021
Short summary
Short summary
We propose an implementation of the material point method using graphical processing units (GPUs) to solve elastoplastic problems in three-dimensional configurations, such as the granular collapse or the slumping mechanics, i.e., landslide. The computational power of GPUs promotes fast code executions, compared to a traditional implementation using central processing units (CPUs). This allows us to study complex three-dimensional problems tackling high spatial resolution.
Emmanuel Wyser, Yury Alkhimenkov, Michel Jaboyedoff, and Yury Y. Podladchikov
Geosci. Model Dev., 13, 6265–6284, https://doi.org/10.5194/gmd-13-6265-2020, https://doi.org/10.5194/gmd-13-6265-2020, 2020
Short summary
Short summary
In this work, we present an efficient and fast material point method (MPM) implementation in MATLAB. We first discuss the vectorization strategies to adapt this numerical method to a MATLAB implementation. We report excellent agreement of the solver compared with classical analysis among the MPM community, such as the cantilever beam problem. The solver achieves a performance gain of 28 compared with a classical iterative implementation.
Yury Alkhimenkov and Yury Y. Podladchikov
Solid Earth, 16, 1227–1247, https://doi.org/10.5194/se-16-1227-2025, https://doi.org/10.5194/se-16-1227-2025, 2025
Short summary
Short summary
We present a thermodynamically consistent derivation of extended Biot poroelasticity, showing that Gassmann, Brown–Korringa, and related models emerge as special cases. Our formulation clarifies the conditions under which Gassmann’s relation holds and extends it by incorporating off-diagonal Hessian terms. Symbolic Maple code with consistency checks ensures full transparency, reproducibility, and accessibility for further research.
Yury Alkhimenkov and Yury Y. Podladchikov
Geosci. Model Dev., 18, 563–583, https://doi.org/10.5194/gmd-18-563-2025, https://doi.org/10.5194/gmd-18-563-2025, 2025
Short summary
Short summary
The accelerated pseudo-transient (APT) method is an efficient way to solve partial differential equations, particularly well-suited for parallel computing. This paper explores the APT method's effectiveness in solving elastic, viscoelastic, and hydromechanical problems, focusing on quasi-static conditions in 1D, 2D, and 3D. The study examines the best numerical settings for fast and accurate solutions. The paper shows how the APT method can handle complex problems in high-resolution models.
Ludovic Räss, Ivan Utkin, Thibault Duretz, Samuel Omlin, and Yuri Y. Podladchikov
Geosci. Model Dev., 15, 5757–5786, https://doi.org/10.5194/gmd-15-5757-2022, https://doi.org/10.5194/gmd-15-5757-2022, 2022
Short summary
Short summary
Continuum mechanics-based modelling of physical processes at large scale requires huge computational resources provided by massively parallel hardware such as graphical processing units. We present a suite of numerical algorithms, implemented using the Julia language, that efficiently leverages the parallelism. We demonstrate that our implementation is efficient, scalable and robust and showcase applications to various geophysical problems.
Emmanuel Wyser, Yury Alkhimenkov, Michel Jaboyedoff, and Yury Y. Podladchikov
Geosci. Model Dev., 14, 7749–7774, https://doi.org/10.5194/gmd-14-7749-2021, https://doi.org/10.5194/gmd-14-7749-2021, 2021
Short summary
Short summary
We propose an implementation of the material point method using graphical processing units (GPUs) to solve elastoplastic problems in three-dimensional configurations, such as the granular collapse or the slumping mechanics, i.e., landslide. The computational power of GPUs promotes fast code executions, compared to a traditional implementation using central processing units (CPUs). This allows us to study complex three-dimensional problems tackling high spatial resolution.
Martin Franz, Michel Jaboyedoff, Ryan P. Mulligan, Yury Podladchikov, and W. Andy Take
Nat. Hazards Earth Syst. Sci., 21, 1229–1245, https://doi.org/10.5194/nhess-21-1229-2021, https://doi.org/10.5194/nhess-21-1229-2021, 2021
Short summary
Short summary
A landslide-generated tsunami is a complex phenomenon that involves landslide dynamics, wave dynamics and their interaction. This phenomenon threatens numerous lives and infrastructures around the world. To assess this natural hazard, we developed an efficient numerical model able to simulate the landslide, the momentum transfer and the wave all at once. The good agreement between the numerical simulations and physical experiments validates our model and its novel momentum transfer approach.
Emmanuel Wyser, Yury Alkhimenkov, Michel Jaboyedoff, and Yury Y. Podladchikov
Geosci. Model Dev., 13, 6265–6284, https://doi.org/10.5194/gmd-13-6265-2020, https://doi.org/10.5194/gmd-13-6265-2020, 2020
Short summary
Short summary
In this work, we present an efficient and fast material point method (MPM) implementation in MATLAB. We first discuss the vectorization strategies to adapt this numerical method to a MATLAB implementation. We report excellent agreement of the solver compared with classical analysis among the MPM community, such as the cantilever beam problem. The solver achieves a performance gain of 28 compared with a classical iterative implementation.
Cited articles
Alkhimenkov, Y. and Podladchikov, Y. Y.: Accelerated pseudo-transient method for elastic, viscoelastic, and coupled hydromechanical problems with applications, Geosci. Model Dev., 18, 563–583, https://doi.org/10.5194/gmd-18-563-2025, 2025. a, b
Alkhimenkov, Y., Räss, L., Khakimova, L., Quintal, B., and Podladchikov, Y.: Resolving wave propagation in anisotropic poroelastic media using graphical processing units (GPUs), Journal of Geophysical Research: Solid Earth, 126, e2020JB021175, https://doi.org/10.1029/2020JB021175, 2021. a
Alkhimenkov, Y., Khakimova, L., Utkin, I., and Podladchikov, Y.: Resolving Strain Localization of Brittle and Ductile Deformation in two-and three-dimensions using Graphical Processing Units (GPUs), arXiv [preprint], https://doi.org/10.48550/arXiv.2305.01701, 2023. a, b, c
Alkhimenkov, Y., Khakimova, L., and Podladchikov, Y.: Shear bands triggered by solitary porosity waves in deforming fluid-saturated porous media, Geophysical Research Letters, 51, e2024GL108789, https://doi.org/10.1029/2024GL108789, 2024a. a
Alkhimenkov, Y., Khakimova, L., Utkin, I., and Podladchikov, Y.: Resolving strain localization in frictional and time-dependent plasticity: Two-and three-dimensional numerical modeling study using graphical processing units (GPUs), Journal of Geophysical Research: Solid Earth, 129, e2023JB028566, https://doi.org/10.1029/2023JB028566, 2024b. a
Alkhimenkov, Y., Khakimova, L., and Podladchikov, Y.: FastLocalization_2D_Stress_Drops, Zenodo [code], https://doi.org/10.5281/zenodo.16412530, 2025. a
Allison, K. L. and Dunham, E. M.: Earthquake cycle simulations with rate-and-state friction and power-law viscoelasticity, Tectonophysics, 733, 232–256, 2018. a
Andrews, D.: Rupture dynamics with energy loss outside the slip zone, Journal of Geophysical Research: Solid Earth, 110, https://doi.org/10.1029/2004JB003191, 2005. a
Buck, W. R.: Effect of lithospheric thickness on the formation of high-and low-angle normal faults, Geology, 21, 933–936, 1993. a
Cundall, P.: Numerical experiments on localization in frictional materials, Ingenieur-archiv, 59, 148–159, 1989. a
Cundall, P.: Numerical modelling of jointed and faulted rock, in: Mechanics of Jointed and Faulted Rock, 1st edn., edited by: Rossmanith, H. P., CRC Press, London, https://doi.org/10.1201/9781003078975, 1990. a
Dal Zilio, L., Lapusta, N., Avouac, J.-P., and Gerya, T.: Subduction earthquake sequences in a non-linear visco-elasto-plastic megathrust, Geophysical Journal International, 229, 1098–1121, https://doi.org/10.1093/gji/ggab521, 2022. a
de Borst, R. and Duretz, T.: On viscoplastic regularisation of strain-softening rocks and soils, International Journal for Numerical and Analytical Methods in Geomechanics, 44, 890–903, 2020. a
De Borst, R., Crisfield, M. A., Remmers, J. J., and Verhoosel, C. V.: Nonlinear finite element analysis of solids and structures, John Wiley & Sons, https://doi.org/10.1002/9781118375938, 2012. a
de Souza Neto, E. A., Peric, D., and Owen, D. R.: Computational methods for plasticity: theory and applications, John Wiley & Sons, https://doi.org/10.1002/9780470694626, 2011. a, b, c
Dieterich, J. H.: Time-dependent friction and the mechanics of stick-slip, Rock friction and earthquake prediction, 790–806, https://doi.org/10.1007/BF00876539, 1978. a
Dieterich, J. H.: Modeling of rock friction: 1. Experimental results and constitutive equations, Journal of Geophysical Research: Solid Earth, 84, 2161–2168, 1979. a
Dormy, E. and Tarantola, A.: Numerical simulation of elastic wave propagation using a finite volume method, Journal of Geophysical Research: Solid Earth, 100, 2123–2133, 1995. a
Drucker, D. C. and Prager, W.: Soil mechanics and plastic analysis or limit design, Quarterly of applied mathematics, 10, 157–165, 1952. a
Duretz, T., de Borst, R., and Yamato, P.: Modeling lithospheric deformation using a compressible visco-elasto-viscoplastic rheology and the effective viscosity approach, Geochemistry, Geophysics, Geosystems, 22, e2021GC009675, https://doi.org/10.1029/2021GC009675, 2021. a
Erickson, B. A., Dunham, E. M., and Khosravifar, A.: A finite difference method for off-fault plasticity throughout the earthquake cycle, Journal of the Mechanics and Physics of Solids, 109, 50–77, 2017. a
Frankel, S. P.: Convergence rates of iterative treatments of partial differential equations, Mathematics of Computation, 4, 65–75, 1950. a
Gabriel, A.-A., Ampuero, J.-P., Dalguer, L., and Mai, P. M.: Source properties of dynamic rupture pulses with off-fault plasticity, Journal of Geophysical Research: Solid Earth, 118, 4117–4126, 2013. a
Heeres, O. M., Suiker, A. S., and de Borst, R.: A comparison between the Perzyna viscoplastic model and the consistency viscoplastic model, European Journal of Mechanics-A/Solids, 21, 1–12, 2002. a
Johnson, P. A., Savage, H., Knuth, M., Gomberg, J., and Marone, C.: Effects of acoustic waves on stick–slip in granular media and implications for earthquakes, Nature, 451, 57–60, 2008. a
Kaneko, Y. and Fialko, Y.: Shallow slip deficit due to large strike-slip earthquakes in dynamic rupture simulations with elasto-plastic off-fault response, Geophysical Journal International, 186, 1389–1403, 2011. a
Lemiale, V., Mühlhaus, H.-B., Moresi, L., and Stafford, J.: Shear banding analysis of plastic models formulated for incompressible viscous flows, Physics of the Earth and Planetary Interiors, 171, 177–186, 2008. a
Ma, S.: A physical model for widespread near-surface and fault zone damage induced by earthquakes, Geochemistry, Geophysics, Geosystems, 9, https://doi.org/10.1029/2008GC002231, 2008. a
Ma, S. and Andrews, D.: Inelastic off-fault response and three-dimensional dynamics of earthquake rupture on a strike-slip fault, Journal of Geophysical Research: Solid Earth, 115, https://doi.org/10.1029/2009JB006382, 2010. a
Moresi, L., Mühlhaus, H.-B., Lemiale, V., and May, D.: Incompressible viscous formulations for deformation and yielding of the lithosphere, Geological Society, London, Special Publications, 282, 457–472, 2007. a
Poliakov, A. and Herrmann, H.: Self-organized criticality of plastic shear bands in rocks, Geophysical Research Letters, 21, 2143–2146, 1994. a
Poliakov, A., Podladchikov, Y., and Talbot, C.: Initiation of salt diapirs with frictional overburdens: numerical experiments, Tectonophysics, 228, 199–210, 1993. a
Pranger, C., Sanan, P., May, D. A., Le Pourhiet, L., and Gabriel, A.-A.: Rate and state friction as a spatially regularized transient viscous flow law, Journal of Geophysical Research: Solid Earth, 127, e2021JB023511, https://doi.org/10.1029/2021JB023511, 2022. a
Preuss, S., Ampuero, J. P., Gerya, T., and van Dinther, Y.: Characteristics of earthquake ruptures and dynamic off-fault deformation on propagating faults, Solid Earth, 11, 1333–1360, https://doi.org/10.5194/se-11-1333-2020, 2020. a
Räss, L., Utkin, I., Duretz, T., Omlin, S., and Podladchikov, Y. Y.: Assessing the robustness and scalability of the accelerated pseudo-transient method, Geosci. Model Dev., 15, 5757–5786, https://doi.org/10.5194/gmd-15-5757-2022, 2022. a
Ruina, A.: Slip instability and state variable friction laws, Journal of Geophysical Research: Solid Earth, 88, 10359–10370, 1983. a
Scuderi, M., Marone, C., Tinti, E., Di Stefano, G., and Collettini, C.: Precursory changes in seismic velocity for the spectrum of earthquake failure modes, Nature geoscience, 9, 695–700, 2016. a
Simpson, G.: Emergence and growth of faults during earthquakes: Insights from a dynamic elasto-plastic continuum model, Tectonophysics, 868, 230089, https://doi.org/10.1016/j.tecto.2023.230089, 2023. a
Templeton, E. L. and Rice, J. R.: Off-fault plasticity and earthquake rupture dynamics: 1. Dry materials or neglect of fluid pressure changes, Journal of Geophysical Research: Solid Earth, 113, https://doi.org/10.1029/2007JB005529, 2008. a, b
Tong, X. and Lavier, L. L.: Simulation of slip transients and earthquakes in finite thickness shear zones with a plastic formulation, Nature communications, 9, 3893, https://doi.org/10.1038/s41467-018-06390-z, 2018. a
Uphoff, C., May, D. A., and Gabriel, A.-A.: A discontinuous Galerkin method for sequences of earthquakes and aseismic slip on multiple faults using unstructured curvilinear grids, Geophysical Journal International, 233, 586–626, 2023. a
Virieux, J.: P-SV wave propagation in heterogeneous media: Velocity-stress finite-difference method, Geophysics, 51, 889–901, 1986. a
Wollherr, S., Gabriel, A.-A., and Uphoff, C.: Off-fault plasticity in three-dimensional dynamic rupture simulations using a modal Discontinuous Galerkin method on unstructured meshes: implementation, verification and application, Geophysical Journal International, 214, 1556–1584, 2018. a
Zienkiewicz, O. C. and Taylor, R. L.: The finite element method for solid and structural mechanics, Elsevier, ISBN 9780080455587, 2005. a
Short summary
This study examines stress drops sequences in elasto-plastic media using 2D simulations, highlighting the importance of high temporal and spatial resolutions in capturing stress evolution and strain fields. Stress drops reflect fault rupture mechanics and emulate earthquake behavior. The non-Gaussian distribution of stress drop amplitudes resembles "solid turbulence." Elasto-plastic models simulate key earthquake processes and could improve seismic hazard assessment.
This study examines stress drops sequences in elasto-plastic media using 2D simulations,...