REVIEW OF THE DEVASTATING SERIES OF EARTHQUAKES IN TÜRKIYE AND SYRIA IN FEBRUARY 2023. EARTHQUAKE ANALISYS USING SAR SATELLITE DATA: PRELIMINARY RESULTS
https://doi.org/10.53656/nat2023-1.05
Резюме. This article is dedicated to the earthquakes of February 6th, 2023 in Türkiye and Syria. The epicenters are located in seismically active region where the East Anatolian Fault is placed. Past and present seismicity is reviewed and together with the aftershock activity is shown as simplified map. The tectonic and geological features of the affected area are discussed in brief. Preliminary assessments of the consequences and damages are presented, according to the data we known up to now. In addition, the main shock generated tsunami waves that were registered by the mareograph network. Using SAR satellite data, we made an estimation of the deformations in the Earth’s crust after the earthquakes. The preliminary results show displacements above 2.5 m.
Ключови думи: earthquake; Türkiye; Syria; SAR satellite data
1. Introduction
Devastating earthquake struck southern and central Türkiye, as well as northern and western Syria on February 6 th, 2023 at 01:17 UTC. The estimated moment magnitude is M w7.8 and the maximum Mercalli intensity of XI, which refers to extreme scale level. The earthquake hypocenter was at a depth of 17.9 km according to USGS 1) and 5 km according to KOERI 2). The epicenter is located 34 km west of Gaziantep and close to the border with Syria. Another strong earthquake (M w7.5) occurred nine hours later near the province Kahramanmaraş. According to GEOSCOPE it had a depth of 13 km3), according to KOERI of 5 km 4). It was followed by two additional mb 6.0 aftershocks. Both of the earthquake epicenters are shown in Figure 1.
The main earthquake is the strongest earthquake in Türkiye since the Erzincan earthquake of 1939, with the same magnitude making it the second strongest in the country‘s history after the North Anatolia earthquake of 1668. The earthquake is the deadliest in Türkiye since the Cilicia earthquake of 1268 and in Syria since the Aleppo earthquake of 1822. It has been felt as far away as Egypt, Israel, Palestine, Lebanon, Cyprus, and Türkiye‘s Black Sea coast.
Three days prior the mainshock, a foreshock of magnitude M4.4 occurred. The earthquakes were followed by more than 3 858 aftershocks 5). The seismic sequence was the result of a shallow strike-slip fault. By February 19, 2023, more than 47 400 deaths had been confirmed: more than 41 000 in Türkiye and more than 6,400 in Syria 6). A severe winter storm hampered rescue efforts, dropping snow on the ruins and lowering temperatures. Because of the freezing temperatures in the area, survivors, especially those trapped under the rubble, were at great risk of hypothermia. Overall, the earthquakes are estimated to have caused $84.1 billion in damage 7), making them the fourth deadliest earthquakes since records began. It is the deadliest natural disaster in Türkiye’s modern history.
For fully understanding and analyzing the consequences of these earthquakes, the presence of satellite data is crucial. The Copernicus Programme was established with the vision on one hand to enable continuous monitoring of the Earth environment and ecosystem for creating added value products, and on the other, to provide information for quick response in case of natural disasters8). It relays on both, in-situ and satellite data to create its products and services. The satellite imagery with its good revisit time and spatial resolution can be implemented in planning and navigating hazard response, but also to evaluate the damage in timely manner and support humanitarian aid.
Radar images acquired by Synthetic Aperture Radar (SAR) can be utilized in estimating the ground motion and thus observing the land cover deformation. The Sentinel-1 satellite is equipped with a SAR antenna, which operates at C-band (approx. 5.5 cm wavelength and 5.405 GHz) and in the Inteferometric Wide (IW) swath mode its spatial resolution is 20 m in azimuth and 5 m in range9).
Information about changes in the terrain topography over an area of interest can be gained from the Interferometric SAR (InSAR) signal. This method makes use of the phase difference between two observations of the same area, obtained from different sensor position. After coregistration of the two images the terrain deformation can be mapped 10).
In this context the InSAR technology can be very important in understanding and evaluating the impact of the devastating earthquake events in southern and central Türkiye and northwestern Syria from 6 th February 2023. In this study we present preliminary results for the deformation of the Earth’s crust based on InSAR image analysis.
Figure 1. Location of the earthquake epicenters from February 6th, 2023. The main shock M7.8 (left panel)11) and the larger aftershock M7.5 (right panel)12)
2. Geology, Seismicity and Tectonic Features
The East Anatolian Fault is a major strike-slip fault zone that runs from eastern to south-central Türkiye. The EAF is a 700 km long left-lateral fault between the Anatolian and Arabian plates and forms the transform tectonic boundary between the Anatolian Plate and the northward-moving Arabian Plate. The difference in the relative motions of the two plates is evident in the left lateral motion along the fault. The East and North Anatolian faults together accommodate the westward motion of the Anatolian plate, which is pushed out by the ongoing collision with the Eurasian plate. Figure 2. Represents the relative motion of the tectonic plates in the region. Several geological and geodetic studies have estimated the slip rates on segments along the EAF and its subparallel segments (McClusky et al. 2000; Reilinger et al. 2006). Geological and geomorphological studies show that the slip rate systematically decreases from ∼10 mm per year-1 at Karlıova to ∼4 mm per year-1 near Türkoğlu. Further west, the slip rate decreases to ∼2.5 mm per year-1 on the main fault and to ∼1 mm per year-1 on subparallel faults. Slip rates estimated from Holocene dislocations are compatible with geodetic fault slip rates (Güvercin 2022).
From 1939 to 1999, a series of earthquakes occurred westward along the North Anatolian Fault. However, since 1998, there has been a series on the East Anatolian Fault. This began with the 1998 Adana-Ceyhan earthquake and includes the 2003 Bingöl earthquake, the 2010 Elâzığ earthquake, the 2020 Elâzığ earthquake, and the 2023 Türkiye-Syria earthquake. Türkiye has been hit by 21 earthquakes of magnitude 7 or higher since 1900 (see Figure 3). Aleppo, in Syria, was devastated several times historically by large earthquakes, though the precise locations and magnitudes of these earthquakes can only be estimated. Aleppo was struck by an estimated magnitude 7.1 earthquake in 1138 and an estimated magnitude 7.0 earthquake in 1822. Fatality estimates of the 1822 earthquake were 20 000 – 60 000.
Figure 4 depicts the seismicity13) in the period between 1st and 13th February 2023. The earthquakes have magnitude equal to or larger than 3.5. The size of the circles is scaled according the magnitude of the earthquakes. The black lines show the seismic faults according the European Database of Seismogenic Faults (Basili et al. 2013).
Figure 2. Relative motion of Anatolian and Arabian plates and major faults14)
Figure 3. Strong earthquakes in the past 100 years in Türkiye15)
Figure 4. Foreshock and aftershock activity in the period 1 Feb. 2023 – 13 Feb. 2023
3. Earthquake Damages
One of the most powerful earthquakes ever recorded in Türkiye killed thousands of people and devastated a densely populated region of southern Türkiye and northern Syria. The area was particularly vulnerable to a large earthquake. Older concrete-framed buildings are common, and thousands were destroyed. The infrastructure in northern Syria was already very fragile after years of airstrikes and bombings during the country’s civil war. While the full extent of the damage is not yet known, collapsed homes, hotels and other buildings have been spotted in dozens of towns near the epicenter on both sides of the border. More than 5 600 buildings were destroyed across Türkiye. In Syria, the earthquakes hit the country’s war-ravaged north. More than a thousand were dead and thousands more were injured, according to the Syrian Health Ministry. As of 9 th February 2023, the casualties in different Turkish provinces are shown in Table 1. As of 17 th February 2023, the total number of live affected and some statistics are collected in Table 2. Figure 5 and Figure 6 depicts destroyed buildings in Türkiye and in Syria. Damage extension according the level of the shake intensity of the main shock if presented in Figure 7.
Table 1. Casualties by Turkish province (as of 9th February 2023)16), 17)
Table 2. Lives affected and statistics as of 17th February 202318)
Figure 5. Destroyed building in Diyarbakir, Türkiye19)
Figure 6. The town of Jandaris, near Afrin, Syria20)
Figure 7. Damage extension in Türkiye and Syria according to the shake intensity of the main earthquake21)
4. Earthquake Analysis Using SAR Satellite Data: Preliminary Results
The satellite data used in our analysis is from the space segment of the Copernicus Programme. The images are freely available for download at the Copernicus Open Access Hub after short registration22). The pre-processing steps which need to be applied are described in Braun & Veci 2021. The images are processed with the software SNAP (Sentinel Applications Toolbox)23). All satellite images are selected to cover as short period as possible before and after a given earthquake. A condition for using the SAR data is that the earthquake has a magnitude greater than 5.5 or 6; presence of vegetation or complex topography may lead to bad results because of the strong interference. The interferometric maps are represented by interferometric fringes (colour-coded from blue to purple) each of them equals to 28 mm, half the wavelength of the instrument. After unwrapping the phase, the deformation map is obtained. It is important to note that the deformation is in Line-of-Sight (LoS) direction of the satellite antenna and not vertical displacement. Both ascending and descending orbits over the given region are used to calculate the vertical displacement, and after calculating the total deformation along the direction of the satellite antenna and using information about the angle of incidence between the radiation and the perpendicular to the base, we can calculate the vertical displacement in the up-down direction and the horizontal displacement in the east-west direction.
Table 3. Main parameters of the selected products
Figure 8. Interferogram formed by Orbit 14 and Orbit 21
Figure 9. Coherence (left panel). Intensity (central panel). DEM (right panel) Orbit 14: 28.01.2023-9.02.2023
Figure 10. Deformation map using images from Orbit 21. The earthquake epicenter is denoted as red star
In total we selected, downloaded and processed 3 pairs of available images. The results are exported as .kml files and visualized via Google Earth Pro software24. The final interferogram composed by different images is presented in Figure 8.
Coherence is used to verify the legitimacy of the derived phase data. When creating the interferogram an important part is to check the coherence of the master and slave images. Usually, the threshold is 0.3. In the case of this huge earthquake for some of the pairs we get a minimum coherence of 0.8. Figure 9 shows the coherence, the intensity and the digital elevation model (DEM) from Orbit 14, master image of 28 January 2023 and slave image of 9 February 2023.
An example of deformation map is presented in Figure 10. We can see that at some areas the displacements reach values greater than 2.5 m. These are preliminary results, thus when more images from different orbits are available, we will combine all the data and create more complex image of the different scenes.
5. Conclusions
An overview of the destructive series of earthquakes in Türkiye and Syria was described in this paper. The number of casualties is not final as the rescuers are still struggling to find many of the people under the collapsed buildings.
The analysis of earthquakes using the data from SAR leads to a better understanding of the generation mechanisms and complements the geophysical and geological information. Using SAR satellite data, we analyzed a series of images and produced interferometric maps. The differences in the values from the different orbits confirm the need, when studying deformations obtained after earthquakes, to prepare interferogram images from both types of orbits in order to get a better idea of them and their subsequent comparison with field measurements or results from other sources.
ACKNOWLEDGEMENTS
This study is conducted within project funded by the National Scientific Program “Young Scientists and Postdoc – 2”, 2022 by Ministry of Education and Science of Republic of Bulgaria.
NOTES
1 United States Geological Survey (USGS): https://earthquake.usgs.gov/ earthquakes/eventpage/us6000jllz/executive
2. Kandilli Observatory KOERI: http://www.koeri.boun.edu.tr/new/en
3. GEOSCOPE: http://geoscope.ipgp.fr/index.php/en/catalog/earthquakedescription?seis=us6000jlqa
4. http://www.koeri.boun.edu.tr/sismo/2/06-subat-2023-ml7-5-ekinozukahramanmaras-depremi/
5. https://www.milliyet.com.tr/galeri/son-depremler-listesi-kandilli-afad-az-oncedeprem-mi-oldu-kac-siddetinde-deprem-bolgesinde-artci-sarsintilar-devamediyor-6902271/1
6. https://www.milliyet.com.tr/gundem/afad-son-durumu-acikladi-depremde-cankaybi-41-bin-20-6906107
7. https://www.bloomberght.com/turkonfed-raporu-depremlerin-mali-hasarinin84-1-milyar-dolar-olmasi-bekleniyor-2325161
8. https://www.copernicus.eu/en
9. https://sentinels.copernicus.eu/web/sentinel/missions/sentinel-1
10. https://sentinels.copernicus.eu/web/sentinel/user-guides/sentinel-1-sar/productoverview/interferometry
11. https://static1.emsc.eu/Images/EVID/121/1218/1218444/1218444.regional.jpg
12. https://static1.emsc.eu/Images/EVID/121/1218/1218771/1218771.regional.jpg
13. http://isc.ac.uk/iscbulletin/search/catalogue/
14.https://en.wikipedia.org/wiki/2023_Turkey%E2%80%93Syria_earthquake#/ media/File:Anatolian_Plate_Vectoral.svg
15. https://www.theguardian.com/world/2023/feb/20/thousands-dead-millionsdisplaced-the-earthquake-fallout-in-turkey-and-syria
16. https://www.sabah.com.tr/galeri/yasam/gaziantep-depremi-olu-ve-yaralisayisi-gaziantepte-kac-bina-yikildi-olu-sayisi-kac-oldu-kac-kisi-yaralandi/7
17. https://www.aa.com.tr/tr/asrin-felaketi/saglik-bakani-koca-10-ilde-17-bin929u-hekim-olmak-uzere-143-bin-829-personelimiz-hizmet-veriyor/2813793#
18. https://www.reuters.com/, accessed on 17.02.2023
19. https://canberraweekly.com.au/australia-to-give-turkey-10-million-afterhorrific-earthquake/
20. Photo credit: Rami al-Sayed/AFP/Getty https://www.nature.com/articles/ d41586-023-00364-y
21. https://www.nytimes.com/interactive/2023/02/06/world/turkey-earthquakedamage.html
22. https://scihub.copernicus.eu/
23. https://step.esa.int/main/
24. https://www.google.com/earth/versions/#earth-pro
REFERENCES
BASILI R., ET AL., 2013. The European Database of Seismogenic Faults (EDSF) compiled in the framework of the Project SHARE. Istituto Nazionale di Geofisica e Vulcanologia (INGV). https://doi.org/10.6092/ ingv.it-share-edsf.
BRAUN, A. & VECI, L., 2021. TOPS Interferometry Tutorial. Sentinel-1 Toolbox.
GÜVERCIN, S. E., ET AL., 2022. Active seismotectonics of the East Anatolian Fault. Geophysical Journal International. 230(1), 50 – 69.
MCCLUSKY, S. ET AL., 2000. Global Positioning System constraints on plate kinematics and dynamics in the eastern Mediterranean and Caucasus. J. Geophys. Res. 105, 5695 – 5719. 10.1029/1999JB900351.
REILINGER, R. ET AL., 2006. GPS constraints on continental deformation in the Africa-Arabia-Eurasia continental collision zone and implications for the dynamics of plate interactions. J. Geophys. Res. 111, doi:10.1029/2005JB004051.