Source code for openquake.hazardlib.gsim.weatherill_2024
# -*- coding: utf-8 -*-# vim: tabstop=4 shiftwidth=4 softtabstop=4## Copyright (C) 2014-2023 GEM Foundation## OpenQuake is free software: you can redistribute it and/or modify it# under the terms of the GNU Affero General Public License as published# by the Free Software Foundation, either version 3 of the License, or# (at your option) any later version.## OpenQuake is distributed in the hope that it will be useful,# but WITHOUT ANY WARRANTY; without even the implied warranty of# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the# GNU Affero General Public License for more details.## You should have received a copy of the GNU Affero General Public License# along with OpenQuake. If not, see <http://www.gnu.org/licenses/>."""Module exports :class:`Weatherill2024ESHM20AvgSA`, :class:`Weatherill2024ESHM20SlopeGeologyAvgSA`, :class:`Weatherill2024ESHM20AvgSAHomoskedastic`"""fromopenquake.hazardlib.imtimportAvgSAfromopenquake.hazardlib.gsim.kotha_2020importKothaEtAl2020ESHM20fromopenquake.hazardlib.gsim.baseimportCoeffsTablefromopenquake.hazardlibimportconst
[docs]classWeatherill2024ESHM20AvgSA(KothaEtAl2020ESHM20):""" This class implements a variation of the Kotha et al (2020; 2022) GMM that was used for the ESHM20, but here the predicted intensity measure is average SA (AvgSA) rather than SA. This is a form of direct AvgSA GMM, which is fit using the same data set as that of KothaEtAl2020 with AvgSA defined according the specifications of (among others) Iacoletti et al. (2023): AvgSA = sqrt(prod([0.2 x T <= T <= 1.5 x T])) where a total of 10 linearly-spaced conditioning periods in the range are used to define the average SA. As the same regression methods were used to fit AvgSA then all of the adjustment terms adopted by the ESHM20 (sigma_mu_epsilon, c3_epsilon, ergodic etc.) can be applied to the AvgSA GMM, which allows the same logic tree to be constructed for the direct AvgSA case. Further details on the compilation and application of the GMM are being developed in the following publication (in preparation): Weatherill, G (2024) "A Regionalised Direct AvgSA Ground Motion Model for Europe", (Journal TBC) As this is in preparation, futue changes to the model are possible so we therefore retain the experimental warning, which will be removed at a future date. """experimental=True#: Set of :mod:`intensity measure types <openquake.hazardlib.imt>`#: this GSIM can calculate. A set should contain classes from module#: :mod:`openquake.hazardlib.imt`.DEFINED_FOR_INTENSITY_MEASURE_TYPES={AvgSA,}#: Supported standard deviation types is are only total std.devDEFINED_FOR_STANDARD_DEVIATION_TYPES={const.StdDev.TOTAL,const.StdDev.INTER_EVENT,const.StdDev.INTRA_EVENT}#: Required site parameters are vs30, vs30measured and the eshm20_regionREQUIRES_SITES_PARAMETERS=set(("region","vs30","vs30measured"))kind="avgsa_ESHM20"# Coefficients obtained direclty from the regression outputsCOEFFS=CoeffsTable(sa_damping=5,table="""\ imt e1 b1 b2 b3 c1 c2 c3 tau_c3 phis2s tau_event_0 tau_l2l phi_0 d0_obs d1_obs phi_s2s_obs d0_inf d1_inf phi_s2s_inf AvgSA(0.050) 4.3468615401 2.1512063226 0.3717451438 0.3840026031 -1.5844962290 0.3056929868 -0.6358197233 0.2763393183 0.6594454507 0.4280092597 0.5380302118 0.4676597188 2.61783357 -0.41984752 0.40314318 1.88712924 -0.29623828 0.54897643 AvgSA(0.100) 4.6536804102 2.1243118148 0.3569931311 0.3628813315 -1.5186700529 0.2841907202 -0.7361755939 0.3101014049 0.6722188832 0.4398426213 0.5568916029 0.4648998708 2.63178976 -0.42260122 0.42319496 1.78617139 -0.28015514 0.57163410 AvgSA(0.150) 4.7945771137 2.1274812721 0.3288504315 0.3795627756 -1.4389896655 0.2501908703 -0.7667357492 0.3157034067 0.6633675533 0.4434411516 0.5202488389 0.4582624492 2.71796680 -0.43603910 0.43472118 1.78935693 -0.28055203 0.57516812 AvgSA(0.200) 4.8268358912 2.1435342480 0.3057719029 0.4094215327 -1.3832602354 0.2227051661 -0.7550607691 0.3110395382 0.6458262629 0.4435212138 0.4841042855 0.4498230477 2.86217015 -0.45823486 0.43749194 1.87918582 -0.29464359 0.56498414 AvgSA(0.250) 4.8051268888 2.1754139829 0.2918562846 0.4426226923 -1.3412962888 0.2019136688 -0.7305447434 0.3035029231 0.6280883485 0.4407924746 0.4529410313 0.4425130460 3.06323423 -0.48950776 0.43044807 2.04792489 -0.32117497 0.54608798 AvgSA(0.300) 4.7620900260 2.2141154082 0.2835145315 0.4665423314 -1.3070686825 0.1860974192 -0.7071940818 0.2957328792 0.6153905644 0.4369210190 0.4322453716 0.4352045167 3.26556277 -0.52103527 0.41182470 2.25496936 -0.35374018 0.52720419 AvgSA(0.400) 4.6332719548 2.2947237222 0.2741828304 0.5159131190 -1.2627673311 0.1619710217 -0.6499258739 0.2794709669 0.5989732123 0.4255359229 0.3951049711 0.4216267430 3.49875383 -0.55740129 0.39503497 2.47642829 -0.38853586 0.50678238 AvgSA(0.500) 4.4778294883 2.3662343667 0.2696192831 0.5660951786 -1.2319346319 0.1446120370 -0.5980583921 0.2666824382 0.5902397100 0.4169401398 0.3687464707 0.4117459145 3.71388275 -0.59068507 0.37900446 2.68871629 -0.42186426 0.49004842 AvgSA(0.600) 4.3185533698 2.4277720222 0.2652148003 0.6175414034 -1.2077156215 0.1318085467 -0.5511719494 0.2536315305 0.5876115938 0.4103320095 0.3515066109 0.4034543724 3.87941606 -0.61580447 0.36505069 2.86525837 -0.44955313 0.47419374 AvgSA(0.700) 4.1672194228 2.4894825880 0.2650749306 0.6679669266 -1.1874512590 0.1228799964 -0.5133671436 0.2433881796 0.5867287159 0.4065674107 0.3347690807 0.3963077735 3.97630461 -0.62992274 0.36262323 2.98664886 -0.46856927 0.46461818 AvgSA(0.800) 4.0288646652 2.5542391769 0.2695180409 0.7085251974 -1.1727356202 0.1173776862 -0.4805689498 0.2315398128 0.5861237737 0.4053007828 0.3200552353 0.3900808395 4.04577748 -0.63992675 0.36701024 3.08292173 -0.48365672 0.45841564 AvgSA(0.900) 3.8947701849 2.6046828029 0.2708607782 0.7424403141 -1.1636251550 0.1140958019 -0.4495322496 0.2191350303 0.5868406504 0.4054880421 0.3128836707 0.3847701262 4.09468998 -0.64675747 0.37430943 3.17474995 -0.49805311 0.45633438 AvgSA(1.000) 3.7694600868 2.6564485947 0.2755580747 0.7773839055 -1.1569230243 0.1119008735 -0.4228857116 0.2077140636 0.5883712385 0.4076373926 0.3052586715 0.3809345790 4.11952534 -0.64980249 0.38194385 3.26430061 -0.51212861 0.45553548 AvgSA(1.250) 3.4706631778 2.7640114111 0.2863971445 0.8515842010 -1.1500139705 0.1098964051 -0.3642458982 0.1878359113 0.5911312275 0.4175841723 0.2877202407 0.3732270011 4.11297148 -0.64828200 0.38649790 3.37305996 -0.52913277 0.45322510 AvgSA(1.500) 3.2190113497 2.8630381148 0.3011679747 0.8833932182 -1.1522344562 0.1111509093 -0.3154533960 0.1759946561 0.5916370848 0.4251689388 0.2791019570 0.3684836574 4.09477274 -0.64520087 0.38407000 3.48254160 -0.54618847 0.45284527 AvgSA(1.750) 3.0052586300 2.9878647999 0.3413194350 0.9057249259 -1.1533431170 0.1156390944 -0.2832581008 0.1664335521 0.5921743466 0.4297625190 0.2848615012 0.3646174735 4.04481293 -0.63766481 0.38178972 3.54860458 -0.55647158 0.45407701 AvgSA(2.000) 2.7977898222 3.0610742969 0.3563412955 0.9382679687 -1.1527918383 0.1203796346 -0.2618124393 0.1609297797 0.5915803874 0.4379178489 0.2782289105 0.3621184180 3.98022776 -0.62800105 0.37429046 3.58857586 -0.56254997 0.45902122 AvgSA(2.500) 2.4526046480 3.2447441319 0.4175895089 0.9935574572 -1.1555045836 0.1319754237 -0.2362881913 0.1606151346 0.5851395458 0.4487655553 0.2727004457 0.3564559472 3.88094598 -0.61300917 0.36664751 3.58945419 -0.56256504 0.46109358 AvgSA(3.000) 2.1846343267 3.3795075185 0.4680279735 1.0273202155 -1.1533536786 0.1423911641 -0.2329857032 0.1478724120 0.5834092888 0.4569015623 0.2736582001 0.3508293944 3.78141847 -0.59784839 0.35646268 3.51803344 -0.55144112 0.45739762 AvgSA(3.500) 1.9409098870 3.4656531050 0.5118696601 1.0856421414 -1.1701293398 0.1553162861 -0.2020706851 0.1410898528 0.5750292172 0.4425943058 0.2958590722 0.3492035567 3.60781997 -0.57117670 0.34930242 3.48405639 -0.54602924 0.45663788 AvgSA(4.000) 1.7015054228 3.5163002844 0.5312213268 1.1456993040 -1.1841200080 0.1711021029 -0.1838747563 0.1437446953 0.5669140794 0.4500152389 0.2825370206 0.3499904629 3.43147290 -0.54394985 0.34725671 3.44930328 -0.54050316 0.45782846 AvgSA(4.500) 1.6213774162 3.7391842699 0.6641371123 1.1288485321 -1.1829236019 0.1847235866 -0.1807585867 0.1213816721 0.5757075561 0.4340712588 0.2587634879 0.3381428158 3.24877975 -0.51560770 0.35330607 3.43813135 -0.53856720 0.46219582 AvgSA(5.000) 1.4308131096 3.7717226445 0.6801268753 1.1816942442 -1.1999994884 0.1988877433 -0.1645299267 0.1239688756 0.5681631639 0.4376824137 0.2528453711 0.3391870321 3.07021475 -0.48775155 0.37025310 3.46440379 -0.54233758 0.47057169 """)
[docs]classWeatherill2024ESHM20AvgSAHomoskedastic(Weatherill2024ESHM20AvgSA):"""Variant of the Weatherill2024ESHM20 direct GMPE for AvgSA with the homoskedastic sigma coming from the original mixed effects regression """experimental=Truekind="avgsa_ESHM20_homoskedastic"