Here we present the results of a kinematic slip model of the 2020 M_w 6.7 Doğanyol-Sivrice, Turkey Earthquake, the most important event in the last 50 years on the East Anatolian Fault zone. Our slip model is constrained by two Sentinel-1 interferograms and by 5 three-component high-rate GNSS recordings close to the earthquake source. We find that most of the slip occurs predominantly in three regions, two of them at between 2 and 10 km depth and a deeper slip region extending down to 20 km depth. We also relocate the first two weeks of aftershocks and find a distribution of events that agrees with these slip features. The HR-GNSS recordings suggest a predominantly unilateral rupture with the effects of a directivity pulse clearly seen in the waveforms and in the measure peak ground velocities. The slip model supports rupture propagation from northeast to southwest at a relatively slow speed of 2.2 km/s and a total source duration of ∼20 s. In the absence of near-source seismic stations, space geodetic data provide the best constraint on the spatial distribution of slip and on its time evolution.

The Corinth Rift in Central Greece has been studied extensively during the past decades, as it is one of the most seismically active regions in Europe. It is characterized by normal faulting and extension rates between 6 and 15 mm yr⁻¹ in an approximately N10E° direction. On 2013 May 21, an earthquake swarm was initiated with a series of small events 4 km southeast of Aigion city. In the next days, the seismic activity became more intense, with outbursts of several stronger events of magnitude between 3.3 and 3.7. The seismicity migrated towards the east during June, followed by a sudden activation of the western part of the swarm on July 15th. More than 1500 events have been detected and manually analysed during the period between 2013 May 21 and August 31, using over 15 local stations in epicentral distances up to 30 km and a local velocity model determined by an error minimization method. Waveform similarity-based analysis was performed, revealing several distinct multiplets within the earthquake swarm. High-resolution relocation was applied using the double-difference algorithm HypoDD, incorporating both catalogue and cross-correlation differential traveltime data, which managed to separate the initial seismic cloud into several smaller, densely concentrated spatial clusters of strongly correlated events. Focal mechanism solutions for over 170 events were determined using P-wave first motion polarities, while regional waveform modelling was applied for the calculation of moment tensors for the 18 largest events of the sequence. Selected events belonging to common spatial groups were considered for the calculation of composite mechanisms to characterize different parts of the swarm. The solutions are mainly in agreement with the regional NNE–SSW extension, representing typical normal faulting on 30–50° north-dipping planes, while a few exhibit slip in an NNE–SSW direction, on a roughly subhorizontal plane. Moment magnitudes were calculated by spectral analysis of S waves, yielding b-values between 1.1 and 1.2 in their frequency–magnitude distribution. The seismic moment release history indicates swarm-like activity during the first phase, which could have acted as a preparatory stage for the second phase (after 12 July) that presented a more typical main-shock–aftershock behaviour. The spatiotemporal analysis reveals that the swarm has occurred in a volume that is likely related with the extension at depth of the NNE-dipping Pirgaki normal fault, outcropping ∼8 km to the south. The slow velocity of eastward migration of the epicentres during June implies triggering by fluids. The situation appears different in the second phase of the sequence, which was probably triggered by a build-up of stress during the first one. The relatively deep hypocentres of the 2013 swarm, compared to the shallower seismic layer within the rift, and their coincidence with the steeply dipping Pirgaki fault, favour an immature rift detachment model. Previous results from instrumental data indicate that approximately the same region had been activated during July–August 1991. The availability of the dense permanent seismological network data thus allowed for a detailed analysis of this crisis, a better understanding of its mechanical context and of the earlier events.

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 efforts 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 first 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.