Earth, Planets and Space (May 2024)

A future scenario earthquake and ground motion hazards for Kathmandu, Nepal

  • Kazuki Koketsu,
  • Hiroe Miyake,
  • Koji Okumura,
  • Haruhiko Suzuki

DOI
https://doi.org/10.1186/s40623-024-02018-3
Journal volume & issue
Vol. 76, no. 1
pp. 1 – 17

Abstract

Read online

Abstract Since Nepal is an earthquake-prone country due to the collision of the Indian and Eurasian plates, it is crucial for its capital of Kathmandu to evaluate ground motion hazards from a future large earthquake. For this purpose, we constructed a realistic scenario earthquake with realistic rupture parameters in a likely location. To obtain the location, a new distribution of the rupture zones of large historical earthquakes along the Main Himalayan Thrust (MHT) was obtained by integrating the distribution from a previous study and the results of trench surveys. Global Navigation Satellite System observations indicated that the plate boundary is strongly coupled from the southern boundary of the MHT to a depth of approximately 10 km and there is almost no lateral change in the coupling. This implies that all regions along the MHT have similar rates of strain increase. Therefore, it is most probable that the rupture zone of the oldest previous event will rupture as a scenario earthquake. In the new distribution, the 1255 earthquake is the oldest. However, large earthquakes occurred in 1934 and 2015 within its rupture zone. Thus, we adopted the area obtained by removing the 1934 and 2015 rupture zones from the western part of the 1255 rupture zone. The relationship between the rupture area size and seismic moment of the 2015 earthquake lies between the scaling formulas for crustal earthquakes and plate-boundary earthquakes, but is closer to the former. Therefore, using this and the scheme for characterized source models, we determined realistic rupture parameters. We then simulated broadband ground motions in Kathmandu using these rupture parameters, our 3-D velocity structure models, and a hybrid method combining the finite-difference method and the stochastic Green’s function method. We obtained the peak ground accelerations (PGAs) of simulated ground motions, and calculated the seismic intensities in the Modified Mercalli intensity scale from the PGAs as indexes of hazards for Kathmandu. Intensities IX coincide with the center of the Kathmandu Valley, and intensities VIII and VII are found in the area surrounded by the sedimentary boundary and the southernmost part of the valley. Graphical abstract

Keywords