Three-dimensional thermomechanical modelling beneath Guerrero, Mexico, and its relation to the occurrence of interplate seismic phenomena

Nobuaki Suenaga1, Shoichi Yoshioka1,2, Vlad Constantin Manea3, Marina Manea3, Erika Jessenia Moreno1 & Yingfeng Ji4,5

1Research Center for Urban Safety and Security, Kobe University, Rokkodai-Cho 1-1, Nada Ward, Kobe, 657-8501, Japan. 2Department of Planetology, Graduate School of Science, Kobe University, Rokkodai-Cho 1-1, Nada Ward, Kobe, 657-8501, Japan. 3Computational Geodynamics Laboratory, Instituto de Geociencias, Universidad Nacional Autónoma de México, Campus Juriquilla, Querétaro, 76230, México. 4Institute of Tibetan Plateau Research, Chinese Academy of Sciences, Lincui Road 16-3, Chaoyang, Beijing, 100101, China. 5University of Chinese Academy of Sciences, Lincui Road 16-3, Chaoyang, Beijing, 100101, China.

Abstract

In this study, we investigated the relationships among the occurrence regions of megathrust earthquakes and slow earthquakes, estimated interplate temperatures, and dehydration processes at the plate interface in the subduction zone of the Cocos plate in Mexico, with a focus on the Guerrero seismic gap. For this purpose, we performed a series of 3-D thermomechanical numerical simulations of subduction based on a combination of multiple unknown model parameters. Comparing the temperature field for the optimal model with the occurrence region of the interplate seismic phenomena, the temperature along the slab surface where megathrust earthquakes tend to occur is estimated to range from 200 to 400 °C. In addition, the occurrence areas of recent long-term slow slip events (hereafter referred to as L-SSEs) with a recurrence interval of ~ 4 years in the range where the amount of slip is 10 cm or more correspond to a temperature range of 350 to 550 °C outside the Guerrero seismic gap region. On the other hand, in the Guerrero seismic gap region, L-SSE slip is observed up to shallower areas, but the updip limit temperature in the region with the largest slip amount (~ 14 cm) is also estimated to be approximately 350 °C. Furthermore, for regions affected by tectonic tremors (hereafter referred to as TTs) within seaside swarms and inland swarms, our temperature estimates range from 500 to 570 °C and from 600 to 700 °C, respectively. Finally, we compared the estimated dehydration processes on the plate interface for the optimal model with interseismic events. These findings suggest that considerable dehydration from MORB within the oceanic crust and ultramafic rock within the slab mantle occurred near long-term slow slip events and tectonic tremor regions.

This is an open-access article provided by Nature, Scientific Reports.

Cite this article

Suenaga, N., Yoshioka, S., Manea, V.C. et al. Three-dimensional thermomechanical modelling beneath Guerrero, Mexico, and its relation to the occurrence of interplate seismic phenomena. Sci Rep 15, 7513 (2025). https://doi.org/10.1038/s41598-025-91257-9