EGU General Assembly 2022
© Author(s) 2022. This work is distributed under
the Creative Commons Attribution 4.0 License.

Sensitivity analysis using the TREMOL code for seismicity forecasting 

Marisol Monterrubio-Velasco and Natalia Zamora
Marisol Monterrubio-Velasco and Natalia Zamora
  • Barcelona Supercompting Center -Centro Nacional de Supercomputación

Forecasting spatio-temporal occurrence of earthquakes is not a trivial step for the seismic and tsunami hazard assessments. Estimating earthquake rates depends on information of a nonlinear system that is poorly known including the source dimensions. Thus, these assessments rely on e.g. seismic catalogues, or geophysical and geological data that could portray the statistical and physical behaviour of the seismogenic zones. In particular, earthquakes could rupture along asperities or areas of the seismogenic zone with high stress accumulation.Those areas have different physical properties than the surrounding area, such as a high frictional strength and larger stress drop (e.g. Madariaga 1979, Corbi et al., 2017). In this work, we apply the TREMOL code (Monterrubio-Velasco et al., 2019), based on the Fiber Bundle Model, to validate it as a tool to reproduce the seismicity occurring by the rupture of large, in some cases, single asperities. We have selected four regions where large earthquakes have occurred: M8.8 Maule 2010 (Chile) earthquakes, M9.1 Tohoku 2011 (Japan), M7.6 Nicoya 2012 (Costa Rica) and M8.3 Coquimbo 2015 (Chile). In these tectonic regions, earthquake sequences are generated based on a discrete model of material failure used in TREMOL. One of the most notable results is that the maximum earthquakes of the real sequences are achieved. Also, in most cases, the magnitude - frequency distribution is similar to those of real data. While the outcomes of TREMOL are given in rupture areas, several area-magnitude scaling laws are employed to obtain moment magnitudes. By carrying out a sensitivity analysis of different scaling laws, we show the bias in the synthetic catalogues which is a critical input in seismic hazard assessment. It is shown that the synthetic seismicity using the Ramirez-Gaytan scaling law (Ramirez-Gaytan et al. 2014) is the best to fit the magnitude of the real series in most of the cases. Following the validation of TREMOL, we provide a new seismic scenario generator of future events to assist e.g. the Probabilistic Seismic/Tsunami Hazard Assessment (PSHA/PTHA) complementing the seismic forecast with other well known statistical tools. 




Corbi, F., Funiciello, F., Brizzi, S., Lallemand, S., and Rosenau, M. (2017). Control of asperities size and spacing on seismic behavior of subduction megathrusts, Geophys. Res. Lett., 44, 8227– 8235, doi:10.1002/2017GL074182.

Madariaga, R. (1979). On the relation between seismic moment and stress drop in the presence of stress and strength heterogeneity, J. Geophys. Res.-Sol. Ea., 84, 2243–2250.

Monterrubio-Velasco et al., (2019). A stochastic rupture earthquake code based on the fiber bundle model (TREMOL v0.1): application to Mexican subduction earthquakes. Geosci. Model Dev., 12, 1809–1831.

Ramírez-Gaytán, A., Aguirre, J., Jaimes, M. A., and Huérfano, V. (2014). Scaling relationships of source parameters of M w 6.9–8.1 earthquakes in the Cocos–Rivera–North American subduction zone, Bulletin of the Seismological Society of America, 104, 840–854.


How to cite: Monterrubio-Velasco, M. and Zamora, N.: Sensitivity analysis using the TREMOL code for seismicity forecasting , EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-7735,, 2022.