gmKriger

gmKriger: A Kriging-based ground motion intensity measure (GMIM) calculator. gmKriger computes GMIMs for past earthquake events given a site's location (latitude and longitude) and the site's Vs30. gmKriger uses Ordinary Kriging interpolation and the spatial correlation models functional forms proposed by Bodenmann et al. (2023) developed using a Bayesian approach and ground motion data.

Installation

Need the following packages:

pip install nvector
pip install torch
pip install geopy
pip install shapely
pip install nvector

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 using this code: here.

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

Examples

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

Acknowledgements

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

Citation

If you use these codes, please cite:

Pretell, R. (2023). RPretellD/gmKriger: Initial release (0.2.0). Zenodo. https://doi.org/10.5281/zenodo.10399419

Pretell, R., Brandenberg, S.J., and Stewart, J.P. "Optimizing ground motion intensity measures for soil liquefaction case histories." (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 v2

Contact

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