Adapting CyberShake for Europe using OpenQuake-Derived Earthquake Rupture Forecasts

Over the past two decades, seismic hazard modeling has advanced along two complementary frontiers: empirical probabilistic frameworks, which systematically capture uncertainty through statistical inference, and physics-based simulation platforms, which directly compute ground motions from the governing equations of wave propagation. This project seeks to unify these two worlds by developing an end-to-end integration between OpenQuake and CyberShake, thereby creating a new generation of seismic hazard models that are globally extensible, probabilistically complete, and physically consistent. CyberShake has been under active development for more than a decade, demonstrating its robustness and scientific maturity through extensive implementations in California. It performs a physics-based probabilistic seismic hazard analysis (PSHA), replacing traditional empirical Ground Motion Prediction Equations (GMPEs) with full 3D numerical simulations of seismic wave propagation. Built upon the UCERF2/3 Earthquake Rupture Forecasts, CyberShake computes hazard curves directly from synthetic seismograms generated via Strain Green’s Tensors and thousands of stochastic rupture variations. This approach enables non-ergodic, site-specific hazard estimation and has set a global benchmark for high-fidelity hazard computation. However, its application has remained geographically limited: both the ERF and 3D velocity models were designed specifically for California, requiring extensive datasets that are rarely available elsewhere. Conversely, OpenQuake, developed by the Global Earthquake Model (GEM) Foundation, provides a fully open-source, Python-based framework for probabilistic seismic hazard and risk analysis. It serves as the computational backbone of large-scale hazard models such as the European Seismic Hazard Model 2020 (ESHM20), which integrates decades of regional expertise into a unified and statistical representation. OpenQuake provides a complete probabilistic framework to build Earthquake Rupture Forecasts (ERFs) that combine declustered catalogs, background seismicity, and multi-branch logic trees, ensuring a balanced and uncertainty-aware representation of regional tectonics. Furthermore, its ecosystem extends seamlessly to vulnerability and exposure modules, enabling the translation of hazard into actionable risk assessments and resilience planning.

This project will establish a direct pipeline from OpenQuake’s event-based results to the generation of an ERF compatible with CyberShake’s simulation framework, ensuring moment–rate consistency. By doing so, it will enable CyberShake simulations to be performed for regions beyond California, extending its use to Europe based on the knowledge contained in the ESHM20. The first pilot region is Istanbul, Turkey, a densely populated metropolis located near the western termination of the North Anatolian Fault. Our initial results show that the workflow is already functioning at the prototype level: we have developed a unified 3D velocity model for the Istanbul region by combining available tomographic models with local datasets; generated preliminary event-based rupture catalogs from ESHM20 using OpenQuake; and demonstrated early convergence behavior in hazard curves, indicating that the rupture sampling strategy is statistically robust. These initial results demonstrate the feasibility of the integration approach and indicate that the essential elements needed for a CyberShake-ready ERF are already in place.