The present paper deals with an inhabited, creeping mountainous landmass with profound surface deformation that affects the local community. The scope of the paper is to gather surficial and subsurface information in order to understand the parameters of this creeping mass, which is usually affected by several parameters, such as its geometry, subsurface water, and shear zone. Therefore, a combined aerial and surface investigation has been conducted. The aerial investigation involves UAV’s LiDAR acquisition for the terrain model and a comparison of historical aerial photographs for land use changes. The multi-technique surface investigation included resistivity (ERT) and seismic (SRT, MASW) measurements and density determination of geological formations. This combination of methods proved to be fruitful since several aspects of the landslide were clarified, such as water flow paths, the internal geological structure of the creeping mass, and its geometrical extent. The depth of the shear zone of the creeping mass is delineated at the first five to ten meters from the surface, especially from the difference in diachronic resistivity change.

Boeotia, located in Central Greece, experiences frequent seismic activity, mainly due to its proximity to the Gulf of Corinth. Significant earthquakes have occurred in the broader study area, such as the ones of Atalanti (Mw = 6.8, 6.9) in 1894, as well as the Alkyonides sequence in 1981 that included three Mw > 6.0 events. In late 2020, a Mw = 4.6 mainshock took place near Thiva, a populated town in Boeotia, followed by the 2021–2022 seismic sequence with three Mw > 4.0 earthquakes. The objective of this study is to perform a Probabilistic Seismic Hazard Assessment (PSHA) for Boeotia through the computation of the Peak Ground Acceleration (PGA) and Peak Ground Velocity (PGV) using two truncation levels (ε = 0 and 3). Moreover, Uniform Hazard Spectra (UHS) are constructed in terms of Spectral acceleration (Sa) for Thiva and Livadia, the capital of Boeotia. To achieve this, three seismotectonic models in the form of area sources are employed in the computational framework. Ground Motion Prediction Equations (GMPEs), using data of the area of Greece, are utilized to estimate PGA and PGV. For each area source, the percentages of normal and non-normal (reverse or strike-slip) fault plane solutions are computed in order to generate minor branches for each GMPE that takes into account the focal mechanism type. This approach introduces variability and reduces uncertainties in PSHA. Additionally, a sensitivity analysis was performed by keeping constant one logic tree, first the GMPE, then the source-model tree, while varying the other, to assess the consistency of individual GMPEs and source models. The findings reveal that western and eastern Boeotia have higher seismic hazard, attributed to the seismotectonics of the study area. Additionally, the seismic hazard level in Thiva is higher compared to Livadia.

The main purpose of this study is the subsurface investigation of two complex geological environments focusing on the improvement of data acquisition and processing parameters regarding electric and seismic tomographic techniques. Two different study areas, in central–east Peloponnese and SE Attica, were selected, where detailed geological mapping and surface geophysical survey were carried out. The applied geophysical survey included the application of electrical resistivity tomography (ERT), seismic refraction tomography (SRT) and ground penetrating radar (GPR). The geoelectrical measurements were acquired with different arrays and electrode configurations. Moreover, various types of seismic sources were used at seventeen shot locations along the seismic arrays. For the processing of geoelectrical data, clustered datasets were created, increasing the depth of investigation and discriminatory capability. The seismic data processing included the following: (a) the creation of synthetic models and seismic records to determine the effectiveness and capabilities of the technique, (b) spectral analysis of the seismic records to determine the optimal seismic source type and (c) inversion of the field data to create representative subsurface velocity models. The results of the two techniques successfully delineated the complex subsurface structure that characterizes these two geological environments. The application of the ERT combined with the SRT are the two dominant, high-resolution techniques for the elucidation of complex subsurface structures.

The Central Ionian Islands of Cephalonia and Ithaca belong to the most seismically active Greek region, mainly due to the presence of the dextral Cephalonia-Lefkada Transform Fault Zone. The study area has experienced strong earthquakes in the 20 th century, including the destructive 1953 sequence with maximum intensity 9.0. The Paliki peninsula, western Cephalonia, hosted two strong earthquakes (M w = 6.1 and 5.8) in 2014, with ground acceleration reaching ~560 cm/s 2 and 735 cm/s 2 , respectively. This study updates the seismic hazard evaluation in Cephalonia and Ithaca using new data and computational techniques to reduce epistemic uncertainties. The probabilistic approach of Cornell and McGuire was used, and the uncertainties are reduced through data variability of the source models, seismicity data, and Ground Motion Prediction Equations using a logic tree approach, sampled by implementing the Latin Hypercube Sampling method. The spatial distribution of Peak Ground Acceleration and Peak Ground Velocity for return periods of 475 and 950 years indicates low variation in the entire study area and that the Paliki peninsula possesses the highest level of seismic hazard. Additionally, site-specific analysis across the three main towns, Lixouri and Argostoli in Cephalonia and Vathi in Ithaca, reveals that Lixouri has the greatest level of seismic hazard, while Vathi the lowest.

The 2020-2021 seismic sequence at the Western Gulf of Corinth, Central Greece, has been thoroughly analyzed using seismological and geodetic data processing. We present a high-resolution dataset of over 4000 relocated earthquakes between June 2020 and February 2021, delineating the activated structures. We examine the evolution of various clusters that were triggered during a three-stage sequence which began on December 23, 2020, with an Mw = 4.6 event near Marathias, migrated eastwards, following an Mw = 5.0 event on January 12, 2021 near Trizonia Island, and culminated on February 17 with an Mw = 5.3 offshore event north of Psatho-pyrgos. Focal mechanisms for the 20 stronger events, determined by moment tensor inversion, and 36 weaker ones, obtained through first motion polarities measurements, revealed dominant normal faulting. The fault plane of the February 17, 2021 major event is modeled through the inversion of geodetic data. Results suggest a very shallow geodetic centroid at 1.5 km depth, consistent with the seismic centroid (3.5 km), but different from the relocated hypocenter (7.5 km). Spatiotemporal analysis reveals seismic migration, following a diffusion law with D values in the range 0.06-0.25 m2s  1, consistent with triggering due to pore-pressure diffusion by fluids intrusion. This enables failure on non-optimally oriented faults, which explains the significant proportion of the observed strike-slip and oblique-normal faulting. Coulomb stress transfer due to the major events shows only weak stress-loading on the fault of the January 12, 2021 earthquake and a higher level of positive stress transfer to the south-dipping than to the north-dipping nodal plane of the 17 February event. The latter could have facilitated the nucleation of the rupture at 7.5 km and its further propagation to shallower depths along a south-dipping plane, where the bulk of its seismic and geodetic moment was released.

Summary The present study aims to present the results of a near-surface geophysical investigation carried out at a construction site in the area of Goudi (Athens) along with its contribution to the determination of the geotechnical characteristics. For this purpose, the Ground Penetrating Radar (GPR), Electrical Resistivity Tomography (ERT), Seismic Refraction Tomography (SRT) and Multichannel Analysis of Surface Waves (MASW) techniques were implemented in the area. The cores of three geotechnical boreholes existing in the area have also been taken into account. The application of the GPR and ERT techniques did not reached a satisfactory investigation depth, however the existence of a possible mechanically degraded zone have been adumbrated at the south part of the excavation area. Through the application of the SRT and MASW techniques, the seismic waves velocities of the investigated lithological formations were calculated and the subsurface structure of the study area was outlined. Additionally, lateral variations in the P-wave seismic velocities, especially in the first investigated seismic layer, provided further evidence for the existence of mechanically unstable zones, which necessitate a more detailed investigation. Finally, from the laboratory determination of the formations’ densities and their seismic wave velocities, their elastic moduli and geotechnical parameters have been calculated.

The Athenet network is the network of the Seismological Laboratory of the National and Kapodistrian University of Athens. We present the geophysical investigation that has been carried out at six seismological stations of the Athenet network for their site characterization. More specifically, at the location of each seismological station, four geophysical methods have been carried out: Seismic Refraction Tomography (SRT), Multichannel Analysis of Surface Waves (MASW), the Horizontal to Vertical Spectral Ratio (HVSR) technique, and Electrical Resistivity Tomography (ERT). The applied geophysical survey provided important information regarding the site characterization at the selected seismological stations, including key parameters such as the fundamental frequency fo, the shear-wave velocity VS, the average shear-wave velocity for the upper 30 m depth (VS30), the seismic bedrock depth, the soil type, and the subsurface geology. Moreover, selected elastic moduli (Poisson’s ratio, shear, bulk, and Young moduli) have been calculated. The site characterization information contributes to the determination of the amplification factors for each site that can lead to more accurate calculation of Peak Ground Acceleration (PGA) or Peak Ground Velocity (PGV) and, therefore, trustworthy Probabilistic and Stochastic Seismic Hazard Assessments. The derived fundamental frequency for the seismological stations of VILL, LOUT, THAL, and EPID have been determined to be equal to 10.4, 2.7, 1.4, and 7.1 Hz and their amplification factors to be 1.9, 3.1, 1.7, and 2.6, respectively. For stations MDRA and ATAL, these parameters could not be determined.

Probabilistic Seismic Hazard Assessment (PSHA) was carried out for the administrative region of Attica (Greece). Peak Ground Acceleration (PGA) and Peak Ground Velocity (PGV) values were calculated for return periods of 475 and 950 years for five sub-areas covering the entire region. PGA hazard curves and Uniform Hazard Spectra (UHS) in terms of spectral acceleration (S-a) values were generated for Athens, Methana, and the capitals of each island of Attica (Salamina, Aegina, Poros, Hydra, Spetses, Kythira, and Antikythira). Area sources were adopted from the Euro-Mediterranean Seismic Hazard Model 2013 (ESHM13) and its update, ESHM20, taking into account both crustal and slab tectonic environments. Ground Motion Prediction Equations (GMPEs) proposed for the Greek territory were ranked for PGA and PGV. Each GMPE was reconstructed as a weighted model, accounting for normal and non-normal focal mechanisms for each area source. PGA, PGV, and S-a values were computed using a logic tree, integrating the seismotectonic models as major branches and sub-logic trees, comprised of multiple ranked GMPEs for each area source, as minor branches. The results showed higher seismic hazard values in sub-areas near the Gulf of Corinth and the slab interface, which could indicate a need to revise the active building code in Attica.

Summary The cave of Alistrati, is located in the Prefecture of Serres, Northern Greece near the foothills of Mount Menoikio, in the area of Petroto. This area is structured by crystalline limestones, where the development of a complex and multilevel karst system is favored. An extensive geomorphological survey was carried out for the accurate mapping of the karst surface above the cave, using UAS. For the investigation of a possible lateral extension of the existing karstic conduit, a detailed surface geophysical investigation was carried out. More specifically, three geophysical techniques were implemented: a) the Electrical Resistivity Tomography (ERT), b) the Ground Penetrating Radar (GPR) technique and c) the Very-Low Frequency (VLF) method. These 13 lines of the three geophysical techniques are fully matched at 4 locations and were join-interpreted, yielding remarkable findings. The comparative results of the above geophysical techniques, as well as their 3D presentation, highlight similar geophysical anomalies, evaluated as different types of karst system structures. Therefore, the combined geophysical survey has indicated the existence and interconnection of the first two karst levels of the area, up to a depth of 50m, as well as the possible extension of the Alistrati karstic conduit to the northeast.

We investigate an earthquake sequence involving an Mw = 4.6 mainshock on 2 December 2020, followed by a seismic swarm in July–October 2021 near Thiva, Central Greece, to identify the activated structures and understand its triggering mechanisms. For this purpose, we employ double-difference relocation to construct a high-resolution earthquake catalogue and examine in detail the distribution of hypocenters and the spatiotemporal evolution of the sequence. Furthermore, we apply instrumental and imaging geodesy to map the local deformation and identify long-term trends or anomalies that could have contributed to stress loading. The 2021 seismic swarm was hosted on a system of conjugate normal faults, including the eastward extension of the Yliki fault, with the main activated structures trending WNW–ESE and dipping south. No pre- or coseismic deformation could be associated with the 2021 swarm, while Coulomb stress transfer due to the Mw = 4.6 mainshock of December 2020 was found to be insufficient to trigger its nucleation. However, the evolution of the swarm is related to stress triggering by its major events and facilitated by pore-fluid pressure diffusion. The re-evaluated seismic history of the area reveals its potential to generate destructive Mw = 6.0 earthquakes; therefore, the continued monitoring of its microseismicity is considered important.

The Gulf of Corinth (Central Greece) is one of the most rapidly extending rifts worldwide, with its western part being the most seismically active, hosting numerous strong (M ≥ 6.0) earthquakes that have caused significant damage. The main objective of this study was the evaluation of seismic hazard through a probabilistic and stochastic methodology. The implementation of three seismotectonic models in the form of area source zones via a logic tree framework revealed the expected level of peak ground acceleration and velocity for return periods of 475 and 950 years. Moreover, PGA values were obtained through the stochastic simulation of strong ground motion by adopting worst-case seismic scenarios of potential earthquake occurrences for known active faults in the area. Site-specific analysis of the most populated urban areas (Patras, Aigion, Nafpaktos) was performed by constructing uniform hazard spectra in terms of spectral acceleration. The relative contribution of each selected fault segment to the seismic hazard characterizing each site was evaluated through response spectra obtained for the adopted scenarios. Almost all parts of the study area were found to exceed the reference value proposed by the current Greek National Building Code; however, the three urban areas are covered by the Eurocode 8 regulations.

In early March 2021, three shallow earthquakes, two mainshocks with M6.3 and M6.0 and one major aftershock with M5.6 impacted both the mountainous Damasi-Tyrnavos region (northern Thessaly, Greece) and the adjacent Plio-Quaternary basin. Each major event was followed by rich aftershock activity recorded by local and regional seismographs and accelerographs. Herein, we present a comprehensive analysis of the seismic sequence, from its foreshock activity starting on 28 February, 2021 and for a period of two months using new high-resolution catalogues of relocated earthquakes and hundreds of focal mechanisms. The results indicate that the aftershocks form a zone that spans ~50 km NW-SE, while focal depths range between 5 and 15 km. More than 400 focal mechanisms, computed for events with M >= 2.5, mainly exhibit normal faulting in a NW-SE direction, while WNW-ESE to E-W normal faulting is also evidenced, in particular after the occurrence of the last major event on 12 March. The stress-field was reconstructed on a local and broader scale by inverting focal mechanism data, revealing a rotation of the sigma 3 axis trend from NNE-SSW, in the Damasi-broader region, to NW-SE northwards, to the region of Kozani-Grevena that hosted an Mw = 6.6 shallow mainshock in 1995. Subcrustal seismicity, present beneath those areas, implies that large-scale tectonics and plate dynamics are likely involved in the deformation of the upper crust. Coulomb stress transfer after the 3 major events of the 2021 Damasi-Tyrnavos sequence reveals that stress-loaded areas include those where most aftershocks were triggered. The analysis provides implications to the seismic hazard of the activated area, as a major NW-SE active normal fault close to Larissa city became stress-loaded, constituting a possible candidate source for significant future earthquakes.

Lefkada and Corfu old towns are located in the western part of Greece, in the Ionian Sea. Their proximity to the Hellenic subduction zone (HSZ) is the reason for their intense seismicity. The main goal of this study was the estimation of the geotechnical characteristics of the subsurface, with the contribution of applied geophysical techniques. Therefore, seismic refraction tomography (SRT) and multichannel analysis of surface waves (MASW) were applied. A total of thirty-three (33) seismic and geoelectrical profiles were performed in both towns in order to evaluate the geotechnical characteristics of the subsurface formations. Additionally, subsurface resistivity distributions were investigated with the application of electrical resistivity tomography (ERT). Some important elastic moduli were calculated through the combination of estimated seismic wave velocities and laboratory density measurements. The horizontal distribution of seismic velocities and mechanical properties (σ, E, K, G) of Corfu town was illustrated in maps, for the depth of 5 m. The geophysical interpretation also revealed that Lefkada’s subsurface consists of only one compact geological formation, with little or no variation of its geophysical-geotechnical characteristics. Beyond that, the ground type classifications for the two towns were determined according to the European Committee for Standardization Eurocode 8, based on VS30 values.

Moderate‐to‐large earthquakes in rifts may occur on leading boundary faults or inner antithetic faults. Here we show a rare case of the 2020–2021 seismic sequence in the Corinth rift, that culminated in the shallow rupture of the antithetic fault, neither preceded nor followed by the leading fault rupture. The hypocenter of the largest shock (Mw 5.3 of 17 February 2021) was located at ∼8 km depth. However, seismic waveform data, supported by satellite‐geodetic and tide gauge measurements, pointed to rupture at shallow depth (∼3 km), where no earthquakes were previously observed. We show that the earthquake most probably ruptured two orthogonal, conjugate fault segments: a weak nucleation phase occurred in the microseismically highly active sub‐horizontal detachment layer, followed – a few seconds later – by a larger, shallow moment release on a high‐angle, south‐dipping normal fault. The latter is the Mornos offshore fault, antithetic to the leading, north‐dipping Psathopyrgos fault. Our study presents the first instrumental/observational evidence of a very shallow Mw 5+ event in this rift – and one of the few reported worldwide. The depth limit of the main shallow slip patch coincides with the expected crossing of the Mornos fault with the Psathopyrgos fault, stressing the importance of fault segmentation and rooting inherited from the rift history. This unusual shallow slip in a depth range with little background seismicity and few aftershocks needs to be further investigated by dynamic modeling as a possible prototype of hazardous events in rift environments.

Coastal areas constitute a very dynamic environment, balancing between numerous natural and anthropogenic processes liable to sometimes hazardous geomorphic phenomena. Especially in tectonically active coastal regions and areas of high economic value, slope failures can have significant impacts and therefore need careful and detailed examination. This work uses Unmanned Aerial System (UAS)-aided photogrammetry and Terrestrial Laser Scanning (TLS) in tectonically active segments of the coastal zone of the Ionian Islands in Greece, to explore how their capabilities can help to improve our understanding of the structural integrity of the slopes. Results show that the two approaches are able to extract large numbers of discontinuity facets, in a more practical, rapid and safe way than conventional methods of rock slope stability analysis extending to unreachable yet important parts of the slope. Through this holistic record of the structural condition of the slope the two applications allow the identification of segments that are more prone to instability and failure. In this way, they improve our understanding of the prioritization of interventions aiming to enhance the prevention of slope failures, mitigating the associated risk and improving local development in these high-value locations.

The Arkalochori village in central Crete was hit by a large earthquake (Mw = 6.0) on 27 September 2021, causing casualties, injuries, and severe damage to the infrastructure. Due to the absence of apparent surface rupture and the initial focal mechanism solution of the seismic event, we initiated complementary, multi-disciplinary research by combining seismological and remote sensing data processing, followed by extensive field validation. Detailed geological mapping, fault surface measuring accompanied with tectonic analysis, fault photorealistic model creation by unmanned aerial system data processing, post-seismic surface deformation analysis by DInSAR image interpretation coupled with accurately relocated epicenters recorded by locally established seismographs have been carried out. The combination of the results obtained from these techniques led to the determination of the contemporary tectonic stress regime that caused the earthquake in central Crete, which was found compatible with extensional processes parallel to the Hellenic arc.

Cephalonia, located in the middle of the central Ionian Islands, has been affected by destructive earthquakes during both the instrumental and the historical period. Despite the fact that it is widely studied from several scientific viewpoints, limited research has been conducted so far regarding the earthquake-triggered landslides (ETL) and the related susceptibility. In the context of the present study, an inventory with 67 ETL from 11 earthquakes that occurred from 1636 to 2014 is presented. Given this record, the study further examines the ETL susceptibility exploiting 10 landslide causal factors in the frame of a GIS-based Analytic Hierarchy Process (AHP). Four factors (i.e., slope, PGA, tectonic structures and lithology) were associated in a higher degree to the locations where ETL occurred on the island. Based on the comparison of the ETL inventory and the landslide susceptibility index (LSI) map, the distribution of ETL in Cephalonia is not random, as their majority (82%) were generated within high to critically high susceptible zones. This fact, along with the AUC values of 80.3%, reveals a fair-to-good accuracy of the landslide susceptibility assessment and indicate that the contribution of the studied variables to the generation of ETL was effectively determined.

Seismicity in the Ionian Sea (W. Greece) is mainly generated along the Cephalonia–Lefkada Transform Fault Zone (CLTFZ) in the central Ionian, and on the northwestern termination of the Hellenic subduction margin in the south. Joint pre-, co- and post-seismic ground deformation and seismological analysis is performed at the broad Ionian area, aiming to homogeneously study the spatiotemporal evolution of the activity prior to and after the occurrence of strong (M > 6) earthquakes during the period of 2014–2018. The 2014 Cephalonia earthquakes (Mw6.1 and Mw5.9) were generated on a faulting system adjacent to CLTFZ, causing local ground deformation. The post-seismic sequence is coupled in space and time with the 2015 Lefkada earthquake (Mw6.4), which occurred on the Lefkada segment of the CLTFZ. Co-seismic displacement was recorded in the broader area. Seismicity was concentrated along the CLTFZ, while its temporal evolution lasted for several months. The 2018 Zakynthos earthquake (Mw6.7) caused regional deformation and alterations on the near-velocity field, with the seismicity rate remaining above background levels until the end of 2021. In the northern Ionian, convergence between the Apulian platform and the Hellenic foreland occurs, exhibiting low seismicity. Seismic hazard assessment revealed high PGA and PGV expected values in the central Ionian.

Messinia is located in SW Peloponnese (Greece), in the vicinity of the Hellenic Arc which is one of the most seismically active areas of Europe. The arc is dominated by reverse faulting, whereas normal faults are mapped onshore, mainly striking N–S. Large earthquakes have occurred in the study area, both in the historical and instrumental periods. In the present study, a Probabilistic Seismic Hazard Assessment (PSHA) is applied to estimate the Peak Ground Acceleration (PGA), Peak Ground Velocity (PGV), Peak Ground Rotational Acceleration (PGRA) and Peak Ground Rotational Velocity (PGRV) for the broader Messinia region. PGRA and PGRV are not often examined in detail in the literature, even though they are useful for the evaluation of the possible damages in structures. The widely used approach proposed by Cornell and McGuire is implemented taking into account: (a) the seismotectonic model proposed by Seismic Hazard Harmonization in Europe (SHARE), (b) an earthquake catalogue for the instrumental period and (c) Ground Motion Prediction Equations (GMPEs) proposed for the Greek territory. The computational grid spacing for Messinia was set to 1 km, in order to accurately calculate the hazard parameters and to reduce the bias of the results through interpolation processes. For PGA and PGV, a logic tree approach is considered, where every branch is a hybrid version of each GMPE considering the percentage of the type (normal or non-normal) of focal mechanisms for all zones of the seismotectonic model. The results show an increase of the hazard values in the NW part of the study area, where the Greek Building Code provides lower PGA values. In addition, hazard curves in terms of PGA for multiple probabilities of exceedance in 50 years are determined for five major towns of Messinia, i.e. Kalamata, Messini, Filiatra, Kyparissia and Pylos. Kyparissia has the higher hazard in all probabilities of exceedance, which is in full agreement with the PGA and PGV results. Pylos and Filiatra have intermediate to high hazard, whereas Kalamata and Messini show intermediate hazard values. Furthermore, the Uniform Hazard Spectrum (UHS) in terms of Spectral Acceleration for the same towns is computed. The latter provides essential information about the design parameters. Lastly, the results are compared to the seismic histories of the five major towns in terms of European Macroseismic Scale (EMS98) intensity, plotted for the last 200 years.

Mediterranean tectonically-active coastal areas are a highly-dynamic environment balancing internal tectonic dynamics with external geomorphic processes, as well as manmade influences. Especially in touristic areas characterized by high built-up pressure and land value, where these dynamics are even more concentrated, the evolution of coastal environments needs careful and high-resolution study to identify localized risk and the processes they derive from. Recently, new advanced remote sensing techniques such as Unmanned Aerial Systems (UAS)- and Terrestrial Laser Scanners (TLS)-aided monitoring have improved our capabilities in understanding the natural processes and the geomorphic risks (i.e. mass movement phenomena). An integrated study comprising Unmanned Aerial Vehicles (UAV) and Light Detection And Ranging (LIDAR) sensors was conducted in coastal areas of the southern Ionian Islands (Western Greece) aiming to the mitigation of earthquake-triggered landslide risk and to responsible coastal development. Located at the northwesternmost part of the Hellenic Arc, this area is characterized by high seismicity and has been affected by destructive earthquakes mainly due to the Cephalonia Transform Fault Zone (CTFZ), which constitutes one of the most seismic active structures in the Eastern Mediterranean region. One of the most common environmental effect triggered by these earthquakes are landslides distributed along fault scarps in developed and highly visited coastal areas. Furthermore, this area is highly susceptible to hydrometeorological hazards inducing intense geomorphic processes, including Medicanes among others. These technologies allow a highly-detailed view of landslide processes, providing insights on the structures and factors controlling and triggering failures along coastal scarps as well as highlighting susceptible zones and high-risk areas with accuracy and mitigating adverse effects with precision and clarity. Overall, by providing a better understanding of the risks the approach used allows a more sustainable development of these coastal segments enhanced by risk mitigation. The study was conducted in the framework of the project “Telemachus - Innovative Operational Seismic Risk Management System of the Ionian Islands”, co-financed by Greece and the European Union (European Regional Development Fund) in Priority Axis “Environmental Protection and Sustainable Development” of the Operational Programme “Ionian Islands 2014–2020”.