Reconstruction of Co-Seismic and Post-Seismic Processes for the February 6, 2023 Earthquake in Turkey from Data of Satellite SAR Interferometry
Authors: Mikhailov V. O. , Babayants I. P. , Volkova M. S. , Timoshkina E. P. , Smirnov V. B. , Tikhotsky S. A. 
Published: April 2024

Abstract—Using different methods for processing SAR images from the Sentinel-1A satellite, the displacement fields were determined in the region of the East Anatolian Fault Zone (EAFZ) and the Sürgü-Çardak faults, as well as a small fault on the continuation of the East Hatay fault zone, which rupture initiated a series of catastrophic earthquakes in Turkey on February 06, 2023. DInSAR and offset methods were applied. The most detailed data on the displacements were obtained by the offset method using images from the descending orbit. When constructing the model from the available SAR data, the data with the maximum signal-to-noise ratio were selected. For the northern part of the region, above 37.4° N, the range displacements obtained by the offset method from a descending orbit were used. South of parallel 37.0° N, we used azimuth displacements from the same descending orbit. The model of the seismic rupture was constructed on the basis of solution of (Pollitz, 1996) of the problem of deformations at the surface of a layered spherical Earth caused by along dip and strike displacements on a rectangular cut located inside the planet. Pollitz (1996) demonstrated that ignoring the spherical layering of the planet leads to errors up to 20%, with the largest errors occurring in the presence of a large strike-slip component. Ignoring sphericity also introduces an error when using the solution in the framework of the idealization of an elastic homogeneous half-space with a flat free boundary (Okada, 1985) which was used when constructing USGS and (Barbot et al., 2023) models. In our model the surfaces of seismic rupture are approximated by 19 rectangular elements along the strike, divided into three levels along the dip. Another element approximated a rupture along the extension of the East Hatay Fault Zone. As in the models of other authors (USGS; Barbot et al., 2023), in our model in the southern part of the EAFZ, the displacements increase from south to north, and are mainly concentrated in the upper part of the Earth’s crust to a depth of 10 km. At the southern end of this rupture, displacements in our model with an amplitude of up to 2 m are obtained at the lower levels of the model, and at its upper level, the displacements were only 0.11 m, and in this area on February 20, 2023 an earthquake of magnitude 6.3 occurred with a hypocentre depth of 11.5 km. The main displacements on the EAFZ are determined on its central segment. Here, the displacements go to a greater depth, their value reaches 10.2 m. On the Sürgü-Çardak fault, significant displacements occurred down to a depth of 20 km; displacements exceeded 9.8 m. In our model, at the northeast end of the seismic rupture along the EAFZ, a displacement area of 6.8 m overlaps with the southwest end of the seismic rupture model of the Doğanyol-Sivrice earthquake of January 24, 2020 with Mw 6.7, published on the USGS website. Therefore, our model does not confirm the hypothesis of the presence of a seismic gap here, which, according to (Barbot et al., 2023), is a zone of a possible nearest earthquake.