# Determination of recent tectonic deformations in the vicinity of Adana–Osmaniye–Hatay–Gaziantep triple junction region by half-space modeling

Selim Serhan Yildiz, Ali Özkan, Hasan Hakan Yavaşoğlu, Frédéric Masson, İbrahim Tiryakioğlu, Mehmet Nurullah Alkan, Serdar Bilgi

<figure><img src="image_1.png" /></figure>

## To cite this version:

Selim Serhan Yildiz, Ali Özkan, Hasan Hakan Yavaşoğlu, Frédéric Masson, İbrahim Tiryakioğlu, et al. Determination of recent tectonic deformations in the vicinity of Adana–Osmaniye–Hatay–Gaziantep triple junction region by half-space modeling. Comptes Rendus. Géoscience, 2020, 352 (3), pp.225-234. 10.5802/crgeos.39. hal-03528813

HAL Id: hal-03528813

https://hal.science/hal-03528813

Submitted on 17 Jan 2022

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Comptes Rendus
## Géoscience
### Sciences de la Planète

Selim Serhan Yıldız, Ali Özkan, Hasan Hakan Yavaşoğlu, Frédéric Masson, İbrahim Tiryakioğlu, Mehmet Nurullah Alkan and Serdar Bilgi

**Determination of recent tectonic deformations in the vicinity of Adana–Osmaniye–Hatay–Gaziantep triple junction region by half-space modeling**

Volume 352, issue 3 (2020), p. 223-232

<https://doi.org/10.5802/crgeos.39>

© Académie des sciences, Paris and the authors, 2020.
Some rights reserved.

<figure><img src="image_3.png" /></figure>

This article is licensed under the
CREATIVE COMMONS ATTRIBUTION 4.0 INTERNATIONAL LICENSE.
http://creativecommons.org/licenses/by/4.0/
Original Article — Tectonics, Tectonophysics

# Determination of recent tectonic deformations in the vicinity of Adana–Osmaniye–Hatay–Gaziantep triple junction region by half-space modeling

Selim Serhan Yıldız<sup>a</sup>, Ali Özkan<sup>*, b</sup>, Hasan Hakan Yavaşoğlu<sup>c</sup>, Frédéric Masson<sup>d</sup>, İbrahim Tiryakioğlu<sup>e</sup>, Mehmet Nurullah Alkan<sup>f</sup> and Serdar Bilgi<sup>c</sup>

**Abstract.** Recent earthquakes on the East Anatolian Fault Zone (EAFZ) and its surroundings in the Eastern Mediterranean region reveal a potential seismic hazard. It is thus needed to constrain deformations within tectonically active zones and clarify kinematics of the triple junction near Adana, Osmaniye, Hatay, and Gaziantep provinces. To understand strain accumulations and interseismic patterns of these active zones, one of the most popular geodetic methods was applied. Slip rates derived from Global Positioning System (GPS) observations by half-space modeling lead us to predict possible magnitudes and even probable locations of future destructive earthquakes. In this context, we form a region-specific GPS network consisting of both campaign observation sites and permanent stations to monitor strain accumulations across major active faults and investigate fault kinematics in detail. The campaign GPS data sets acquired in 2009, 2010, 2011, and 2019 were merged with data from different local permanent GPS networks for an integrated analysis of the determination of recent tectonic deformations. The TDEFNODE software was used to construct a kinematic model in an elastic and homogeneous half-space. The modeling approach needs inputs such as block and fault geometries, GPS site velocities, and locking depth assumptions. Our kinematic model indicates that an almost purely sinistral strike-slip rate of 7.5 mm/yr without any significant normal or reverse component is dominant on the main branch of the EAFZ. This estimated slip rate supports a prediction for a

<footnote><sup>a</sup> Department of Geomatics Engineering, Osmaniye Korkut Ata University, Osmaniye, Turkey</footnote>

<footnote><sup>b</sup> Osmaniye Vocational School, Osmaniye Korkut Ata University, Osmaniye, Turkey</footnote>

<footnote><sup>c</sup> Department of Geomatics Engineering, Istanbul Technical University, Istanbul, Turkey</footnote>

<footnote><sup>d</sup> Institut de Physique du Globe, UMR 7516, Université de Strasbourg/EOST CNRS, Strasbourg, France</footnote>

<footnote><sup>e</sup> Department of Geomatics Engineering, Afyon Kocatepe University, Afyonkarahisar, Turkey</footnote>

<footnote><sup>f</sup> Osmancık Vocational School, Hitit University, Corum, Turkey</footnote>

<footnote>*E-mails*: serhan@osmaniye.edu.tr (S. Serhan Yıldız), aliozkan@osmaniye.edu.tr (A. Özkan), yavasoglu@itu.edu.tr (H. Hakan Yavaşoğlu), frederic.masson@unistra.fr (F. Masson), itiryakioglu@gmail.com (I. Tiryakioğlu), nurullahalkan@hitit.edu.tr (M. Nurullah Alkan), bilgi@itu.edu.tr (S. Bilgi)</footnote>

<footnote>* Corresponding author.</footnote>
magnitude margin of 7.2–7.6 for the next probable large earthquake on the Turkoglu–Golbasi segment. The Karatas–Osmaniye Fault has a 3.4 mm/yr dextral strike-slip rate together with a reverse-slip rate of 3.1 mm/yr, which corresponds to a possible magnitude of 6.8–7.1 for an earthquake on the south-west end of the East Anatolian Fault. The Karasu Fault extends from the Turkoglu Triple Junction on the EAFZ in the north to the Hatay Triple Junction in the south. This critical fault segment has a sinistral 4.4–5.4 mm/yr slip rate with a normal component slip rate of 3.0 mm/yr, supporting earthquake prediction with a possible magnitude of 6.8–7.2. As a result, the determination of such kinematic indications in active tectonic zones by GPS observations plays a key role in seismic hazard analysis for the Eastern Mediterranean region.

**Keywords.** Active tectonics, GPS, East Anatolian Fault Zone, Slip rate, Fault locking depth.

Manuscript received 28th May 2020, revised 6th October 2020, accepted 13th October 2020.

# 1. Introduction

The East Anatolian Fault Zone (EAFZ), one of the most important neotectonic and seismogenic fault systems in Turkey, is a sinistral strike-slip active fault zone extending from Karliova in the east of Erzincan to the Hatay region, where it reaches the Dead Sea Fault Zone further south [Arpat and Saroglu, 1972, Westaway and Arger, 1996]. The EAFZ has intense seismic activity depending on the northward relative motions of Arabian and African plates, the westward motion of the Anatolian plate, and continental collision occurring along the Bitlis–Zagros Fold and Thrust Belt [Saroglu and Yilmaz, 1990].

In terms of seismicity, the Karasu Fault (KF; Amanous Fault) segment of the EAFZ and Toprakkale, Karatas, and Yumurtalık Faults are considerable fault segments within the complex structure of the northeast Mediterranean region. These faults were evaluated as active faults on the Active Fault Map of Turkey updated recently [Emre et al., 2018]. On the other hand, the fact that few geodetic studies conducted to investigate the kinematics of these active faults makes it difficult to understand the structure of these faults and analyze earthquake hazards [Aktuğ et al., 2016, Mahmoud et al., 2012, Nocquet, 2012, Meghraoui et al., 2011, Reilinger et al., 2006]. To estimate the potential hazard of a probable earthquake, there is a need to study the slip rates along these active faults and also fault locking parameters using geodetic methods. However, the use of permanent Global Positioning System (GPS) stations together with campaign GPS observations is the most efficient method for this purpose, but the geometry and distinctive features of relevant faults are important criteria for geodetic networks.

This study was conceived with the aim of monitoring active faults with a high probability of gen-

erating earthquakes around the southern end of the EAFZ. Accordingly, we intend to estimate the recent slip rates along these faults and fault locking parameters with high accuracy and spatial resolution by modeling the observed GPS data. For this purpose, campaign observation sites were used along profiles perpendicular to the related faults. A wide geodetic network consisting of campaign observation sites together with permanent GPS stations was observed during the execution of the above proposed study. Raw data obtained from previous studies related to geodesy, paleoseismology, and active tectonics were analyzed together with the latest GPS observations to estimate current strain accumulations for the Eastern Mediterranean region. By analyzing all the obtained data, slip rates along the faults and fault locking depths were studied spatially and temporally in detail.

During the past decades, strain accumulations on the EAFZ have sometimes caused life-threatening earthquakes in eastern Turkey such as Elazığ–Karakocan (Mw 6.1, 2010) and Bingöl (Mw 6.4, 2003) earthquakes. The last destructive earthquake on the EAFZ occurred with an epicenter close to the town of Sivrice in Elazığ province on the Puturge segment of the EAFZ in January 2020. The earthquake solution published by the United States Geological Survey (USGS) is a magnitude of Mw 6.7 and a depth of 10.0 km. However, seismic hazard analysis may prevent substantial loss of life and property associated with such relatively major earthquakes. In this context, the potential moment magnitude for a possible earthquake within the modeled segment of an active fault can be predicted using empirical relationships. These involve either total displacement considering historical and instrumental seismic data [Wells and Coppersmith, 1994] or fault slip rate together with the length of fault rupture [Anderson et al., 1996].
Since the predicted moment magnitude is converted into average displacement with the help of the approach proposed by Wells and Coppersmith [1994], the recurrence interval of an earthquake can be calculated. From this perspective, the determination of tectonic deformations in active fault zones within the Adana, Osmaniye, Hatay, and Gaziantep provinces is very crucial in terms of drawing attention to seismic hazards in the Eastern Mediterranean region.

## 2. Kinematic setting

Since the motions of tectonic plates in the Earth's lithosphere are defined in accordance with plate tectonics theory, it is believed that the Arabian, Anatolian, and Sinai plates interact around the region of the Eastern Mediterranean. According to this hypothesis, the Arabian plate moves to the north at an approximate rate of 20 mm/yr and the Sinai plate converges to a subduction zone along the Cyprus Arc at approximately 10 mm/yr. Furthermore, the westward motion of the Anatolian plate is accommodated by a sinistral strike slip along the EAFZ and a dextral strike slip along the North Anatolian Fault Zone at rates of 10 mm/yr and 20–25 mm/yr, respectively [Reilinger et al., 1997, McClusky et al., 2000, Reilinger et al., 2006, Mahmoud et al., 2012, Aktuğ et al., 2016].

Considering the role of geodetic techniques in active tectonics, it is now possible to estimate the slip rates more precisely along active faults. The slip rates for EAFZ segments are estimated in different geodetic studies as $15 \pm 3$ mm/yr [Reilinger et al., 1997] for the whole EAFZ and $10.0 \pm 0.3$ mm/yr for the Palu–Hazar Lake segment and $9.9 \pm 0.2$ mm/yr for the Hazar Lake–Turkoglu segment along the EAFZ [Reilinger et al., 2006]. The slip rates are calculated to be $9.7 \pm 0.9$ mm/yr for EAFZ and $5.5 \pm 1.5$ mm/yr for the Karatas–Osmaniye Fault (KOF) segment [Bertrand et al., 2006].

More recently, one of the GPS block modeling studies about slip rate estimations has reported the results as $8.8 \pm 0.3$ mm/yr for the EAFZ and $3.6 \pm 0.6$ mm/yr for the KOF segment [Meghraoui et al., 2011]. However, Aktuğ et al. [2016] obtained slip rates ranging from $10.3 \pm 0.7$ mm/yr to $13.5 \pm 1.3$ mm/yr for the segments between Turkoglu and Karliova. Additionally, the rates were indicated to be $4.5 \pm 1.1$ mm/yr for the KF and $2.7 \pm 1.4$ mm/yr along the KOF in the same study.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** GPS network in the study area. Abbreviations for faults are EAF (East Anatolian Fault), DSF (Dead Sea Fault), KF (Karasu Fault), KOF (Karatas–Osmaniye Fault), and CA (Cyprus Arc). Fault mapping from Meghraoui et al. [2011]. Orange squares show provinces in the region.</figcaption></figure>

Tectonic processes in the Eastern Mediterranean region trigger seismic activities around the area. This is why active fault zones generating seismic hazards need to be monitored continuously using geodetic techniques (Figure 1). Some well-documented destructive historical seismic activities indicate potential risks in seismic zones, namely between the Turkoglu–Golbasi segment of the EAFZ and the KOF (1114 AD earthquake, Ms > 7.8), northern segment of the Dead Sea Fault (DSF) (1408 AD earthquake, M 7.4), KOF (1513 AD earthquake, M ~ 7.4), and KF (1822 AD earthquake, M 7.0 and 1872 AD earthquake, M 7.2) [Ambraseys and Jackson, 1998, Sbeinati et al., 2005, Ambraseys, 2004].
<table><caption>**Table 1.** Transformation results for GPS velocity field solutions</caption><thead><tr><th>Solution</th><th>wrms (mm/yr)</th><th>No. of common observation sites</th></tr></thead><tbody><tr><td>This study</td><td>0.76</td><td>97</td></tr><tr><td>Aktuğ et al. [2013]</td><td>0.67</td><td>5</td></tr><tr><td>Mahmoud et al. [2012]</td><td>0.57</td><td>6</td></tr><tr><td>Al-Tarazi et al. [2011]</td><td>0.72</td><td>20</td></tr><tr><td>Alchalbi et al. [2010]</td><td>0.74</td><td>10</td></tr><tr><td>Le Béon et al. [2008]</td><td>0.71</td><td>32</td></tr><tr><td>Gomez et al. [2007]</td><td>0.77</td><td>5</td></tr></tbody><tfoot><tr><td colspan="3">Reference solution: Reilinger et al. [2006]</td></tr></tfoot></table>

### 3. GPS data set and analysis

Using the GPS network presented in Section 2, campaign GPS observations were carried out at 24 campaign sites for four campaign sets in 2009, 2010, 2011, and 2019. Since 4 of the 24 campaign sites were destroyed, the last campaign in 2019 was achieved with the remaining 20 sites. The strategy adopted during campaign GPS observations is based on data logging in two sessions on consecutive days. In addition, each session covers an observation span of at least 8 hours. All campaign observations were performed within the same time periods of years to mitigate seasonal effects on GPS observations.

The GPS data analysis was performed in the GAMIT/GLOBK software [King et al., 2009] using validated strategies [Reilinger et al., 2006] for both campaign GPS observations and for also data obtained from municipalities, Turkish National Permanent GNSS Network-Active, and IGS permanent GNSS stations. First, in the GAMIT stage, station coordinates, zenith delay of the atmosphere at each station, and orbital and Earth orientation parameters (EOPs) were estimated from daily GPS phase observations using a weighted least squares algorithm. Second, in the GLOBK stage, Kalman filtering was applied to loosely constrained estimates of station coordinates, orbits, and EOPs and their covariances for a reference frame definition with a consistent set of coordinates and velocities. Eighteen stations selected from the IGS global networks were incorporated both to tie our local network with a global network and to estimate the orbital and earth rotation parameters with high accuracy. Apart from IGS stations, a total of 30 permanent GPS stations from local

municipalities and TUSAGA-Aktif networks were introduced into the daily GAMIT processing. Loosely constrained daily solutions from GAMIT processing were investigated on expected normalized root mean square (nrms) values between 0.15 and 0.25. After discarding outliers, coordinate repeatabilities were generated by controlling weighted root mean square (wrms) parameters approximately 2–4 mm for horizontal components and 10–15 mm for vertical components. Since the monumentation of GPS stations fluctuates [Langbein and Johnson, 1997], a random walk of 2 mm/√yr in horizontal positions of the stations was applied to obtain more realistic uncertainties for GPS velocities [McClusky et al., 2000].

Following the procedures for the definition of the reference frame for both coordinates and velocities [Aktuğ et al., 2009], transformation parameters to ITRF2008 coordinates of selected IGS stations were estimated for the stabilization frame. The stabilization was realized for the Eurasia-fixed reference frame by minimizing the adjustments to the horizontal velocities of the 17 stations where these velocities were obtained with a post root mean square (rms) value of 0.31 mm/yr.

The GPS velocity field obtained from this study was combined individually with published velocities while introducing the velocities by Reilinger et al. [2006] as the reference solution. The transformation accuracy of the velocity field solutions was expressed by the wrms parameter depending on the common observation sites (Table 1). The combined velocity field obtained from the transformation model was referenced to the ITRF2000 Eurasia-fixed reference frame (Figure 2).
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Eastern Mediterranean GPS velocity field in the ITRF2000 Eurasia-fixed reference frame with 95% confidence ellipses (fault mapping from GEM GAF-DB).</figcaption></figure>

## 4. Kinematic modeling

Kinematic modeling relying on GPS velocity fields instead of interpreting site velocities individually strengthens tectonic analysis in the region. Blocks bounded by faults have different characteristics inside the mass and along their boundaries. Although the block motions are systematic, the resistance due to the friction of rocks prevents free movement at their boundaries. This velocity difference, known as “slip deficit”, causes strain accumulations generally ending in an earthquake. Slip deficit describes the divergence of surface velocities from expected velocities of the blocks depending on geological factors.

The block modeling approach on the basis of inversion of GPS velocities provides the estimation of block rotation rates, fault locking parameters, and internal strain rates in an elastic and homogeneous half-space [Okada, 1985]. The components of the strain-rate tensor aim to reveal internal deformations within the blocks that are caused by unsettled faults. Apart from GPS site velocities, block and fault geometries formed in a three-dimensional (3D) elastic half-space, euler poles for blocks and locking depths are considered as model inputs for parameter estimation in the inversion. However, defects in block modeling emerge from poor density of observation sites. In this case, slip rate estimations become more susceptible to introduced block geometries and a priori locking depths. The TDEFNODE software [McCaffrey, 2009] based on the simulated annealing method [Press et al., 1986] was used for block modeling by minimizing the residuals between observed and model velocities. According to the approach in McCaffrey [2002], minimizing the misfit of the velocities allows the simultaneous estimation of block rotations using Euler poles and coupling fractions along block boundaries. The data misfit minimization is achieved by reduced chi-square statistic tests.

The locking depth in a seismogenic zone roughly reflects the characteristic mechanism of the fault. The state of being locked or creeping for a fault segment is controlled by coupling fraction. This parameter denoted by $\phi$ ($\varphi$) depends on the short-term slip rate ($V_c$) and the plate velocity ($V$), where $\phi = 0$ corresponds to creeping and $\varphi = 1$ to a fully locked fault while partially locked faults have values between 0 and 1 as per (1).

$$ \varphi = 1 - \frac{V_c}{V}. \qquad (1) $$

Before inversion, fault geometries are considered along the strike and dip directions as all the faults introduced in the kinematic model reach a depth of 30 km downward in almost vertical planes. The nodes along strike and down dip surrounding the blocks in the model define the fault surfaces in 3D with horizontal positions in degrees and depth in kilometers. Since the fault surface is composed of patches (quadrilaterals) between nodes, the kinematic model herein assembles patches bounded by $2 \times 2$ km along the strike and dip directions on the faults to estimate strain accumulations in the region.
The kinematic model consists of creeping, fully or partially locked segments [Wang et al., 2003], while the segments in the tectonic structure accommodate purely strike-slip or along with dip-slip components in the Eastern Mediterranean region [Yavaşoğlu, 2009, Tiryakioğlu, 2012]. The kinematic model in our study is assumed to be locked from the surface down to a locking depth fixed at 7 km in the crust. This approach constrains the fault from creeping since a thicker zone greater than 7 km in depth is recommended when exploring continental transform faults [Genrich et al., 2000, Le Béon et al., 2008, Walters et al., 2011]. In a deeper layer below the locked zone, the fault is assumed to be locked partially down to the depth below where full creeping begins. In our model, the faults were introduced as not creeping with a constraint of $\varphi = 1$ from 0 to 7 km in the zone close to the surface. Following this, $\varphi$ parameters were estimated between 0 and 1 in the inversion for depths ranging from 7 km to the fully creeping limit. Below the partially locked zone, the faults creep fully with $\varphi = 0$. Between locked and fully creeping zones, which is called the “effective transition zone” by Wang et al. [2003], the down dip transition was estimated using a factor of $\lambda = 0.2$ based on an exponential function, which indicates the shape of the slip distribution. During the inversion process, along-strike smoothing was introduced using a penalty function with a smoothing factor referred to above. No special data weighting strategy was followed since all the velocity solutions in modeling were equally weighted. However, the velocity field published by Mahmoud et al. [2012] was not incorporated into the kinematic model since the site velocity uncertainties on average are greater than those in each independent solution.

This study was initiated from a point of view comprising several blocks and main faults, namely Arabia (arab-reference block), Anatolia (anat), and Sinai (sina) blocks and East Anatolian Fault (EAF), KF, Cyprus Arc (CA), and DSF. The east of the Amanous (E-AMNS) block boundary equates to the KF, while the western boundary of the Amanous block was assumed to be free slipping (nonlocked) without any elastic strain accumulated. As discussed in Mahmoud [2012] and Mahmoud et al. [2012], block and fault geometries in the kinematic model vary from the simplest to the most complex related to the tectonic structures in the region. The different

kinematic models from Mahmoud [2012] were considered. However, evaluating the statistical results of the inversion, the geometry of the proposed kinematic model fits better with observations. The simplest model with the blocks of arab, sina, and anat comprises EAF, KF, CA, and DSF. The reduced chi-square parameter of this simplest model was estimated as 4.816 with a wrms of 1.476 mm and an nrms of 2.124. The more complex model consists of one additional block and one fault segment compared to the simplest model. The addition of the Iskenderun Block (iskd) and the KOF improves the reduced chi-square parameter to 3.869 with a wrms

<figure><img src="image_6.png" /></figure>

**Figure 3.** Fault geometry for Eastern Mediterranean region (red and black vectors represent the model and observed velocities with 95% confidence ellipses, respectively). Purple lines are fault segments and blue circles are nodes connecting fault segments to each other. Orange boxes denote block names defined in the model. EAF (East Anatolian Fault), KOF (Karatas–Osmaniye Fault), DSF (Dead Sea Fault), CA (Cyprus Arc), N-AMNS (North Boundary of Amanous Block), E-AMNS (East Boundary of Amanous Block), and S-AMNS (South Boundary of Amanous Block).
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Residuals between the model and observed velocities with 95% confidence ellipses.</figcaption></figure>

of 1.122 and an rms of 1.943. On the other hand, the Amanous block (amns) was introduced additionally into the model in the last stage. Eventually, the GPS observations fit well with the kinematic model formed with five blocks and five fault segments (Figure 3). The significant nonsystematic residuals are noted with 1.35 mm overall data rms between the model and observations (Figure 4). Depending on the block and fault geometries introduced, the reduced chi-square parameter of the model converges to 1.814.

In light of the findings from this research, a left-lateral strike slip at a rate of 7.5 mm/yr on the EAF allows significant stress transfer to the KF at 4.4–5.4 mm/yr sinistral strike-slip rates (Figure 5). Moreover, it contributes to a sinistral strike slip at a rate of 3.4 mm/yr and a reverse slip at a rate of 3.1 mm/yr along the KOF. Reverse faulting on the KOF indicates contraction in agreement with the result by Aktuğ et al. [2016]. Normal faulting on the KF has a slip at a rate approximately 3 mm/yr indicating extension, which is consistent with the result by Reilinger et al. [2006]. The internal block strain-rate tensors, which are shown by blue double-headed

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Green boxes indicate the strike-slip rates in which negative values correspond to left-lateral motion. White boxes indicate dip-slip rates corresponding to extension, while yellow boxes denote contraction. Red vectors demonstrate GPS-derived block velocities and black rotated arrows denote block rotation rates with estimated values on the right side when the Arabian block is assumed as the reference block. Blue double-headed arrows indicate the estimated internal block strain-rate axes. Orange boxes denote block names defined in the model.</figcaption></figure>

arrows in the same figure, clearly point out a compression of approximately 30 nanostrain/year within the Amanous Block. However, the greatest strain accumulation was loaded on the Iskenderun Block at approximately 40 nanostrain/year, inducing extension parallel to the main faults and compression perpendicular to them. It is also possible to see a similar pattern on the Sina Block with a strain rate almost half of the Iskenderun Block. The slowest straining region seems to be the Anatolian Block with an internal deformation of approximately 10 nanostrain/year. Furthermore, a locking depth of less than 10 km was estimated for the EAF. However, it is partially locked
<figure><img src="image_9.png" /></figure>

**Figure 6.** Distribution of fault locking depths. Phi represents fully locked fault zones in red, while purple stands for free slipping. Color tones between red and purple demonstrate partially locked zones.

down to a 25 km depth. In contrast, just like the KOF, the KF has a locking depth of 20 km (Figure 6).

**5. Discussion and conclusion**

In this study, we explored the strain accumula-
tions on major tectonic features around the Eastern
Mediterranean region by the inversion of GPS veloc-
ity estimations to obtain block rotation rates, fault
locking parameters, and internal strain rates in a 3D
elastic dislocation model. For this purpose, a region-
specific GPS network, composed of campaign obser-
vation sites and permanent stations, was formed to
clarify kinematic characteristics of the tectonically
active zones around the Hatay Triple Junction in
the region. Initially, former campaign data sets from
2009, 2010, and 2011 acquired in joint campaign ob-
servation sites from previous studies were collected.
The campaign data sets from previous years provide
a historical kinematic background for the region.
Moreover, it was required to determine recent tec-
tonic deformations. Ongoing strain accumulations
lead this research to review the estimations for slip
distributions and locking depths using new GPS ob-
servations across major faults. Therefore, GPS ob-
servations performed in 2019 on the same network,
which was the fourth time after 2009, 2010, and 2011,

have great importance for understanding the present
kinematics of the Eastern Mediterranean region. The
campaign GPS observations in 2019 helped improve
the data quality in estimating velocity uncertainties.
In addition to the new campaign data set acquired in
2019, permanent GPS stations from local networks
operated by municipalities were merged into our re-
gional network to densify observation sites spatially
and temporally. Both the new campaign data set and
permanent GPS stations, regarded as the true value
of this study, constitute unique data sets that were
not available in the past.

Accordingly, slip rates for the faults mentioned
above were published in several studies recently.
Earlier, Aktuğ et al. [2016] and Mahmoud et al.
[2012] estimated 10.5 mm/yr and 9.0 mm/yr strike-
slip rates, respectively, while Reilinger et al. [2006]
reported 10.0 mm/yr slip rates for the left-lateral
EAFZ. However, there seems to be no consensus on
normal/reverse fault components of the EAF since
Reilinger et al. [2006] argued for extension at a
5.1 mm/yr dip-slip rate, while Aktuğ et al. [2016]
proposed contraction at 2.4–6.3 mm/yr reverse-slip
rates. On the other hand, Mahmoud et al. [2012]
and our study assert no significant normal/reverse
fault components for the EAF. Similarly, the KF was
also studied many times recently so that the esti-
mations for strike-slip rates are 4.5 mm/yr by Ak-
tuğ et al. [2016], 4.0 mm/yr by Mahmoud et al.
[2012], 6.8 mm/yr by Reilinger et al. [2006], and 4.4–
5.4 mm/yr from this study. There are also some de-
bates about normal/reverse faulting on the KF. In
contrast, the KOF needs to be studied further in
detail.

In conclusion, a sinistral slip rate of 7.5 mm/yr
estimated from our model was initialized to accu-
mulate a strain of approximately 6.8 m in total af-
ter the 1114 AD earthquake (Ms > 7.8) between the
Turkoglu–Golbasi segment of the EAFZ and the KOF.
The total strain accumulation within the Golbasi–
Turkoglu segment supports the prediction of the next
probable major earthquake with a magnitude of 7.2–
7.6 if it ruptures entirely over its full 90 km length.
Similarly, a 3.4 mm/yr dextral slip rate for the KOF
corresponds to a 1.7 m total strain since the 1513 AD
earthquake (M ~ 7.4) and makes M 7.1 possible for
the next probable major earthquake. However, Mw
6.8–7.0 is predicted in the case of two segments with a
rupture length of approximately 40 km in each. Con-
sidering the 1822 AD (M 7.0) and 1872 AD (M 7.2) earthquakes on the KF, a total strain accumulation of 0.65–0.80 m at a 4.4–5.4 mm/yr sinistral slip rate confirms the prediction of a possible M 6.8–6.9 earthquake along 150 km based on Wells and Coppersmith [1994]. As per the approach by Anderson et al. [1996], another prediction is a magnitude of Mw 7.2 if one of the 75 km segments ruptures along the KF.

## Acknowledgments

The authors would like to thank all participants in this project who helped during the field work. This article was supported financially by the Coordinator of Scientific Research Projects (BAP) of Hitit University (project no. ODMYO19001.19.002) and Istanbul Technical University (project no. MGA-2019-42243). For GPS data processing, the GAMIT/GLOBK software developed by Massachusetts Institute of Technology (MIT) was used. The TDEFNODE software [McCaffrey, 2009] was used for block modeling. The GMT software [Wessel et al., 2013] was utilized to plot maps in the paper. The global data set of active fault lines was obtained from the database of the GEM Global Active Faults project.

## References

* Aktuğ, B., Kiliçoğlu, A., Lenk, O., and Gürdal, M. A. (2009). Establishment of regional reference frames for quantifying active deformation areas in Anatolia. *Stud. Geophys. Geod.*, 53:169–183.

* Aktuğ, B., Özener, H., Doğru, A., Sabuncu, A., Turgut, B., Halıcıoğlu, K., Yılmaz, O., and Havazlı, E. (2016). Slip rates and seismic potential on the east anatolian fault system using an improved GPS velocity field. *J. Geodyn.*, 94:95–112.

* Aktuğ, B., Parmaksız, E., Kurt, M., Lenk, O., Kılıçoğlu, A., Gürdal, M. A., and Özdemir, S. (2013). Deformation of central anatolia: GPS implications. *J. Geodynamics*, 67:78–96.

* Al-Tarazi, E., Rajab, J. A., Gomez, F., Cochran, W., JAAFAR, R., and Ferry, M. (2011). GPS measurements of near-field deformation along the southern Dead Sea Fault System. *Geochem. Geophys. Geosyst.*, 12.

* Alchalbi, A., Daoud, M., Gomez, F., McClusky, S., Reilinger, R., Romeyeh, M. A., Alsouod, A., Yassiminh, R., Ballani, B., and Darawcheh, R. (2010). Crustal deformation in northwestern Arabia from

GPS measurements in Syria: slow slip rate along the northern Dead Sea Fault. *Geophys. J. Intl*, 180:125–135.

* Ambraseys, N. N. (2004). The 12th century seismic paroxysm in the Middle East: a historical perspective. *Ann. Geophys.*, 47:733–758.

* Ambraseys, N. N. and Jackson, J. A. (1998). Faulting associated with historical and recent earthquakes in the Eastern Mediterranean region. *Geophys. J. Intl*, 133:390–406.

* Anderson, J. G., Wesnousky, S. G., and Stirling, M. W. (1996). Earthquake size as a function of fault slip rate. *Bull. Seismol. Soc. Am.*, 86:683–690.

* Arpat, E. and Saroglu, F. (1972). The East Anatolian fault system: thoughts on its development. *Bull. Miner. Res. Explor. Inst. Turk*, 78:33–39.

* Bertrand, S., Meghraoui, M., Mcclusky, S., Altunel, E., Ergintav, S., and Reilinger, R. (2006). Present-day crustal motions at the triple junction between the Dead Sea Fault, the East Anatolian Fault and the Cyprus Arc (SE Turkey). *Eur. Geosci. Union Geophys. Res. Abstr.*, 8.

* Emre, Ö., Duman, T. Y., Özalp, S., Şaroğlu, F., and Olgun, Ş. (2018). Active fault database of Turkey. *Bull. Earthq. Eng.*, 16(8):3229–3275.

* Genrich, J. F., Bock, Y., McCaffrey, R., Prawirodirdjo, L., Stevens, C. W., Puntodewo, S. S. O., Subarya, C., and Wdowinski, S. (2000). Distribution of slip at the northern Sumatran fault system. *J. Geophys. Res.*, 105(B12):28327–28341.

* Gomez, F., Karam, G., Khawlie, M., McClusky, S., Vervant, P., Reilinger, R., Jaafar, R., Tabet, C., Khair, K., and Barazangi, M. (2007). Global Positioning System measurements of strain accumulation and slip transfer through the restraining bend along the Dead Sea fault system in Lebanon. *Geophys. J. Int.*, 168:1021–1028.

* King, R. W., McClusky, S., and Herring, T. (2009). *Documentation of the GAMIT GPS Analysis Software, Release 10.3*. Massachusetts Institute of Technology, Cambridge.

* Langbein, J. and Johnson, H. (1997). Correlated errors in geodetic time series: implications for time-dependent deformation. *J. Geophys. Res.*, 102:591–603.

* Le Béon, M., Klinger, Y., Amrat, A. Q., Agnon, A., Dorbath, L., Baer, G., Ruegg, J. C., Charade, O., and Mayyas, O. (2008). Slip rate and locking depth from GPS profiles across the southern Dead Sea Trans-
form. *J. Geophys. Res.*, 113:1–19.

Mahmoud, Y. (2012). Caractérisation géodésique de la déformation active du point triple d'Hatay (Syrie-Turquie), Université de Strasbourg, Ecole et Observatoire des Science de la Terre, Thèse de doctorat de l'Université de Strasbourg.

Mahmoud, Y., Masson, F., Meghraoui, M., Cakir, Z., Alchalbi, A., Yavasoglu, H., Yonlu, O., Daoud, M., Ergintav, S., and Inan, S. (2012). Kinematic study at the junction of the East Anatolian fault and the Dead Sea fault from GPS measurements. *J. Geodyn.*, 67:30–39.

McCaffrey, R. (2002). Crustal block rotations and plate coupling. *Plate Bound. Zones, Geodyn. Ser.*, 30:101–122.

McCaffrey, R. (2009). Time-dependent inversion of three-component continuous GPS for steady and transient sources in northern Cascadia. *Geophys. Res. Lett.*, 36.

McClusky, S., Balassanian, S., Barka, A., Demir, C., Ergintav, S., Georgiev, I., Gurkan, O., Hamburger, M., Hurst, K., and Kahle, H. (2000). Global Positioning System constraints on plate kinematics and dynamics in the eastern Mediterranean and Caucasus. *J. Geophys. Res.*, 105:5695–5719.

Meghraoui, M., Cakir, Z., Masson, F., Mahmoud, Y., Ergintav, S., Alchalbi, A., Inan, S., Daoud, M., Yonlu, O., and Altunel, E. (2011). Kinematic modelling at the triple junction between the Anatolian, Arabian, African plates (NW Syria and in SE Turkey). *Geo-phys. Res. Abstr*, 13.

Nocquet, J. M. (2012). Present-day kinematics of the Mediterranean: a comprehensive overview of GPS results. *Tectonophysics*, 579:220–242.

Okada, Y. (1985). Surface deformation due to shear and tensile faults in a half-space. *Bull. Seismol. Soc. Am.*, 75:1135–1154.

Press, W. H., Flannery, B. P., Teukolsky, S. A., and Vetterling, W. T. (1986). *Numerical Recipes*, Vol. 547. Cambridge University Press, New York, USA.

Reilinger, R., McClusky, S., Vernant, P., Lawrence, S., Ergintav, S., Cakmak, R., Ozener, H., Kadirov, F., Guliev, I., and Stepanyan, R. (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.

Reilinger, R. E., McClusky, S. C., Oral, M. B., King, R. W., Toksoz, M. N., Barka, A. A., Kinik, I., Lenk, O., and Sanli, I. (1997). Global Positioning System measurements of present-day crustal movements in the Arabia–Africa–Eurasia plate collision zone. *J. Geophys. Res.*, 102:9983–9999.

Saroglu, F. and Yilmaz, Y. (1990). Tectonics of the Karliova triple junction. *Bull. Tech. Univ. Istanb.*, 44:475–493.

Sbeinati, M. R., Darawcheh, R., and Mouty, M. (2005). The historical earthquakes of Syria: an analysis of large and moderate earthquakes from 1365 BC to 1900 AD. *Ann. Geophys.*, 48(3):347–435.

Tiryakioğlu, İ. (2012). GNSS ölçüleri ile Güneybatı Anadolu'daki blok hareketleri ve gerinim alanlarının belirlenmesi. Yıldız Teknik Üniversitesi Fen Bilimleri Enstitüsü, Doktora Tezi, YTÜ, İstanbul.

Walters, R. J., Holley, R. J., Parsons, B., and Wright, T. J. (2011). Interseismic strain accumulation across the North Anatolian Fault from Envisat InSAR measurements. *Geophy. Res. Lett.*, 38.

Wang, K., Wells, R., Mazzotti, S., Hyndman, R. D., and Sagiya, T. (2003). A revised dislocation model of interseismic deformation of the Cascadia subduction zone. *J. Geophys. Res.*, 108(B1):1–13.

Wells, D. L. and Coppersmith, K. J. (1994). New empirical relationships among magnitude, rupture length, rupture width, rupture area, and surface displacement. *Bull. Seismol. Soc. Am.*, 84:974–1002.

Wessel, P., Smith, W. H. F., Scharroo, R., Luis, J. F., and Wobbe, F. (2013). Generic mapping tools: improved version released. *EOS Trans. AGU*, 94:409–410.

Westaway, R. O. B. and Arger, J. A. N. (1996). The Golbasi basin, southeastern Turkey: a complex discontinuity in a major strike-slip fault zone. *J. Geol. Soc.*, 153:729–744.

Yavaşoğlu, H. (2009). *Kuzey Anadolu Fayının Orta Anadolu bölümündeki güncel tektonik aktivitenin jeodezik yöntemler ve elastik yarı uzay modelleme ile belirlenmesi*. İstanbul Teknik Üniversitesi Fen Bilimleri Enstitüsü, Doktora Tezi, İTÜ, İstanbul.