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.

A dense microseismic network was installed in Northwestern Greece for a period of eleven months. A total of 1368 events were recorded and located using a 1D model. These events were also used to derive a 3D velocity model for the area. This work presents results from further processing of the data using (a) simple location method of events in a 1D medium through Hypo71 standard procedure; (b) location via the probabilistic, non-linear earthquake location method in 3D medium; (c) relocation of the events using the Double - Difference method in 1D medium; and (d) the same relocation  procedure  invoking  3D  medium.  The  application  of  different  location methodologies results in slightly different locations, which are evaluated using as criterion the compactness of hypocenter distribution. The three point method was used in order to derive linear characteristics from the hypocenter distribution and the final results were compared against the focal mechanisms of the events as computed using the polarity method and the 3D velocity model. The combination of accurately computed hypocenters and focal mechanisms provides important information for the seismotectonics of Epirus

Following the installation of a temporary seismological network in western Greece north of the Gulf of Patras, we determined the quality of the sites of each of the 10 stations in the network. For this, we used the horizontal-to-vertical spectral ratio (HVSR) method and calculated an average curve over randomly selected days between 0 Hz and 10 Hz. The daily HVSR curve is determined by the HVSR 12-hour calculation (one hour every two) without distinction between seismic ambient noise and earthquake signal. The HVSR curves obtained can be classified in three categories: flat curves without amplification, curves with a amplification peaks covering a large frequency range, and curves with one or more narrow peaks. In this third category C3, one station has one peak, two have two, and one has three. On the contrary of what it is commonly assumed, the amplitudes and the resonance frequencies of these narrow peaks are not stable over time in C3. We determined the maximum of the amplitude of each peak with the corresponding central frequency for each day during 2.5 years. Except for the station with three peaks, which finally appears stable within the uncertainties, the principal peak exhibits a seasonal variation, with a maximum in winter and a minimum in summer, the observations being more dispersed during winter. The second peak, when it exists, varies in the same way except at one station where it varies oppositely. These variations are clearly correlated with the loading and unloading cycle of the underlying aquifers as shown by the comparison with water level and yield measurements from wells located close to the stations. Moreover, they are also correlated with the vertical surface displacements observed at continuously recording GPS stations. The dispersion of the observed maximum amplitude in winter is probably related to the rainfall and the soil moisture modifying the S-wave velocity as revealed by other studies. From this study, we would like to emphasize that the use the HVSR method to constrain the S-wave velocity and the thickness of the sediment layer over the bedrock in the basin, has to be done with caution. Upon further confirmation of its robustness, the HVSR methodology presented here could be a good and easy-to-use tool for a qualitative survey of the aquifer backdrop and its seasonal behaviour, and of the soil moisture conditions.

Climate change constitutes a serious global challenge with consequences that are directly affecting infrastructure. Thus, there is a great need to develop reliable cost-effective systems, which integrate remote sensing data, in situ measurements and advanced methods for infrastructure monitoring. In this framework, the European Union and the Hellenic government are financially supporting a R&D project, named "PROIΟΝ". The purpose of the project is the development of a platform for the continuous monitoring of high-priority infrastructures (public infrastructure, dams, bridges, etc.) which are located in a particularly active area in terms of tectonics and seismicity. Monitoring is based on the combination of instrumental and remote sensing measurements along with fuzzy logic networks methods and machine learning algorithms in order to generate an innovative decision-making and decision-support tool. Specifically, measurements derived from three-axis accelerometers, Global Navigation Satellite System (GNSS) receivers and Persistent Scatterer Interferometry are imported into the platform. The measurements will be validated using high-accuracy reference representations arising from data acquired by Terrestrial Laser Scanning (TLS) surveys and Unmanned Aerial Vehicles (UAV) campaigns and subsequently, deformation maps will be generated. Intelligent data analysis methods will contribute to making decisions on the current as well as the future state of the infrastructure. At this initial stage of the project, the proposed monitoring system is described in detail.