Figure 2
Figure 1
Figure 2
Figure 7
Figure 1
Figure 2
Figure 7
Figure 2
Figure 5
Figure 5
Figure 5
Figure 8
Figure 4a
Figure 6
Table 1
Table 2
The 1000-km-long left-lateral Dead Sea fault is a major tectonic structure of the oriental Mediterranean basin, bounding the Arabian Plate to the west. The fault is located in a region with an exceptionally long and rich historical record, allowing to document historical seismicity catalogues with unprecedented level of details. However, if the earthquake time series is well documented, location and lateral extent of past earthquakes remain often difficult to establish, if only based on historical testimonies. We excavated a palaeoseismic trench in a site located in a kilometre-size extensional jog, south of the Dead Sea, in the Wadi Araba. Based on the stratigraphy exposed in the trench, we present evidence for nine earthquakes that produced surface ruptures during a time period spanning 5000 yr. Abundance of datable material allows us to tie the five most recent events to historical earthquakes with little ambiguities, and to constrain the possible location of these historical earthquakes. The events identified at our site are the 1458 C.E., 1212 C.E., 1068 C.E., one event during the 8th century crisis, and the 363 C.E. earthquake. Four other events are also identified, which correlation with historical events remains more speculative. The magnitude of earthquakes is difficult to assess based on evidence at one site only. The deformation observed in the excavation, however, allows discriminating between two classes of events that produced vertical deformation with one order of amplitude difference, suggesting that we could distinguish earthquakes that started or stopped at our site from earthquakes that potentially ruptured most of the Wadi Araba fault. The time distribution of earthquakes during the past 5000 yr is uneven. The early period shows little activity with return interval of ~500 yr or longer. It is followed by a ~1500-yr-long period with more frequent events, about every 200 yr. Then, for the past ~550 yr, the fault has switched back to a quieter mode with no significant earthquake along the entire southern part of the Dead Sea fault, between the Dead Sea and the Gulf of Aqaba. We computed the Coefficient of Variation for our site and three other sites along the Dead Sea fault, south of Lebanon, to compare time distribution of earthquakes at different locations along the fault. With one exception at a site located next to Lake Tiberias, the three other sites are consistent to show some temporal clustering at the scale of few thousands years.
The 1000-km-long left-lateral Dead Sea fault (Fig. 1) marks the boundary between the Arabian Plate and the Sinai block (Quennell 1959; Freund et al. 1968; Garfunkel et al. 1981; Klinger et al. 2000a). The fault is mostly linear, trending N5◦E–N10◦E, with only a few noticeable extensional steps such as the Hula basin (Heimann & Ron 1993) and the Dead Sea basin (Garfunkel 1981; Ben Avraham 1992), and one major bend in its central section, in Lebanon. There, the Dead Sea fault bends eastward about 20◦, to form a major transpressional step with the motion partitioned between the on-shore Yammouneh strikeslip fault and the off-shore thrust faults (Daeron ¨ et al. 2005; Daeron ¨ et al. 2007; Elias et al. 2007). Mount Lebanon, resulting from this bend, tops at 3090 m. Along the central and southern sections, the slip rate has been estimated at several places, and over various timescales, from few years to several hundred thousands of years. Results converge around an average value of 5 ± 2 mm yr-1 (Klinger et al. 2000a; Niemi et al. 2001; Daeron ¨ et al. 2004; Marco et al. 2005; Gomez et al. 2007; Le Beon ´ et al. 2008, 2010, 2012; Makovsky et al. 2008; Ferry et al. 2011). North of Lebanon, however, slip rates estimated over different timescales, from present-day using GPS to historical times, diverge significantly, with values ranging from 2 to 6 mm yr-1 (Meghraoui et al. 2003; Alchalbi et al. 2010).
The trenching site (29◦ 48' 46.4" N, 35◦ 04' 26.3" E) is located in Jordan, 30 km north of the town of Aqaba, close to a village named Qatar after which the site is named. It is located at the southern tip of the Yotvata extensional relay-zone (Figs 1 and 2). The relay is about 1.5 km wide and 5 km long. Sand dunes have invaded part of the jog, obscuring the geometry of the northern tip of the basin and the exact continuation of the fault northward. The central part of the extensive jog presents the typical morphology of a dry lake in arid environment, with a flat bottom covered by a thin crust of salt-rich sediments (Figs 2 and 3). South of the Yotvata relay-zone, the fault clearly disrupts several alluvial surfaces (Le Beon ´ et al. 2012). Where the fault enters the jog, it is marked by a series of aligned pressure ridges (Fig. 3), about 1 m high and 20 m wide, which progressively disappear when getting into the basin. Further north into the basin, the location of the fault can be inferred from the location of two pressure ridges, which top about 50 cm above the flat bottom of the playa (Fig. 3), and from the offset of palaeoshorelines associated with the lake (Fig. 2). Southwest of the main fault strand entering the jog, the lake is possibly limited by normal faults with a minor vertical throw.
The trench is about 25 m long and 4 m wide. The depth, limited by the water table, is about 3 m. The azimuth of the trench, perpendicular to the local azimuth of the fault, is N110◦E. Beneath a salty crust, the stratigraphy is characterized by a series of individualized layers, white to dark brown in colour (Figs 4 and S1). The different strata consist in medium to coarse sand with very few clasts and in silts. The thickness of the layers varies from a couple of centimetres to few decimetres. Units tend to be thicker and more homogeneous in the lower part of the trench. Conversely, the upper part of the trench is characterized by thinner and more discontinuous deposits, resulting from units interfingering. Several units can be recognized across a large part of the trench, which are labelled A to H (Fig. 5). Note that E and H correspond to packages of several related units rather than to a single unit. These units are used as reference to discuss event horizons in the following sections.
The identification of earthquake ruptures is mostly based on consistent upward terminations of sets of cracks that vertically offset sedimentary units (McCalpin 1996). Horizontal motion is more difficult to recognize when sedimentary unit thickness is similar over a large area, as expected in playa deposits. In addition, in our trench the medium to coarse grain size of many sandy units makes it more difficult to definitely identify the exact location of some of the cracks termination (Fig. S1). Hence, our analysis focuses only on the large events that have extensively ruptured the ground surface at the trench site, and are responsible for primary deformation and/or cracking. Significant ground ruptures involve that the magnitude associated to these events is most probably Mw ∼ 6.5 or above (McCalpin 1996), when the entire thickness of the brittle crust would rupture. In this respect, as we benefitted from rich stratigraphy and do not have massive homogeneous layer in excess of 10–15 cm thickness, at least in the main fault zone, it is very unlikely that we could miss any major event. Smaller magnitude earthquakes responsible for limited surface ruptures, or secondary faulting related to distant earthquakes, might also be visible in the trench. At this stage, however, the lack of congruent observations in the trench prevents us from discussing them further in details.
Time constraints were derived from accelerator mass spectrometry (AMS) radiocarbon dating of detrital charcoals. The trench exposure proved to be unexpectedly rich in datable material and a total of 39 samples were dated (Tables 1 and 2). Only four samples had to be discarded, as they likely correspond to reworked samples yielding ages significantly older than other ages, relative to their stratigraphic position. The regular distribution of ages down to unit G (Fig. 6) confirms, in support to the logged section, that no major hiatus affects the stratigraphy and that the sedimentological record seems complete, especially for the last 2000 yr. For the earlier time period, because the contact between units H and G is erosional, up to 700 yr may be missing in the stratigraphy, making the earthquake record likely incomplete. The layers at the base of the trench date back to nearly 5 kyr BP.
Overall, the Qatar record displays evidence for nine unambiguous paleoearthquakes and two additional events, more uncertain, over a period extending back to about 5 kyr BP. Within this time span, we have reasonable time constraints on an ‘historical’ sequence for the last 1700 yr, and fewer indications about the earthquake activity during the earlier ‘prehistoric’ sequence. Hence, dates for individual events E1–E3, E6 and E7, are well constrained with several radiocarbon samples (14C) between each individual event horizon. The other identified events (E4 and E5) and (E8 and E9) interestingly appear to occur somehow by pair. The age of each pair is well constrained by several samples. If such pairing is a bias due to fewer samples in the lower section, or is real remains unclear. The latter possibility is explored in the discussion section.
In the Middle East, during the last 2000 yr, most of the significant earthquakes have been well recorded in historical accounts. However, the use of macroseismic data only has proved to be quite insufficient to locate specific earthquakes accurately, and additional paleoseismic data have often been decisive (Daeron ¨ et al. 2005, 2007). The situation of the southern section of the Dead Sea fault, in this respect, is even more difficult as the number of cities that could potentially report damages was always limited due to very arid climate. The new paleoseismological data exposed in our trench allow to identify several of the classic historical earthquakes and to constrain their area of rupture (Fig. 7). Historical earthquakes of 1458 C.E., 1212 C.E., 1068 C.E., one major earthquake during the 8th century seismic crisis of 746 C.E.–757 C.E., and the 363 C.E. event all ruptured the southern section of the Dead Sea fault, in the Wadi Araba. Interesting as well, some historical earthquakes have not been found at our trench site, meaning that they most probably ruptured other sections of the Dead Sea fault. The 1834 C.E. earthquake was widely felt in Palestine and caused some destruction in Jerusalem. Destructions are also reported East of the Dead Sea. Interestingly, it has been reported that large amount of asphalt was cast on the shore of the Dead Sea, suggesting that the earthquake could have been located in the Dead Sea or in direct vicinity (Ambraseys 2009; Agnon 2014). A ground rupture at the site of Tilah could also be correlated to this event although this correlation is far from being univocal (Haynes et al. 2006). Hence, for the 1834 event the location remains loosely constrained (Ambraseys 2009; Agnon 2014). The 1588 C.E. earthquake, which mostly affected the city of Aila and produced damage as far as Cairo in Egypt (Ambraseys 2009) was most likely limited to the Gulf of Aqaba (Agnon 2014), with no rupture propagating along the Wadi Araba section, similarly to the 1995 Nuweiba earthquake. The 1546 C.E. earthquake is also not found in our trench. This event is classically described as an event that had affected the northern part of the Dead Sea basin. The southern extent of the rupture, however, is difficult to ascertain, as damage reports seem to have been grossly exaggerated (Ambraseys & Karcz 1992; Ambraseys 2009). Haynes et al. (2006) reported a ground rupture about 15 km south of the Dead Sea, at the Tilah site, that they associated to the earthquake of 1546 C.E. If this were to be confirmed, it would probably correspond to the southern end of the rupture. If conversely, as suggested by (Agnon 2014), the last rupture at Tilah is 1834 C.E., then the 1546 C.E. rupture would be limited to the Dead Sea basin only with no extension in the Wadi Araba. Finally, the 1293 C.E. earthquake is not found either in our trench. This is consistent with the fact that damages associated with this earthquake are mostly located north of the Wadi Araba, with the citadel of Al-Karak severely destroyed. More surprisingly, the 1293 C.E. rupture is not identified by Haynes et al. (2006), suggesting that this event would have been only limited to the Dead Sea basin, with no rupture along the Wadi Araba segments.
The number of strike-slip faults where the earthquake record is long enough to be able to study fault seismic behaviour remains very small. Although attempts have been numerous, sites that preserved a long record of past earthquakes and enough material suitable for dating are difficult to find. Here, we present a site located at the southern end of the on-land part of the Dead Sea fault where we could identify up to nine events grouped in eight bursts of activity, over a period of about 5000 yr. Several of these earthquakes can be tied to historical earthquakes in chronicles, but for which location was remaining uncertain until now. The series of earthquakes recorded at our site, along with earthquakes recorded at three other sites further north along the Dead Sea fault, show some variability in the time distribution of earthquakes, with period of intense activity over few centuries and period of seismic quiescence. In addition, the fact that we observe several pairs of earthquakes, which occurred only few years to few decades apart, suggests that during periods of seismic activity, rupture of one section of the Dead Sea fault might trigger rupture on other nearby fault sections, leading to the total rupture of the Dead Sea fault in a relatively short time of few centuries. Such pattern of temporal clustering is not unique to the Dead Sea fault and it has been unambiguously evidenced along other large strike-slip fault such as the North Anatolian Fault, in Turkey (Barka 1996; Stein et al. 1997), and the San Jacinto fault, in United States (Rockwell et al. 2014). In contrast, however, other records along the Alpine fault, in New Zealand (Berryman et al. 2012) and the San Andreas fault, in United States (Scharer et al. 2010), show that on these faults the seismicity might be quasi-periodic in time. Hence, time clustering does not seem to be an intrinsic behaviour of strike-slip faults. Because one could not do any obvious correlation between a specific earthquake behaviour and a simple property of the system, such as fault geometry or isolation of the fault system relatively to other active faults, we suggest instead that it is more likely that the key to temporal clustering might depend on some relation between slip-rate and the existence of a time- dependent rheology (a rheology that would keep memory of earthquake history) in the lower crust and the upper mantle, as suggested in several model-driven approaches (Ben-Zion et al. 1999; Lyakhovsky et al. 2001; Kenner & Simons 2005; DiCaprio et al. 2008).