The New Seismotectonic Atlas of Greece (v1.0) and Its Implementation

: Knowledge and visualization of the present-day relationship between earthquakes, active tectonics and crustal deformation is a key to understanding geodynamic processes, and is also essential for risk mitigation and the management of geo-reservoirs for energy and waste. The study of the complexity of the Greek tectonics has been the subject of intense e ﬀ orts of our working group, employing multidisciplinary methodologies that include detailed geological mapping, geophysical and seismological data processing using innovative methods and geodetic data processing, involved in surveying at various scales. The data and results from these studies are merged with existing or updated datasets to compose the new Seismotectonic Atlas of Greece. The main objective of the Atlas is to harmonize and integrate the most recent seismological, geological, tectonic, geophysical and geodetic data in an interactive, online GIS environment. To demonstrate the wealth of information available in the end product, herein, we present thematic layers of important seismotectonic and geophysical content, which facilitates the comprehensive visualization and ﬁrst order insight into seismic and other risks of the Greek territories. The future prospect of the Atlas is the incorporation of tools and algorithms for joint analysis and appraisal of these datasets, so as to enable rapid seismotectonic analysis and scenario-based seismic risk assessment.


Introduction
From a geodynamics point of view, Greece and its surrounding areas constitute the most active region of the eastern Mediterranean and, indeed, of Europe [1][2][3][4][5][6][7][8][9][10][11]. The main geotectonic features entail: (a) Continental convergence, comprising the subduction of the oceanic part of the North African plate beneath the European plate, associated with intense crustal shortening and an uplift rate of a few mm/yr along the Hellenic Arc due to accretion of sediments of the African plate beneath the overriding Aegean plate; (b) widespread, high-rate extension in the back-arc region due to the rollback of the subducting African slab; (c) significant right-lateral strike-slip motion along the North Aegean Trough (NAT), and the Cephalonia-Lefkada Transform Zone (CTFZ), which is due to the offset between the oceanic-continental convergence in the west and the westward propagation of the Anatolian plate in the east (Figure 1). Notably, the combination of all the above factors leads to significant left-lateral motion along the southeastern front of the Aegean-African interface. Blue and red arrows indicate trend and plunge of the major (S 1 ) and minor (S 3 ) principal stress axes [12]. Focal mechanisms with M w ≥ 5.0 from the database of NKUA-SL are presented, with colour corresponding to faulting type according to Zoback's classification criteria [13]. CTFZ Inasmuch as earthquake hazards crown the list of possible natural threats in Greece, seismotectonic studies and maps are necessary for the prevention and mitigation of pertinent disasters. It has been over 30 years since the first Seismotectonic Map of Greece was published by the Institute of Geological and Mining Research in 1989 [14]. Therefore, it is imperative to update the map, given that since an extensive body of geological, tectonic, seismological, geophysical, geodetic and volcanological information has been accumulated and dozens of significant earthquakes have occurred. This wealth of information has improved our understanding of the geodynamics of Greece, and could be used to improve our capacity to forecast earthquake hazards and plan the mitigation of disasters in the context of rapidly expanding industrial and economic activities.
The complexity of the contemporary geodynamics setting of Greece has been subjected to continuous and intense study by our research group during the last two decades, in which multidisciplinary state-of-the-art methodologies were applied to the analysis of tectonic, seismological, geophysical and volcanological problems at various scales. The products of these studies are the pieces of a puzzle that we attempt to integrate and superimpose on elevation, land use and other relevant data, in order to compose the digital edition of a new Seismotectonic Atlas of Greece.
The configuration of the tectonic, seismicity and deformation styles of the Greek territories, the spatial configuration of stress and strain, the locations of active faults and their slip-rates, the patterns of geophysical potential fields and a wealth of other information are much better known today, than they were 30 years ago. GIS tools can be used to illustrate such features, typically at hectokilometric to decakilometric scales and, exceptionally, at kilometric scales, wherever there is sufficient information. Accordingly, the new digital Seismotectonic Atlas should be able to generate comprehensive multi-thematic images of up-to-date parameters related to the dynamic behavior of the Earth's crust, which can then be used by earth scientists, higher education students, civil engineers and surveyors, public and private agency operatives and others, in applications ranging from academic research to practical geological hazard and risk assessment, industrial planning, disaster prevention/mitigation etc.
Herein, we present the first edition of the new Seismotectonic Atlas of Greece, which is composed of a set of thematic maps to be used in earthquake and volcanic hazard assessment. Therefore, we introduce the datasets used in the compilation of the new Atlas and their curation to suit its purposes, we provide basic information about the respective thematic layers and, finally, we discuss the outlook and prospect the project. The thematic layers, together with tables containing quantified descriptions of critical parameters, are accessible at http://www.geophysics.geol.uoa.gr/atlas.html.

Design and Implementation of the New Seismotectonic Atlas
The data incorporated in the new Seismotectonic Atlas of Greece are composed of sets of active faults, historical and instrumental era seismicity, focal mechanisms, stress and strain-rate field distributions, active volcanoes, historical and instrumental era tsunamis, gravity and magnetic anomalies. Each of these datasets has been constructed as a GIS layer, augmented with interactive access to enable the user to inspect the respective attributes. Each and any combination of these layers can be superimposed on user-selectable digital elevation models of Greece and surrounding areas. For the layers that include data publicly available by third parties, their associated metadata can be accessed through the Atlas web-page (with due reference to the source). Atlas has been developed with respect to the INSPIRE 2007/2/EK EU Directive, which urges stakeholders to ensure easy access to detailed data and metadata for all public and private sectors, free of charge, or at a cost covering the maintenance of data sets and related data services.

GIS Mapping
The various thematic layers were constructed using the ArcGIS Pro v2.6.2 software. Datasets for seismicity, volcanism, stress-field, strain-rate and tsunami events were collected into MS-Excel sheets and then converted into Point Feature shapefiles. Some data were already available in shapefile format (e.g., infrastructure data) and were imported as layers in ArcGIS Pro. Layers, such as stress-shape, geomagnetic and gravity anomalies, were constructed from the point data after interpolation and rasterization, and were saved into GeoTIFF images. To convey the important information for each layer, several symbology methods were applied. For seismicity we used "Graduated Symbols" with different size classes for specific magnitude ranges, while a secondary classification was applied to the fill color of each symbol upon the focal depth of each event. "Unique values" symbology is used for layers such as Volcanism and Tsunami events, with different symbols and colours assigned according to the type/age of the data.
For vector field representations, the data were divided into multiple layers to facilitate the usage of different symbols for different attributes. More specifically, for the crustal stress-field, one layer depicts the S 1 (blue, double arrows, converging to the center) and another the S 3 principal stress-axes (red, double arrows, diverging from the center). In addition, different classes for the symbol size were applied, proportional to the cosine of the respective axes plunge, i.e., a larger symbol for the more horizontal axes and smaller for the subvertical axes. Secondary symbology is applied to vary rotation according to the axes trend (azimuth).
Similarly, for the surface strain-rate field, a Graduated Symbols primary symbology was applied to relate the symbol size with the magnitude of the strain-rate (in nstrain/yr), while secondary symbology depicts the rotation according to the vectors azimuth. Four layers were used, depending on the sign of each of the two principal strain-rate axes, e 1 , e 2 , with positive values (extension) being depicted by red, divergent double-arrows while negative values (shortening) is represented by blue, converging double-arrows. Simpler symbology was used for layers, such as the active faults, which are line features,

Active Faults: The NOA Database (NOAFaults)
The first edition of the "NOAFaults" database of active faults was published by [15]. The new Seismotectonic Atlas incorporates an updated second edition of the database [16,17]. The active fault database was initially intended to support the monitoring of seismicity by the National Observatory of Athens (NOA); it was constructed with material extracted from at least 110 scientific papers published in international journals as of 1972 and has proven to be useful to many researchers or practitioners of the civil engineering community and public authorities.
The active faults were digitized from the original fault maps that were published in peer-review journals at various cartographic scales. Inside the database there is a field in the fault attribute table that is called "location reliability" with three (3) classifications (precise, approximate and inferred) depending on the cartographic scale used in the digitization (1:50,000, 1:50,000-1:250,000 and above 1:250,000, respectively). In relation to the fault trace geometry, first the line is digitized on-screen so that its average azimuth is in agreement with the right-hand rule in seismology (i.e., a fault strike of 110 • indicates a fault striking ESE-WNW and dipping to the south). In all cases, we inserted in our database the fault attributes provided in the original papers; where those are missing, we put a value of 60 • for a typical normal fault and a rake angle of −90 • for the orientation of the slip vector, etc. We noted that the "NOAFaults" database does not currently consider the seismogenic part of the subducting African (oceanic) plate where most reverse-slip earthquakes originate (see blue beachballs in Figure 1).
Each entry (fault) in the database contains information on the geometry of the fault, its strike, dip direction and dip angle, the type of rupture and, wherever possible, the kinematics of the rupture. The fault database is shown in Figure 2 and can be imported to GIS programs or used as input for other applications and queries. For well-studied faults [18][19][20][21][22][23][24][25][26], the database also reports the long-term slip rates (or its creep) and/or the co-seismic slip of the last earthquake to have taken place. In this first version of the Atlas, we focus on the active faults of the Aegean and Eurasian plates and the back-arc region of the Hellenic Arc (2437 entries). Of those, 95% exhibit normal, 3% strike-slip and only 1% reverse faulting. The published literature information has permitted the correlation between epicenters of strong seismic events with the location of active faults, thus identifying the seismic fault for 57 strong events. The magnitude of those events ranges from M w 5.0 to 7.4. The fault database shows that almost 50% of active faults is longer than 5 km and, therefore, comprise a moderate to high risk level when these faults cross or lie in the vicinity of urban areas. strike-slip and only 1% reverse faulting. The published literature information has permitted the correlation between epicenters of strong seismic events with the location of active faults, thus identifying the seismic fault for 57 strong events. The magnitude of those events ranges from Mw 5.0 to 7.4. The fault database shows that almost 50% of active faults is longer than 5 km and, therefore, comprise a moderate to high risk level when these faults cross or lie in the vicinity of urban areas.

Figure 2.
The database of active faults in the broader Greek region [16,17]. The fault length and rupture type data can be used for the estimation of the maximum expected magnitude and average co-seismic displacement, using appropriate empirical relationships [27][28][29]. This type of information is very useful in planning the use of land and in designing infrastructures and critical facilities close to, or across active faults. Finally, the fault database, in combination with the seismicity, focal mechanism, stress/strain and other databases of the Atlas, could be used for rupture and hazard forecasting, with the incorporation of additional information and with appropriate geoprocessing tools [30][31][32]. Such capabilities will be included in future editions of the Atlas, as discussed in Section 3 below.

Volcanism and Plutonism
Volcanic activity in Greece is not very frequent. Nevertheless, given historical destructive eruptive episodes, it has to be taken into account in modern multi-hazard assessment since numerous Quaternary to Holocene volcanic centers exist throughout the Greek territory ( Figure 3). Several authors posit that this is due to south-westward migration of the volcanic centers caused by the rollback of the subducting African plate [33][34][35]. The most significant features related to volcanic activity are shown in Figure 3 as follows: The fault length and rupture type data can be used for the estimation of the maximum expected magnitude and average co-seismic displacement, using appropriate empirical relationships [27][28][29]. This type of information is very useful in planning the use of land and in designing infrastructures and critical facilities close to, or across active faults. Finally, the fault database, in combination with the seismicity, focal mechanism, stress/strain and other databases of the Atlas, could be used for rupture and hazard forecasting, with the incorporation of additional information and with appropriate geoprocessing tools [30][31][32]. Such capabilities will be included in future editions of the Atlas, as discussed in Section 3 below.

Volcanism and Plutonism
Volcanic activity in Greece is not very frequent. Nevertheless, given historical destructive eruptive episodes, it has to be taken into account in modern multi-hazard assessment since numerous Quaternary to Holocene volcanic centers exist throughout the Greek territory ( Figure 3). Several authors posit that this is due to south-westward migration of the volcanic centers caused by the rollback of the subducting African plate [33][34][35]. The most significant features related to volcanic activity are shown in Figure 3 as follows: 1.
The Hellenic Volcanic Arc, a.k.a. South Aegean active Volcanic Arc (SAVA) of Plio-Quaternary age. Extensive and detailed information can be found in the topical volume edited by [36]. It extends from Sousaki/Methana in the west, to Milos and Santorini in the center and Nisyros to the east, and is characterized by calc-alkaline and high-K calc-alkaline activity ( Figure 3).
The western terminus of the SAVA is a broad volcanic field, consisting of the Sousaki, Aegina, Poros and Methana volcanic centers (Figure 3), all attributed to Pliocene calc-alkaline volcanism. Recent geochemical and geophysical studies have shown that magmas resulted from the upwards migration of fluids in the asthenospheric mantle wedge [43][44][45]. At present, Sousaki is a low-intensity solfatara, whereas Aegina and Poros are practically extinct. The only remaining major active element of the field is the Methana volcanic terrain. Activity may have begun during the late Pliocene [46,47] and built the volcanic edifice through the Pleistocene and Quaternary in several phases (seven according to [47]). The last significant eruption took place circa 238 BCE and effused lava that flowed up to 500 m north-northwest of the coastline; it was mentioned by Pausanias (who also described the post-eruptive thermal activity), Strabo and Ovid. The most recent noteworthy paroxysmic activity was a small-scale submarine effusion in 1700 CE, approximately 2 km N-NW of the 238 BCE eruption (at the so-called Pausanias submarine volcano). It has been argued [47] that tectonic activity is closely knit with the inception and evolution of the volcano in a top-down sense. This has been confirmed [48], and it has also been demonstrated that all recent (eruptive and non-eruptive) activity is completely controlled by contemporary tectonics.
The central SAVA is composed of the Santorini Volcanic Complex, Milos Island and Christiana (islets SW of Santorini) volcanic fields, all associated with tholeiitic and calc-alkaline magmatism. Featuring a high temperature geothermal field, Milos is a presently dormant stratovolcano; the most recent significant activity there was related to a series of phreatic eruptions between 380 and 90 Ka ago [49,50], although other studies propose younger ages [51]. Phreatic explosions, commonly producing overlapping craters of usually sub-kilometric diameter, continued from the late-Pleistocene to recent times. A lahar deposit in south-east Mílos, buried the walls of a Roman harbor town and is thought to have originated from a small phreatic explosion.
The most significant volcanic field of SAVA is the Santorini Volcanic Complex (SVC), which came into being in the Middle Pleistocene. At least twelve eruptions of Plinian intensity have taken place during the last 360 Ka. Each of those discharged material volumes of a few to several cubic kilometres and all together formed pyroclastic deposits, with a thickness of 200 m, containing relics of at least five large shield volcanoes. The intervals between Plinian eruptions average to 30 Ka. The explosive activity triggered at least four caldera collapses and resulted in the formation of the present-day composite caldera structure [52]. The last caldera-forming explosion was the renowned Minoan eruption of the late Bronze Age (1645-1500 BCE) [52][53][54]. Several significant submarine and sub-aerial paroxysms have occurred since, the last one taking place in 1950. The most recent activity occurred during 2011-2012 and was accompanied by crustal dilation and intense seismic activity [55,56]. The contemporary active tectonic modes of the SVC were recently clarified by [57], who illustrated their complete top-down control of the volcanism, and argued that it is physically shaped by tectonic rather than volcanic activity, including the form of the caldera.
The eastern terminus of the SAVA is characterized by several volcanic centers of Late Miocene-Pliocene age, the most important of which are the Nisyros and Kos islands. Kos is currently characterized by the complete absence of volcanic activity, while thermal manifestations are sparse and limited to a few low-temperature thermal springs. Nisyros, for which very extensive and thorough information can be found in the topical volume edited by [58], is still active. The volcano came into being at approximately 300 Ka BCE, with submarine extrusive activity that developed into subaerial dacitic and rhyodacitic after about 161 Ka BCE. A series of Plinian intensity eruptions accompanied by central caldera formation occurred between 44-24 Ka BCE. The post-caldera era began with the extrusion of the rhyolitic-dacitic magmas that cover the western part of the caldera (24 Ka BCE), followed by intermittent low intensity volcanic and continuous high intensity hydrothermal activity. At least eleven phreatic explosions have been reported between 1422 CE and 1888 CE, with the last of these episodes (1871-1873) accompanied by intense seismic activity and damage. The latest minor episode of unrest has taken place in 1996-1997. As evident from [58], local and regional tectonics has had a principal role in the inception and evolution of the volcano. Finally, it has been argued [59] that the contemporary thermal activity is completely controlled and guided by active faulting.
At the north of Evia Island and along an 80 km E-W stretch to the west, there are andesitic outcrops (Lichades islets, offshore the westernmost tip of Evia Island) and several major thermal springs which include the famous Thermopylae. All these are due to mainly plutonic calc-alkaline magmatism clearly identified by geophysical (magnetic anomaly) data and related to the SAVA system [43,[58][59][60].
came into being at approximately 300 Ka BCE, with submarine extrusive activity that developed into subaerial dacitic and rhyodacitic after about 161 Ka BCE. A series of Plinian intensity eruptions accompanied by central caldera formation occurred between 44-24 Ka BCE. The post-caldera era began with the extrusion of the rhyolitic-dacitic magmas that cover the western part of the caldera (24 Ka BCE), followed by intermittent low intensity volcanic and continuous high intensity hydrothermal activity. At least eleven phreatic explosions have been reported between 1422 CE and 1888 CE, with the last of these episodes (1871-1873) accompanied by intense seismic activity and damage. The latest minor episode of unrest has taken place in 1996-1997. As evident from [58], local and regional tectonics has had a principal role in the inception and evolution of the volcano. Finally, it has been argued [59] that the contemporary thermal activity is completely controlled and guided by active faulting.
At the north of Evia Island and along an 80 km E-W stretch to the west, there are andesitic outcrops (Lichades islets, offshore the westernmost tip of Evia Island) and several major thermal springs which include the famous Thermopylae. All these are due to mainly plutonic calc-alkaline magmatism clearly identified by geophysical (magnetic anomaly) data and related to the SAVA system [43,[58][59][60].

Seismicity
Given that the probabilistic and semi-probabilistic (hybrid) assessment of earthquake threat is heavily dependent on the existence of a reliable seismic record, earthquake catalogues are essential in the appraisal of seismic hazard and the mitigation of earthquake disasters. Accordingly, we consider that the inclusion of this type of information in the new Seismotectonic Atlas of Greece, together with tools, currently under development, to assist in the analysis of seismic hazard, structural vulnerability and risk, is absolutely essential. In response, we compiled a unified earthquake catalogue of Greece and surrounding areas on the basis of seven different source catalogues spanning the period 1000 to June 2020 CE and available through different national and international agencies. This allowed us to assemble all the events that have been recorded during the late historical and instrumental eras and may have been included in some catalogues but not in others. The source catalogues are: Having converted all magnitudes to a unified M w , we merged the catalogues, eliminated multiple records of the same events and added from the NKUA catalogue, all M W ≥ 3 events that occurred as of 1/1/2018. The priority (from high to low) in which we retained instrumental events from multiple records was: SHEEC, MKK, ISC-GEM, ISC-EHB and ISC bulletin, while the applied criteria for their detection was an origin time difference of up to 20 s, epicentral distance less than 50 km and magnitude difference (in unified M w ) smaller than 0.5. Finally, we added all M L ≥ 3 from the NOA (January-May 2018) and NKUA-SL catalogues (June 2018 onwards) and also converted the magnitude to M w . We generated a comprehensive seismicity catalogue covering the area of N33.0 • -42.5 • , E18 • -33 • . This is the Atlas Earthquake Catalogue and contains 116,268 events with M W = 2.7-8.3 at a focal depth range 0-263 km. The attributes of the seismicity database, include the date, time, location, depth and unified M w (average of the converted M w magnitudes, in case the single or preferred source reports more than one type of magnitude). Generally speaking, the completeness of the catalogue is 5.8 for the period 1000-1899, drops to 5.2 between 1900 and 1963, then to approximately 4 during 1964-1989 and, finally, to 3.6 during the last 30 years (in the unified M w scale). It is worth noting that seismicity has become spatially denser in particular after 2008 [70], when the various Greek seismological networks were integrated into the Hellenic Unified Seismological Network (HUSN) of broadband receivers.
As can be seen in Figure 4 and the respective online Atlas application, the largest part of seismicity is crustal, at focal depths of less than 30 km, concentrated along major active seismic zones such as the Hellenic Arc, the N. Aegean Trough (NAT) and the Corinth Rift (CR). The NE and SE Aegean exhibit intense activity. Earthquake foci become progressively deeper as one moves from the front of the Hellenic Arc to the Aegean, demonstrating the activity of the north-easterly dipping Benioff zone of the Hellenic Subduction System (HSS). The central Aegean and northern Greece exhibit low seismicity rates.
Geosciences 2019, 9, x FOR PEER REVIEW 9 of 29 when the various Greek seismological networks were integrated into the Hellenic Unified Seismological Network (HUSN) of broadband receivers. As can be seen in Figure 4 and the respective online Atlas application, the largest part of seismicity is crustal, at focal depths of less than 30 km, concentrated along major active seismic zones such as the Hellenic Arc, the N. Aegean Trough (NAT) and the Corinth Rift (CR). The NE and SE Aegean exhibit intense activity. Earthquake foci become progressively deeper as one moves from the front of the Hellenic Arc to the Aegean, demonstrating the activity of the north-easterly dipping Benioff zone of the Hellenic Subduction System (HSS). The central Aegean and northern Greece exhibit low seismicity rates. Highest seismicity rates are observed along the CTFZ and the Hellenic Arc. The western part of the Hellenic Arc presents the major part of seismic energy release across the Greek territory. In particular the area of the CTFZ and offshore NW Peloponnese, which according to the European Seismic Hazard Model (ESHM-13) [71] is the most seismically prone area of Europe, has been often activated recently, producing strong and occasionally damaging earthquakes [72][73][74][75][76][77][78]. An interesting feature is a NE-SW boundary of seismicity between the CTFZ and Zakynthos Island, marking an area that has been suggested to host a Subduction-Transform-Edge-Propagator (STEP) fault, due to a trench-perpendicular tearing of the slab [79][80][81][82][83], or a ramp-type structure [84,85] beneath the CTFZ, bridging NAT with the Hellenic Trench [86], possibly the margin between the continental and oceanic slab of the Hellenic Subduction System (HSS) [43,84,87]. This area is well-known for transpressional tectonics [12,88], while the recent Mw6.7 earthquake on 25th October 2018 provided strong indications for seismic coupling in the western region of the Hellenic Arc with most of co-seismic deformation taken up by the upper (Aegean) plate [77], with a potentially significant effect for the regional seismic hazard. Highest seismicity rates are observed along the CTFZ and the Hellenic Arc. The western part of the Hellenic Arc presents the major part of seismic energy release across the Greek territory. In particular the area of the CTFZ and offshore NW Peloponnese, which according to the European Seismic Hazard Model (ESHM-13) [71] is the most seismically prone area of Europe, has been often activated recently, producing strong and occasionally damaging earthquakes [72][73][74][75][76][77][78]. An interesting feature is a NE-SW boundary of seismicity between the CTFZ and Zakynthos Island, marking an area that has been suggested to host a Subduction-Transform-Edge-Propagator (STEP) fault, due to a trench-perpendicular tearing of the slab [79][80][81][82][83], or a ramp-type structure [84,85] beneath the CTFZ, bridging NAT with the Hellenic Trench [86], possibly the margin between the continental and oceanic slab of the Hellenic Subduction System (HSS) [43,84,87]. This area is well-known for transpressional tectonics [12,88], while the recent M w 6.7 earthquake on 25th October 2018 provided strong indications for seismic coupling in the western region of the Hellenic Arc with most of co-seismic deformation taken up by the upper (Aegean) plate [77], with a potentially significant effect for the regional seismic hazard.
Gaps of seismicity are aligned orthogonally to the western Hellenic Arc (N-E to N-S), extending from the south of CTFZ to central Crete, related to~N-S normal faulting of the upper (Aegean) plate. Seismicity gaps are observed only in the part of the Hellenic Arc where N-S normal faulting occurs. The easternmost seismicity gap occurs in an area where faulting rotates from~N-S to~E-W.
These seismicity gaps could be explained by slab tears perpendicular to the western Hellenic Arc [89], allowing extensive subcrustal convection/flux from the mantle wedge [43,84]. However, the seismic gap, located offshore south Peloponnese and western Crete, is probably due to the low earthquake detectability of the seismological networks.
East of central Crete, seismicity appears continuous, distributed in a WNW-ESE direction, related to E-W and NE-SW normal faults. East of Crete, a NNW-SSE structure, perpendicular to the Hellenic Arc bounds dense shallow seismicity to its west and sparser one, mainly subcrustal to the east. This area is suggested to be related to a STEP fault beneath the western coast of the south Anatolian plate, separating the Hellenic Subduction Zone (HSZ) from the western Cyprus Subduction Zone (CSZ) [80,[90][91][92][93], or horizontal tearing of the slab that likely explains the presence of a NW dipping slab [94]. This NNW-oriented seismicity boundary is in agreement with [87], who suggest slab segmentation east of Crete, with the eastern segment dipping with a steeper angle compared to the western one. Despite the importance of the pattern for the regional geodynamics, the new seismicity catalogue indicates that deep seismic sources mainly control the seismic hazard of this area. East of HSZ, seismicity is sparser, clustered, distributed in a NE-SW direction, along the southern margin of the hydrothermal field of west Anatolia [61][62][63].
Instrumental subcrustal seismicity, located at depths greater than 100 km, highlights the amphitheatrically-shaped SAVA north of the convergence front of North Africa and Aegean. Subcrustal earthquakes coexist with crustal seismicity in the SAVA, with the major seismicity rates observed east of the SVC, mainly due to the recent unrest of Santorini volcano during 2011-2012 [54], and the M6. 6 Kos earthquake on 20 July 2017 [95]. A prominent feature highlighted from seismicity is a 55 km long concentration of shallow events along the Amorgos fault that produced a catastrophic tsunamigenic event in 1956 with M7.7, the largest magnitude observed in the Aegean during the instrumental era [96,97]. To the West of the SVC, seismicity is sparse. However, an increase in seismic activity is observed at the westernmost tip of SAVA, in the region of Methana volcano and Saronikos Gulf, monitored after April 2016 [98]. This phenomenon should be taken under consideration with respect to a likely increase of seismic hazard in the Athens metropolitan area, the center of which lies less than 50 km away from the activated region. In this respect, it should be noted that Athens was recently struck by an M w = 5.1 event on 19 July 2019 [99,100], which reactivated part of the structure that had hosted the M w = 6.0 earthquake of 7 September 1999 [101].
Central Aegean appears aseismic, while E-W structures offshore west Anatolia and along the Corinth Rift are delineated by seismicity at these latitudes, signature of the Late Cenozoic evolution of the Aegean [102,103]. In North Aegean, seismicity indicates major E-W to NE-SW structures related to the westward extrusion of North Anatolia in the Aegean. This area has been often activated during the last years, producing strong earthquakes at times damaging for the nearby islands [103][104][105][106][107][108][109][110][111], with the most recent destructive one having occurred offshore southern Lesvos Island on 12th June 2017 [111,112]. The epicenters distribution of the compiled catalogue does not imply continuation of these strike-slip structures in continental Greece, where the most active structure is the Corinth Rift, and in particular its western part. On the contrary, major active NW-SE-oriented structures, such the one that produced the 2001 M6.4 Skyros earthquake [103][104][105], Evia Island and the Corinth Rift likely constitute barriers for the western propagation of NAT.
West of Corinth Rift, the aseismic Cephalonia-Lefkada-Acarnania Block (CLAB) is evidenced, a tectonically complex area bounded by transpressional and transtensional structures [12,88,113,114]. North of CLAB the main part of seismicity is distributed beneath NW Greece and Albania, mostly associated with thrust faulting [12,[114][115][116]. East of Amvrakikos Gulf, seismicity is related to normal faulting, therefore, a continuation of NAT within continental Greece is likely not supported by seismic activity. Clustered seismicity in continental Greece is mainly related to seismic sequences and particularly those monitored by local temporary networks.

Focal Mechanisms, Stress and Strain-Rate Field
Focal mechanisms are essential in understanding the seismotectonic setting of an area and, together with seismicity, are necessary for the definition of seismic sources to be used in seismic hazard assessment. In this section, we present the catalogue of double-couple earthquake focal mechanisms, used in the new Seismotectonic Atlas, together with the stress field derived therefrom by formal inversion.
In cases of duplicate events, a "preferred" solution was chosen by expert judgment, based on the compliance of a focal mechanism with known local and regional geodynamics and the precision of its hypocentral location [114]. The majority of post-2008 focal mechanisms correspond to events located using recordings of the HUSN. However, a large number of events have been relocated [114] using different local and regional crustal velocity models, each optimized so as to yield improved hypocentral locations for earthquakes recorded (primarily) by local networks, as well as by the HUSN [114]. Figure 5 presents the distribution of the focal mechanisms, color-coded according to the type of rupture following the classification criteria of [13]. The type of rupture according to the focal mechanisms is generally consistent with the type of rupture registered in the active faults database (see Section 2.1; Figure 2). Specifically, normal faulting predominates in the Greek mainland and along the eastern Aegean island and the western coast of Turkey, thrust faulting predominates in the Hellenic Arc and the Apulian collision front (NW Greece-W Albania). Finally, strike-slip faulting is the main type of rupture along the CTFZ and NW Peloponnese, the N. Aegean Trough and southeastern boundary of the Aegean and African plates (Strabo and Pliny trenches; see Figure 2). Focal mechanisms indicate transpressional tectonics along the Hellenic Arc and transtensional across the overriding Aegean plate. The attributes of the focal mechanism database [12] include the location, magnitude (M w ), strike, dip and rake for both nodal planes, the trend and plunge of the principal moment tensor axes, the slip vectors, the faulting type according to [13] and the bibliographic source of each entry. In addition, the Atlas currently includes a layer with focal mechanisms of significant earthquakes at both shallow and intermediate focal depths for the period 1995-June 2020 from the database of NKUA-SL.  [13]. Data from [12].
The determination of the stress field is important for the assessment of a region's seismotectonic setting and associated seismic hazard. The inversion of focal mechanisms of earthquakes is a robust procedure, which in combination with strain-rates deduced from satellite geodesy and tectonic field observations, can provide constraints on the three-dimensional (3D) deformation of the crust [114]. Using the focal mechanism database and implementing a least-squares approach, based on the Wallace-Bott hypothesis [125] that faults slip in the direction of the maximum shear stress on the fault plane, the contemporary stress-field of Greek territories and surrounding areas has been computed [12] in a grid of 0.25° × 0.25° cells (Figures 6-8). As described in [12] the inversion was performed using a damped least-squares approach, with the assumption of stress uniformity within each individual cell while, simultaneously, minimizing differences between neighboring cells. In summary: • The distribution of the stress field is consistent with the faulting style observed for major faults.

•
Transtensional stress is observed in the Aegean, the Corinth Rift and W. Greece.

•
Transpressional stress predominates in NW Greece, the Hellenic Arc and Crete.

•
In the Aegean, the SAVA roughly marks a zone of transition, across which the stress field rotates by 90°: The S1 axis changes from almost E-W to the north of approximately 36.5° N to almost N-S to the south of that parallel. It also coincides with a sharp contrast in the crustal stress-shape ( Figure 7) and a deflection in the orientation of SKS and Rayleigh wave anisotropy [43; references therein]. An analogous transition can be observed in the Peloponnese relative to approximately 37° N, with the difference that field there is mainly N-S extensional to the north of this parallel and mainly E-W extensional to the south; this effect is highlighted by high stress-shape values R = 1 − (σ2 − σ3)/(σ1 − σ3), where σ1, σ2, σ3 are the magnitudes of the principal stress axes, with σ1 > σ2 > σ3).

•
These abovementioned observations are in good agreement with the strain-rate field of the Global Strain Rate Model [126]. Surface Strain Rate (SSR; Figure 8) appears compatible with seismicity ( Figure 4). High shortening rate is observed along the Hellenic Arc, whereas high extensional The determination of the stress field is important for the assessment of a region's seismotectonic setting and associated seismic hazard. The inversion of focal mechanisms of earthquakes is a robust procedure, which in combination with strain-rates deduced from satellite geodesy and tectonic field observations, can provide constraints on the three-dimensional (3D) deformation of the crust [114]. Using the focal mechanism database and implementing a least-squares approach, based on the Wallace-Bott hypothesis [125] that faults slip in the direction of the maximum shear stress on the fault plane, the contemporary stress-field of Greek territories and surrounding areas has been computed [12] in a grid of 0.25 • × 0.25 • cells (Figures 6-8). As described in [12] the inversion was performed using a damped least-squares approach, with the assumption of stress uniformity within each individual cell while, simultaneously, minimizing differences between neighboring cells. In summary: • The distribution of the stress field is consistent with the faulting style observed for major faults.

•
Transtensional stress is observed in the Aegean, the Corinth Rift and W. Greece.

•
Transpressional stress predominates in NW Greece, the Hellenic Arc and Crete.

•
In the Aegean, the SAVA roughly marks a zone of transition, across which the stress field rotates by 90 • : The S 1 axis changes from almost E-W to the north of approximately 36.5 • N to almost N-S to the south of that parallel. It also coincides with a sharp contrast in the crustal stress-shape ( Figure 7) and a deflection in the orientation of SKS and Rayleigh wave anisotropy [43; references therein]. An analogous transition can be observed in the Peloponnese relative to approximately 37 • N, with the difference that field there is mainly N-S extensional to the north of this parallel and mainly E-W extensional to the south; this effect is highlighted by high stress-shape values R = 1 − (σ 2 − σ 3 )/(σ 1 − σ 3 ), where σ 1 , σ 2 , σ 3 are the magnitudes of the principal stress axes, with σ 1 > σ 2 > σ 3 ).

•
These abovementioned observations are in good agreement with the strain-rate field of the Global Strain Rate Model [126]. Surface Strain Rate (SSR; Figure 8) appears compatible with seismicity ( Figure 4). High shortening rate is observed along the Hellenic Arc, whereas high extensional deformation rate is observed along a line that connects CTF and NW Peloponnese with NAT, an area that has been suggested to constitute the boundary between the continental and oceanic slab of the HSS [43,84]. An interesting implication is the consistency between stress and strain-rate, regarding the inferred faulting type, particularly apparent in areas of complex tectonics, i.e., transpressional or transtensional. Low strain-rate can be observed at regions where the seismicity is also low, i.e., in central Aegean and northern Greece (see also Figure 4).
Geosciences 2019, 9, x FOR PEER REVIEW 13 of 29 deformation rate is observed along a line that connects CTF and NW Peloponnese with NAT, an area that has been suggested to constitute the boundary between the continental and oceanic slab of the HSS [43,84]. An interesting implication is the consistency between stress and strain-rate, regarding the inferred faulting type, particularly apparent in areas of complex tectonics, i.e., transpressional or transtensional. Low strain-rate can be observed at regions where the seismicity is also low, i.e., in central Aegean and northern Greece (see also Figure 4). Figure 6. Crustal stress field derived from the inversion of earthquake focal mechanisms [12]. Blue and red arrows indicate major (S1), and minor principal stress axes (S3), respectively. Black lines are active faults [16,17]. Figure 6. Crustal stress field derived from the inversion of earthquake focal mechanisms [12]. Blue and red arrows indicate major (S 1 ), and minor principal stress axes (S 3 ), respectively. Black lines are active faults [16,17].   σ 3 ), derived from the inversion of earthquake focal mechanisms [12]. Black lines are active faults [16,17]. Figure 7. Crustal stress shape, R = 1 − (σ2 − σ3)/(σ1 − σ3), derived from the inversion of earthquake focal mechanisms [12]. Black lines are active faults [16,17].

Tsunamis
The new Seismotectonic Atlas of Greece contains the tsunami catalogue of [127]. This includes 162 events, starting with the late Bronze Age (Minoan eruption of Santorini volcano) and ending Figure 8. Surface Strain-Rate (SSR) [126]. Black lines are active faults [16,17].

Tsunamis
The new Seismotectonic Atlas of Greece contains the tsunami catalogue of [127]. This includes 162 events, starting with the late Bronze Age (Minoan eruption of Santorini volcano) and ending with a sea-level oscillation observed at Heraklion, Crete, in year 2000 CE. Information for each entry of the catalogue is derived from original historical documents, previous tsunami catalogues, scientific reports, studies and books and, in a few cases, from field observations made by [127]. Of this dataset, we chose to include 160 events that have been associated with geographic coordinates (Figure 9; also see Table 3 for the general characteristics). The attributes of each event in the original tsunami database [127], include the ID number, the time of occurrence, the Reliability Index (RI) of the time of occurrence, the region (M1 = Greek and adjacent regions in the European Tsunami Catalogue), the originating event, the location of the originating event in geographic coordinates, the RI of that location, a short description of the event, the intensity (I 0 ) in the Modified Mercalli Intensity scale, the surface-wave magnitude (M s ) and the focal depth of the originating event (if the tsunami was caused by an earthquake), the Volcanic Explosivity Index of the originating event (if the tsunami was caused by a volcanic eruption), the maximum vertical tsunami run-up (in cm), the Tsunami Intensity (TI) in the scale of [128], its respective RI and, finally, an index showing if the tsunami parameters were revised (Y) or not (N) in reference to previous source catalogues.  Table 3 for explanation). The year of occurrence is also indicated. Data from [127]. Black lines are active faults [16,17].
As evident in Figure 9, 23 tsunami events are located inside the Sea of Marmara and the Dardanelles, 19 in the gulfs of Corinth and Patras, 15 around Rhodes and SW Turkey, 19 across the western Peloponnese and the central Ionian Sea, and 13 in the south Adriatic Sea; these are associated with shallow, strong (M6+) earthquakes. Recent events, such as the M6. 5 Lefkada earthquake of 2015 [77], the M6.6 Kos earthquake of 2017 [95] and the M6.7 Zakynthos event of 2018 [125], confirmed the recurrence of tsunamis in the Dodecanese and the central Ionian Sea. Overall, it appears that significant tsunamis are frequent in Greece and its surrounding areas, and that the accurate mapping of tsunami occurrence is important for both academic research and for tsunami disaster management [130].

Magnetic Field
The onshore elements of the unified and homogeneous magnetic anomaly map of Greece are based on the 1:50,000 aeromagnetic map series of IGME, produced by ABEM AB and Hunting Geology and Geophysics Ltd. Details can be found in [131,132]. The ABEM maps [133] were produced with measurements at a nominal constant ground clearance (CGC) of 275 ± 75 m above ground level (AGL). Flight line (track) direction was NE-SW with a mean distance (spacing) of 1 km between tracks. Connection (tie) lines for crossover adjustment were also flown in a NE-SW direction, with a mean inter-track distance of 10 km. Hunting Geology and Geophysics Ltd. has generally used a measurement spacing of 200-250 m and has flown in a NE-SW direction but at different CGC and track densities. These would be 150 m AGL at 400 m spacing in areas of known massive ultramafic rock outcrops and 300 m AGL at 800 m spacing otherwise. The high  Table 2 for explanation). The year of occurrence is also indicated. Data from [127]. Black lines are active faults [16,17]. Table 3. Characteristics of the tsunami database included in the Atlas. RI: Reliability Index (RI-Rel_2 field in the Atlas database). Greece and adjacent areas 4 45 In terms of the reliability of the tsunami events, a modified version of the tsunami reliability scale of [129] was adopted (0 = very improbable tsunami, 1 = improbable, 2 = questionable, 3 = probable, 4 = definite). Events assigned with a reliability index of 0 were not included in the database. Some others, described as sea disturbances or sea-quakes, were included with a reliability of 1 or 2, if it was not certain that the event was not a tsunami.

Number of Events 160 RI Counts
Tsunami origins are classified in three main classes that include sub-classes: Class 1 includes earthquakes with subclasses ER (submarine earthquake), EA (earthquake associated), EL (earthquake-caused landslides), ES (earthquake-caused submarine landslides). Class 2 includes volcano paroxysms with subclasses VO (submarine eruption) VA (volcano-associated), VL (volcano-related landslides), VS (volcano-related submarine landslides). Class 3 includes only slumps, i.e., gravitational submarine slides (GS). The database contains 143 tsunamis related to earthquakes, 5 related to volcanic activity and 3 related to slumps (Table 2 and Figure 9).
As evident in Figure 9, 23 tsunami events are located inside the Sea of Marmara and the Dardanelles, 19 in the gulfs of Corinth and Patras, 15 around Rhodes and SW Turkey, 19 across the western Peloponnese and the central Ionian Sea, and 13 in the south Adriatic Sea; these are associated with shallow, strong (M6+) earthquakes. Recent events, such as the M6. 5 Lefkada earthquake of 2015 [77], the M6.6 Kos earthquake of 2017 [95] and the M6.7 Zakynthos event of 2018 [125], confirmed the recurrence of tsunamis in the Dodecanese and the central Ionian Sea. Overall, it appears that significant tsunamis are frequent in Greece and its surrounding areas, and that the accurate mapping of tsunami occurrence is important for both academic research and for tsunami disaster management [130].

Magnetic Field
The onshore elements of the unified and homogeneous magnetic anomaly map of Greece are based on the 1:50,000 aeromagnetic map series of IGME, produced by ABEM AB and Hunting Geology and Geophysics Ltd. Details can be found in [131,132]. The ABEM maps [133] were produced with measurements at a nominal constant ground clearance (CGC) of 275 ± 75 m above ground level (AGL). Flight line (track) direction was NE-SW with a mean distance (spacing) of 1 km between tracks. Connection (tie) lines for crossover adjustment were also flown in a NE-SW direction, with a mean inter-track distance of 10 km. Hunting Geology and Geophysics Ltd. has generally used a measurement spacing of 200-250 m and has flown in a NE-SW direction but at different CGC and track densities. These would be 150 m AGL at 400 m spacing in areas of known massive ultramafic rock outcrops and 300 m AGL at 800 m spacing otherwise. The high mountainous zones Olympus in Thessaly/Central Macedonia and Taygetus in the Southern Peloponnese were surveyed at constant elevations of 3000 m and 2300 m above mean sea level (AMSL) respectively, with 1000 m track spacing in both cases. In all cases, connection lines were flown in a NE-SW direction and a spacing of 10 km. The IGRF (International Geomagnetic Reference Field) correction was based on the IGRF model for the epoch 1977.3. A constant value of 150 nT was added to the magnetic anomaly values before plotting the contour maps. The offshore elements were extracted from the magnetic anomaly maps of [134], spanning the north and south Aegean, Mirtoan Sea and large parts of the Cretan seas.
The original map sheets were converted to high-resolution raster images and digitized to vector form in image coordinates. Using the corners of the map sheets as control points, the image coordinates were geo-located and transformed to Cartesian in the UTM (Universal Transverse Mercator) projection. The digitized data was finally interpolated to a rectangular grid with uniform spacing of 250 m and collated to generate Digital Magnetic Anomaly Models, one for each survey area and data collection scheme. For Hunting data, the arbitrary 150 nT value was removed from the grids. In the case of onshore data, the IGRF model of epochs 1977.3 (Hunting survey) and 1966.5 (ABEM survey) was re-calculated at 0 m elevation and added to the grids. Thus, the original total intensity measurements were reproduced as representatively as possible. Next, the Digital Magnetic Anomaly Models were re-computed, this time by subtracting the Definite IGRF for the corresponding epochs and CGC's of each magnetic anomaly model; calculations were performed on the basis of the 3 arcsec Shuttle Radar Topography Mission (SRTM) digital elevation model. The improvement in the accuracy of the magnetic anomalies afforded by this reprocessing scheme is significant and ranges from a few nT to several tens of nT at low elevation areas [132]. In a final step, the mosaic of magnetic anomaly models was "homogenized" by referral to the common CGC of 1300m AGL/AMSL. This entailed field continuation between arbitrary surfaces and was carried out with the equivalent sources method of [135]. The details of the procedure can be found in [132]. The continued Magnetic Anomaly Models were subsequently collated and reduced to the pole (RTP) to generate the map shown in Figure 10.
the mosaic of magnetic anomaly models was "homogenized" by referral to the common CGC of 1300m AGL/AMSL. This entailed field continuation between arbitrary surfaces and was carried out with the equivalent sources method of [135]. The details of the procedure can be found in [132]. The continued Magnetic Anomaly Models were subsequently collated and reduced to the pole (RTP) to generate the map shown in Figure 10.  The principal sources of high-amplitude magnetic responses in mainland Greece are ophiolite sequences that are clearly observed in Macedonia and Thessaly (zonal NW-SE trending anomalies). In the rest of the mainland, magnetic anomalies are more or less local and generated by large near-surface ophiolitic, olistolites and olistostromes embedded in tectonic mélange, or by dismembered ophiolite complexes emplaced by alpine thin-skinned tectonics. Accordingly, they do not stand out at the elevation and scale of Figure 10. The most remarkable observation, however, is the abundance and intensity of marine magnetic anomalies, the most reasonable explanation of which is widespread (calc-alkaline) plutonic magmatism. The alignments of such anomalies can clearly be observed along the North Aegean Trough (NAT), in the Central Aegean, across the NE Peloponnese (Mirtoan Sea-Argolic Gulf), along the Ikaria-Naxos-Paros-Milos line in the South Aegean Sea and along the northern margin of the Cretan Sea. In all cases, magmatism develops either directly along, or adjacent to crustal-scale dislocation surfaces that can be strike-slip, oblique-normal or normal fault zones. A number of studies along the SAVA have also demonstrated that magma emplacement and thermal fluid circulation is directly associated with contemporary active faulting, and specifically at the intersection of regional-scale fault zones [42,46,55]. This would imply that earthquake and volcanic hazards are causally related and in the grand scheme of things, it might be possible to improve the appraisal of volcanic risks by detailed monitoring of tectonic deformation in the broader areas of major volcanic fields.

Gravity Field
The Gravity Anomalies Map of the region is based on the data included in the Gravity Data Bank (UoA-GDB), compiled at the Section of Geophysics & Geothermics by [136]. The number of gravity stations in the UoA-GDB was originally 33,000 points; these have been augmented with data collected from a variety of sources, so that the present count is approximately 290,000. All these have undergone a procedure of homogenization that reduced the datasets from different sources to a common reference (datum), and included recalculation of the free-air, Bouguer and terrain corrections on the basis of common digital elevation models. The calculations were carried out using vertical prisms with the algorithm originally described by [137,138], as modified by [139] to calculate the gravity effect of topography on a spherical Earth. After applying the necessary reductions, the standard Bouguer anomaly was calculated using a density 2.67 gr/cm 3 . The terrain corrections were calculated to a radius of 167 km around each gravity station. For land data, a detailed 100 m DEM was used for radii up to 22 km around the station and a coarse 1 km DEM for radii between 22 and 167 km. On the basis of the final homogeneous data bank, and depending on the density of gravity stations, local and regional gravity anomaly maps can be produced by interpolation, to the desired level of detail: the regional complete Bouguer anomaly map shown in Figure 11 is based on a 1 × 1 km grid. The map clearly shows the mode of crustal thickening as one moves from the abyssal plains of the oceanic African crust (Ionian Sea, Libyan Sea and Eastern Mediterranean) into the continental crust of the Aegean plate, and the subsequent thinning of the backarc region (Cretan sea-Aegean Sea). In the Southern Aegean area the gravity signal of the dipping African Slab is also detectable.

Isostatic Anomalies
The theory of Isostasy explains the response of the Earth's crust to the gravitational load exerted by changes in topography (topographic load). The quantification of this response is based on isostatic models, and its application to the complete Bouguer anomaly (isostatic correction) generates the isostatic anomaly. At regional (hectokilometric) scales, isostatic anomalies expose the presence (or lack thereof) uncompensated or overcompensated topographic loads. At local (kilometric/decakilometric) scales, isostatic anomalies enhance the expression of density contrasts due to geology and tectonics. In a tectonically active realm, both results are useful in the context of a seismotectonic Atlas. Regional isostatic anomalies provide information related to the geodynamics, as well as to "deep" structures connected with crustal-scale tectonics and seismicity. Local isostatic anomalies illustrate the geometric relationship between the shape of geological structures and active tectonics, as the latter is expressed by active fault lineaments and/or seismicity. Previous studies [136][137][138][139] indicated that the best fitting Airy model computed on the basis of our data is one with a mean crustal depth of 25 km. Subsequent unpublished validation tests have led to similar conclusions. Accordingly, the application of isostatic corrections based on this model to the Bouguer anomaly map of Figure 11, generated the isostatic anomaly map of Figure 12, which is the version Figure 11. The complete Bouguer Anomaly Map. Fault lines are after [16,17]. AG: Argolic Gulf, NAT: North Aegean Trough.

Isostatic Anomalies
The theory of Isostasy explains the response of the Earth's crust to the gravitational load exerted by changes in topography (topographic load). The quantification of this response is based on isostatic models, and its application to the complete Bouguer anomaly (isostatic correction) generates the isostatic anomaly. At regional (hectokilometric) scales, isostatic anomalies expose the presence (or lack thereof) uncompensated or overcompensated topographic loads. At local (kilometric/decakilometric) scales, isostatic anomalies enhance the expression of density contrasts due to geology and tectonics. In a tectonically active realm, both results are useful in the context of a seismotectonic Atlas. Regional isostatic anomalies provide information related to the geodynamics, as well as to "deep" structures connected with crustal-scale tectonics and seismicity. Local isostatic anomalies illustrate the geometric relationship between the shape of geological structures and active tectonics, as the latter is expressed by active fault lineaments and/or seismicity. Previous studies [136][137][138][139] indicated that the best fitting Airy model computed on the basis of our data is one with a mean crustal depth of 25 km. Subsequent unpublished validation tests have led to similar conclusions. Accordingly, the application of isostatic corrections based on this model to the Bouguer anomaly map of Figure 11, generated the isostatic anomaly map of Figure 12, which is the version incorporated in the new Seismotectonic Atlas.

Summary and Conclusions
We developed the first version of the new Seismotectonic Atlas of Greece, in order to build a comprehensive inventory and reference source of up-to-date data, which will be useful in a variety of applications and will benefit the theoretical and applied geosciences and to the society at large. Specifically, the new Seismotectonic Atlas is a fully digital, GIS-based platform, comprising a collection of up-to-date elevation, neotectonic (active faults), volcanological, seismological (historical/instrumental seismicity and focal mechanisms), stress and strain-rate (deformation), geophysical (gravity and magnetic) and tsunami data. These datasets were assembled from a variety of sources, including public international databases, published scientific and unpublished grey literature and seismicity/earthquake focal mechanism catalogues, compiled by the Seismological Laboratory of the NKUA (SL-NKUA).
These datasets are kept in separate databases and constitute independent thematic layers that can be displayed in various combinations. This enables the creation of detailed multi-thematic maps of different types of geodynamic and anthropogenic activity at decakilometric to hectokilometric scales. These maps are useful to several geo-scientific applications, both at the academic and practical levels of studying and understanding the earthquake risks of Greek territories. Single-and multi-thematic maps of this type are accessible at the URL http://www.geophysics.geol.uoa.gr/atlas.html. At present, the user can view the layers in the open access online platform choose which one to make visible, click on point and line features to acquire specific information, or collectively view the attribute table of feature-type data layers and import layers to proprietary or free GIS applications.
While already rich in information, the new Seismotectonic Atlas is planned to be extended with several additional thematic layers. Specifically, future plans include:

•
Introduction of geological and geotechnical data (physical and mechanical rock properties, VS30 estimates wherever available, etc.). This may enable the rapid classification of soils and could be used in the realistic quantitative estimation of seismic hazard and risk. VS30 and site amplification data have already been obtained by SL-NKUA through microtremor and

Summary and Conclusions
We developed the first version of the new Seismotectonic Atlas of Greece, in order to build a comprehensive inventory and reference source of up-to-date data, which will be useful in a variety of applications and will benefit the theoretical and applied geosciences and to the society at large. Specifically, the new Seismotectonic Atlas is a fully digital, GIS-based platform, comprising a collection of up-to-date elevation, neotectonic (active faults), volcanological, seismological (historical/instrumental seismicity and focal mechanisms), stress and strain-rate (deformation), geophysical (gravity and magnetic) and tsunami data. These datasets were assembled from a variety of sources, including public international databases, published scientific and unpublished grey literature and seismicity/earthquake focal mechanism catalogues, compiled by the Seismological Laboratory of the NKUA (SL-NKUA).
These datasets are kept in separate databases and constitute independent thematic layers that can be displayed in various combinations. This enables the creation of detailed multi-thematic maps of different types of geodynamic and anthropogenic activity at decakilometric to hectokilometric scales. These maps are useful to several geo-scientific applications, both at the academic and practical levels of studying and understanding the earthquake risks of Greek territories. Single-and multi-thematic maps of this type are accessible at the URL http://www.geophysics.geol.uoa.gr/atlas.html. At present, the user can view the layers in the open access online platform choose which one to make visible, click on point and line features to acquire specific information, or collectively view the attribute table of feature-type data layers and import layers to proprietary or free GIS applications.
While already rich in information, the new Seismotectonic Atlas is planned to be extended with several additional thematic layers. Specifically, future plans include:

•
Introduction of geological and geotechnical data (physical and mechanical rock properties, V S30 estimates wherever available, etc.). This may enable the rapid classification of soils and With a view to improve the appraisal of earthquake hazard, risk and damage, the introduction of GIS and other tools will enable the extraction of data from any combination of databases (thematic layers) and their joint statistical analysis and interpretation. Deterministic and/or probabilistic (or hybrid) methods of Seismic Hazard Assessment (SHA) at regional and local scale will be facilitated in terms of data detail and completeness. For example, HELPOS-Risk [143] is a GIS tool currently implemented by SL-NKUA in a semi-automatic manner, by applying Deterministic SHA (DSHA) while taking into account empirical seismic vulnerability and site effects. The Atlas will soon incorporate urban risk maps on the basis of earthquake scenarios in seismically prone areas. The connection of the tool to the SL-NKUA real-time SeisComP system, where it will receive alerts and information about the seismic source, will facilitate the automatic HELPOS-Risk execution. This will allow rapid damage assessments. It will support timely disaster management decisions and more accurate first responses. A simplified version of Atlas maps accompanied with a tutorial on the terms of usage will be developed for educational purposes of students and information of civilians and stakeholders.
Last but not least, it is important that the databases be updated on a regular basis for the new Seismotectonic Atlas of Greece to remain state-of-the-art. This could likely be ensured by systematic updates of seismicity layers by the SL-NKUA and by collaboration with third parties. As an example, we refer to a planned incorporation of Atlas within cutting-edge multi-hazard platforms that combine geohazard data with emergency management. Ideally, this would include integration of real-time shake maps produced by dense ground motion recordings and calibrated crowdsourcing information on the severity of an earthquake, as for instance LastQuake (https://www.emsc-csem.org/service/ Browser-extension/Lastquake/) which could largely improve the response during crises.