The seismogenic structure of March 2021 Tyrnavos (central Greece) doublet (M_w 6.3 and M_w 6.0), constrained by aftershock locations and geodetic data

The Northern Thessaly Basin in central Greece ranks amongst the most well pronounced extensional (graben) basins in the backarc Aegean Sea region, with well-mapped faults having an & SIM;E-W orientation, compatible with the ongoing predominant & SIM;N-S extension. The southern margin of the basin is bounded by major faults associated with strong (M6 to M7) earthquakes, whereas along its northern margin, strong events are more scarce, in the documented catalogues. Along this northern margin, a weak, albeit persisting foreshock activity, culminated within 3 d, to an M-w 6.3 earthquake on 3 March 2021 associated with a 15-km-long NE dipping fault segment. It was followed the next day, by the second M-w 6.0 main shock associated with a 13-km-long NE dipping fault segment and 9 d later by an M-w 5.5 earthquake associated with an 8-km-long SW dipping fault segment, with its aligned epicentres, showcasing the cascade type activation of adjacent fault segments. The sequence, evolved to be very productive, with aftershocks extending & SIM;50 km along a & SIM;NW-SE trending narrow seismic zone. All events indicate pure normal faulting, with an NNE-SSW oriented extensional axis, oblique to our previous consensus of the prevalence of & SIM;N-S extension. This observation documents that inherited fault fabric can be reactivated within the modern tectonic stress field. We use high-quality seismological data, alongside Interferometric Synthetic Aperture Radar (InSAR) methodology and Global Navigation Satellite System (GNSS) data, to study the temporal and spatial evolution of the sequence, and to provide inferred kinematic models that describe the complexity of the seismic process, in terms of heterogeneous slip distribution, activated fault planes, fault geometry and displacement field. Cross-sections show that the activity defines the crustal seismogenic layer at depths between 5 and 10 km, associated with low-angle fault segments dipping to the NE. Other faults, both antithetic and secondary ones, appear active and accommodated aftershocks clusters. Using our preferred finite fault source model, we calculated the changes of Coulomb failure stress on the neighbouring faults.