# Potentially large post-1505 AD earthquakes in western Nepal revealed by a lake sediment record

Z. Ghazoui, S. Bertrand, K. Vanneste, Y. Yokoyama, J. Nomade, P. A. van

Der Beek

► To cite this version:

Z. Ghazoui, S. Bertrand, K. Vanneste, Y. Yokoyama, J. Nomade, et al. Potentially large post-1505 AD earthquakes in western Nepal revealed by a lake sediment record. Nature Communications, 2019, 10, pp.621-636. 10.1038/s41467-019-10093-4. insu-03595785

HAL Id: insu-03595785

https://insu.hal.science/insu-03595785

Submitted on 3 Mar 2022

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
<figure><img src="image_2.png" /></figure>

ARTICLE

https://doi.org/10.1038/s41467-019-10093-4

OPEN

# Potentially large post-1505 AD earthquakes in western Nepal revealed by a lake sediment record

Z. Ghazoui ID <fnref n="1" />,², S. Bertrand ID ², K. Vanneste³, Y. Yokoyama ID ⁴, J. Nomade ID <fnref n="1" />, A.P. Gajurel⁵ & P.A. van der Beek<fnref n="1" />

According to paleoseismological studies, the last earthquake that ruptured the Main Frontal Thrust in western Nepal occurred in 1505 AD. No evidence of large earthquakes has been documented since, giving rise to the concept of a seismic gap in the central Himalaya. Here, we report on a new record of earthquake-triggered turbidites from Lake Rara, western Nepal. Our lake-sediment record contains eight possibly moderate-to-large earthquake-triggered turbidites during the last 800 years, three of which overlap in age with previously reported $M_w \ge 7$ events in western Nepal. Shaking intensity modelling, together with instrumental records, suggests that near-field earthquakes ($\le 15$ km) should have a minimum $M_w$ 5.6, and regional earthquakes ($\le 80$ km) a $M_w$ > ~ 6.5, to trigger turbidites. We present a likely scenario that western Nepal may be as seismically active as central Nepal; however, more data are needed to revaluate the seismic risk in the central Himalaya.

<footnote n="1">¹ISTerre, Institut des Sciences de la Terre, Université Grenoble Alpes, CNRS, 38058 Grenoble, France. ²RCMG, Renard Centre of Marine Geology, Ghent University, 9000 Ghent, Belgium. ³Department of Seismology and Gravimetry, Royal Observatory of Belgium, 1180 Brussels, Belgium. ⁴Atmosphere and Ocean Research Institute, The University of Tokyo, 277-0882 Chiba, Japan. ⁵Department of Geology, Tribhuvan University, 44618 Kathmandu, Nepal. Correspondence and requests for materials should be addressed to Z.G. (email: zakaria.ghazoui@ugent.be)</footnote>
The Himalayan collision, in which India underthrusts below Tibet and the Himalaya along a major crustal detachment known as the Main Himalayan Thrust (MHT), regularly produces major destructive earthquakes, as elastic deformation accumulated during underthrusting of the Indian plate is released periodically by slip along the MHT fault plane¹⁻⁴. The destructive 2015 M<sub>w</sub> 7.8 Gorkha earthquake⁵,⁶ represented an intermediate-size event in this process as it ruptured only the lower, northern part of the MHT without breaking through to the surface⁴,⁵.

Evaluation of seismic hazard in the Himalaya has been based on the comparison of geodetic strain rates with the amount of seismic moment release, measured by instrumental seismicity and inferred from paleoseismology¹,⁷. Whereas historical seismicity and trenching studies record at least four, and possibly up to eight, major earthquakes over the last 800 years in central and eastern Nepal⁸, the last known major event to have affected western Nepal and northern India, rupturing a long portion of the Main Frontal Thrust (MFT), was the M<sub>s</sub> ~8.2 earthquake of 1505 AD⁹,¹⁰ (Fig. 1). The intervening 500 years have resulted in the accumulation of a >10 m slip deficit along this segment of the MHT, leading to the concept of a seismic gap in western Nepal

and adjacent areas in northern India, which could potentially
trigger a great earthquake in the near future¹,⁷,¹¹.

However, both historical and paleoseismic records of earthquakes are inherently incomplete<sup>8,12,13</sup> and the age, extent, and correlation of surface ruptures inferred from paleoseismic trenching studies in the central Himalaya are the subject of significant controversy<sup>6,8,14–17</sup>.

Historical written archives in western Nepal are limited due to
several factors. For the period prior to the 20th century, references
to earthquakes in western Nepal and the northern Himalaya are
scattered throughout Tibetan literature and, until the mid-20th
century, access to Tibet and Tibetan documents was extremely
restricted to foreign researchers<sup>12</sup>. Moreover, translation and
interpretation of traditional references to earthquakes is compli-
cated by the religious undertone in the Tibetan hagiographic
literature<sup>12</sup>. In addition, western Nepal has always been much less
densely populated than central Nepal, lacking major urban cen-
tres (such as Kathmandu, Pokhara or Gorkha) and it was virtually
inaccessible to foreign researchers during the Nepal civil war of
1996–2006. Furthermore, seismic activity based on written
archives is exclusively evaluated in terms of macroseismic effects

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

**Fig. 1** Geomorphic and seismotectonic setting of Lake Rara. **a** ASTER Global Digital Elevation model (GDEM) of the central Himalaya, showing active faults³⁷,⁵⁷, instrumental seismicity (M<sub>w</sub> ≥ 6) from 1908-2017 (2017 update of ISC-GEM³⁴), and inferred rupture lengths of historical earthquakes (dates as indicated). Squares indicate paleoseismological trenching sites, in red for the great 1505 AD earthquake, in black for the Tibrikot earthquake and in blue for central Nepal earthquakes (see main text for references). Green dots are locations for which historical chronicles record destruction during the 1505 AD earthquake⁵⁸. Instrumental seismicity is coloured according to magnitude, earthquakes discussed in the text are indicated by date. Active faults are: KF—Karakoram Fault; MFT—Main Frontal Thrust; TG—Thakkola Graben; WNFS—Western Nepal Fault System (composed in part by: DF—Dharma Fault; GMHF—Gurla Mandhata-Humla Fault system; TF—Tibrikot Fault). The red rectangle shows the location of Lake Rara (inset **b**). **b** Bathymetric map of Lake Rara (from ref. ⁴⁰ and our own measurements) showing the sampling sites A and B, overlain on a satellite image of the lake catchment (outlined by yellow dashed line, the full catchment is shown in Supplementary Fig. 1; Landsat-7 image courtesy of the U.S. Geological Survey). ASTER GDEM is a product of NASA and METI
and is subject to misinterpretations¹², as the vulnerability of buildings exposed to earthquakes varies greatly in the Himalayan region depending on the type of architecture. In particular, traditional timber-laced buildings common to the western Himalaya (including western but not central Nepal)¹⁸,¹⁹ have proven their ability to resist earthquake loading much better than modern structures, but their exceptional resistance is not taken into account in any of the intensity scales¹²,¹十八,¹⁹.

Trench-based paleoseismic records are rendered equivocal by the strong vegetation cover and erosional activity of the Himalayan front, leading to poor access and preservation potential of fault scarps. Furthermore, the unknown and variable age inheritance in the charcoals used for radiocarbon dating of observed surface ruptures, i.e., the time lapse between charcoal formation and its incorporation in the sediments from which it is sampled²⁰, may lead to significant overestimation of the ages of seismic events⁸.

Finally, trenches on the Himalayan front only record surface-breaking earthquakes. However, not all earthquakes produce slip that reaches the surface (e.g., the 2015 M<sub>w</sub> 7.8 Gorkha earthquake<sup>5,6)</sup> and, therefore, would be detected from trenches. Moreover, while some ruptures reach the surface at the MFT, producing fault escarpments or fault-related folds, others (e.g., the 2005 M<sub>w</sub> 7.6 Kashmir earthquake)⁴ reach the surface through out-of-sequence thrusting.

Compared to the historical and terrestrial archives discussed above, lake sediments may provide a complementary and more continuous paleoseismic record, as earthquake-triggered slope failures and surficial sediment remobilization form turbidite deposits<sup>21-27</sup>. In spite of their potential to complement the paleoseismic inventory, lake records have hitherto not been investigated in the Himalaya.

In order to reconstruct past earthquake activity in the inferred seismic gap, we collected three short sediment cores from Lake Rara in western Nepal (29°32'N, 82°05'E; Fig. 1; Supplementary Fig. 1). We cored two sites within the lake, selected on the basis of a bathymetry survey (See Methods), in the deepest basin at water depths of 168 m (site A) and in the shallower northeastern arm of the main basin at 120 m depth (site B), using a gravity corer operated from an inflatable dingy (Fig. 1). The cores were analysed using X-ray Computerized Tomography (CT), X-ray fluorescence (XRF) core scanning, logging of physical properties, bulk organic geochemistry, and high-resolution grain-size measurements (Fig. 2; Supplementary Figs. 2, 3; see Methods). Chronology was established on cores RA14-SC05 and RA14-SC06 by combining radionuclide (²¹⁰Pb and ¹³⁷Cs) and radiocarbon (¹⁴C) dating (Supplementary Fig. 4, Supplementary Tables 1–3). Radiocarbon dating was performed on terrestrial leaf material and thus does not suffer from age inheritance. Our lake-sediment record contains eight earthquake-triggered turbidites during the last 800 years, five of which were previously undocumented and may represent large seismic events. These results suggest that western Nepal may be as seismically active as central Nepal and call for a revaluation of the risk of a major earthquake affecting western Nepal and northern India.

## Results

**Lake Rara turbidites.** The cored sediments are mainly composed of mud, interrupted by 10–25 mm thick, dense, magnetic and Ti-rich micaceous sandy silt layers (Fig. 2), which are characterized and differentiated on the basis of variations in mean grain size, magnetic susceptibility, Ti concentration, C/N ratios and radio-density (see Methods). As their distinctive signatures on the magnetic susceptibility and XRF profiles clearly contrast with the homogeneous muddy background, these fining-upward deposits were readily identified as turbidites (see Supplementary text). Up

to eight of these deposits were identified in the cores; they are most clearly represented in core RA14-SC05 (Fig. 2; Supplementary Fig. 2). Age-depth models (Supplementary Fig. 4; Methods) show that sedimentation rates are ~0.4 mm/yr at site A, located in the deepest part of the lake (168 m), and approximately half that at site B. Both the lake physiography and the sedimentation rates at site A favour the ability of this site to record a complete event history<sup>23,24</sup>; we, therefore, focus here on the interpretation of core RA14-SC05.

**Origin of turbidites.** Turbidites within lake sediments can be triggered by various processes such as floods, landslides, spontaneous slope failures, or earthquakes (see Supplementary text). The limited hydrographic system (Fig. 1, Supplementary Fig. 1), with short (the largest stream length is ~4 km from source to lake) and small (streams entering the lake are typically less than 30 cm wide and 15 cm deep) streams, and the relatively low relief of the Lake Rara catchment, render floods and landslides unlikely triggering mechanisms. This inference is consistent with the very low sedimentation rates measured at the shallowest site RA14-SC06 (~0.2 mm/yr), which reflects the low hydrodynamic activity in the catchment of Lake Rara, and render spontaneous slope failures also unlikely. Hence, the turbidites identified in the sediments of Lake Rara most likely represent either earthquake-triggered slope failures or earthquake-triggered remobilization of surficial lake sediments. This interpretation is supported by their C/N ratios, which range between 13 and 17 (Fig. 2), reflecting a mixture of organic matter of aquatic (C/N < 8) and terrestrial (C/N > 20) origin<sup>22</sup>. Although the absolute values depend on grain size (Supplementary Fig. 3), these results suggest that the turbidites originated from reworking of sediment previously deposited on shallower parts of the slopes within the lake, where C/N values are slightly higher than in the centre of the lake due to continuous input of terrestrial organic matter from the catchment. Flood or landslide-triggered turbidites, in contrast, would have a purely terrestrial C/N signature<sup>22</sup>.

The strongest argument for attributing an earthquake origin to lake turbidites is a temporal correlation with known historical events<sup>22,23,26</sup>. A good starting point for this correlation is the great 1505 AD earthquake. Historical records describe widespread devastation in western Nepal, northern India and southern Tibet (Fig. 1), from which a magnitude M<sub>s</sub> 8.2 has been estimated<sup>12</sup>. Surface ruptures on the MFT immediately south and southwest of Lake Rara have been attributed to the 1505 event<sup>9,10</sup>. This earthquake, therefore, seems an obvious candidate to have triggered slope failures within Lake Rara. The age range of turbidite T7 in core RA14-SC05 (1399–1570 AD) overlaps the timing of the 1505 AD earthquake. This is the most prominent event deposit; it can be correlated to turbidite TA in core RA14-SC06 from site B based on the independent age models (Supplementary Figs. 4, 5), and it is present in all three cores (Supplementary Fig. 2). We suggest that these events have been caused by the great 1505 AD earthquake.

Another known earthquake that affected the region of Lake Rara is the 1916 M<sub>w</sub> ~7.0–7.2 Dharchula earthquake<sup>28</sup>. The ages of the two topmost turbidites in our record (T1 and T2; Fig. 2; Supplementary Fig. 4) both overlap with this earthquake, with the weighted-mean age of T2 (1917) being very close to it. Therefore, it appears likely that the 1916 earthquake triggered either T1 or T2, but the age overlap between these does not allow us to pick one or the other (Supplementary Fig. 4). The association of these two turbidites with known earthquakes supports the seismic nature of the lake turbidites recorded in Lake Rara. The others (Fig. 2) thus most likely represent previously undocumented earthquakes.
<figure><img src="image_4.png" /><figcaption>**Fig. 2** Identification of turbidites in sediment core RA14-SC05. From left to right are shown: a synthetic log, 3D coronal view of computed X-ray tomography images in pseudo-colour, mean grain size, magnetic susceptibility, Itrax XRF Ti profile (grey curve: raw data; orange curve: 5 pt-weighted average) and carbon to nitrogen (C/N) atomic ratio. These five parameters were used to identify turbidites T1 to T8. Their ages were extracted from the age-depth model (see Methods)</figcaption></figure>

**Earthquake turbidite-triggering threshold.** Observations in lakes worldwide suggest that earthquake-triggered turbidites can originate either from subaqueous slope failures or from surficial sediment remobilisation, when local shaking exceeds Modified Mercalli Intensities (MMI) 6–7 and 6, respectively²¹⁻²⁴,²⁶,²⁷. In order to predict shaking intensity and constrain the local earthquake turbidite-triggering threshold (E<sub>QTT</sub>)²⁴ for Lake Rara, we produced MMI shaking maps for several historical and instrumental events (Fig. 3) based on a set of intensity-prediction equations²⁸⁻³³ (IPEs; see Methods).

The E<sub>QTT</sub> for Lake Rara was estimated using two instrumental earthquakes; the 1980 M<sub>w</sub> 6.5 Bajhang and 1916 M<sub>w</sub> ~7.0–7.2 Dharchula earthquakes. The most recent update of the ISC-GEM catalogue³⁴ reassessed the location and magnitude of the Dharchula earthquake at an epicentral distance of ~130 km from Lake Rara (29.730°N 80.745°E; Supplementary Fig. 8). As discussed above, this earthquake can be correlated to T1 or T2 in the sediments of Lake Rara. The 1980 M<sub>w</sub> 6.5 Bajhang earthquake³⁴ (80-km epicentral distance from Lake Rara) constitutes the largest among six M<sub>w</sub> 6.2–6.5 instrumental earthquakes recorded in the 1974–2017 USGS catalogue within an 80-km radius of Lake Rara (Supplementary Fig. 8). Since their ages are much younger than the most recent turbidite T1 in Lake Rara (Fig. 2), we infer that those six M<sub>w</sub> 6.2–6.5 earthquakes did not trigger slope failures or surficial sediment remobilization in the lake. The E<sub>QTT</sub> for Lake Rara should, therefore, lie between the shaking intensities felt at the lake for these two large earthquakes. Mean shaking intensities obtained with our IPE selection suggest that the E<sub>QTT</sub> is situated between MMI ~4.5 and MMI ~5.5 (Fig. 3a, b). However, comparison with observed intensities for the Bajhang earthquake²⁸ suggests that our modelled near-field intensities may underestimate actual intensities by a half to one unit (see Methods; Fig. 3c, d). From the IPEs in our selection, it appears that the IPE calibrated for India and the Himalaya²⁸ shows the best match with the observed intensities²⁸ in the distance range of interest to this study. Hence, with the available data and the associated uncertainties, we

cannot tightly constrain the E<sub>QTT</sub> for Lake Rara, although a value close to the lower bound of the range reported in the literature (MMI 6) is likely. In the following, we will, therefore, use an indicative estimate of the local E<sub>QTT</sub> in the range of MMI 5.3–5.7, based on the Himalayan IPE²⁸.

**Possible correlation with other historical earthquakes.** Attributing T1 or T2 to the M<sub>w</sub> ~7.0–7.2 1916 Dharchula earthquake and turbidite T7 to the 1505 M<sub>s</sub> ~8.2 earthquake implies that five other post-1505 AD events are recorded in Lake Rara sediments (Fig. 2), some of which may correspond to known earthquakes that occurred in Nepal or northern India during the last two centuries. In particular, the 1833/08/26 and 1934/01/15 earthquakes affecting central and eastern Nepal (Supplementary Fig. 6) are well documented both geologically and macroseismically (Supplementary Table 5). Both earthquakes ruptured the MHT east of Kathmandu at epicentral distances of ~365 km and ~545 km from Lake Rara, respectively, and have been attributed moment magnitudes M<sub>w</sub> ~7.3–7.7 and ~8.1–8.4, respectively⁸,²⁸,³⁰. Surface ruptures on the MFT were attributed to the latter event⁸,¹³,¹⁵,³⁰, although this interpretation has recently been questioned¹⁷. Published isoseismals for these earthquakes vary significantly but in most of these, Lake Rara lies outside the MMI = 6 isoseismal¹⁵,²⁸,³⁰. We have modelled the isoseismals for both earthquakes using the rupture planes and magnitudes inferred by ref. ³⁵ (Fig. 3, Supplementary Fig. 6). For the 1934 earthquake, modelled intensities at the lake are MMI < 5, whereas they are 5 < MMI < 6 for the 1833 earthquake. For this far-field earthquake, our predicted intensities correspond much better to the observed intensities than for the near-field earthquakes discussed above. We thus conclude that these modelled intensities are reasonable and that the 1934 event is not likely to have generated a turbidite in the lake, even though it is consistent with the age of turbidite T1. Either T1 or T2 thus corresponds to a previously unknown late 19th–early 20th-century event, which could have gone unnoticed if it was not felt in India. Whereas modern documentation of earthquakes (including instrumental
<figure><img src="image_5.png" /><figcaption>**Fig. 3** Modified Mercalli Intensity isoseismals modelled for historical and instrumental earthquakes in central/western Nepal, and plots of observed and predicted intensity vs. distance. Modelled isoseismals and observed intensities for instrumental earthquakes in western Nepal; **a** 1980 (MHT?) Bhajang earthquake (modelled *M*<sub>w</sub> 6.5; intensity observations<sup>28</sup> in map region are indicated), **b** 1916 (MHT?) Dharchula earthquake (modelled *M*<sub>w</sub> 7.2; no intensity observations in map region). **c, d** Plots of observed and predicted intensity vs. distance for the 1980 Bhajang (**c**) and 1916 Dharchula (**d**) earthquakes, respectively. The vertical blue line corresponds to the epicentral distance of Lake Rara; the dashed grey lines indicate the lower (**c**) and upper (**d**) bounds of our approximate local earthquake turbidite-triggering threshold. Red dots are observed intensities<sup>28</sup>, coloured lines show predicted intensity vs. distance according to the different IPEs<sup>28,29,31-33</sup> and their average. Also shown is USGS Shakemap fit for the 1980 Bhajang earthquake (see Methods for discussion). **e–g** Modelled isoseismals for historical earthquakes in central/western Nepal: **e** 1934 MHT earthquake (modelled *M*<sub>w</sub> 8.4), **f** 1833 MHT earthquake (modelled *M*<sub>w</sub> 7.7) and **g** 1165–1400 AD Tibrikot Fault (part of the WNFS) earthquake (modelled *M*<sub>w</sub> 7.9) The surface projection of the rupture plane is represented by the grey rectangle; the top edge of the rupture is shown by a solid white line. Rupture scenarios for **e** and **f** are from ref. <sup>35</sup>; for **g** from ref. <sup>37</sup>. The location of Lake Rara is indicated</figcaption></figure>

seismicity, newspaper reports, geological accounts) goes back to
this time in India, such documentation for western Nepal only
exists since the 1970's. The 1833 event could possibly have trig-
gered a turbidite in lake Rara, but none of the turbidites overlaps
with this age. Likewise, the 1803/09/01 Kumaon earthquake, at an
epicentral distance of ~385 km, has been attributed a magnitude
*M*<sub>w</sub> ~7.3–7.5<sup>28,30</sup>. Published isoseismals (Supplementary Fig. 7)
suggest an intensity of IV (MSK) for Lake Rara<sup>30</sup>, which is clearly
insufficient to trigger a turbidite. We did not attempt to model
isoseismals for this earthquake, given the expected low intensities
at Lake Rara.

The Lake Rara sediment record also includes at least one
medieval earthquake (T8 in core RA14-SC05; Fig. 2,

Supplementary Fig. 4). Although two or three great medieval
earthquakes were inferred from paleoseismic trenches in Nepal,
of which the historic 1255 earthquake has attracted the most
attention, their exact magnitude and extent of slip remain
controversial<sup>8,14–16</sup>. Surface ruptures attributed to the 1255
earthquake have recently been documented from locations to
the south and southwest of Kathmandu<sup>36</sup>. However, given the
apparent similarity of the 1255 and 1934 events<sup>8,35</sup>, we consider it
unlikely that T8 was triggered by the 1255 AD earthquake.

The Western Nepal Fault System (WNFS) has been recognised
as an active trans-tensional fault system that accommodates
oblique convergence in the western central Himalaya and
connects to the Karakorum Fault<sup>37,38</sup> (Fig. 1). Although the
<figure><img src="image_6.png" /><figcaption>**Fig. 4** Overview of historical and paleoseismic data along the Himalayan front. The map is modified and updated from ref. ⁵⁹. a Digital terrain model of the Himalayan region. Yellow areas correspond to inferred ruptures of historical earthquakes (the surface ruptures of the 1255 AD and 1833 AD earthquakes overlap that of the 1934 earthquake and are not represented); their timing and magnitude are indicated. Red dots labelled A to N represent paleoseismic trenches along the MFT; references are given in parentheses below plot (b). b Ages and extents of historical Himalayan $M_w$ > ~7.5 earthquakes (horizontal yellow bars). Red vertical bars reflect ¹⁴C ages that bracket the age of surface ruptures in calendar years AD, with 2σ uncertainties. Blue vertical bars represent the turbidite ages extracted from the age-depth model (2σ; this study); blue dots indicate the weighted mean age. Data from refs. 8-10,14,15,17,43,60,61</figcaption></figure>

activity of this system cannot be resolved from geodetic data³⁹, it presents features that are diagnostic of active faulting. In particular, a surface-rupturing event has been documented on the Tibrikot Fault (Fig. 1), tentatively dated between 1165 and 1400 AD and attributed to an $M_w$ ~7.9 earthquake³⁷. Our modelling shows that Lake Rara is within the MMI = 7 isoseismal for this rupture scenario (Fig. 3e). We, therefore, suggest that turbidite T8 (1135–1303 AD) records seismic shaking associated with the Tibrikot earthquake³⁷.

**Significance of previously unknown events.** Overall, five turbidites post-dating the great 1505 earthquake cannot be attributed to any previously described earthquake (Fig. 4). To evaluate whether or not background seismicity ($M_w$ < 5.6) could have

triggered these turbidites, we consider all earthquakes reported by the USGS within a radius of 20 km of Lake Rara between 1974 and 2017, the period during which the USGS catalogue is considered complete. Our search yielded 11 earthquakes with $M_w$ ranging from 4 to 5.6, none of which is recorded as a turbidite in Lake Rara. These results suggest that background seismicity, i.e., events $M_w$ < ~6, do not significantly remobilize surficial sediments to generate turbidites in Lake Rara. In addition, we tested the sensitivity of the lake site to near-field moderate to large earthquakes. In this test (see Methods), we used the aforementioned IPEs to infer the magnitude required to reach the inferred approximate (MMI 5.3–5.7) threshold intensity at Lake Rara for hypothetical mid- to upper-crustal moderate-magnitude events on the closest known faults (Supplementary Fig. 8): one on the MHT directly below the lake (at a depth of 26 km³⁹), and two
others on the Dharma (DF) and Hula (HF) Fault segments of the WNFS, at epicentral distances of <15 km and 80 km, respectively. The results (Supplementary Fig. 8) indicate that, in order to be recorded in Lake Rara sediments, nearby earthquakes should have a minimum magnitude $M_w \ge 5.3-5.6$ for the Rara-MHT and DF scenarios, with the required magnitude increasing rapidly with epicentral distance, and $M_w \ge 6.5-6.9$ for the HF scenario. Notwithstanding the uncertainty associated with the IPE's discussed previously, the upper range of these magnitudes is consistent with the observation that none of the reported near-field earthquakes with magnitudes $M_w \le 5.6$ produced a turbidite in Lake Rara, suggesting that the turbidites record earthquakes with magnitudes above the regional background seismicity during the last decades. The rapid increase in required magnitude with epicentral distance is by the fact that the 1980 $M_w$ 6.5 earthquake at an epicentral distance of ~100 km, and six instrumental $M_w \ge 6.2$ earthquakes in western Nepal and northern India, failed to trigger turbidites in the lake during the last 40 years.

**Implications for seismic activity in western Nepal.** The above findings suggest that, on average every 50–100 years since 1505, earthquakes with a minimum magnitude $M_w \ge 5.3-5.6$ have produced turbidites in Lake Rara. These minimum magnitudes are for near-field earthquakes occurring in the immediate vicinity (<15 km) of the lake. Regional earthquakes (≥80 km epicentral distance from the lake), should have magnitudes $M_w > 6.5$ in order to trigger turbidites in the lake. Such earthquakes could occur on either the MHT or the WNFS. In both cases, the earthquake slip would contribute to reducing the stored seismic moment in western Nepal, since the WNFS accommodates oblique convergence in the western-central Himalaya through slip partitioning$^{37,38}$. Our data suggest that western Nepal may be as seismically active as central Nepal (Fig. 4) and call for reconsideration of the notion of a seismic gap in western Nepal and northern India. Further testing the seismic-gap hypothesis will require studies involving longer time series and additional lacustrine paleoseismic investigations in this region. Future studies of other Himalayan lakes should allow correlating these events over a broader region. We also note the critical need for IPE calibration in this region in order to better quantify the most likely threshold intensity and magnitude to trigger lake turbidites. By using lake-sediment records as a paleo-seismometer for the first time in the Himalaya, this study complements the record obtained from paleoseismological and historical archives, demonstrating the importance of a holistic and diversified approach in paleoseismology to improve seismic hazard assessment.

## Methods

**Sediment core collection and analysis.** In order to reconstruct the past earthquake activity of the Lake Rara area, we collected three short sediment cores with an average length of 40 cm. The coring sites were selected after a preliminary bathymetric survey that complemented existing depth information$^{40}$. The coring sites are located in two different areas at water depths of 168 m and 120 m, respectively (Supplementary Figs. 1, 2). The cores were obtained in October 2014 with an Uwitec gravity corer operated from an inflatable dingy.

Cores were split lengthwise using a Geotek core splitter and described macroscopically. Their physical properties, including γ-ray attenuation density and magnetic susceptibility, were obtained using a Geotek multi-sensor core logger at 2-mm resolution. High-resolution grain-size analysis was performed on core RA14-SC05 with a step of 2.5 mm using a Malvern Mastersizer 3000 after removing organic matter, calcium carbonate and biogenic silica$^{41}$. Grain-size distribution parameters were obtained using Gradistat V8$^{42}$ and were calculated according to ref. 43. The relative concentration of major elements was obtained by X-ray fluorescence spectrometry (Itrax XRF core scanner) at 500-μm resolution at the Stockholm University Slamlab (Sweden). Bulk organic geochemistry (TOC, TN) was measured with an elemental analyser (PDZ Europa ANCA-GSL) at the UC Davis Stable Isotope Facility. Before analysis, samples were placed in silver

capsules and decarbonated using 1 N sulphurous acid$^{44}$. The amount of sediment was optimized based on preliminary LOI$_{550}$ measurements.

To reconstruct the 3D structure of the split cores we made use of a Siemens Somatom Definition Medical X-ray Computerized Tomography (CT) scanner$^{45}$, which produces CT-slice images composed of voxels (volume elements) with a resolution of 0.15 × 0.15 × 0.60 mm. The 3D reconstruction was obtained by processing the contiguous set of CT slices with VGStudio v2.1. Radio-density values were extracted from the grey levels of the CT-slice images and represented using a colour chart to highlight variations in density within the cores. The CT grey levels correspond to X-ray attenuation, reflecting the proportion of X-rays absorbed or scattered as they pass through each voxel, which is primarily a function of X-ray energy and the density and composition of the material being analysed.

**Age models.** A chronology for the Lake Rara sediment record was established on cores RA14-SC05 (site A) and RA14-SC06 (site B) by combining radiocarbon (Supplementary Table 1) and radionuclide ($^{210}$Pb and $^{137}$Cs; Supplementary Tables 2, 3) dating. Samples for radiocarbon dating were picked outside of the turbidites, as these are considered to be instantaneous deposits.

Radionuclides $^{210}$Pb and $^{137}$Cs were measured respectively on 11 and 9 bulk-sediment samples from core RA14-SC05 (Supplementary Fig. 4) by Flett Research Ltd. (Winnipeg, Canada) following the methods of refs. 46,47. Radiocarbon ages for cores RA14-SC05 and RA14-SC06 were obtained by dating four and five leaf samples, respectively. After acid-alkali-acid pre-treatment, the samples were converted to graphite following the procedure of ref. 48. Isotopic analysis of the graphite targets was performed using accelerator mass spectrometry (AMS) at the University of Tokyo, Japan$^{49}$. All ages were calibrated using the calibration curve for Northern Hemisphere terrestrial $^{14}$C dates IntCal13$^{50}$.

Age-depth models of cores RA14-SC05 and RA14-SC06 (Supplementary Fig. 2) were produced using Bacon 2.2 software$^{51}$ after removal of the turbidites. The ages indicated in Fig. 2 and Supplementary Fig. 4 correspond to the base of each turbidite

**Modelling shaking intensity.** In order to evaluate the potential impact of historical earthquakes on the Lake Rara sediment record, we used custom software built on top of the core Python library of the open-source seismic hazard engine OpenQuake$^{52}$ to compute shaking intensities based on the rupture parameters inferred for these earthquakes and on a set of intensity-prediction equations (IPEs). We evaluated five events: a surface-rupturing earthquake ($M_W \sim 7.9$) on the Tibrikot fault (Western Nepal Fault System; WNFS) between AD 1165 and 1400$^{37}$, the ruptures inferred in ref. 35 for the 1833 ($M_w \sim 7.3-7.7$) and 1934 ($M_w \sim 8.1-8.4$) earthquakes on the Main Himalayan Thrust (MHT), and the instrumental 1916 Dharchula ($M_w \sim 7.0-7.2$)$^{28,34}$ and 1980 Bajhang ($M_w 6.5$)$^{34}$ earthquakes inferred to involve the MHT. The rupture parameters are summarized in Supplementary Table 5. A large number of IPEs is available in the literature, predicting shaking intensity as a function of source and path parameters (mainly magnitude and distance). They differ in the intensity scale (MMI, EMS-98, MSK), magnitude scale ($M_W, M_S$) and distance metric (epicentral, hypocentral or rupture distance) that is considered, in the geographic area or tectonic region for which they are representative, and in the number and quality of input macroseismic data points. As it is currently not possible to identify a single IPE that is best suited to model shaking intensities for fault ruptures in western Nepal, we applied a mix of five different IPEs: we selected three IPEs$^{29,31,32}$ that were developed for California and other tectonically active regions globally, and performed best in an evaluation for application in the Global Shakemap$^{53}$, supplemented with two IPEs developed specifically for the Himalaya$^{28,33}$. The main parameters of these IPEs are summarized in Supplementary Table 4. Unfortunately, the Himalayan IPEs use a different intensity scale (EMS-98 and MSK, respectively, vs. MMI), but the difference with MMI is probably less than the variability among different IPEs, and one of these studies$^{28}$ made a direct comparison with the two Californian IPEs$^{29,31}$, concluding that there is good agreement in intensity attenuation between the two regions. Using the selected IPEs, we computed maps of mean shaking intensity (Fig. 3a–g) as the arithmetic average of the mean intensity predicted by the different IPEs. We note that, whereas the predicted and observed intensities for the 1934 and 1833 central Nepal earthquakes (Fig. 3e, f) appear to overlap reasonably well (Supplementary Fig. 6), significant discrepancies appear between predicted and observed intensities for the 1980 Bajhang earthquake (Fig. 3a). These discrepancies result both from poorly constrained locations (hypocentre, epicentre, depth) and magnitudes of the modelled events$^{34}$, as well as from the use of IPEs$^{28,29,31–33}$ that are not well calibrated for the region and possibly underestimate intensities (Fig. 3c, d). In particular for the 1916 Dharchula earthquake, the closest recording station used for constraining the source parameters was in Calcutta, at a distance of ~10° (>1000 km)$^{34}$. Our intensity map for the 1980 Bajhang earthquake, based on the IPE mix$^{28,29,31–33}$ appears to underestimate the observed near-field intensities$^{28}$ (Fig. 3c, f). In the USGS Shakemap solution, based on one of the IPEs in our mix$^{29}$, it appears that a value of 0.76 has been added to the event magnitude to obtain predicted intensities matching with the observations$^{28}$, explaining why the published Shakemap (https://earthquake.usgs.gov/ earthquakes/eventpage/usp0001959#shakemap) shows higher intensities at the site of Lake Rara than our solution. From the IPEs in our selection, it appears that the Himalayan IPE$^{28}$ shows the best match with the observed intensities in the distance
range of interest to this study (up to about 150 km), although this is not surprising as ref. 28 fitted their IPE to the same observations, which are in EMS-98. Using the IPE from ref. 28, the estimated earthquake-recording threshold of Lake Rara would lie between MMI = 5.3–5.7 (Fig. 3).

**Modelling sensitivity to near-field background seismicity.** To assess the sensitivity of Lake Rara to moderate-magnitude near-field earthquakes, we estimated the magnitude necessary to cause shaking of intensity MMI ≥ 5.3–5.7 at the lake for different scenarios using the same IPE mix as above. We computed intensities (mean ± 1 standard deviation) for a range of possible magnitudes, and interpolated the magnitude corresponding to MMI = 5.3–5.7 from the obtained magnitude-intensity curves (Supplementary Fig. 8). We considered three different rupture scenarios on the faults that are closest to Lake Rara: one on the MHT directly below Lake Rara (Rara-MHT), and two others on the Dharma (DF) and Humla (HF) Faults, which are both part of the Western Nepal Fault System (WNFS; Fig. 1, Supplementary Fig. 8). Instrumental seismicity in western Nepal mainly occurs around the MHT and is concentrated in two bands north and south of Lake Rara⁵⁴. These scenarios are modelled by selecting a fixed hypocentral location, which is the midpoint of ruptures that increase in size with the considered magnitude. In the Rara-MHT scenario, the hypocentre is directly below Lake Rara at a depth of 26 km³⁹. The DF and HF scenarios were modelled at locations corresponding to recent earthquakes, which occurred respectively 2008/12/08 (M<sub>w</sub> 5.3), and 1980/06/22 (M<sub>w</sub> 5.1). Their source parameters were taken from the Global Centroid Moment Tensor database⁵⁵,⁵⁶, and we selected the nodal plane that agrees best with the corresponding faults and style of faulting³⁸,⁵⁷. The rupture parameters are summarized in Supplementary Table 6. Because epicentral distance is zero in the Rara-MHT scenario, we decided to leave out IPEs that use epicentral distance, as this may lead to unrealistically high or undefined intensities for this case. In the other scenarios, the distance between the lake and the rupture decreases with increasing magnitude. The magnitude vs. intensity plots obtained for the three scenarios are shown in Supplementary Fig. 8.

**Data availability**

All the relevant data from this study are available from the corresponding author.

Received: 8 August 2017 Accepted: 3 April 2019

Published online: 21 May 2019

**References**

1. Bilham, R., Gaur, V. K. & Molnar, P. Himalayan seismic hazard. *Science* **293**, 1442–1444 (2001).

2. Avouac, J. P. Mountain building, erosion, and the seismic cycle in the Nepal Himalaya. *Adv. Geophys.* **46**, 1–80 (2003).

3. Feldl, N. & Bilham, R. Great Himalayan earthquakes and the Tibetan plateau. *Nature* **444**, 165–170 (2006).

4. Mugnier, J. L. et al. Structural interpretation of the great earthquakes of the last millennium in the central Himalaya. *Earth Sci. Rev.* **127**, 30–47 (2013).

5. Avouac, J.-P., Meng, L., Wei, S., Wang, T. & Ampuero, J.-P. Lower edge of locked Main Himalayan Thrust unzipped by the 2015 Gorkha earthquake. *Nat. Geosci.* **8**, 708–711 (2015).

6. Hayes, G. P. et al. Rapid characterization of the 2015 Mw 7.8 Gorkha, Nepal, earthquake sequence and its seismotectonic context. *Seism. Res. Lett.* **86**, 1557–1567 (2015).

7. Stevens, V. L. & Avouac, J. P. Millenary Mw > 9.0 earthquakes required by geodetic strain in the Himalaya. *Geophys. Res. Lett.* **43**, 1118–1123 (2016).

8. Bollinger, L., Tapponnier, P., Sapkota, S. N. & Klinger, Y. Slip deficit in central Nepal: omen for a repeat of the 1344 AD earthquake? *Earth, Plan. Space* **68**, 1–12 (2016).

9. Kumar, S. et al. Paleoseismic evidence of great surface rupture earthquakes along the Indian Himalaya. *J. Geophys. Res.* **111**, B03304–B03319 (2006).

10. Yule, D., Dawson, S., Lavé, J., Sapkota, S. & Tiwari, D. Possible evidence for surface rupture of the Main Frontal Thrust during the great 1505 Himalayan earthquake, far-western Nepal. Am. Geophys. Union Fall Meeting Abstracts #S33C–05 (2006).

11. Rajendran, C. P., John, B. & Rajendran, K. Medieval pulse of great earthquakes in the central Himalaya: Viewing past activities on the frontal thrust. *J. Geophys. Res.* **120**, 1623–1641 (2015).

12. Ambraseys, N. & Jackson, D. A note on early earthquakes in northern India and southern Tibet. *Curr. Sci.* **84**, 570–582 (2003).

13. Bilham, R. Earthquakes in India and the Himalaya: tectonics, geodesy and history. *Ann. Geophys.* **47**, 839–858 (2004).

14. Lavé, J., Yule, D., Sapkota, S., Basant, K. & Madden, C. Evidence for a great Medieval earthquake (~1100 AD) in the central Himalayas, Nepal. *Science* **307**, 1302–1305 (2005).

15. Sapkota, S. N. et al. Primary surface ruptures of the great Himalayan earthquakes in 1934 and 1255. *Nat. Geosci.* **6**, 71–76 (2012).

16. Pierce, I. & Wesnousky, S. G. On a flawed conclusion that the 1255 A.D. earthquake ruptured 800 km of the Himalayan Frontal Thrust east of Kathmandu. *Geophys. Res. Lett.* **43**, 9026–9029 (2016).

17. Wesnousky, S. G., Kumahara, Y., Nakata, T., Chamlagain, D. & Neupane, P. New observations disagree with previous interpretations of surface rupture along the Himalayan Frontal Thrust during the great 1934 Bihar-Nepal earthquake. *Geophys. Res. Lett.* **45**, 2652–2658 (2018).

18. Rai, D. C. & Murty, C. V. R. Effects of the 2005 Muzaffarabad (Kashmir) earthquake on built environment. *Curr. Sci.* **90**, 1066–1070 (2006).

19. Langenbach, R. in *Nonconventional and Vernacular Construction Materials* 3–26 (Woodhead Publishing, 2016).

20. Frueh, W. T. & Lancaster, S. T. Correction of deposit ages for inherited ages of charcoal: implications for sediment dynamics inferred from random sampling of deposits on headwater valley floors. *Quat. Sci. Rev.* **88**, 110–124 (2014).

21. Strasser, M., Anselmetti, F. S., Fäh, D., Giardini, D. & Schnellmann, M. Magnitudes and source areas of large prehistoric northern Alpine earthquakes revealed by slope failures in lakes. *Geology* **34**, 1005–1008 (2006).

22. Howarth, J. D., Fitzsimons, S. J., Norris, R. J. & Jacobsen, G. E. Lake sediments record high intensity shaking that provides insight into the location and rupture length of large earthquakes on the Alpine Fault, New Zealand. *Earth Planet. Sci. Lett.* **403**, 340–351 (2014).

23. Moernaut, J. et al. Lacustrine turbidites as a tool for quantitative earthquake reconstruction: New evidence for a variable rupture mode in south central Chile. *J. Geophys. Res.* **119**, 1607–1633 (2014).

24. Wilhelm, B. et al. Quantified sensitivity of small lake sediments to record historic earthquakes: Implications for paleoseismology. *J. Geophys. Res.* **121**, 2–16 (2016).

25. Kremer, K. et al. Lake-sediment based paleoseismology: Limitations and perspectives from the Swiss Alps. *Quat. Sci. Rev.* **168**, 1–18 (2017).

26. Monecke, K., Anselmetti, F. S., Becker, A., Sturm, M. & Giardini, D. The record of historic earthquakes in lake sediments of Central Switzerland. *Tectonophysics* **394**, 21–40 (2004).

27. Moernaut, J. et al. Lacustrine turbidites produced by surficial slope sediment remobilization: A mechanism for continuous and sensitive turbidite paleoseismic records. *Mar. Geol.* **384**, 159–176 (2017).

28. Szeliga, W., Hough, S., Martin, S. & Bilham, R. Intensity, magnitude, location, and attenuation in India for felt earthquakes since 1762. *Bull. Seism. Soc. Am.* **100**, 570–584 (2010).

29. Allen, T. I., Wald, D. J. & Worden, C. B. Intensity attenuation for active crustal regions. *J. Seismol.* **16**, 409–433 (2012).

30. Ambraseys, N. N. & Douglas, J. Magnitude calibration of north Indian earthquakes. *Geophys. J. Int.* **159**, 165–206 (2004).

31. Bakun, W. H. & Wentworth, C. M. Estimating earthquake location and magnitude from seismic intensity data. *Bull. Seism. Soc. Am.* **87**, 1502–1521 (1997).

32. Atkinson, G. M. & Wald, D. J. 'Did you feel it?' intensity data: A surprisingly good measure of earthquake ground motion. *Seism. Res. Lett.* **78**, 362–368 (2007).

33. Ghosh, G. K. & Mahajan, A. K. Intensity attenuation relation at Chamba-Garhwal area in north-west Himalaya with epicentral distance and magnitude. *J. Earth Syst. Sci.* **122**, 107–122 (2013).

34. Storchak, D. A. et al. Public release of the ISC-GEM global instrumental earthquake catalogue (1900-2009). *Seism. Res. Lett.* **84**, 810–815 (2013).

35. Hubbard, J. et al. Structural segmentation controlled the 2015 Mw 7.8 Gorkha earthquake rupture in Nepal. *Geology* **44**, 639–642 (2016).

36. Wesnousky, S. G. et al. Geological observations on large earthquakes along the Himalayan frontal fault near Kathmandu, Nepal. *Earth Planet. Sci. Lett.* **457**, 366–375 (2017).

37. Murphy, M. A. et al. Limit of strain partitioning in the Himalaya marked by large earthquakes in western Nepal. *Nat. Geosci.* **7**, 38–42 (2014).

38. Silver, C., Murphy, M. A., Taylor, M. H., Gosse, J. & Baltz, T. Neotectonics of the Western Nepal Fault System: implications for Himalayan strain partitioning. *Tectonics* **34**, 2494–2513 (2015).

39. Jouanne, F., Mugnier, J.-L., Sapkota, S. N., Bascou, P. & Pêcher, A. Estimation of coupling along the Main Himalayan Thrust in the central Himalaya. *J. Asian Earth Sci.* **133**, 62–71 (2017).

40. Okino, T. & Satoh, Y. Morphology, physics, chemistry and biology of Lake Rara in West Nepal. *Hydrobiologia* **140**, 125–133 (1986).

41. McGregor, H. V., Dupont, L., Stuut, J.-B. W. & Kuhlmann, H. Vegetation change, goats, and religion: a 2000-year history of land use in southern Morocco. *Quat. Sci. Rev.* **28**, 1434–1448 (2009).

42. Blott, S. J. & Pye, K. GRADISTAT: a grain size distribution and statistics package for the analysis of unconsolidated sediments. *Earth Surf. Proc. Landf.* **26**, 1237–1248 (2001).

43. Folk, R. L. & Ward, W. C. Brazos River bar: A study in the significance of grain size parameters. *J. Sedim. Res.* **27**, 3–26 (1957).
44. Verardo, D. J., Froelich, P. N. & McIntyre, A. Determination of organic carbon and nitrogen in marine sediments using the Carlo Erba NA-1500 Analyzer. *Deep-Sea Res.* **37**, 157–165 (1990).

45. Kak, A. C. & Slaney, M. Principles of computerized tomographic imaging (IEEE Press, 1988).

46. Eakins, J. D. & Morrison, R. T. A new procedure for the determination of lead-210 in lake and marine sediments. *Int. J. Appl. Radiat. Isot.* **29**, 531–536 (1978).

47. Mathieu, G. G., Biscaye, P. E., Lupton, R. A. & Hammond, D. E. System for measurement of 222Rn at low levels in natural waters. *Health Phys.* **55**, 989–992 (1988).

48. Yokoyama, Y., Miyairi, Y., Matsuzaki, H. & Tsunomori, F. Relation between acid dissolution time in the vacuum test tube and time required for graphitization for AMS target preparation. *Nucl. Instrum. Methods Phys. Res.* **259**, 330–334 (2007).

49. Yokoyama, Y., Koizumi, M., Matsuzaki, H., Miyairi, Y. & Ohkouchi, N. Developing ultra small-scale radiocarbon sample measurement at the University of Tokyo. *Radiocarbon* **52**, 310–318 (2010).

50. Reimer, P. J. et al. IntCal13 and Marine13 radiocarbon age calibration curves 0–50,000 years cal BP. *Radiocarbon* **55**, 1869–1887 (2013).

51. Blaauw, M. & Christen, J. A. Flexible paleoclimate age-depth models using an autoregressive gamma process. *Bayesian Anal.* **6**, 457–474 (2011).

52. Pagani, M. et al. OpenQuake Engine: An open hazard (and risk) software for the Global Earthquake Model. *Seismol. Res. Lett.* **85**, 692–702 (2014).

53. Cua, G. et al. "Best practices" for using macroseismic intensity and ground motion intensity conversion equations for hazard and loss models in GEM1 (GEM Technical Report 2010-4, GEM Foundation, Pavia, Italy, 2010).

54. Ader, T. et al. Convergence rate across the Nepal Himalaya and interseismic coupling on the Main Himalayan Thrust: Implications for seismic hazard. *J. Geophys. Res.* **117**, B04403 (2012).

55. Dziewonski, A. M., Chou, T.-A. & Woodhouse, J. H. Determination of earthquake source parameters from waveform data for studies of global and regional seismicity. *J. Geophys. Res.* **86**, 2825–2852 (1981).

56. Ekström, G., Nettles, M. & Dziewonski, A. M. The global CMT project 2004–2010: Centroid-moment tensors for 13,017 earthquakes. *Phys. Earth Planet. Inter.* **200-201**, 1–9 (2012).

57. Taylor, M. & Yin, A. Active structures of the Himalayan-Tibetan orogen and their relationships to earthquake distribution, contemporary strain field, and Cenozoic volcanism. *Geosphere* **5**, 199–214 (2009).

58. Jackson, D. in *Tibet, Past and Present* 147–159 (Brill, Leiden, 2002).

59. Kumar, S. et al. Paleoseismological evidence of surface faulting along the northeastern Himalayan front, India: Timing, size, and spatial extent of great earthquakes. *J. Geophys. Res.* **115**, B12422–20 (2010).

60. Nakata, T., Kumura, K. & Rockwell, T. First successful paleoseismic trench study on active faults in the Himalaya. *EOS. Trans. Am. Geophys. Union* **45**, F615 (1998).

61. Le Roux-Mallouf, R. et al. First paleoseismic evidence for great surface-rupturing earthquakes in the Bhutan Himalayas. *J. Geophys. Res.* **121**, 7271–7283 (2016).

## Acknowledgements

This research was supported by Marie Curie Initial Training Network iTECC funded by the EU REA under the FP7 implementation of the Marie Curie Action, grant agreement # 316966. Malin Kylander (Stockholm University) provided access to the Itrax XRF core scanner and analytical support. Y. Miyairi, C. Sawada and T. Aze assisted with the AMS measurements, which were funded by grants JSPS KAKENHI JP15KK0151 and JP17H1168 to YY. Lorenzo Gemignani, Arnaud Watlet, Guillaume Lippinois and Rham Kumar provided scientific and logistical support during and after the field campaign. Thomas Vandorpe and Maarten Van Daele contributed to the CT analysis of the sediment cores. We thank Thomas Lecocq, Tom Roland, Arnaud Watlet and Thomas Lemaire, as well as Jamie Howarth and an anonymous reviewer, for insightful comments on earlier drafts of this manuscript. iSTerre is part of the Labex OSUG@2020 (ANR10 LABX56).

## Author contributions

Z.G. and S.B. conceived the study. Z.G., S.B. and P.v.d.B interpreted the results. K.V. implemented the IPEs and rupture scenarios in OpenQuake, and performed all shaking intensity computations. Y.Y. performed the radiocarbon measurements. J.N. contributed to the sampling/coring. A.P.G. facilitated the logistics of the field expedition. Z.G., S.B., P. v.d.B. and K.V. wrote the paper with inputs from all co-authors.

## Additional information

Supplementary Information accompanies this paper at https://doi.org/10.1038/s41467-019-10093-4.

**Competing interests:** The authors declare no competing interests.

**Reprints and permission** information is available online at http://npg.nature.com/ reprintsandpermissions/

**Journal Peer Review Information:** *Nature Communications* thanks Jamie Howarth and the other anonymous reviewers for their contribution to the peer review of this work.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2019