This is an implementation in OpenQuake of the probabilistic seismic hazard analysis (PSHA) model of Nath & Thingbaijam (2012).
Files specific to the first version of the model are subscripted _v0.
This model is described in the report by Ackerley (2016).
A revised version, subscripted _v1 and incorporating clarifications from the authors of
Nath & Thingbaijam (2012), was included in the initial release of the Global Seismic Seismic Hazard Map, v. 2018.1.
Work is in progress to develop a _v2 model to address outstanding issues.
Zipped job files can be downloaded for various releases.
Smoothed-gridded model files in Natural hazards’ Risk Markup Language (NRML) format are too big to fit in a github repository but they can be regenerated from smaller input data files and scripts, all included in the repository.
The OpenQuake NRML .xml model files are generated from .csv (source model) and .tsv (logic tree) files using scripts that depend on the OpenQuake Hazard Modeler's Toolikt (HMTK). This is done using python .py and Jupyter Notebook scripts in .ipynb format.
Ensure that the prerequisites are installed:
Clone the repository:
cd ~/src
git clone https://github.com/nackerley/indian-subcontinent-psha.gitCreate an environment with the necessary python packages:
cd ~/src/indian-subcontinent-psha
conda env create --file oq.yml
source activate oqMake the (oq) kernel accessible to jupyter in the root
environment by running the following:
source activate oq
python -m ipykernel install --user --name oqTo open jupyter notebooks which use this environment, run the following in your root environment:
jupyter notebookand select the appropriate files from the browser window which opens.
Add ~/src/indian-subcontinent-psha/utilities to your Python path.
One way is to append the following to your ~/.bashrc:
export PYTHONPATH=$PYTHONPATH:~/src/indian-subcontinent-psha/utilitiesVerify your installation by trying:
import source_model_tools as smtTo regenerate the logic tree XML, open and run the following jupyter notebook:
"Logic Trees/logic_trees_nt2012.ipynb"Note that the logic tree model XML files are small, and already included in the repository, so this is only necessary if changes are made to logic tree TSV description files.
To regenerate the source model XML, in your (oq) environment run:
cd ~/src/indian-subcontinent-psha/Source\ Models/
python write_source_models_nt2012.py --version VERSION --fullThe smooothed-gridded models in particualr are large, so only "thinned"
versions are included in the repository for quality control purposes. The version
number only refers to the major version. Without the --full argument
only a "thinned" version of the smoothed-gridded model is generated; this is
quicker and useful for testing.
The following jupyter notebooks are useful for investigating and visualizing the source models:
"Source Models/areal_source_models_nt2012.ipynb""Source Models/smoothed_source_models_nt2012.ipynb""Source Models/collapsed_source_models_nt2012.ipynb"Job configuration job.ini files are created manually rather than
automatically generated. Symbolic links to data files are include in each
folder, to avoid duplication of files.
The key job variants supported are:
Jobs/cities_collapsed_v0/Jobs/cities_collapsed_v1/Jobs/cities_full_disaggregation_v1/Jobs/map_collapsed_v1/Note:
collapsed vs. full refers to whether or not frequency-magnitude
distributions were collapsed prior to hazard calculation.
collapsed source models will give correct results for mean hazard but not
for hazard quantiles or deaggregation.cities jobs are site-specific analyses, for the 18 cities listed
in Table 3 of Nath & Thingbaijam (2012).map refers to the 0.2° grid of 8102 points used to generate the data in the
Figure 7 and the electronic supplement of Nath & Thingbaijam (2012).The script tar_all.py in Jobs zipped versions of all jobs found in
subdirectories.
Files for generating maps using QGIS are stored at Maps/.
Maps incorporate data which can be downloaded separately from:
This project draws heavily on published work. Portions of some of these publications are included in this repository.
This model incorporates data from the following publications:
Data/nath2012probabilistic/:
tables, digitized hazard curves and electronic supplement of
Nath & Thingbaijam (2012)Data/nath2011peak/:
tables of
Nath & Thingbaijam (2011)Data/thingbaijam2011seismogenic/:
Table 1 from "Thingbaijam & Nath (2011) A Seismogenic Source Framework for the Indian subcontinent" (unpublished)Catalogue/ for model validation:
Nath, Thingbaijam & Ghosh (2010) Earthquake catalogue of South Asia - a generic MW scale frameworkA set of jupyter notebooks in ~/src/indian-subcontinent-psha/GMPEs were
used to generate validation data for the following modules of
openquake.hazardlib.gsim:
GMPEs/gupta2010response/:
openquake.hazardlib.gsim.gupta_2010GMPEs/kanno2006new/:
openquake.hazardlib.gsim.kanno_2006GMPEs/nath2012ground/:
openquake.hazardlib.gsim.nath_2012GMPEs/raghukanth2007estimation/:
openquake.hazardlib.gsim.raghukanth_iyengar_2007GMPEs/sharma2009ground/:
openquake.hazardlib.gsim.sharma_2009GMPEs/atkinson2003empirical/:
openquake.hazardlib.gsim.atkinson_boore_2003.AtkinsonBoore2003SSlabJapanPull requests are welcome (though not expected, given the nature of the repository).
Useful tools for development:
The Indian Subcontinent PSHA model is released under the GNU Affero Public License 3.