Eastern Turkey Earthquake Research Articles

From Interseismic Deformation With Near‐Repeating Earthquakes to Co‐Seismic Rupture: A Unified View of the 2020 Mw6.8 Sivrice (Elazığ) Eastern Turkey Earthquake

AbstractThe East Anatolian Fault (EAF) is a left‐lateral transform fault accommodating the relative motion between the Anatolian and Arabian plates. On January 24, 2020, Mw6.8 Sivrice (Elazığ) earthquake is the largest event that occurred along the EAF since the nineteenth century. The earthquake provides a unique opportunity to capture a critical stage of the seismic cycle from the interseismic deformation to co‐seismic rupture. In this study, we examine the relationship between the interseismic fault activity and co‐seismic behavior of the earthquake. A kinematic model of the earthquake obtained from strong‐motion, GNSS and teleseismic waveforms along with static displacements from GNSS and InSAR data shows that the mainshock ruptured only 45 km of the 95 km long Sivrice‐Pütürge segment. Rupture initiated adjacent to the interseismically weakly coupled northeastern section and propagated unilaterally toward southwest with a rupture velocity of ∼2.5 km/s, stopping ∼30 km before the southwestern segment boundary. The earthquake did not generate any surface offsets. We identified 4 long‐term near‐repeating earthquake clusters beneath the highest co‐seismic slip zone adjacent to the northeastern creeping section. The mainshock was dynamically triggered by a M ∼ 5.4 foreshock located within the zone of near‐repeating earthquakes. We suggest that creep along the weakly coupled section of the fault loaded the neighboring locked section leading to the repeating earthquakes below the locked zone. The earthquake partially ruptured a fault segment characterized by high rate of diffuse seismicity, structural complexities and heterogeneous coupling, leading to a source characterized by a complex source time function with relatively slow rupture velocity.

Rupture process of 2011 Mw7.1 Van, Eastern Turkey earthquake from joint inversion of strong-motion, high-rate GPS, teleseismic, and GPS data

We analyzed the rupture process of the 2011 Mw7.1 Van, Eastern Turkey earthquake using teleseismic, strong-motion, 1-Hz GPS waveforms and static GPS displacement measurements. We performed data sensitivity analyses using four different rupture scenarios. Overall, when geodetic and seismic datasets are modeled jointly, slip distribution is well-constrained and rupture velocity can be obtained. The 2011 Van earthquake is a case where none of the available datasets are sufficient to constrain the slip distribution and the rupture kinematics on their own. This study confirms that rather than fitting one dataset perfectly, using multiple datasets jointly leads to a better-constrained slip distribution. The kinematic model obtained from the joint inversion of all the available data shows a 45-km-long bilateral rupture with two sub-events; one larger slip patch propagating up-dip and toward west of the hypocenter and a smaller second slip patch toward the east. The highest slip is to the west of the hypocenter with a peak value of 4.5 m. The slip is confined to the depths of 7.5 to 20 km and the shallower part of the fault remains unbroken. The average rupture velocity is around 3 km/s, close the Rayleigh wave velocity. The rupture is faster with shorter rise times for the larger sub-event toward the west, while the rise times are longer for the smaller sub-event to the east. This difference in seismic behavior might be related to the segmentation of the fault.

Rupture pattern of the Oct 23, 2011 Van-Merkez, Eastern Turkey earthquake

High-frequency rupture process of the Oct 23, 2011 Van-Merkez earthquake is imaged by back-projection method using high-quality teleseismic P wave data from the US Array, and prestack Kirchhoff migration using P wave data from a subarray of global seismic networks. The rupture model with two asperities is confirmed by previous two methods. In low-frequency imaging, a large asperity derived from the migration method corresponds to the second one from the high-frequency P waves. The consistency of the locations of asperities from datasets with different frequency bands indicates that there is possible insignificance of the frequency-dependent feature for the earthquake. The resultant images illustrate the spatial and temporal evolution of the rupture, which mainly propagated WSW over a length of 33 km during the first 18 s, accompanying with bursts of two asperities at 3 and 11–13 s. The rupture direction is confirmed by the S wave corner frequency variations of strong ground accelerations. The rupture fronts are mainly located at the updip of the causative fault. Based on polarities of the P waveforms and focal mechanisms of the mainshock and aftershocks, the failure of these two asperities is determined to have occurred on a reverse fault with a dip angle of 47°. Hence, the rupture pattern of the 2011 Van-Merkez earthquake was dominated by a unilateral rupture toward the west-southwest direction.

Ground motion simulations for the 23 October 2011 Van, Eastern Turkey earthquake using stochastic finite fault approach

The 23 October 2011 Van (Mw 7.1) earthquake that occurred in Eastern Turkey resulted in heavy damage particularly in the city of Van and town of Ercis. This paper presents ground motion simulations of Van earthquake by using stochastic finite fault method (EXSIM, Motazedian and Atkinson in Bull Seismol Soc Am 95:995–1010, 2005; Boore in Bull Seismol Soc Am 99:3202–3216, 2009) that provides a simple and effective tool to generate high frequency strong motion. The input parameters related to source, path, and site effects are calibrated on the basis of minimizing the error functions between simulations and observations both in time and frequency domain. Validated model parameters are used to produce synthetics in regional extent with the aim of understanding the level and distribution of the ground shaking particularly in the near fault region where no recordings are available within the 40 km of the epicenter. This paper evaluates the effect of two different slip models on ground motion intensity measures over the area of interest and addresses the variability in the near fault region associated with the source effect. The synthetics are compared with the corresponding estimations of ground motion prediction equations by Boore and Atkinson (Earthq Spectra 24:99–138, 2008), Akkar and Bommer (Seismol Res Lett 81:195–206, 2010) and Akkar and Cagnan (Bull Seismol Soc Am 100:2978–2995, 2010). Our results indicate that despite the limitation of the method for incorporating the directivity effect and inadequate representation of the soil conditions at the individual stations, a satisfactory match between synthetics and observations are obtained both in time and frequency domain. Spatial distributions of the synthetics in regional level also show reasonable correlation with ground motion prediction equations and damage observations.

Site attenuation and source parameters on the North Anatolian Fault zone, eastern Turkey estimated from the aftershocks of 13 March 1992 Erzincan earthquake

Site attenuation and source characteristics of 45 aftershocks of the 13 March 1992 Erzincan, eastern Turkey earthquake have been determined from SH-wave spectra using a least-squares best-fit method. Although the most of the seismograms were recorded on the ophiolitic rock sites and the average regional attenuation correction of Q(f) = 35 f0.83 was applied, the high fall-off rates γ ranging between 3.0 and 4.5 beyond the corner frequencies of 8–13 Hz were observed on the SH spectra. The site attenuation value, κ, has been calculated from the slope of the high frequency part of the SH spectra. It was found that κ varied in the range of 0.0124–0.0364 s and the average was 0.0246 s with a standard deviation of 0.0047 s. The high fall-off rates of observed spectra have decreased considerably as a result of site attenuation corrections and converged to an average value of 2.3. We concluded that the high fall-off rates of SH wave spectra are mainly controlled by highly deformed Miocene ophiolitic formations which covers a wide area in the North Anatolian Fault zone in the vicinity of Erzincan region. Using the spectra that were corrected for regional and site attenuation, and assuming a Brune's source model; seismic moment, source radius, and stress-drops of the aftershocks were computed. We found that stress-drops for some of the aftershocks in Erzincan area have slightly decreased after removal of site attenuation over SH wave observed spectra.