Here we present a joint analysis of the geodetic, seismological and geological data of the March 2021 Northern Thessaly seismic sequence, that were gathered and processed as of April 30, 2021. First, we relocated seismicity data from regional and local networks and inferred the dip-direction (NE) and dip-angle (38°) of the March 3, 2021 rupture plane. Furthermore, we used ascending and descending SAR images acquired by the Sentinel-1 satellites to map the co-seismic displacement field. Our results indicate that the March 3, 2021 Mw=6.3 rupture occurred on a NE-dipping, 39° normal fault located between the villages Zarko (Trikala) and Damasi (Larissa). The event of March 4, 2021 occurred northwest of Damasi, along a fault oriented WNW-ESE and produced less deformation than the event of the previous day. The third event occurred on March 12, 2021 along a south-dipping normal fault. We computed 22 focal mechanisms of aftershocks with M≥4.0 using P-wave first motion polarities. Nearly all focal mechanisms exhibit normal kinematics or have a dominant normal dip-slip component. The use of InSAR was crucial to differentiate the ground deformation between the ruptures. The majority of deformation occurs in the vertical component, with a maximum of 0.39 m of subsidence over the Mw=6.3 rupture plane, south and west of Damasi. A total amount of 0.3 m horizontal displacement (E-W) was measured. We also used GNSS data (at 30-s sampling interval) from twelve permanent stations near the epicentres to obtain 3D seismic offsets of station positions. Only the first event produces significant displacement at the GNSS stations (as predicted by the fault models, themselves very well constrained by InSAR). We calculated several post-seismic interferograms, yet we have observed that there is almost no post-seismic deformation, except in the footwall area (Zarkos mountain). This post-seismic deformation is below the 7 mm level (quarter of a fringe) in the near field and below the 1 mm level at the GNSS sites. The cascading activation of the three events in a SE to NW direction points to a pattern of domino-style earthquakes, along neighbouring fault segments. The kinematics of the ruptures point to a counter-clockwise change in the extension direction of the upper crust (from NE-SW near Damasi to N-S towards northwest, near Verdikoussa).

Physical laws governing friction on shallow faults in the Earth and spatial heterogeneity of parameters are critical to our understanding of earthquake physics and the assessment of earthquake hazards. Here we use a laboratory-derived fault-friction law and high-quality strong-motion seismic recordings of the 2020 Elazığ earthquake, Turkey, to reveal the complex rupture dynamics. We discover an initial Mw 5.8 rupture stage and explain how cascading behavior of the event, involving at least three episodes, each of M > 6, caused it to evolve into a large earthquake, contrarily to other M5+ events on this part of the East Anatolian Fault. Although the dynamic stress transfer during the rupture did not overcome the strength of the uppermost ~5 kilometers, surface ruptures during future earthquakes cannot be ruled out. We foresee that future, routine dynamic inversions will improve understanding of earthquake rupture parameters, an essential component of modern, physics-based earthquake hazard assessment.

The National Observatory of Athens data center for the European Integrated Data Archive (EIDA@NOA) is the national and regional node that supports International Federation of Digital Seismograph Networks and related webservices for seismic waveform data coming from the southeastern Mediterranean and the Balkans. At present, it serves data from eight permanent broadband and strong-motion networks from Greece and Cyprus, individual stations from the Balkans, temporary networks and aftershock deployments, and earthquake engineering experimental facilities. EIDA@NOA provides open and unlimited access from redundant node end points, intended mainly for research purposes (see Data and Resources). Analysis and quality control of the complete seismic data archive is performed initially by calculating waveform metrics and data availability. Seismic ambient noise metrics are estimated based on power spectral densities, and an assessment of each station’s statistical mode is achieved within each network and across networks. Moreover, the minimum ambient noise level expected for strong-motion installations is defined. Sensor orientation is estimated using surface-wave polarization methods to detect stations with misalignment on particular epochs. A single data center that hosts the complete seismic data archives with their respective metadata from networks covering similar geographical areas allows coordination between network operators and facilitates the adhesion to widely used best practices regarding station installation, data curation, and metadata definition. The overall achievement is harmonization among all contributing networks and a wider usage of all data archives, ultimately strengthening seismological research efforts in the region.

No abstract

We studied the kinematic behaviour of active landslides at several localities in the area of Panachaikon Mountain, Achaia (Peloponnese, Greece) using Sentinel (C-band) InSAR time series analysis. We processed LiCSAR interferograms using the SBAS tool, and we obtained average displacement maps for the period 2016–2021. We found that the maximum displacement rate of each landslide is located at about the center of it. The average E-W velocity of the Krini landslide is ~3 cm/year (toward the east) and 0.6 cm/year downward. The line-of-sight (LOS) velocity of the landslide (descending orbit) compares well to a co-located GNSS station within (±) 3 mm/yr. Our results also suggest a correlation between rainfall and landslide motion. For the Krini landslide, a cross-correlation analysis of our data suggests that the mean time lag was 13.5 days between the maximum seasonal rainfall and the change in the LOS displacement rate. We also found that the amount of total seasonal rainfall controls the increase in the displacement rate, as 40–550% changes in the displacement rate of the Krini landslide were detected, following to a seasonal maximum of rainfall values at the nearby meteorological station of Kato Vlassia. According to our results, it seems that large part of this mountainous region of Achaia suffers from slope instability that is manifested in various degrees of ground displacement greatly affecting its morphological features and inhabited areas.