Package: gmKriger
gmKrigergmKriger: A Kriging-based ground motion intensity measure (GMIM) calculator.gmKrigercomputes 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.InstallationNeed the following packages:pipinstallnvectorpipinstalltorchpipinstallgeopypipinstallshapelypipinstallnvectorpipinstallgeostatspipinstallgmmspipinstallpygmmpipinstallgmKrigerHow to useInputssite: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').RunimportgmKrigersite=['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)SiteLat (deg)Lon (deg)PGA (g)sigma_PGA (ln)Alameda Naval Air Station37.7857-122.3090.1893120.386172Treasure Island37.8261-122.3710.1353610.273654Alameda Bay Farm Island37.7338-122.2500.1519780.429026Farris Farm36.9103-121.7440.4584040.494209POO737.8052-122.3390.1544220.350264Run this example using this code:here.Models availableGet them usingthis code.The spatial correlation models for the events and ground motion intensity measures below are accessible viaDesignSafe(Pretell et al. 2023).EarthquakePGAPGVIaCAVCAVdpPSA(0.100)PSA(0.300)PSA(0.600)PSA(1.000)PSA(3.000)PSA(6.000)PSA(10.000)1968 M8.2 Tokachi-OkiNoNoYesYesNoNoNoNoNoNoNoNo1971 M6.6 San FernandoYesYesYesYesNoNoNoNoNoNoNoNo1978 M7.7 Miyagiken-OkiNoNoYesYesNoNoNoNoNoNoNoNo1979 M6.5 Imperial ValleyYesYesYesYesNoNoNoNoNoNoNoNo1980 M6.3 VictoriaYesYesYesYesNoNoNoNoNoNoNoNo1981 M5.9 WestmorlandYesYesYesYesNoNoNoNoNoNoNoNo1983 M7.7 Nihonkai-ChubuNoNoYesYesNoNoNoNoNoNoNoNo1983 M6.8 Nihonkai-ChubuNoNoYesYesNoNoNoNoNoNoNoNo1987 M6.5 Superstition HillsYesYesYesYesNoNoNoNoNoNoNoNo1989 M6.9 Loma PrietaYesYesYesYesNoNoNoNoNoNoNoNo1993 M7.6 Kushiro-Oki HokkaidoNoNoYesYesNoNoNoNoNoNoNoNo1994 M6.7 NorthridgeYesYesYesYesNoNoNoNoNoNoNoNo1994 M8.3 Toho-Oki HokkaidoNoNoYesYesNoNoNoNoNoNoNoNo1995 M6.9 KobeYesYesYesYesNoNoNoNoNoNoNoNo1999 M7.5 KocaeliYesYesYesYesNoNoNoNoNoNoNoNo1999 M7.6 Chi-ChiYesYesYesYesNoNoNoNoNoNoNoNo2000 M6.6 TottoriYesYesYesYesNoNoNoNoNoNoNoNo2002 M5.0 Au Sable ForksNoNoNoNoNoNoNoNoNoNoNoNo2003 M8.3 TokachiNoNoYesYesNoNoNoNoNoNoNoNo2007 M6.8 Chuetsu-okiYesYesYesYesNoNoNoNoNoNoNoNo2010 M7.2 El Mayor-CucapahYesYesYesYesNoNoNoNoNoNoNoNo2010 M7.0 DarfieldYesYesYesYesNoNoNoNoNoNoNoNo2010 M8.8 MauleNoNoYesYesNoNoNoNoNoNoNoNo2011 M6.2 ChristchurchYesYesYesYesNoNoNoNoNoNoNoNo2011 M5.0 ChristchurchYesNoNoNoNoNoNoNoNoNoNoNo2011 M6.0 ChristchurchYesYesYesYesNoNoNoNoNoNoNoNo2011 M5.9 LyttletonYesYesYesYesNoNoNoNoNoNoNoNo2011 M9.1 Tohoku-OkiNoNoYesYesNoNoNoNoNoNoNoNo2012 M6.1 EmiliaYesYesYesYesNoNoNoNoNoNoNoNo2012 M6.0 EmiliaYesYesYesYesNoNoNoNoNoNoNoNo2019 M7.06 RidgecrestYesYesYesYesNoNoNoNoNoNoNoNo2019 M6.48 RidgecrestYesYesYesYesNoNoNoNoNoNoNoNo2020 M7.0 SamosYesYesYesYesNoNoNoNoNoNoNoNo2023 M7.81 PazarcikYesYesYesYesYesYesYesYesYesYesYesYes2023 M7.74 KahramanmarasYesYesYesYesYesYesYesYesYesYesYesYes2023 M6.81 NurdagiYesYesYesYesYesYesYesYesYesYesYesYes2023 M6.37 YayladagiYesYesYesYesYesYesYesYesYesYesYesYesGround motion intensity measures (GMIMs) and unitsGMIMKeyUnitPeak ground acelerationPGAgPeak ground velocityPGVcm/sArias intensityIam/sCumulative absolute velocityCAVm/sDamage-potential CAVCAVdpg-sPseudo-spectral acceleration @ 0.1 sPSA(0.100)gPseudo-spectral acceleration @ 0.3 sPSA(0.300)gPseudo-spectral acceleration @ 0.6 sPSA(0.600)gPseudo-spectral acceleration @ 1.0 sPSA(1.000)gPseudo-spectral acceleration @ 3.0 sPSA(3.000)gPseudo-spectral acceleration @ 6.0 sPSA(6.000)gPseudo-spectral acceleration @ 10.0 sPSA(10.000)gExamplesExample 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.AcknowledgementsThe implementation of gmKriger greatly benefitted from discussions with Scott J. Brandenberg and Jonathan P. Stewart.CitationIf you use these codes, please cite:Pretell, R. (2023). RPretellD/gmKriger: Initial release (0.2.0). Zenodo.https://doi.org/10.5281/zenodo.10399419Pretell, 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-t096v2ContactFor any questions or comments, contact Renmin Pretell atrpretell@unr.edu.