# Resolving the slip-rate inconsistency of the northern Dead Sea fault

Xing Li, Sigurjón Jónsson, Shaozhuo Liu, Zhangfeng Ma, Nicolás

Castro-Perdomo, Simone Cesca, Frédéric Masson, Yann Klinger

► To cite this version:

Xing Li, Sigurjón Jónsson, Shaozhuo Liu, Zhangfeng Ma, Nicolás Castro-Perdomo, et al. Resolving the slip-rate inconsistency of the northern Dead Sea fault. Science Advances, 2024, 10 (11), 10.1126/sciadv.Adj8408. [https://doi.org/10.1126/sciadv.Adj8408](https://doi.org/10.1126/sciadv.Adj8408). [arXiv:2406.08548](https://arxiv.org/abs/2406.08548)

HAL Id: hal-04723588

https://hal.science/hal-04723588v1

Submitted on 8 Oct 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished 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_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
GEOLOGY

# Resolving the slip-rate inconsistency of the northern Dead Sea fault

Xing Li<fnref n="1" />, Sigurjón Jónsson<fnref n="1" />*, Shaozhuo Liu<fnref n="1" />, Zhangfeng Ma², Nicolás Castro-Perdomo<fnref n="1" />,³, Simone Cesca⁴, Frédéric Masson⁵, Yann Klinger⁶

Reported fault slip rates, a key quantity for earthquake hazard and risk analyses, have been inconsistent for the northern Dead Sea fault (DSF). Studies of offset geological and archeological structures suggest a slip rate of 4 to 6 millimeters per year, consistent with the southern DSF, whereas geodetic slip-rate estimates are only 2 to 3 millimeters per year. To resolve this inconsistency and overcome limited access to the northern DSF in Syria, we here use burst-overlap interferometric time-series analysis of satellite radar images to provide an independent slip-rate estimate of ~2.8 millimeters per year. We also show that the high geologic slip rate could, by chance, be inflated by earthquake clustering and suggest that the slip-rate decrease from the southern to northern DSF can be explained by splay faults and diffuse offshore deformation. These results suggest a microplate west of the northern DSF and a lower earthquake hazard for that part of the fault.

Copyright © 2024 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution License 4.0 (CC BY).

INTRODUCTION

The Dead Sea fault (DSF), the ~1000-km-long left-lateral transform plate boundary between the Sinai and Arabian plates, has been extensively studied since the 1950s (1, 2). The over 2000-year-long historical earthquake catalog of the DSF and numerous geological investigations have made the fault a test bed for a range of hypotheses related to earthquake occurrence patterns and clustering, characteristic earthquakes, fault segmentation, and earthquake interactions (3–7). Structurally, the DSF can be divided into three main sections from south to north (Fig. 1A). The southern section extends from the Red Sea and the transtensional Gulf of Aqaba along Wadi Arabah and the Jordan Valley, including the extensional step of the Dead Sea basin. The central section consists of the Lebanese restraining bend, which has several splay faults, with the Yammouneh fault as the main through-going structure. Further north, the northern DSF section that runs through the coastal ranges in northwestern Syria consists of the Misyaf fault segment and the faults bounding the Ghab basin and connects to the East Anatolian Fault at the Hatay triple junction in southern Türkiye (8), near the southern fault-rupture termination of the first [moment magnitude (Mw) 7.8] mainshock of the 6 February 2023 Kharamanmaraş earthquake doublet (9–11).

The geodetic and geologic slip-rate estimates are largely consistent for the southern and central DSF. In the south, geodetic slip rates of ~5 mm/year are similar to those from offset geological markers spanning a range of periods since the formation of the DSF (2, 12–23). Also, the geodetic slip rate of 3 to 4 mm/year at the Yammouneh fault of the Lebanese restraining bend agrees with the lower bound of the Late Pleistocene/Holocene slip rate of 3.8 to 6.4 mm/year (24–27). These results, therefore, show no evidence of a secular change in slip rate along the DSF.

Estimated slip rates along the northern DSF from offsets of several archeological sites as well as of geologic and geomorphological features indicate the range of ~4 to 6 mm/year (28–34), i.e., similar to geologic and geodetic slip rates on the southern DSF. Modeling of GNSS (Global Navigation Satellite System) velocities, on the other hand, indicates that the geodetic slip rate for the northern DSF is only 2.2 to 2.5 mm/year (35, 36), notably lower than estimates of slip rate on other sections of the DSF. Resolving this slip-rate inconsistency on the northern DSF has, however, remained challenging, in part due to the inaccessibility of the fault in Syria during the past decade.

Remote sensing observations, such as from interferometric synthetic aperture radar (InSAR), can often help to overcome inaccessibility challenges. In this case, however, the north-trending strike of the DSF poses a problem, as conventional InSAR inherently has limited sensitivity to north-south displacements (37, 38). While along-track pixel offsets (39, 40) or multiple aperture interferometry (41, 42) provide information about north-south displacements from SAR images, these methods only resolve displacements larger than tens of centimeter and are not capable of resolving millimeter per year interseismic velocities. Burst-overlap interferometry (BOI) from Sentinel-1 radar data, on the other hand, can provide more precise along-track observations, albeit only along narrow swaths where radar bursts overlap (43, 44). Time-series analysis of a large number of BOI data has been shown to allow retrieval of interseismic velocities as low as 5 mm/year (20, 45, 46); thus, here we extend our work from the Gulf of Aqaba (20) and use BOI time-series analysis (see Materials and Methods) for the entire DSF with data from 2014 to 2021 to provide critical constraints on the interseismic deformation and geodetic slip rates along the DSF.

RESULTS

Fault-parallel velocities from burst-overlap interferometry

We use 310 ascending and 330 descending radar images from the Sentinel-1A/B satellites, acquired between October 2014 and March 2021, to derive the horizontal velocities along the DSF. The descending along-track velocities closely approximate the fault-parallel motion as most of the DSF is almost parallel to the flight direction of the

<footnote n="1">¹King Abdullah University of Science and Technology (KAUST), Thuwal 23955, Saudi Arabia. ²Earth Observatory of Singapore, Singapore 639798, Singapore. ³Department of Earth and Atmospheric Sciences, Indiana University, Bloomington, IN 47405, USA. ⁴GeoForschungZentrum (GFZ), 14473 Potsdam, Germany. ⁵ITES, EOST Université de Strasbourg, CNRS, 67084 Strasbourg, France. ⁶Université Paris Cité, Institut de Physique du Globe de Paris, CNRS, F-75005 Paris, France.</footnote>

<footnote>*Corresponding author. Email: sigurjon.jonsson@kaust.edu.sa</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Seismicity and geodetic observations along the DSF. The 1000-km-long Dead Sea left-lateral transform fault system extends from the Red Sea in the south to the Hatay triple junction (HTJ) and East Anatolian Fault (EAF) in the north, here shown with (A) earthquake locations (magnitude, >2) from 2000 to 2022 (Euro-Mediterranean Seismological Centre, EMSC), (B) fault-parallel velocities with respect to the Arabian plate (negative values indicate southward motion) estimated from time-series analysis of burst-overlap interferometric (BOI) observations, and (C) Global Navigation Satellite System (GNSS) velocities from combining multiple sources, again with respect to the Arabian plate. NDSF, northern Dead Sea fault; YF, Yammouneh fault; JV, Jordan Valley; DS, Dead Sea; WA, Wadi Arabah; NGA, northern Gulf of Aqaba; SGA, southern Gulf of Aqaba.</figcaption></figure>

descending orbit. The horizontal fault-parallel velocities from BOI (i.e., the descending-orbit along-track velocities) attest to subtle left-lateral motion in a northward direction across the DSF (Fig. 1B). In the south, across the Gulf of Aqaba and Wadi Arabah, our BOI velocities (in Fig. 1B) clearly show the ~5 mm/year relative plate motion. Across the Jordan Valley and in Lebanon, the relative motion is also clear in most of the BOI velocities, but the signals are somewhat smaller and noisier than in the south. Further north, across the northern DSF, the fault parallel motion is noisy and less pronounced, due to lower coherence. Similar results are found from time-series analysis of ascending-orbit BOI data (see fig. S1),

although the orientation of the ascending-orbit BOI profiles is less favorable to assess the fault-parallel motion.

The northward reduction in fault parallel velocities in the BOI results is in accord with what GNSS observations have shown (36). We compiled multiple GNSS results into one solution (12, 14, 15, 35, 36, 47) and when displayed in the Arabian plate reference frame (Fig. 1C), the relative motion between Arabia and the coastal ranges of Syria is small, clearly smaller than the relative motion across the DSF further south. Both the GNSS and BOI results thus show that the southern Sinai plate (west of the southern DSF) is moving faster to the south (relative to Arabia) than the northern Sinai plate (west

Downloaded from https://www.science.org on October 08, 2024
of the northern DSF), indicating substantial relative motion between the northern and southern parts of the Sinai plate (36) (Fig. 1, B and C).

To reduce noise and get more robust fault-parallel velocities, we stack all neighboring BOI profiles across each fault segment (as indicated by black lines in Fig. 1B), down-sample the observations and select high-quality points, and assess the uncertainty (see Materials and Methods). The stacked BOI profiles clearly show gradients in the fault-parallel velocities, diagnostic of interseismic strain accumulation, which are consistent with those obtained from GNSS observations (Fig. 2, A to D). Also, the level of noise is comparable between the GNSS and BOI results, but many more BOI observations ensure a more complete mapping of the fault parallel velocities across the different parts of the DSF. As for the BOI and GNSS observations in map view (Fig. 1, B and C), the profiles show the same clear reduction of fault parallel velocities from about 5 mm/year in the south to 2 to 3 mm/year across the northern DSF (Fig. 2, A to D).

For each of the profile data, we estimate the geodetic slip rate (48) that confirms the trend of decreasing slip rate from south to north. The best-fitting slip rates for the southern half of the DSF are ~5 mm/year but then decrease to 3.8, 3.4, and 2.8 mm/year for the Jordan Valley, Yammouneh fault, and northern DSF, respectively (Fig. 2E and fig. S2).

The Jordan Valley slip rate decreases from the southern part, just north of the Dead Sea, to the northern part, south of the Sea of Galilee, from 4.3 to 3.5 mm/year (fig. S3). This suggests that a part of the slip rate is transferred from the DSF to the Carmel-Gilboa fault, which has been previously reported in GNSS studies estimating its slip rate in the range of 0.8 to 1.2 mm/year (14, 36). Within the restraining bend, the BOI deformation is not clearly associated with a single fault and we thus presume that the Yammouneh fault accounts for most of the observed deformation. Within uncertainties, the result is consistent with the reported block model rate of 3.8 mm/year (26) and with the lower bound of the Late Pleistocene/Holocene slip rate of 3.8 to 6.4 mm/year (24).

Together, the descending and ascending BOI results, the GNSS velocities, as well as our modeling results are consistent with the geologic slip rates (2, 21–23) for the southern and central parts of the ~1000-km DSF. They are also consistent with the relative plate motion rates according to the latest estimate of the Sinai-Arabia Euler pole (49). The exception is the northern DSF, where the geodetic results clearly point to lower slip rates that are inconsistent with the higher estimated geological rates and with what plate motion models predict (Fig. 2E).

<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Fault-parallel velocities and estimated geodetic slip rates. (A to D) Profiles across the DSF at different locations showing the fault-parallel velocity (into direction N12°E) from stacking several BOI profiles (colored circles) in comparison with GNSS velocities (gray); see Fig. 1B for locations and fig. S2 for further profile results. (E) Estimated geodetic slip rates and locking depths along the DSF (bold squares, 1-sigma confidence limits; see table S1) compared with other geodetic estimates (other symbols) and the predicted Sinai-Arabia plate motion prediction (49) (dashed line). For display purposes, some of the symbols have been slightly shifted vertically from their exact locations of estimation.</figcaption></figure>
DISCUSSION

Possible earthquake clustering

Although most estimates of geological slip rates on the northern DSF are in the range of ~4 to 6 mm/year, there is considerable variability in the reported slip rates. A relatively low slip rate of 1.4 to 4.5 mm/year was reported by Searle *et al.* (34), a result criticized later by Westaway (50) because of inaccurate dating of lava flow units offset by the fault and because of mistaken alignment of piercing points. A faster slip rate of 4.9 to 6.0 mm/year was estimated in the Amik basin based on offset alluvial fans and three archeological sites (29, 31). These rates, however, are not representative of the slip rate on the northern DSF, as the Amik basin is within the Hatay triple junction and thus influenced by the East Anatolian Fault and onland extensions of faults related to the Cyprus arc. The main geologic slip-rate results of the northern DSF come from the Misyaf segment of the fault, where paleoseismic trenching and archaeoseismic studies of the faulted Al-Harif Roman aqueduct in Syria by four earthquakes during the past 3500 years, each having an estimated offset of ~4 m, yield a slip rate of 4.6 mm/year, or even higher in the range of 4.9 to 6.3 mm/year (28, 33), depending on how the slip rate is calculated.

We postulate that the high estimated geological slip rate for the Misyaf segment could be elevated due to earthquake clustering during the past 3500 years. Earthquake clustering is widely reported in

paleoseismological studies for strike-slip faults where long enough records have been retrieved (4, 5, 51, 52). To test this conjecture, we constructed a synthetic earthquake catalog based on the unique 220,000-year-long earthquake record from the central Dead Sea (53), which provides a distribution of inter-event times between major earthquakes. Normalizing the inter-event distribution and then scaling it such that characteristic 4-m fault-offset Misyaf mainshocks (33) yield a slip rate of 2.8 mm/year, as our geodetic results indicate (Fig. 2A), we compute the probability of fast apparent slip rates within 3500-year time intervals in the synthetic catalog (Fig. 3). About 12% of the time windows include four or more mainshocks, resulting in elevated slip-rates greater than 4.6 mm/year. These results show that it is possible that the estimated 4.6 to 6.3 mm/year fault slip rate on the Misyaf segment of the northern DSF is due to earthquake clustering in the past 3500 years.

There are not many geological slip rate estimates of the northern DSF due to a lack of accurate piercing points and in situ field observations have not been possible for over the past decade. However, one geomorphological study based on a digital elevation model from InSAR yielded a post-Miocene slip rate of 3.3 to 4 mm/year along the northern DSF (30), i.e., only slightly higher than the geodetic slip rates. Taking together the geodetic results, the lower bound of this secular slip-rate estimate, and our reasoning above that the high rate on the Misyaf segment could be due to earthquake

<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Synthetic earthquake catalog and clustering. Possible earthquake clustering along the northern DSF, assessed from a synthetic earthquake catalog based on the inter-event time statistics observed in the central Dead Sea (53). (A) Normalized distribution of Dead Sea inter-event times and scaled such that for characteristic 4-m fault-offset mainshocks [as found at Misyaf (28, 35)], it yields a slip rate of 2.8 mm/year (as our geodetic results indicate; Fig. 2A), with an average recurrence interval of ~1430 years. (B and C) The number of mainshocks is most often one to three in any 3500-year-long time interval of the synthetic catalog, which was generated from randomly drawing 2000 inter-event time samples, yielding ~3-million-year-long synthetic catalog with an accumulated 8-km fault offset. (D) Example of estimated apparent slip rate in a 3500-year-long window of the catalog. (E) Histogram of ~15,000 apparent slip-rate estimates from sliding the 3500-year window in 200-year-long steps throughout the synthetic catalog. The apparent slip rate is high when, by chance, there are three or more mainshocks in the 3500-year-long window, but low if they are less. In ~12% of the cases, there are four or more mainshocks, resulting in elevated slip rates greater than 4.6 mm/year. Using smaller (3 m) or larger (5 m) characteristic earthquake offsets yield similar results of ~10 and ~14%, respectively (see fig. S4).</figcaption></figure>
clustering, the slip-rate inconsistency between geological and geodetic slip rates on the northern DSF has been resolved.

## Tectonic block configuration in the eastern Mediterranean

The question remains about where the ~2 mm/year difference between the low geodetic slip rate on the northern DSF and the predicted Sinai-Arabia plate motion is accommodated. We explored whether the low geodetic rate could be related to interactions with the offshore Cyprus arc or due to viscoelastic earthquake cycle effects, but neither possibility is likely (see fig. S5). Crustal shortening (1 mm/year) within the Arabian plate over the Palmyrides fold and thrust belt has been suggested (35), along with north-south extension (1.5 mm/year) just north of the Lebanese restraining bend. There is, however, scant evidence to support this hypothesis in the GNSS or BOI data (Figs. 1B and 2A). A north-south block boundary just offshore the Syrian coastal ranges has also been proposed in GNSS block modeling (36), representing strike-slip motion parallel to the northern DSF (blue dashed line in Fig. 4A). However, results of extensive offshore seismic imaging show no evidence for a strike-slip fault at this location (54).

Given the above, it appears that the northern Sinai plate moves separately from the main Sinai plate to the south. The Carmel-Gilboa and Roum faults splay off to the northwest from the DSF and into the Mediterranean, but no clear indication has been found for an offshore plate boundary extending from the coast to the Cyprus arc. Maps of the regional seismicity show distributed seismicity in the area from the coast of northern Israel and Lebanon extending to the Cyprus arc that is markedly stronger than the neighboring offshore seismicity north of Egypt and west of Syria (Fig. 4A). This seismicity suggests a broad zone of diffuse deformation in the eastern Mediterranean, similar to what has been found for some other diffuse oceanic plate boundaries (55, 56). We estimated deviatoric

moment tensor solutions for magnitude 3.3 to 5.0 events in this zone and they support (fig. S8 and table S3) that this zone accommodates the relative SSW-NNE extension between the main Sinai plate and its part west of the northern DSF.

Our results thus indicate that the northern Sinai plate is a separate microplate, here named the Latakia-Tartus microplate. It is bounded by the Cyprus arc to the west, the Hatay triple junction to the north, northern DSF to the east, and the Lebanon restraining bend and the offshore diffuse seismicity zone to the south. The GNSS velocities in southern Hatay, southwest of the Amik basin, and in northwestern Syria, clearly show a systematic pattern of velocities different from the velocities to the north (northern Hatay and Anatolia), east (Arabian plate), and south (Fig. 4A). Together, the results resolve the slip-rate inconsistency along the northern DSF and show that the seismic hazard associated with this fault is notably lower than if it were moving at 4 to 6 mm/year. Still, with the last major earthquake on the northern DSF in 1170 (32), a slip deficit amounting to about 2.4 m of fault slip has accumulated. While this might appear to indicate that the northern DSF is ripe for another major earthquake, the recent clustering of characteristic 4-m offset events suggests that it might not be the case.

## MATERIALS AND METHODS

We first geometrically aligned the reference and secondary radar images using precise orbits (https://dataspace.copernicus.eu/) for orbital and topographic corrections and formed interferograms (standard InSAR observations) after refining the image co-registration process with the enhanced spectral diversity method (57). For the narrow areas of overlapping radar bursts, we subtracted forward- and backward-looking interferograms and generated a redundant network of burst-overlap interferograms with a temporal threshold

<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Seismicity and proposed tectonic block configuration in the eastern Mediterranean. (A) Seismicity (from EMSC) showing elevated offshore seismicity between Lebanon/Israel and the Cyprus arc. This diffuse seismicity zone appears to bound the Latakia-Tartus microplate to the southwest, with the Cyprus arc bounding it to the west, Hatay triple junction in the north, and the northern DSF in the east. (B) Green vectors are GNSS velocities as in Fig. 1C, with respect to stable Arabia, and the blue dashed line in (A) marks the block boundary proposed by Gomez *et al.* (36).</figcaption></figure>
of 60 days. To reduce ionospheric delays and decorrelation noise, we filtered all the burst-overlap interferograms using a local low-pass filter. In addition, to further reduce noise in the descending-orbit data, we applied a threshold of 0.1 rad (2 cm) in each reconstructed interferogram. We then combined multiple BOI data with a Small Baseline Subset (SBAS) time-series analysis to retrieve the horizontal displacement along the entire DSF. For more details, see the Supplementary Materials.

We inverted for the along-track BOI ground velocity using conventional time-series analysis based on the reconstructed interferograms. However, due to low coherence along the central and northern DSF, the accuracy of each point could not meet the desired level. Therefore, we selected only the high-quality points after velocity stacking. To retain as many coherent observations as possible for velocity stacking, we did not apply any point selection to each overlap velocity. Instead, we applied point selection based on the SD after stacking. The smaller the SD, the better, as the slip rate is assumed to be similar at the same distance from the fault.

To perform the profile stacking, we followed these steps: (i) We down-sampled the observations at 3-km intervals perpendicular to the fault and included all possible overlap areas within each fault segment (as labeled in Fig. 1B). (ii) Using the down-sampled points, we calculated the mean and SD of the points at the same distance from the fault, assuming that the overlap areas are independent of each other. To obtain each median value and estimated 1σ uncertainty (at fault-perpendicular distance x), we included all the along-track BOI observations within the same fault-perpendicular distance bin (x ±0.6 km). (iii) We applied a threshold to each stacked profile, retaining only the most coherent points.

## Supplementary Materials

**This PDF file includes:**

Supplementary Text

Figs. S1 to S8

Tables S1 to S3

References

## REFERENCES AND NOTES

1. L. Dubertret, Problemes de la geologie du Levant. *Bull. Soc. Géol. Fr.* **5**, 3–31 (1947).

2. A. M. Quennell, The structural and geomorphic evolution of the Dead Sea rift. *Quart. J. Geol. Soc.* **114**, 1–24 (1958).

3. M. Ferry, M. Meghraoui, N. Abou Karaki, M. Al-Taj, L. Khalil, Episodic behavior of the Jordan Valley section of the Dead Sea fault inferred from a 14-ka-long integrated catalog of large earthquakes. *Bull. Seismol. Soc. Am.* **101**, 39–67 (2011).

4. N. Wechsler, T. K. Rockwell, Y. Klinger, Variable slip-rate and slip-per-event on a plate boundary fault: The Dead Sea fault in northern Israel. *Tectonophysics* **722**, 210–226 (2018).

5. M. Lefevre, Y. Klinger, M. Al-Qaryouti, M. Le Beon, K. Moumani, Slip deficit and temporal clustering along the Dead Sea fault from paleoseismological investigations. *Sci. Rep.* **8**, 4511 (2018).

6. M. Meghraoui, Paleoseismic history of the Dead Sea fault zone. *J. Earthq. Eng.*, 1–20 (2015).

7. Y. Klinger, M. Le Béon, M. Al-Qaryouti, 5000 yr of paleoseismicity along the southern Dead Sea fault. *Geophys. J. Int.* **202**, 313–327 (2015).

8. G. Brew, J. Lupa, M. Barazangi, T. Sawaf, A. Al-Imam, T. Zaza, Structure and tectonic development of the Ghab basin and the Dead Sea fault system, Syria. *J. Geol. Soc.* **158**, 665–674 (2001).

9. P. M. Mai, T. Aspiotis, T. A. Aquib, E. V. Cano, D. Castro-Cruz, A. Espindola-Carmona, B. Li, X. Li, J. Liu, R. Matrau, A. Nobile, K. H. Palgunadi, M. Ribot, L. Parisi, C. Suhendi, Y. Tang, B. Yalcin, U. Avşar, Y. Klinger, S. Jónsson, The destructive earthquake doublet of 6 February 2023 in South-Central Türkiye and Northwestern Syria: Initial observations and analyses. *Seism. Rec.* **3**, 105–115 (2023).

10. Z. Jia, Z. Jin, M. Marchandon, T. Ulrich, A.-A. Gabriel, W. Fan, P. Shearer, X. Zou, J. Rekoske, F. Bulut, The complex dynamics of the 2023 Kahramanmaraş, Turkey, M<sub>w</sub> 7.8-7.7 earthquake doublet. *Science* **381**, 985–990 (2023).

11. C. Liu, T. Lay, R. Wang, T. Taymaz, Z. Xie, X. Xiong, T. S. Irmak, M. Kahraman, C. Erman, Complex multi-fault rupture and triggering during the 2023 earthquake doublet in southeastern Türkiye. *Nat. Commun.* **14**, 5564 (2023).

12. N. Castro-Perdomo, R. Viltres, F. Masson, Y. Klinger, S. Z. Liu, M. Dhahry, P. Ulrich, J. D. Bernard, R. Matrau, A. Alothman, H. Zahran, R. Reilinger, P. M. Mai, S. Jonsson, Interseismic deformation in the Gulf of Aqaba from GPS measurements. *Geophys. J. Int.* **228**, 477–492 (2022).

13. M. Sadeh, Y. Hamiel, A. Ziv, Y. Bock, P. Fang, S. Wdowinski, Crustal deformation along the Dead Sea Transform and the Carmel Fault inferred from 12 years of GPS measurements. *J. Geophys. Res. Solid Earth* **117**, 2012JB009241 (2012).

14. Y. Hamiel, O. Piatibratova, Spatial variations of slip and creep rates along the Southern and Central Dead Sea fault and the Carmel-Gilboa fault system. *J. Geophys. Res. Solid Earth* **126**, e2020JB021585 (2021).

15. S. McClusky, R. Reilinger, S. Mahmoud, D. Ben Sari, A. Tealeb, GPS constraints on Africa (Nubia) and Arabia plate motions. *Geophys. J. Int.* **155**, 126–138 (2003).

16. S. Wdowinski, Y. Bock, G. Baer, L. Prawirodirdjo, N. Bechor, S. Naaman, R. Knafo, Y. Forrai, Y. Melzer, GPS measurements of current crustal movements along the Dead Sea fault. *J. Geophys. Res. Solid Earth* **109**, 2003JB002640 (2004).

17. M. Le Beon, Y. Klinger, A. Q. Amrat, A. Agnon, L. Dorbath, G. Baer, J. C. Ruegg, O. Charade, O. Mayyas, Slip rate and locking depth from GPS profiles across the southern Dead Sea Transform. *J. Geophys. Res. Solid Earth* **113**, 2007JB005280 (2008).

18. F. Masson, Y. Hamiel, A. Agnon, Y. Klinger, A. Deprez, Variable behavior of the Dead Sea fault along the southern Arava segment from GPS measurements. *C. R. Geosci.* **347**, 161–169 (2015).

19. Y. Klinger, J. P. Avouac, N. Abou Karaki, L. Dorbath, D. Bourles, J. L. Reyss, Slip rate on the Dead Sea transform fault in northern Araba valley (Jordan). *Geophys. J. Int.* **142**, 755–768 (2000).

20. X. Li, S. Jónsson, Y. M. Cao, Interseismic deformation from Sentinel-1 burst-overlap Interferometry: Application to the Southern Dead Sea fault. *Geophys. Res. Lett.* **48**, e2021GL093481 (2021).

21. R. Freund, Z. Garfunkel, I. Zak, M. Goldberg, T. Weissbrod, B. Derin, F. Bender, F. Wellings, R. Girdler, The shear along the Dead Sea rift. *Philos. Trans. R. Soc. Lond. A Math. Phys. Sci.* **267**, 107–130 (1970).

22. Z. Garfunkel, Internal structure of the Dead Sea Leaky Transform (Rift) in relation to plate kinematics. *Tectonophysics* **80**, 81–108 (1981).

23. S. Joffe, Z. Garfunkel, Plate kinematics of the circum Red Sea—A re-evaluation. *Tectonophysics* **141**, 5–22 (1987).

24. M. Daëron, L. Benedetti, P. Tapponnier, A. Sursock, R. C. Finkel, Constraints on the post-25-ka slip rate of the Yammoûneh fault (Lebanon) using in situ cosmogenic 36Cl dating of offset limestone-clast fans. *Earth Planet. Sci. Lett.* **227**, 105–119 (2004).

25. S. Mahmoud, R. Reilinger, S. McClusky, P. Vernant, A. Tealeb, GPS evidence for northward motion of the Sinai Block: Implications for E. Mediterranean tectonics. *Earth Planet. Sci. Lett.* **238**, 217–224 (2005).

26. F. Gomez, G. Karam, M. Khawlie, S. McClusky, P. Vernant, R. Reilinger, R. Jaafar, C. Tabet, K. Khair, M. Barazangi, 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 (2007).

27. F. Gomez, T. Nemer, C. Tabet, M. Khawlie, M. Meghraoui, M. Barazangi, Strain partitioning of active transpression within the Lebanese restraining bend of the Dead Sea fault (Lebanon and SW Syria). *Spec. Publ. Geol. Soc. Lond.* **290**, 285–303 (2007).

28. M. Meghraoui, F. Gomez, R. Sbeinati, J. Van der Woerd, M. Mouty, A. N. Darkal, Y. Radwan, I. Layrous, H. Al Najjar, R. Darawcheh, F. Hijazi, R. Al-Ghazzi, M. Barazangi, Evidence for 830 years of seismic quiescence from palaeoseismology, archaeoseismology and historical seismicity along the Dead Sea fault in Syria. *Earth Planet. Sci. Lett.* **210**, 35–52 (2003).

29. E. Altunel, M. Meghraoui, V. Karabacak, S. H. Akyuz, M. Ferry, C. Yalcmer, M. Munschy, Archaeological sites (tell and road) offset by the Dead Sea fault in the Amik Basin, Southern Turkey. *Geophys. J. Int.* **179**, 1313–1329 (2009).

30. F. Gomez, M. Khawlie, C. Tabet, A. N. Darkal, K. Khair, M. Barazangi, Late Cenozoic uplift along the northern Dead Sea transform in Lebanon and Syria. *Earth Planet. Sci. Lett.* **241**, 913–931 (2006).

31. V. Karabacak, E. Altunel, M. Meghraoui, H. S. Akyuz, Field evidences from northern Dead Sea fault zone (South Turkey): New findings for the initiation age and slip rate. *Tectonophysics* **480**, 172–182 (2010).

32. M. Meghraoui, R. Toussaint, M. Ferry, P. Nguema-Edzang, in *EGU General Assembly Conference Abstracts* (2015), p. 5218.

33. M. R. Sbeinati, M. Meghraoui, G. Suleyman, F. Gomez, P. Grootes, M.-J. Nadeau, H. Al Najjar, Timing of earthquake ruptures at the Al Harif Roman aqueduct (Dead Sea fault, Syria) from archaeoseismology and paleoseismology. *Ancient Earthquakes* **243**, (2010).

34. M. P. Searle, S. L. Chung, C. H. Lo, Geological offsets and age constraints along the northern Dead Sea fault, Syria. *J. Geol. Soc.* **167**, 1001–1008 (2010).

35. A. Alchalbi, M. Daoud, F. Gomez, S. McClusky, R. Reilinger, M. Abu Romeyh, A. Alsouod, R. Yassminh, B. Ballani, R. Darawcheh, R. Sbeinati, Y. Radwan, R. Al Masri, M. Bayerly,
R. Al Ghazzi, M. Barazangi, Crustal deformation in northwestern Arabia from GPS measurements in Syria: Slow slip rate along the northern Dead Sea fault. *Geophys. J. Int.* **180**, 125–135 (2010).

36. F. Gomez, W. J. Cochran, R. Yassminh, R. Jaafar, R. Reilinger, M. Floyd, R. W. King, M. Barazangi, Fragmentation of the Sinai plate indicated by spatial variation in present-day slip rate along the Dead Sea Fault System. *Geophys. J. Int.* **221**, 1913–1940 (2020).

37. J. H. Liu, J. Hu, Z. W. Li, Z. F. Ma, L. X. Wu, P. P. Jiang, G. C. Feng, J. J. Zhu, Complete three-dimensional coseismic displacements due to the 2021 Maduo earthquake in Qinghai Province, China from Sentinel-1 and ALOS-2 SAR images. *Sci. China Earth Sci.* **65**, 687–697 (2022).

38. Y. Fialko, M. Simons, D. Agnew, The complete (3-D) surface displacement field in the epicentral area of the 1999*Mw*7.1 Hector Mine earthquake, California, from space geodetic observations. *Geophys. Res. Lett.* **28**, 3063–3066 (2001).

39. T. Wang, S. Jónsson, Improved SAR amplitude image offset measurements for deriving three-dimensional coseismic displacements. *IEEE J-Stars* **8**, 3271–3278 (2015).

40. R. Michel, J. P. Avouac, J. Taboury, Measuring ground displacements from SAR amplitude images: Application to the Landers earthquake. *Geophys. Res. Lett.* **26**, 875–878 (1999).

41. H. S. Jung, W. J. Lee, L. Zhang, Theoretical accuracy of along-track displacement measurements from multiple-aperture interferometry (MAI). *Sensors* **14**, 17703–17724 (2014).

42. N. B. Bechor, H. A. Zebker, Measuring two-dimensional movements using a single InSAR pair. *Geophys. Res. Lett.* **33**, 2006GL026883 (2006).

43. H. Jiang, G. Feng, T. Wang, R. Bürgmann, Toward full exploitation of coherent and incoherent information in Sentinel-1 TOPS data for retrieving surface displacement: Application to the 2016 Kumamoto (Japan) earthquake. *Geophys. Res. Lett.* **44**, 1758–1767 (2017).

44. R. Grandin, E. Klein, M. Métois, C. Vigny, Three-dimensional displacement field of the 2015 *Mw*8.3 Illapel earthquake (Chile) from across- and along-track Sentinel-1 TOPS interferometry. *Geophys. Res. Lett.* **43**, 2552–2561 (2016).

45. Z. F. Ma, S. J. Wei, X. Li, Y. Aoki, J. H. Liu, X. J. Liu, W. F. Mao, S. Yang, N. X. Wang, Q. H. Huang, T. Huang, S. H. Yun, Challenges and prospects to time series Burst Overlap Interferometry (BOI): Some insights from a new BOI algorithm test over the Chaman fault. *IEEE Trans. Geosci. Remote Sens.* **60**, 5118219 (2022).

46. M. Lazecký, A. Hooper, P. Piromthong, InSAR-derived horizontal velocities in a global reference frame. *Geophys. Res. Lett.* **50**, e2022GL011173 (2023).

47. R. Reilinger, S. McClusky, P. Vernant, S. Lawrence, S. Ergintav, R. Cakmak, H. Ozener, F. Kadirov, I. Guliev, R. Stepanyan, M. Nadariya, G. Hahubia, S. Mahmoud, K. Sakr, A. ArRajehi, D. Paradissis, A. Al-Aydrus, M. Prilepin, T. Guseva, E. Evren, A. Dmitrotsa, S. V. Filikov, F. Gomez, R. Al-Ghazzi, G. Karam, GPS constraints on continental deformation in the Africa-Arabia-Eurasia continental collision zone and implications for the dynamics of plate interactions. *J. Geophys. Res. Solid Earth* **111**, 2005JB004051 (2006).

48. J. C. Savage, R. O. Burford, Geodetic determination of relative plate motion in central California. *J. Geophys. Res.* **78**, 832–845 (1973).

49. R. Viltres, S. Jónsson, A. O. Alothman, S. Z. Liu, S. Leroy, F. Masson, C. Doubre, R. Reilinger, Present-day motion of the Arabian plate. *Tectonics* **41**, e2021TC007013 (2022).

50. R. Westaway, Discussion on 'geological offsets and age constraints along the northern Dead Sea fault, Syria'. *J. Geol. Soc. London* **168**, 621–623 (2011).

51. S. Marco, A. Agnon, High-resolution stratigraphy reveals repeated earthquake faulting in the Masada Fault Zone, Dead Sea Transform. *Tectonophysics* **408**, 101–112 (2005).

52. R. Weldon, K. Scharer, T. Fumal, G. Biasi, Wrightwood and the earthquake cycle: What a long recurrence record tells us about how faults work. *GSA Today* **14**, 4–10 (2004).

53. Y. Lu, N. Wetzler, N. Waldmann, A. Agnon, G. P. Biasi, S. Marco, A 220,000-year-long continuous large earthquake record on a slow-slipping plate boundary. *Sci. Adv.* **6**, ab4170 (2020).

54. S. A. Bowman, Regional seismic interpretation of the hydrocarbon prospectivity of offshore Syria. *Geoarabia* **16**, 95–124 (2011).

55. J. A. Wiens, C. S. Crawford, J. R. Gosz, Boundary dynamics: A conceptual-framework for studying landscape ecosystems. *Oikos* **45**, 421–427 (1985).

56. R. G. Gordon, The plate tectonic approximation: Plate nonrigidity, diffuse plate boundaries, and global plate reconstructions. *Annu. Rev. Earth Planet Sci.* **26**, 615–642 (1998).

57. P. Prats-Iraola, R. Scheiber, L. Marotti, S. Wollstadt, A. Reigber, TOPS interferometry with TerraSAR-X. *IEEE Trans. Geosci. Remote Sens.* **50**, 3179–3188 (2012).

58. H. Fattahi, M. Simons, P. Agram, InSAR time-series estimation of the ionospheric phase delay: An extension of the split range-spectrum technique. *IEEE Trans. Geosci. Remote Sens.* **55**, 5984–5996 (2017).

59. C. Liang, E. J. Fielding, Measuring azimuth deformation with L-Band ALOS-2 ScanSAR interferometry. *IEEE Trans. Geosci. Remote Sens.* **55**, 2725–2738 (2017).

60. C. R. Liang, P. Agram, M. Simons, E. J. Fielding, Ionospheric correction of InSAR time series analysis of C-band Sentinel-1 TOPS data. *IEEE Trans. Geosci. Remote Sens.* **57**, 6755–6773 (2019).

61. H. S. Jung, D. T. Lee, Z. Lu, J. S. Won, Ionospheric correction of SAR Interferograms by multiple-aperture interferometry. *IEEE Trans. Geosci. Remote Sens.* **51**, 3191–3199 (2013).

62. A. M. Guarnieri, S. Mancon, S. Tebaldini, in FRINGE Workshop (2015), pp. 1–6.

63. C. Fernández, "Sentinels POD service file format specifications" (Tech. Rep. GMESGSEG-EOPG-FS-10-0075, European Space Agency, 2011).

64. J. H. Liu, J. Hu, Z. W. Li, J. J. Zhu, Q. Sun, J. Gan, A method for measuring 3-D surface deformations with InSAR based on strain model and variance component estimation. *IEEE Trans. Geosci. Remote Sens.* **56**, 239–250 (2018).

65. K. Mosegaard, A. Tarantola, Monte Carlo sampling of solutions to inverse problems. *J. Geophys. Res. Solid Earth* **100**, 12431–12447 (1995).

66. J. Savage, W. Prescott, Asthenosphere readjustment and the earthquake cycle. *J. Geophys. Res.* **83**, 3369–3376 (1978).

67. F. Q. Diao, X. Xiong, R. J. Wang, T. R. Walter, Y. B. Wang, K. Wang, Slip rate variation along the Kunlun Fault (Tibet): Results from new GPS observations and a viscoelastic earthquake-cycle deformation model. *Geophys. Res. Lett.* **46**, 2524–2533 (2019).

68. J. S. Tong, T. Ma, K. R. Shen, H. Y. Zhang, S. P. Wu, A criterion of asphalt pavement rutting based on the thermal-visco-elastic-plastic model. *Int. J. Pavement Eng.* **23**, 1134–1144 (2022).

69. K. L. Wang, Y. J. Zhu, E. Nissen, Z. K. Shen, On the relevance of geodetic deformation rates to earthquake potential. *Geophys. Res. Lett.* **48**, e2021GL093231 (2021).

70. S. Heimann, M. Isken, D. Kühn, H. Sudhaus, A. Steinberg, S. Daout, S. Cesca, H. Bathke, T. Dahm, Grond: A probabilistic earthquake source inversion framework (2018).

71. GEOFON Data Centre, "GEOFON seismic network" [dataset] (Deutsches GeoForschungsZentrum GFZ, 1993).

72. G. Asch, A. Mohsen, B. Braeuer, R. Hofstetter, D. Jaser, R. El-Kelani, M. Weber, "DESIRE - Seismology (2006/2008)" [dataset] (GFZ Data Services, 2006).

73. E. Sandvol, "Continental dynamics/central anatolian tectonics: Surface to mantle dynamics during collision to escape" [dataset] (International Federation of Digital Seismograph Networks, 2013).

74. I. Bastow, C. Ogden, S. Pilidou, I. Dimitriados, P. Iosif, C. "Constantinou, TROODOS: Tomography and receiver function observations of an Ophiolite using data obtained from seismology" (International Federation of Digital Seismograph Networks, 2016).

75. "KO: Kandilli Observatory And Earthquake Research Institute (KOERI)" [dataset] (International Federation of Digital Seismograph Networks, 1971).

76. R. J. Wang, A simple orthonormalization method for stable and efficient computation of Green's functions. *Bull. Seismol. Soc. Am.* **89**, 733–741 (1999).

77. C. Bassin, The current limits of resolution for surface wave tomography in North America. *Eos. Trans. AGU* **81**, F897 (2000).

**Acknowledgments:** We thank two anonymous reviewers for helpful comments that improved the manuscript. **Funding:** This study was supported by King Abdullah University of Science and Technology (KAUST) awards BAS/1/1353-01-01 and OSR-CRG2020-4335. **Author contributions:** X.L. processed the BOI data, interpreted the results, and drafted the manuscript. S.J. initiated and supervised the paper, interpreted the results, and edited the manuscript. S.L. contributed to the modeling section. Z.M. contributed to the BOI data processing. N.C.-P. compiled the GNSS velocities. S.C. estimated the moment tensor solutions. F.M. and Y.K. contributed to the editing of the manuscript. All authors discussed the results and contributed to the final manuscript. **Competing interests:** The authors declare that they have no competing interests. **Data and materials availability:** All original Sentinel-1 SAR data are available at the Copernicus Open Access Hub (https://dataspace.copernicus.eu/) and GNSS velocities are from published studies (12, 14, 15, 35, 36, 47). The seismic data are available from EMSC (https://www.emsc-csem.org/). All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. The geodetic datasets and the synthetic earthquake catalog can be downloaded (DOI: 10.5281/zenodo.10158386).

Submitted 19 July 2023

Accepted 9 February 2024

Published 15 March 2024

10.1126/sciadv.adj8408