gmKriger: Kriging-based ground motion intensity measure calculator

gmKriger computes ground motion intensity measures for past earthquake events given a site's location (latitude and longitude) and Vs30.

Installation

Need the following packages:

pip install geostats
pip install gmms
pip install pygmm
pip install gmKriger

How to use

Inputs

site:
Site ID(s) or site name(s).

latitude:
Site's latitude(s).

longitude:
Site's longitude(s).

Vs30:
Time-average shear-wave velocity in the top 30 m for the site(s).

earthquake:
Event from the available models (e.g., '1989 M6.9 Loma Prieta').

model:

• realizations: To use 1000 spatial correlation models.
  
• MAP: To use the maximum a posteriori spatial correlation model.
  
• all: To use all the 1000 and the maximum a posteriori spatial correlation models.
  

gmim:
Ground motion intensity measure from the available models (e.g., 'PGA').

Run

import gmKriger

site      = ['Alameda Naval Air Station', 'Treasure Island', 'Alameda Bay Farm Island', 'Farris Farm', 'POO7']
latitude  = [37.785748,37.8261394,37.73380567,36.91026828,37.805242]
longitude = [-122.309346,-122.3712351,-122.250101,-121.7437891,-122.339702]
Vs30      = [186.2,181.1,230.7,209.5,223]

earthquake = '1989 M6.9 Loma Prieta'
model      = 'realizations'
gmim       = 'PGA'

gmKriger.get_Kgmim(site,latitude,longitude,Vs30,earthquake,model,gmim)

RUN this example YOURSELF using this Jupyter Notebook.

Models available

Get them using this code.

The spatial correlation models for the events and ground motion intensity measures below are accessible via DesignSafe (Pretell et al. 2023).

Ground motion intensity measures (GMIMs) and units

Example

• Example 1: Compute PGA for the 1989 Loma Prieta Earthquake using all the 1000 spatial correlation models.
• Example 2: Compute several ground motion intensity measures for the 2023 M7.8 Pazarcik Earthquake using the maximum a posteriori spatial correlation model.
• Example 3: Compute several ground motion intensity measures for the 1987 M6.5 Superstition Hills Earthquake all the 1000 and the maximum a posteriori spatial correlation model.

Acknowledgements

• The implementation of these codes greatly benefitted from discussions with Scott J. Brandenberg and Jonathan P. Stewart.

Citation

If you use these codes, please cite:

Pretell, R. (2023). gmKriger: Kriging-based ground motion intensity measure calculator (0.1.0). Zenodo. https://doi.org/10.5281/zenodo.10399419

Pretell, R., Brandenberg, S.J., and Stewart, J.P. "Consistently computed ground motion intensity measures at liquefaction case history sites". (In preparation).

Pretell, R., Brandenberg, S., and Stewart, J. (2023) "Consistently computed ground motion intensity measures for liquefaction triggering assessment." DesignSafe-CI. https://doi.org/10.17603/ds2-6vj1-t096

Contact

For any questions or comments, contact Renmin Pretell at rpretell@unr.edu.