India (IND)


N. Ackerley, K.K.S. Thingbaijam, S.K. Nath


Coverage of the Indian subcontinent is with the hazard model developed by Nath and Thingbaijam (2012). This model covers India, Bangladesh, Bhutan and Nepal. The model has been updated and translated from its original format into the OpenQuake (OQ) engine in collaboration with Natural Resources Canada. Additional information, material and documentation about the implementation of the model can be found at:

Information about the OQ model versions and input files can be found on the Results and Dissemination page.

Interactive Viewer

The viewer below depicts the seismic sources and hazard results in terms of PGA for a return period of 475 years. Click on the menu in the upper right corner to select the layer.


Tectonic overview

The Indian subcontinent moves northward with respect to Eurasia and is colliding with the southern Asian margin at 35-45 mm/yr; this collision has contributed greatly to the uplift of the great mountain ranges of central and eastern Asia, including the Himalaya, the Tien Shan, the Pamir, and the Tibetan plateau. The Indian subcontinent is much stronger than the Asian continental crust with which it collides, though, so most of the deformation and seismicity that results from the plate collision is located in and north of the Himalaya outside of India. Exceptions include in northwest and northeast India where the Himalaya are within the national borders, and the Shillong region of eastern India, where some component of the plate convergence is accommodated on intraplate faults. The faults in the Himalayan belt take up at least half of the total plate convergence and represent the greatest source of seismic hazard to India; earthquakes on the Himalayan thrusts may be large enough that those in Nepal, Bhutan and Pakistan may still produce dangerous seismic shaking in densely-populated northern India.

Nonetheless, the strong Indian crust is capable of transmitting compressive stresses over great distances, and large, damaging earthquakes have occurred well within the Indian borders as a result. A prime example of this is the 2001 Mw 7.6 Bhuj earthquake in Gujarat, which killed over 2000 people (Bodin and Horton, 2004, Bulletin of the Seismological Society of America)

Basic Datasets

See Nath and Thingbaijam (2012) for a description of the datasets used for developing the hazard model.

Hazard Model

Seismic Source Characterisation

Nath and Thingbaijam (2012) proposed three source models: a single set of areal seismogenic source zones, and two smoothed-gridded point source models. Epistemic uncertainty is considered through alternative values of b and Mmax in each source zone.

Ground Motion Characterisation

A wide range of tectonic regions are considered, and epistemic uncertainty is accounted for by using multiple GMPEs per tectonic region.

active shallow crust normal Weight
AkkarBommer2010 0.25
BooreAtkinson2008 0.25
CampbellBozorgnia2008 0.25
Kanno2006Shallow 0.25
active shallow crust strike-slip reverse Weight
AkkarBommer2010 0.2
BooreAtkinson2008 0.2
CampbellBozorgnia2008 0.2
Kanno2006Shallow 0.2
SharmaEtAl2009 0.2
intraplate margin lower Weight
AtkinsonBoore2006 0.25
ToroEtAl2002SHARE 0.25
SharmaEtAl2009 0.25
NathEtAl2012Lower 0.25
intraplate margin upper Weight
AtkinsonBoore2006 0.25
ToroEtAl2002SHARE 0.25
SharmaEtAl2009 0.25
NathEtAl2012Upper 0.25
stable shallow crust Weight
AtkinsonBoore2006 0.25
ToroEtAl2002SHARE 0.25
Campbell2003MwNSHMP2008 0.25
RaghukanthIyengar2007 0.25
subduction interface Weight
AtkinsonBoore2003SInter 0.334
ZhaoEtAl2006SInter 0.333
Kanno2006Deep 0.333
subduction interface megathrust Weight
AtkinsonBoore2003SInter 0.25
AtkinsonMacias2009 0.25
ZhaoEtAl2006SInter 0.25
Kanno2006Deep 0.25
subduction intraslab Himalayas Weight
AtkinsonBoore2003SSlabJapan 0.25
YoungsEtAl1997SSlab 0.25
ZhaoEtAl2006SSlab 0.25
LinLee2008SSlab 0.25
subduction intraslab Weight
AtkinsonBoore2003SSlabCascadia 0.25
YoungsEtAl1997SSlab 0.25
ZhaoEtAl2006SSlab 0.25
Gupta2010SSlab 0.25


Hazard curves were computed with the OQ engine for peak ground acceleration (PGA) and spectral acceleration (SA) at 0.2s, 0.5s, 1.0s, and 2s. The computation was performed on a grid of 55038 sites (spaced at approximately 10 km) with reference soil conditions with shear wave velocity in the upper 30 meters (Vs30) of 760-800 m/s.

The hazard map for PGA corresponding to a 10% probability of exceedance in 50 years (475 year return period), can be seen using the interactive viewer. For a more comprehensive set of hazard and risk results, please see the GEM Visualization Tools.


Nath, S. K. and Thingbaijam, K. K. S. (2012). Probabilistic seismic hazard assessment of India. Seismological Research Letters, 83(1):135–149.