## 5 Conclusions  
A holistic seismic hazard assessment, including both probabilistic and scenario-based methods, was conducted for the highly active Western Gulf of Corinth area in Central Greece. Initially, seismic hazard was evaluated using a probabilistic approach, taking into account all the M\({}_{\text{w}}\geq 4.0\) earthquakes recorded during 1900-2019, in order to determine values for PGA, PGV, and S\({}_{\text{a}}\) in a dense grid. For this purpose, the Cornell-McGuire method was utilized. Aiming to incorporate a multitude of seismotectonic models in order to qualitatively cope with epistemic uncertainties, three models--namely, ESHM13, ESHM20, and VAM16--were considered, and GMPEs were introduced for each zone. The PGA and PGV results were obtained through an equal-weight logic tree approach in which each major branch was a seismotectonic model and each minor branch was a modified GMPE. The logic tree procedure incorporated with the five GMPEs that were developed using Greek data qualitatively reduced epistemic uncertainties, which simple PSHA approaches may carry. In addition, PGA values were computed for the most populated urban areas (Patras, Nafpaktos, and Aigion) for various exceedance probabilities in 50 years. The results obtained herein highlight the great significance of seismic hazard assessment using a combined approach that takes into account not only the evaluation of ground parameters in terms of probabilities of occurrence in a given time frame but also the anticipated effects of deterministic worst-case scenarios. Future work could include disaggregation of PSHA results so that the parameter pair of magnitude and epicentral distance that contributes most to seismic hazard can be identified. In this way, scenarios could be generated based on this result. Furthermore, a study of the Peak Ground Rotational Acceleration (PGRA) and Velocity (PGRV) values could be undertaken, as it has been proven that their results provide aid for engineering purposes [58; 77; 120]. In addition, the incorporation of a comprehensive microzonation scheme could provide valuable insight into the impact of future strong earthquakes on the major cities of the study area by identifying the possible amplification trends with respect to the structural response. Lastly, the results of the present seismic hazard study can be exploited in the future to assess seismic risk at urban centers in the WGoC area after incorporating structural vulnerability data. The following supporting information can be downloaded at [[https://www.mdpi.com/article/10.3390/app12211152/s1](https://www.mdpi.com/article/10.3390/app12211152/s1)]([https://www.mdpi.com/article/10.3390/app12211152/s1](https://www.mdpi.com/article/10.3390/app12211152/s1)), Figure S1. The three zonation models that were implemented for the WGoC, (a) the ESHM13 model [62; 63], (b) the ESHM20 model [64] and (c) the model VAM16 from [PERSON] et al. [65] proposed for the Greek territory; Table S1: The seismicity and geometry data that was computed for each Area Source Zone (ASZ). Conceptualization, G.K.; Data curation, G.K. and A.Z. ; Formal analysis, G.K. and A.Z. ; Funding acquisition, G.K.; Investigation, G.K., A.Z., P.B. and V.K. ; Methodology, G.K.; Resources, N.V.; Software, A.Z. and P.B. ; Supervision, G.K. and N.V.; Validation, G.K. and A.Z. ; Visualization, P.B. and V.K. ; Writing--original draft, G.K., A.Z., P.B. and V.K. ; Writing--review & editing, G.K., A.Z., V.K. and N.V. All authors have read and agreed to the published version of the manuscript. This research received no external funding. Not applicable. We express our gratitude to the three anonymous reviewers for their constructive comments that helped improve the manuscript, as well as to the editor and the guest editors for providing us with the opportunity to publish this study.