We examine the 6 February 2023 Türkiye–Syria earthquakes using an extensive SAR dataset, addressing some limitations of previous studies. Large surface displacements caused significant loss of coherence in the Sentinel-1 Differential SAR Interferometry (DInSAR) results and prior analyses using Pixel Offset Tracking (POT) were limited by the poor azimuthal resolution of the available Sentinel-1 and ALOS-2 SAR images. For the first time, we present high-azimuth-resolution displacement measurements obtained thanks to the SAOCOM-1 sensors. Azimuth information is important considering that the main movements occurred were horizontal and, in some areas, with an important N-S component. By exploiting multi-frequency Sentinel-1, ALOS-2, and SAOCOM-1 SAR data and applying the DInSAR and POT techniques, where appropriate, we derived a detailed displacement field and retrieved an elaborated fault model comprising 22 segments; this model accurately characterizes the geometry and kinematics of the two main faults. Maximum slip reaches ∼15 m for both faults, and the total seismic moment corresponds to Mw 7.9. A finite-fault ShakeMap generated from this source model shows improved agreement with near-field ground motions relative to point-source formulations, while on-fault static stress changes identify low-slip areas that remained unbroken during rupture. The three-dimensional displacement field reveals a broad uplifted region and opposing horizontal motions between the two main ruptures, indicating distributed deformation within an interfault block that accommodates part of the Arabia–Anatolia convergence. This off-fault deformation has not been documented previously for the 2023 sequence and provides new constraints on strain partitioning and future seismic hazard along the East Anatolian Fault Zone.
In 2023, seismic activity of considerable magnitude occurred along the Türkiye-Syria border, characterised by an Mw 7.8 earthquake on the 6th of February and was followed by an Mw 7.5 event, nine hours later. These earthquakes, which are the strongest recorded in recent years, resulted in over 50,000 casualties and are related with the activity of the East Anatolian Fault Zone —a 600 km-long plate boundary where the Arabian and Anatolian plates meet. To analyse these seismic events, we leveraged data from diverse satellites, including SAOCOM-1, Sentinel-1, and ALOS-2. Employing InSAR techniques, such as conventional interferometry and Pixel Offset tracking, we assessed surface deformations caused by the events. The high-resolution Synthetic Aperture Radar displacement results underwent non-linear and linear inversions, enabling the creation of detailed variable slip fault models. A meticulous multiscale sampling approach was applied, that facilitated a comprehensive examination of the tectonic structures triggering these events. The fault zone exhibited a pronounced left-lateral strike-slip character, with components of dip-slip movements observed in specific segments. Additionally, we capitalised the detailed slip models, to estimate the distribution of the intensity of ground motions in the affected region.
In 2023, seismic activity of considerable magnitude occurred along the Türkiye-Syria border, characterised by an Mw 7.8 earthquake on the 6th of February and was followed by an Mw 7.5 event, nine hours later. These earthquakes, which are the strongest recorded in recent years, resulted in over 50,000 casualties and are related with the activity of the East Anatolian Fault Zone —a 600 km-long plate boundary where the Arabian and Anatolian plates meet. To analyse these seismic events, we leveraged data from diverse satellites, including SAOCOM-1, Sentinel-1, and ALOS-2. Employing InSAR techniques, such as conventional interferometry and Pixel Offset tracking, we assessed surface deformations caused by the events. The high-resolution Synthetic Aperture Radar displacement results underwent non-linear and linear inversions, enabling the creation of detailed variable slip fault models. A meticulous multiscale sampling approach was applied, that facilitated a comprehensive examination of the tectonic structures triggering these events. The fault zone exhibited a pronounced left-lateral strike-slip character, with components of dip-slip movements observed in specific segments. Additionally, we capitalised the detailed slip models, to estimate the distribution of the intensity of ground motions in the affected region.
ABSTRACT In March 2021, a series of three moderate events with moment magnitudes Mw 6.3, 6.0, and 5.5 occurred within a span of 10 days in northern Thessaly, Greece, resulting in a stop–start pattern. The moment tensors (MTs) obtained from the events suggested normal faulting along distinct yet adjacent southeast–northwest-trending faults, indicating a sequential triggering process. We applied two methods to investigate the spatial and temporal characteristics of the source process of the triplet. This approach includes linear slip inversion of regional seismic waveforms and static Global Positioning System offsets, along with backprojection of regional waveforms using the Source-Scanning Algorithm technique. The detailed modeling suggests that three adjacent, subparallel (∼300°), and low-angle, detachment-type structures were responsible for the earthquakes. The stress axes orientation was determined by inverting a dataset of recomputed MTs combined with published solutions of the sequence. The results indicated a clockwise rotation of the local extensional axis by ∼25° compared to the regional ∼north–south extension. This rotation is consistent with the orientation of the modeled seismogenic structures. Inversions for both nodal planes and grid-searching fault geometry and orientation showed that for the first-two strongest events, the activated fault planes dipped to the northeast, whereas an antithetic (southwest dipping) fault was activated during the last weaker event. For the Mw 6.3 event, both the independent analyses indicated an overall rupture duration of ∼10 s, energy bursts above and near the hypocenter, up-dip (toward the southwest) rupture propagation in the early stages of the source process, and a bilateral rupture at later stages.
ABSTRACT We explore a hypothetical zero-latency earthquake early warning (EEW) system in Greece, aiming to provide alerts before warning thresholds of the intensity of ground motion are exceeded. Within the seismotectonic context of Greece, both shallow- and intermediate-depth earthquakes (along the Hellenic subduction zone) are plausible and, thus, examined. Using regionally applicable attenuation relations, we combine and adjust the methodologies of Minson et al. (2018) and Hoshiba (2020) to examine what are the minimum magnitudes required to invoke the warning thresholds at the user site. With simple modeling, we examine how fast an alert can be issued and what is the available warning time when taking into account delays due to finite-fault rupture propagation, alongside other delays. These computations are merged with delays introduced due to the present-day configuration of the Greek national monitoring network (varying spatial density of permanent monitoring stations). This approach serves as a tool to assess the feasibility of an EEW system at specific sites and to redesign the national permanent monitoring network to serve such a system more effectively (we provide results for four sites.). Warning times for on-land crustal earthquakes are found to be shorter, whereas for intermediate-depth earthquakes in Greece an EEW system is feasible (provides warning times of several tens of seconds at large cities, e.g., on Crete Island) even with the current configuration of the national monitoring network, which is quite sparse in the southern part of the country. The current network configuration also provides sufficient early warning (e.g., of the order of 10 s for a warning threshold of 0.05g) at the center of Athens from earthquakes of the eastern Gulf of Corinth—a zone posing elevated hazard in the broader area of the Greek capital. Several additional assumptions and factors affecting the operability of an EEW system in Greece (i.e., source process complexity and uncertainty in attenuation laws) are also discussed.
We carry out a probabilistic seismic hazard assessment (PSHA) for Lesvos Island, in the northeastern Aegean Sea. Being the most populated island in the northern Aegean Sea and hosting the capital of the prefecture, its seismic potential has significant social-economic meaning. For the seismic hazard estimation, the newest version of the R-CRISIS module, which has high efficiency and flexibility in model selection, is used. We incorporate into the calculations eight (8) ground motion prediction equations (GMPEs). The measures used are peak ground acceleration, (PGA), peak ground velocity, (PGV), and spectral acceleration, (SA), at T=0.2 sec representative of the building stock. We calculate hazard curves for selected sites on the island, sampling the southern and northern parts: Mytilene, the capital, the village of Vrisa, Mithymna and Sigri. Hazard maps are also presented in terms of all three intensity measures, for a mean return period of 475 years (or 10% probability of exceedance in 50 years), assuming a Poisson process. Our results are comparable to the predictions of on-going EU hazard models, but higher than the provisions of the Greek Seismic Code. Finally, we perform disaggregation of hazard to depict the relative contribution of different earthquake sources and magnitudes to the results.