the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
| 22 Dec 2025
Consistency of seismic hazard estimates from a physics-based earthquake simulator: a case study in south-eastern Spain
Abstract. Estimating seismic hazard is crucial for enhancing societal resilience and risk mitigation strategies. Probabilistic Seismic Hazard Analysis (PSHA) is the current standard framework, traditionally relying on empirical earthquake rupture forecasts (ERFs) and ground-motion models. In this framework, physics-based earthquake cycle simulators are emerging as powerful tools in PSHA, capable of replicating observed seismicity and seismic hazard statistics. Here, we present a quantitative consistency evaluation of physics-based PSHA models at the Eastern Betic Shear Zone in south-eastern Spain against both historical macroseismic intensity data and instrumental ground-shaking records. We use synthetic catalogues from RSQSim earthquake-cycle simulations to construct two physics-based ERFs that we pipeline into a PSHA calculation. Results indicate that the physics-based ERFs derived from the best-performing simulation model, previously ranked against empirical benchmarks, achieve the best overall agreement with observed macroseismic intensities and acceleration records at 10 sites, outperforming both the lower-performing simulation and a traditional area-source model. Our findings highlight that the incorporation of physics-based models into a PSHA framework is reliable, enabling the inclusion of rupture behaviour and fault-system interactions, key challenges in PSHA. We also advocate for the complementarity of physics-based models with traditional approaches in PSHA to better capture epistemic uncertainties in the hazard representation.
- Preprint
(2460 KB) - Metadata XML
-
Supplement
(1720 KB) - BibTeX
- EndNote
Status: open (until 07 Feb 2026)
- RC1: 'Comment on egusphere-2025-5485', Anonymous Referee #1, 07 Jan 2026 reply
-
RC2: 'Comment on egusphere-2025-5485', Anonymous Referee #2, 09 Jan 2026
reply
Overall, the manuscript is well written and technically sound, addressing an important and timely topic on the integration and evaluation of physics-based earthquake simulators within a PSHA framework. The methodology is rigorous, the workflow is clearly designed, and the consistency tests against both macroseismic and instrumental observations are carefully implemented, making the study suitable for publication after revision. That said, the introduction would benefit from a clearer hierarchical structure and stronger logical progression, particularly by more explicitly distinguishing the limitations of traditional PSHA, the advantages of physics-based approaches, and the specific scientific gap this work fills, as well as by summarizing the main objectives and novel contributions more concisely toward the end of the section. In addition, while the discussion is generally solid, its depth and breadth could be further enhanced by providing a more integrative synthesis of the results, clarifying the broader implications for seismic hazard practice in low-to-moderate strain regions, and more explicitly discussing the transferability of the proposed workflow to other tectonic settings and its role within future hybrid or logic-tree PSHA frameworks. Addressing these points would improve the clarity, balance, and overall impact of the manuscript without requiring major changes to the core analysis or results.
Citation: https://doi.org/10.5194/egusphere-2025-5485-RC2
Data sets
Earthquake cycle seismic hazard models for the Eastern Betic Shear Zone Octavi Gómez-Novell et al. https://doi.org/10.5281/zenodo.17482084
Model code and software
Simulator to Openquake Octavi Gómez-Novell https://github.com/octavigomez/simulator_to_openquake
Viewed
| HTML | XML | Total | Supplement | BibTeX | EndNote | |
|---|---|---|---|---|---|---|
| 148 | 51 | 8 | 207 | 27 | 17 | 25 |
- HTML: 148
- PDF: 51
- XML: 8
- Total: 207
- Supplement: 27
- BibTeX: 17
- EndNote: 25
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1
This manuscript addresses an important topic in seismic hazard assessment: the integration of physics-based earthquake simulators into probabilistic seismic hazard analysis (PSHA) frameworks. The work demonstrates that properly calibrated physics-based models can provide consistent hazard estimates in low-to-moderate strain regions, which represents a valuable contribution to the field. However, some issues require clarification.
MAJOR COMMENTS
Lines 220-235: The paper presents a consistency testing framework, but there are opportunities to employ more rigorous statistical validation methods. In statistical practice, methods exist to estimate out-of-sample performance and predictive power of models, such as cross-validation. While applying such methods to an entire PSHA model may be technically challenging, components of the PSHA (particularly the occurrence models) can be tested using these approaches. Since this paper focuses on physics-based synthetic catalogs, testing the catalog or the occurrence models based on the synthetic catalogs before integrating them into the full PSHA would strengthen the methodology significantly. If this is out of the scope of your work, or if that would lead to a lengthy article, please discuss this briefly and state these issues should be tested in a companion report or paper. If such tests have been conducted in other studies, they should be cited here. Additionally, the adopted testing procedure evaluates the joint model (occurrence rates combined with GMPE), which should be explicitly stated and briefly discussed.
Lines 250-254: Some clarifications are needed to rule out methodological inconsistency. The synthetic catalogs are tested against the stationary Poisson process hypothesis, which is typically used in PSHA for declustered events (mainshocks only). The manuscript suggests that there may be some foreshocks and aftershocks in the synthetic catalogs, which do not significantly affect the PSHA results, which is acceptable. However, readers may be confused and think the paper performs PSHA for full non-declustered seismicity. This distinction must be clarified explicitly. Please discuss briefly whether the presence of clustering in the catalogs affects the interpretation of results and whether a sensitivity analysis excluding clustered events would alter conclusions.
Line 208: The authors use "Classical PSHA" in OpenQuake, but the OpenQuake Engine also allows event-based PSHA using synthetic catalogs directly. This approach would avoid fitting occurrence laws, a procedure that potentially loses information contained in the catalogs. Moreover, OpenQuake allows scenario calculations that could be run for every event in the synthetic catalogs to calculate exceedances directly. These alternatives should be discussed, including why the classical approach was chosen over these potentially more direct methods. This discussion would clarify whether methodological choices limit the ability to fully leverage the physics-based catalog information.
Line 137: While the authors state that Cat-21 and Cat-18 represent the best and worst performing catalogs based on benchmarking, the rationale for explicitly including the worst-performing catalogue requires explanation. What insights does Cat-18 provide that justify its inclusion in the hazard analysis? Is it meant to demonstrate the importance of proper model selection, or does it represent a plausible epistemic uncertainty branch? This clarification would help readers understand whether poor-performing models should ever be included in operational hazard assessments or logic trees.
Line 23 (Abstract) and throughout: The manuscript advocates for combining physics-based and traditional approaches, describing this as "complementarity." This concept is well-established across many scientific fields as hybrid modeling, which systematically combines physics-based models with data-driven approaches. The authors should acknowledge this broader context and cite relevant literature on hybrid models in seismic hazard or related fields (e.g., hydrology, climate science). This would strengthen the theoretical foundation and help position the work within established methodological frameworks.
MINOR COMMENTS
ABSTRACT
Line 21: The phrase "both the lower-performing simulation and" could be omitted for brevity without loss of meaning.
Line 22: The term "reliable" should be verified to ensure it is fully supported by the results presented. Given the consistency testing (rather than validation) performed for some components, this wording may overstate the conclusions.
1. INTRODUCTION
Line 28: The statement that PSHA was "formalized by Cornell (1968)" is not 100% fair. Please see McGuire (2008), Probabilistic seismic hazard analysis: Early history, Earthquake Engng Struct. Dyn. 2008; 37:329–338. DOI: 10.1002/eqe.765
Line 86: The reference to Ellingwood and Wen (2005) does not support the statement about high-impact, low-probability events as written. This citation should be removed or replaced with more appropriate references.
End of Introduction: A paragraph should be added that clearly states the objectives of this paper and provides a roadmap of the manuscript structure. This would help readers understand the overall contribution and organization.
2. DATA AND METHODS
Line 175: Please clarify whether each rupture in the catalogs has only a single occurrence, or whether repeated similar ruptures can occur.
Line 262: The statement that macroseismic records at close distances are mostly related to faults requires explanation.
Line 290: The statement that PGV "is the most likely linked to damage" is too strong and not generally true. Please see the literature on seismic fragility curves, e.g.
Luco N., Cornell C.A. (2007) Structure-Specific Scalar Intensity Measures for Near-Source and Ordinary Earthquake Ground Motions, Earthquake Spectra, Volume 23, No. 2, pages 357–392. https://doi.org/10.1193/1.2723158
Figure 5c: The jump in the curve at the last point for intensity MI=12 requires explanation. Is this a computational artifact, a feature of the GMICE, or physically meaningful?
Line 324: Please provide the justification for taking the average p-value across the four cases and then computing its logarithm.
3. RESULTS
Line 340: The choice of 2% probability of exceedance in 50 years should be justified. While this is a valid choice, the most common selection for design purposes is 10% in 50 years. Was this choice made for specific reasons related to the EBSZ, or to match existing hazard maps? This should be stated explicitly.
Figure 6: Please provide specific commentary on the differences between the hazard maps in the area of the city of Vera.
Figure 7: The current presentation makes comparison difficult. It would be more effective to have one subplot per city/station showing all three hazard curves (Cat-21, Cat-18, and area source) overlaid. This would facilitate direct comparison of model performance.
Line 455: The summing of LogP values to rank models requires justification.
EDITORIAL COMMENTS
Line 234: The term "consistency check" is introduced but not formally defined until later in the text. Provide a brief definition at first use.
Line 590: "areas model" should be "area source model" throughout.