# -*- coding: utf-8 -*-
# vim: tabstop=4 shiftwidth=4 softtabstop=4
#
# Copyright (C) 2012-2018 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:`ZhaoEtAl2016Asc`,
:class:`ZhaoEtAl2016AscSiteSigma`,
:class:`ZhaoEtAl2016UpperMantle`,
:class:`ZhaoEtAl2016UpperMantleSiteSigma`,
:class:`ZhaoEtAl2016SInter`,
:class:`ZhaoEtAl2016SInterSiteSigma`,
:class:`ZhaoEtAl2016SSlab`,
:class:`ZhaoEtAl2016SSlabSiteSigma`
"""
import numpy as np
from openquake.hazardlib.gsim.base import GMPE, CoeffsTable
from openquake.hazardlib import const
from openquake.hazardlib.imt import PGA, SA
[docs]class ZhaoEtAl2016Asc(GMPE):
"""
Implements the GMPE of Zhao et al (2016a) for shallow crustal and upper
mantle events from Japan. Only the shallow crustal version is implemented
here.
Zhao, J. X., Zhou, S., Zhou, J., Zhao, C., Zhang, H., Zhang, Y., Gao, P.,
Lan, X., Rhoades, D. A., Fukushima, Y., Somerville, P., Irikura, K.,
(2016c), "Ground-Motion Prediction Equations for Shallow Crustal and
Uppe-Mantle Earthquakes in Japan Using Site Class and Simple Geometric
Attenuation Functions", Bulletin of the Seismological Society of America,
106(4), 1518-1534
Main version with standard deviations independent of site term
"""
#: Supported tectonic region type is active shallow crust
DEFINED_FOR_TECTONIC_REGION_TYPE = const.TRT.ACTIVE_SHALLOW_CRUST
#: Supported intensity measure types are spectral acceleration,
#: and peak ground acceleration
DEFINED_FOR_INTENSITY_MEASURE_TYPES = set([
PGA,
SA
])
#: Supported intensity measure component is geometric mean
#: of two horizontal components :
DEFINED_FOR_INTENSITY_MEASURE_COMPONENT = const.IMC.AVERAGE_HORIZONTAL
#: Supported standard deviation types are inter-event, intra-event
#: and total
DEFINED_FOR_STANDARD_DEVIATION_TYPES = set([
const.StdDev.TOTAL,
const.StdDev.INTER_EVENT,
const.StdDev.INTRA_EVENT
])
#: Required site parameters is Vs30 (converted to site class)
REQUIRES_SITES_PARAMETERS = set(('vs30', ))
#: Required rupture parameters are magnitude, top-of-rupture depth and
#: style of faulting (rake)
REQUIRES_RUPTURE_PARAMETERS = set(('mag', 'rake', 'ztor'))
#: Required distance measure is Rrup and Rvolc
REQUIRES_DISTANCES = set(('rrup', 'rvolc'))
[docs] def get_mean_and_stddevs(self, sites, rup, dists, imt, stddev_types):
"""
See :meth:`superclass method
<.base.GroundShakingIntensityModel.get_mean_and_stddevs>`
for spec of input and result values.
"""
# extracting dictionary of coefficients specific to required
# intensity measure type.
C = self.COEFFS[imt]
C_SITE = self.SITE_COEFFS[imt]
s_c, idx = self._get_site_classification(sites.vs30)
sa_rock = (self.get_magnitude_scaling_term(C, rup) +
self.get_sof_term(C, rup) +
self.get_depth_term(C, rup) +
self.get_distance_term(C, dists, rup))
sa_soil = self.add_site_amplification(C, C_SITE, sites,
sa_rock, idx, rup)
stddevs = self.get_stddevs(C, sites.vs30.shape, idx, stddev_types)
return sa_soil, stddevs
[docs] def get_magnitude_scaling_term(self, C, rup):
"""
Returns the magnitude scaling term in equations 1 and 2
"""
if rup.mag <= self.CONSTANTS["m_c"]:
return C["ccr"] * rup.mag
else:
return (C["ccr"] * self.CONSTANTS["m_c"]) +\
(C["dcr"] * (rup.mag - self.CONSTANTS["m_c"]))
[docs] def get_sof_term(self, C, rup):
"""
Shallow crustal faults have a style-of-faulting dependence as
normal faulting is found to produce higher ground motion (equation 1)
"""
if rup.rake <= -45.0 and rup.rake >= -135.0:
# Normal faulting
return C["FN_CR"]
else:
# No adjustment for strike-slip or reverse faulting
return 0.0
[docs] def get_depth_term(self, C, rup):
"""
Returns the top-of-rupture depth scaling (equation 1)
"""
return C["bcr"] * rup.ztor
[docs] def get_distance_term(self, C, dists, rup):
"""
Returns the distance scaling term defined in equation 3
"""
x_ij = dists.rrup
gn_exp = np.exp(C["c1"] + 6.5 * C["c2"])
# Geometric attenuation scaling described in equation 6
g_n = C["gcrN"] * np.log(self.CONSTANTS["xcro"] + 30. + gn_exp) *\
np.ones_like(x_ij)
idx = x_ij <= 30.0
if np.any(idx):
g_n[idx] = C["gcrN"] * np.log(self.CONSTANTS["xcro"] +
x_ij[idx] + gn_exp)
# equation 5
c_m = min(rup.mag, self.CONSTANTS["m_c"])
# equation 4
r_ij = self.CONSTANTS["xcro"] + x_ij + np.exp(C["c1"] + C["c2"] * c_m)
return C["gcr"] * np.log(r_ij) + C["gcrL"] * np.log(x_ij + 200.0) +\
g_n + C["ecr"] * x_ij + C["ecrV"] * dists.rvolc + C["gamma_S"]
[docs] def add_site_amplification(self, C, C_SITE, sites, sa_rock, idx, rup):
"""
Applies the site amplification scaling defined in equations from 10
to 15
"""
n_sites = sites.vs30.shape
# Convert from reference rock to hard rock
hard_rock_sa = sa_rock - C["lnSC1AM"]
# Gets the elastic site amplification ratio
ln_a_n_max = self._get_ln_a_n_max(C, n_sites, idx, rup)
# Retrieves coefficients needed to determine smr
sreff, sreffc, f_sr = self._get_smr_coeffs(C, C_SITE, idx, n_sites,
hard_rock_sa)
snc = np.zeros(n_sites)
alpha = self.CONSTANTS["alpha"]
beta = self.CONSTANTS["beta"]
smr = np.zeros(n_sites)
sa_soil = hard_rock_sa + ln_a_n_max
# Get lnSF
ln_sf = self._get_ln_sf(C, C_SITE, idx, n_sites, rup)
lnamax_idx = np.exp(ln_a_n_max) < 1.25
not_lnamax_idx = np.logical_not(lnamax_idx)
for i in range(1, 5):
idx_i = idx[i]
if not np.any(idx_i):
# No sites of the given site class
continue
idx2 = np.logical_and(lnamax_idx, idx_i)
if np.any(idx2):
# Use the approximate method for SRC and SNC
c_a = C_SITE["LnAmax1D{:g}".format(i)] /\
(np.log(beta) - np.log(sreffc[idx2] ** alpha + beta))
c_b = -c_a * np.log(sreffc[idx2] ** alpha + beta)
snc[idx2] = np.exp((c_a * (alpha - 1.) *
np.log(beta) * np.log(10.0 * beta) -
np.log(10.0) * (c_b + ln_sf[idx2])) /
(c_a * (alpha * np.log(10.0 * beta) -
np.log(beta))))
# For the cases when ln_a_n_max >= 1.25
idx2 = np.logical_and(not_lnamax_idx, idx_i)
if np.any(idx2):
snc[idx2] = (np.exp((ln_a_n_max[idx2] *
np.log(sreffc[idx2] ** alpha + beta) -
ln_sf[idx2] * np.log(beta)) /
C_SITE["LnAmax1D{:g}".format(i)]) - beta) **\
(1.0 / alpha)
smr[idx_i] = sreff[idx_i] * (snc[idx_i] / sreffc[idx_i]) *\
f_sr[idx_i]
# For the cases when site class = i and SMR != 0
idx2 = np.logical_and(idx_i, np.fabs(smr) > 0.0)
if np.any(idx2):
sa_soil[idx2] += (-C_SITE["LnAmax1D{:g}".format(i)] *
(np.log(smr[idx2] ** alpha + beta) -
np.log(beta)) /
(np.log(sreffc[idx2] ** alpha + beta) -
np.log(beta)))
return sa_soil
def _get_smr_coeffs(self, C, C_SITE, idx, n_sites, sa_rock):
"""
Returns the SReff and SReffC terms needed for equation 14 and 15
"""
# Get SR
sreff = np.zeros(n_sites)
sreffc = np.zeros(n_sites)
f_sr = np.zeros(n_sites)
for i in range(1, 5):
sreff[idx[i]] += (np.exp(sa_rock[idx[i]]) * self.IMF[i])
sreffc[idx[i]] += (C_SITE["Src1D{:g}".format(i)] * self.IMF[i])
# Get f_SR
f_sr[idx[i]] += C_SITE["fsr{:g}".format(i)]
return sreff, sreffc, f_sr
def _get_ln_a_n_max(self, C, n_sites, idx, rup):
"""
Defines the rock site amplification defined in equations 10a and 10b
"""
ln_a_n_max = C["lnSC1AM"] * np.ones(n_sites)
for i in [2, 3, 4]:
if np.any(idx[i]):
ln_a_n_max[idx[i]] += C["S{:g}".format(i)]
return ln_a_n_max
def _get_ln_sf(self, C, C_SITE, idx, n_sites, rup):
"""
Returns the log SF term required for equation 12
"""
ln_sf = np.zeros(n_sites)
for i in range(1, 5):
ln_sf_i = (C["lnSC1AM"] - C_SITE["LnAmax1D{:g}".format(i)])
if i > 1:
ln_sf_i += C["S{:g}".format(i)]
ln_sf[idx[i]] += ln_sf_i
return ln_sf
def _get_site_classification(self, vs30):
"""
Define the site class categories based on Vs30. Returns a
vector of site class values and a dictionary containing logical
vectors for each of the site classes
"""
site_class = np.ones(vs30.shape, dtype=int)
idx = {}
idx[1] = vs30 > 600.
idx[2] = np.logical_and(vs30 > 300., vs30 <= 600.)
idx[3] = np.logical_and(vs30 > 200., vs30 <= 300.)
idx[4] = vs30 <= 200.
for i in [2, 3, 4]:
site_class[idx[i]] = i
return site_class, idx
[docs] def get_stddevs(self, C, n_sites, idx, stddev_types):
"""
Retuns the standard deviation
"""
stddevs = []
phi = C["sigma"] + np.zeros(n_sites)
tau = C["tau"] + np.zeros(n_sites)
for stddev_type in stddev_types:
assert stddev_type in self.DEFINED_FOR_STANDARD_DEVIATION_TYPES
if stddev_type == const.StdDev.TOTAL:
stddevs.append(np.sqrt(phi ** 2. + tau ** 2.))
elif stddev_type == const.StdDev.INTRA_EVENT:
stddevs.append(phi)
elif stddev_type == const.StdDev.INTER_EVENT:
stddevs.append(tau)
return stddevs
# Coefficients taken from Excel spreadsheet provided by the author
COEFFS = CoeffsTable(sa_damping=5, table="""\
imt c1 c2 ccr cum dcr FN_CR FRV_UM FN_UM FUM bcr gcr gUM gcrN gcrL ecr eum ecrV gamma_S lnSC1AM S2 S3 S4 sigma tau sigma_T sc1_sigma_S sc1_tau_S sc1_sigma_ST sc2_sigma_S sc2_tau_S sc2_sigma_ST sc3_sigma_S sc3_tau_S sc3_sigma_ST sc4_sigma_S sc4_tau_S sc4_sigma_ST
pga -3.224 0.900 1.07312 1.07312 0.20000 0.312820000000000 -0.20236 0.251935430429416 0.0 0.009070000000000 -1.26034 -1.09985 -0.49919 1.265640000000000 -0.00794 -0.010830000000000 -0.00628 -9.087242726439600 0.323 0.28877 0.12210 0.20813 0.556 0.391 0.680 0.4153 0.3581 0.5484 0.4363 0.3620 0.5669 0.4235 0.3486 0.5485 0.4426 0.2778 0.5225
0.005 -3.224 0.900 1.07312 1.07312 0.20000 0.312820000000000 -0.20236 0.251935430429416 0.0 0.009070000000000 -1.26034 -1.09985 -0.49919 1.265640000000000 -0.00794 -0.010830000000000 -0.00628 -9.087242726439600 0.323 0.28877 0.12210 0.20813 0.556 0.391 0.680 0.4153 0.3581 0.5484 0.4363 0.3620 0.5669 0.4235 0.3486 0.5485 0.4426 0.2778 0.5225
0.010 -3.357 0.909 1.07850 1.07850 0.20000 0.315680000000000 -0.21425 0.258372637504816 0.0 0.009070000000000 -1.26949 -1.10720 -0.48684 1.241487880676370 -0.00772 -0.010580000000000 -0.00629 -9.052058022318450 0.205 0.29991 0.10321 0.21925 0.556 0.390 0.679 0.4142 0.3581 0.5475 0.4356 0.3627 0.5668 0.4230 0.3482 0.5479 0.4419 0.2804 0.5233
0.020 -3.552 0.927 1.07254 1.07254 0.20000 0.318540000000000 -0.22125 0.260697166781076 0.0 0.009070000000000 -1.28775 -1.11766 -0.44645 1.198941940195550 -0.00756 -0.010505216376958 -0.00634 -8.857750749461700 0.083 0.29778 0.11109 0.21713 0.555 0.396 0.682 0.4123 0.3581 0.5460 0.4344 0.3646 0.5671 0.4226 0.3483 0.5477 0.4418 0.2823 0.5243
0.030 -3.640 0.937 1.05736 1.05736 0.20000 0.320220000000000 -0.22306 0.260150407731831 0.0 0.009070000000000 -1.29619 -1.11392 -0.42175 1.186537103020870 -0.00788 -0.011005243836384 -0.00647 -8.627942448779750 0.041 0.22535 0.08972 0.15948 0.553 0.408 0.687 0.4089 0.3625 0.5465 0.4307 0.3580 0.5601 0.4207 0.3544 0.5501 0.4415 0.2857 0.5259
0.040 -3.758 0.944 1.03568 1.03568 0.20000 0.321410000000000 -0.22334 0.258909018258196 0.0 0.009070000000000 -1.25147 -1.07971 -0.37622 1.142133802990220 -0.00863 -0.011487892527436 -0.00681 -8.398219253745690 0.034 0.15865 0.04037 0.07029 0.558 0.438 0.710 0.4072 0.3832 0.5591 0.4275 0.3666 0.5632 0.4208 0.3544 0.5501 0.4435 0.2881 0.5289
0.050 -3.826 0.948 1.00578 1.00578 0.20000 0.321850000000000 -0.22297 0.257458255516038 0.0 0.009070000000000 -1.14724 -0.98606 -0.43576 1.141373265957230 -0.00988 -0.012291799210266 -0.00710 -8.153070000000000 0.046 0.08257 -0.05934 -0.03529 0.564 0.460 0.728 0.4059 0.4052 0.5735 0.4294 0.3760 0.5708 0.4168 0.3508 0.5448 0.4446 0.2945 0.5332
0.060 -3.890 0.956 0.98413 0.98413 0.20000 0.324710000000000 -0.22229 0.255956496918906 0.0 0.009070000000000 -1.09126 -0.93901 -0.46789 1.157491249527060 -0.01075 -0.012739417109939 -0.00724 -8.003270000000000 0.069 0.05542 -0.10956 -0.09017 0.577 0.481 0.751 0.4066 0.4341 0.5948 0.4304 0.3911 0.5816 0.4156 0.3453 0.5404 0.4415 0.3082 0.5384
0.070 -3.965 0.967 0.98059 0.98059 0.20000 0.328600000000000 -0.22145 0.254464785542036 0.0 0.009070000000000 -1.04676 -0.90036 -0.50260 1.186791451212470 -0.01130 -0.013080135792042 -0.00734 -8.047830000000000 0.098 0.04177 -0.13262 -0.09583 0.599 0.488 0.772 0.4171 0.4517 0.6148 0.4329 0.4057 0.5933 0.4215 0.3605 0.5546 0.4407 0.3378 0.5552
0.080 -4.055 0.980 0.98634 0.98634 0.20000 0.329844812501483 -0.22053 0.253008191936176 0.0 0.009070000000000 -1.01364 -0.87064 -0.53802 1.223378486073060 -0.01164 -0.013312834482819 -0.00741 -8.196910000000000 0.132 0.06040 -0.13537 -0.07321 0.616 0.481 0.782 0.4255 0.4745 0.6374 0.4361 0.4247 0.6088 0.4220 0.3821 0.5693 0.4351 0.3557 0.5619
0.090 -4.153 0.995 0.99121 0.99121 0.20000 0.336373431249984 -0.21956 0.251596694854290 0.0 0.009070000000000 -0.99232 -0.85080 -0.57301 1.264040913700510 -0.01182 -0.013454810030999 -0.00746 -8.350790000000000 0.169 0.09304 -0.12513 -0.03629 0.629 0.477 0.789 0.4315 0.4884 0.6517 0.4417 0.4445 0.6266 0.4161 0.3958 0.5743 0.4374 0.3784 0.5784
0.100 -4.255 1.009 1.00033 1.00033 0.20000 0.340964949991461 -0.21858 0.250233442665255 0.0 0.009070000000000 -0.98315 -0.84175 -0.60839 1.305257984164180 -0.01184 -0.013478869951054 -0.00749 -8.524300000000000 0.208 0.15210 -0.07090 0.01110 0.641 0.460 0.789 0.4371 0.4932 0.6590 0.4478 0.4485 0.6338 0.4241 0.4229 0.5989 0.4402 0.4182 0.6072
0.120 -4.466 1.040 1.03440 1.03440 0.20000 0.345623419543235 -0.21661 0.247649908306841 0.0 0.009579251858235 -0.97311 -0.82925 -0.67012 1.392835221341830 -0.01174 -0.013487990193394 -0.00751 -9.054780000000000 0.288 0.26161 0.00261 0.10298 0.657 0.437 0.789 0.4515 0.4927 0.6683 0.4476 0.4641 0.6448 0.4332 0.4768 0.6442 0.4594 0.4216 0.6236
0.140 -4.677 1.070 1.08388 1.08388 0.20000 0.345964793275106 -0.21469 0.245243568556044 0.0 0.010548033973543 -0.98324 -0.83484 -0.72524 1.477695244453800 -0.01142 -0.013335243001700 -0.00748 -9.676400000000000 0.370 0.35467 0.11202 0.21434 0.663 0.414 0.782 0.4559 0.4906 0.6697 0.4553 0.4758 0.6586 0.4347 0.4592 0.6323 0.4697 0.3869 0.6085
0.150 -4.781 1.085 1.10602 1.10602 0.20000 0.344999481637875 -0.21374 0.244100676565134 0.0 0.011218017682426 -0.99261 -0.84153 -0.74977 1.518814263594270 -0.01123 -0.013231344388198 -0.00746 -9.965950000000000 0.412 0.39410 0.15670 0.25574 0.666 0.402 0.778 0.4594 0.4815 0.6655 0.4547 0.4817 0.6624 0.4380 0.4519 0.6294 0.4707 0.3867 0.6092
0.160 -4.883 1.100 1.12674 1.12674 0.20000 0.343465417412490 -0.21282 0.242994687815017 0.0 0.011700559401692 -1.00423 -0.85033 -0.77240 1.558838786529320 -0.01101 -0.013114569632387 -0.00743 -10.241600000000000 0.453 0.42523 0.19678 0.29811 0.671 0.384 0.773 0.4647 0.4725 0.6627 0.4607 0.4820 0.6668 0.4419 0.4617 0.6391 0.4744 0.3788 0.6071
0.180 -5.085 1.129 1.16455 1.16455 0.20000 0.339107164829684 -0.21102 0.240884343952180 0.0 0.012332909878059 -1.03299 -0.87328 -0.81285 1.634757732666460 -0.01054 -0.012840476069927 -0.00735 -10.752000000000000 0.535 0.47642 0.27414 0.36360 0.680 0.380 0.779 0.4721 0.4640 0.6619 0.4697 0.4830 0.6737 0.4514 0.4660 0.6488 0.4784 0.3903 0.6174
0.200 -5.233 1.151 1.19837 1.19837 0.20000 0.333545099152249 -0.20929 0.238895952015333 0.0 0.013457084876985 -1.06493 -0.89945 -0.84625 1.706045133394880 -0.01007 -0.012558176324852 -0.00725 -11.224920000000000 0.606 0.51169 0.33785 0.42568 0.692 0.359 0.780 0.4722 0.4708 0.6668 0.4800 0.4932 0.6882 0.4493 0.4452 0.6325 0.4825 0.4010 0.6274
0.250 -5.229 1.151 1.27000 1.27000 0.20000 0.316942763583834 -0.20526 0.234370987754789 0.0 0.016170000000000 -1.15136 -0.97276 -0.90575 1.862468560694100 -0.00890 -0.011816876971470 -0.00696 -12.262750000000000 0.670 0.55177 0.47868 0.55643 0.694 0.340 0.773 0.4780 0.4650 0.6669 0.4882 0.5436 0.7306 0.4626 0.4192 0.6243 0.4843 0.3942 0.6244
0.300 -5.226 1.151 1.32852 1.32852 0.20000 0.299124734648411 -0.20159 0.230357423625167 0.0 0.018310000000000 -1.23723 -1.04809 -0.93933 1.991624592628190 -0.00784 -0.011099993928613 -0.00661 -13.141810000000000 0.710 0.55325 0.57391 0.65844 0.688 0.344 0.770 0.4906 0.4595 0.6722 0.4991 0.5246 0.7241 0.4833 0.4356 0.6507 0.4930 0.3725 0.6179
0.350 -5.223 1.151 1.37801 1.37801 0.20000 0.281669557686284 -0.19822 0.226742010930494 0.0 0.019800000000000 -1.31736 -1.12034 -0.95403 2.098225191698050 -0.00691 -0.010436248715154 -0.00624 -13.902540000000000 0.719 0.53555 0.63824 0.73536 0.675 0.353 0.762 0.4918 0.4432 0.6621 0.5043 0.4872 0.7012 0.4783 0.4689 0.6698 0.4914 0.3421 0.5988
0.400 -5.221 1.151 1.42087 1.42087 0.20000 0.265209162746323 -0.19511 0.223445783566271 0.0 0.020780000000000 -1.38989 -1.18735 -0.95484 2.186300573524260 -0.00610 -0.009830759490494 -0.00586 -14.573550000000000 0.706 0.50855 0.67945 0.79334 0.667 0.363 0.759 0.4951 0.4296 0.6555 0.5074 0.4722 0.6931 0.4720 0.4821 0.6747 0.4783 0.3466 0.5907
0.450 -5.218 1.151 1.45868 1.45868 0.20000 0.249966529741596 -0.19221 0.220411614509675 0.0 0.021380000000000 -1.45461 -1.24859 -0.94513 2.259336928898360 -0.00541 -0.009285306710044 -0.00548 -15.173840000000000 0.693 0.47913 0.70517 0.83853 0.665 0.359 0.756 0.4927 0.4310 0.6546 0.5074 0.4651 0.6883 0.4944 0.4638 0.6779 0.4787 0.3594 0.5986
0.500 -5.216 1.151 1.49250 1.49250 0.19000 0.235980866829892 -0.18950 0.217596832776964 0.0 0.021680000000000 -1.51219 -1.30439 -0.92771 2.319726842585720 -0.00482 -0.008788091918160 -0.00511 -15.715670000000000 0.681 0.45094 0.72065 0.87505 0.664 0.361 0.756 0.4890 0.4323 0.6527 0.5149 0.4493 0.6834 0.4949 0.4327 0.6574 0.4764 0.3998 0.6220
0.600 -5.213 1.151 1.55102 1.55102 0.17800 0.211571038601388 -0.18454 0.212501492886060 0.0 0.021610000000000 -1.60834 -1.40097 -0.87661 2.410718513737190 -0.00387 -0.007921609069809 -0.00439 -16.662190000000000 0.658 0.39849 0.72820 0.92621 0.669 0.362 0.761 0.4792 0.4402 0.6507 0.5177 0.4472 0.6841 0.4848 0.4577 0.6668 0.4854 0.4454 0.6588
0.700 -5.210 1.151 1.60051 1.60051 0.16200 0.191334420072367 -0.18008 0.207971451817460 0.0 0.020940000000000 -1.68374 -1.48074 -0.81164 2.471787807358160 -0.00315 -0.007193263566848 -0.00374 -17.467280000000000 0.637 0.35254 0.71582 0.95627 0.670 0.370 0.765 0.4665 0.4390 0.6406 0.5044 0.4515 0.6770 0.4816 0.4712 0.6737 0.4707 0.4982 0.6854
0.800 -5.208 1.151 1.64337 1.64337 0.14800 0.174557580328734 -0.17602 0.203882937318070 0.0 0.019870000000000 -1.74306 -1.54711 -0.73894 2.511003526994950 -0.00261 -0.006567782184498 -0.00315 -18.164530000000000 0.618 0.31707 0.69660 0.97619 0.674 0.378 0.773 0.4600 0.4491 0.6429 0.4954 0.4504 0.6696 0.4751 0.4622 0.6628 0.4837 0.5183 0.7089
0.900 -5.206 1.151 1.68118 1.68118 0.13600 0.160616866070260 -0.17229 0.200149921480523 0.0 0.018540000000000 -1.78967 -1.60263 -0.66214 2.534178788957430 -0.00220 -0.006025164356083 -0.00261 -18.777980000000000 0.600 0.28896 0.67366 0.98819 0.673 0.386 0.775 0.4541 0.4579 0.6449 0.4834 0.4559 0.6645 0.4554 0.4781 0.6603 0.4729 0.5358 0.7146
1.000 -5.204 1.151 1.71500 1.71500 0.12500 0.149000000000001 -0.16883 0.196710000000000 0.0 0.017040000000000 -1.82603 -1.64914 -0.58347 2.545383253429940 -0.00189 -0.005550701319875 -0.00214 -19.325620000000000 0.583 0.26685 0.64985 0.99485 0.670 0.390 0.776 0.4475 0.4663 0.6463 0.4720 0.4593 0.6586 0.4415 0.4743 0.6480 0.4709 0.5240 0.7045
1.250 -5.200 1.151 1.78663 1.78663 0.10100 0.127608506221474 -0.16113 0.189110834221069 0.0 0.012920000000000 -1.88672 -1.73810 -0.40237 2.547066045238380 -0.00139 -0.004579681056724 -0.00119 -20.467050000000000 0.545 0.22889 0.59394 0.99639 0.660 0.399 0.771 0.4404 0.4644 0.6400 0.4595 0.4606 0.6506 0.4380 0.4495 0.6276 0.4627 0.4932 0.6763
1.500 -5.196 1.151 1.84515 1.84515 0.08300 0.113816631375633 -0.15447 0.182585464558957 0.0 0.008630000000000 -1.91705 -1.79802 -0.24053 2.525806060977580 -0.00116 -0.003841842878586 -0.00054 -21.388220000000000 0.511 0.20778 0.55120 0.98814 0.656 0.397 0.767 0.4284 0.4546 0.6247 0.4499 0.4688 0.6497 0.4219 0.4633 0.6266 0.4564 0.4942 0.6727
2.000 -5.191 1.151 1.93750 1.93750 0.05300 0.099176091595634 -0.14326 0.171710489423896 0.0 0.000420000000000 -1.93178 -1.86912 -0.03590 2.497512890223570 -0.00109 -0.002838451894609 0.00000 -22.811700000000000 0.454 0.18675 0.50247 0.96004 0.631 0.381 0.737 0.4148 0.4190 0.5896 0.4374 0.4530 0.6297 0.4172 0.4718 0.6298 0.4323 0.5195 0.6758
2.500 -5.187 1.151 2.00913 2.00913 0.02967 0.093289958829313 -0.13399 0.162787337116202 0.0 -0.006870000000000 -1.91787 -1.90393 0.08743 2.478890119830560 -0.00125 -0.002226291691253 0.00000 -23.905120000000000 0.406 0.17749 0.48879 0.92903 0.605 0.376 0.713 0.4003 0.4084 0.5718 0.4343 0.4310 0.6119 0.3955 0.4530 0.6014 0.4128 0.5059 0.6529
3.000 -5.183 1.151 2.06765 2.06765 0.01078 0.091055639803074 -0.12604 0.155180191692808 0.0 -0.013150000000000 -1.89632 -1.92241 0.17408 2.464735244149810 -0.00144 -0.001844374770841 0.00000 -24.794480000000000 0.364 0.16807 0.48245 0.89351 0.591 0.363 0.694 0.4033 0.4052 0.5717 0.4328 0.4152 0.5997 0.3995 0.4030 0.5674 0.4170 0.4886 0.6424
3.500 -5.181 1.151 2.12357 2.12357 0.00000 0.090028774676100 -0.11905 0.148526449306405 0.0 -0.018490000000000 -1.87695 -1.93474 0.24501 2.452674190639510 -0.00159 -0.001603938312670 0.00000 -25.588720000000000 0.327 0.15988 0.47612 0.85603 0.578 0.363 0.682 0.4026 0.3861 0.5578 0.4214 0.3953 0.5778 0.4131 0.3803 0.5615 0.4265 0.4928 0.6517
4.000 -5.178 1.151 2.13734 2.13734 0.00000 0.088965758728213 -0.11280 0.142598301048652 0.0 -0.022980000000000 -1.86168 -1.94360 0.31173 2.442705939276230 -0.00171 -0.001487903801684 0.00000 -26.046350000000000 0.294 0.15433 0.46978 0.81676 0.556 0.376 0.671 0.3804 0.3756 0.5346 0.4264 0.3681 0.5633 0.4060 0.3228 0.5187 0.4259 0.4510 0.6203
4.500 -5.176 1.151 2.13734 2.13734 0.00000 0.087216625027152 -0.10713 0.137242613236394 0.0 -0.026720000000000 -1.85415 -1.95343 0.37949 2.433909230947000 -0.00176 -0.001451647632538 0.00000 -26.370080000000000 0.263 0.15116 0.46344 0.77379 0.542 0.377 0.661 0.3736 0.3595 0.5185 0.4290 0.3589 0.5593 0.3921 0.3069 0.4979 0.4256 0.4386 0.6111
5.000 -5.174 1.151 2.13734 2.13734 0.00000 0.084444226520033 -0.10194 0.132351162212196 0.0 -0.029800000000000 -1.85292 -1.96352 0.45204 2.426728340621800 -0.00177 -0.001498616447348 0.00000 -26.677970000000000 0.235 0.15275 0.45711 0.72802 0.538 0.395 0.667 0.3831 0.3594 0.5253 0.4378 0.3298 0.5481 0.3779 0.3076 0.4873 0.4195 0.3774 0.5643
""")
# Coefficients specific to the site amplification
SITE_COEFFS = CoeffsTable(sa_damping=5, table="""\
imt LnAmax1D1 LnAmax1D2 LnAmax1D3 LnAmax1D4 Src1D1 Src1D2 Src1D3 Src1D4 fsr1 fsr2 fsr3 fsr4
pga 0.650220 0.709730 0.644340 0.404280 8.429000 1.913680 1.117140 0.836440 1.000000 1.000000 1.000000 1.000000
0.005 0.650220 0.709730 0.644340 0.404280 8.429000 1.913680 1.117140 0.836440 1.000000 1.000000 1.000000 1.000000
0.010 0.651810 0.706790 0.646240 0.404280 8.090000 1.882560 1.114440 0.836440 1.000000 1.200000 1.000000 1.000000
0.020 0.653620 0.694650 0.638650 0.387890 6.992000 1.778610 1.124370 0.830000 1.000000 1.300000 1.000000 0.949000
0.030 0.654670 0.687550 0.634210 0.378300 6.350000 1.717810 1.130170 0.826240 1.000000 1.253000 1.000000 0.550000
0.040 0.652850 0.698920 0.606040 0.317370 4.883000 2.052340 1.150800 0.767580 1.000000 1.064000 1.000000 0.477000
0.050 0.672640 0.701370 0.617160 0.309340 5.043000 2.387130 1.239710 0.786320 1.000000 1.120000 1.000000 0.492000
0.060 0.699660 0.724450 0.637970 0.325300 6.271000 2.833990 1.348190 0.837750 1.000000 1.207000 1.000000 0.531000
0.070 0.717130 0.743430 0.654370 0.354120 7.667000 3.294470 1.451810 0.926160 1.000000 1.238000 1.000000 0.613000
0.080 0.716030 0.785980 0.680190 0.392820 9.034000 3.990910 1.583150 1.022280 1.000000 1.360000 1.000000 0.693000
0.090 0.725610 0.797210 0.708890 0.421840 11.251000 4.465760 1.732920 1.118020 1.000000 1.355000 1.000000 0.780000
0.100 0.742000 0.816680 0.718810 0.437360 14.817000 5.045610 1.841340 1.165780 1.000000 1.341000 1.080000 0.816000
0.120 0.762360 0.845230 0.725810 0.472080 14.817000 5.899600 2.030290 1.285510 0.000000 1.195000 1.093000 0.998000
0.140 0.752150 0.782960 0.745250 0.512780 14.817000 5.053530 2.281330 1.398080 0.000000 0.835000 0.948000 0.954000
0.150 0.738190 0.794800 0.761030 0.534320 14.817000 5.204900 2.444130 1.443270 0.000000 0.781000 0.908000 0.942000
0.160 0.719110 0.808610 0.768130 0.550220 14.817000 5.386940 2.580170 1.471770 0.000000 0.738000 0.862000 0.927000
0.180 0.654080 0.843310 0.756900 0.572790 14.817000 5.871650 2.741610 1.546940 0.000000 0.684000 0.745000 0.927000
0.200 0.583950 0.877700 0.717850 0.596740 14.817000 6.573910 2.825870 1.644010 0.000000 0.654000 0.623000 0.949000
0.250 0.583950 0.937670 0.654700 0.611360 14.817000 8.500000 2.718930 1.790130 0.000000 0.683000 0.436000 0.970000
0.300 0.583950 0.950000 0.696190 0.626380 14.817000 10.670300 2.417590 1.823450 0.000000 0.691000 0.400000 0.961000
0.350 0.583950 1.000000 0.779070 0.630120 14.817000 10.670300 2.303750 1.790370 0.000000 0.800000 0.433000 0.948000
0.400 0.583950 1.000000 0.827760 0.647730 14.817000 10.670300 2.236250 1.768440 0.000000 0.876000 0.460000 0.965000
0.450 0.583950 1.000000 0.876450 0.641520 14.817000 10.670300 2.216780 1.675390 0.000000 0.966000 0.496000 0.958000
0.500 0.583950 1.000000 0.925140 0.655820 14.817000 10.670300 2.243380 1.625390 0.000000 1.034000 0.529000 0.984000
0.600 0.583950 1.000000 0.973830 0.686680 14.817000 10.670300 2.805350 1.524530 0.000000 1.206000 0.578000 1.055000
0.700 0.583950 1.000000 1.022520 0.705600 14.817000 10.670300 6.658390 1.397240 0.000000 1.314000 0.578000 1.114000
0.800 0.583950 1.000000 1.071220 0.714290 14.817000 10.670300 30.000000 1.320290 0.000000 1.357000 0.578000 1.175000
0.900 0.583950 1.000000 1.119910 0.703880 14.817000 10.670300 30.000000 1.266370 0.000000 1.357000 0.578000 1.216000
1.000 0.583950 1.000000 1.168600 0.678130 14.817000 10.670300 30.000000 1.226800 0.000000 1.357000 0.000000 1.230000
1.250 0.583950 1.000000 1.217290 0.611190 14.817000 10.670300 30.000000 1.220650 0.000000 0.000000 0.000000 1.192000
1.500 0.583950 1.000000 1.265980 0.547360 14.817000 10.670300 30.000000 1.318050 0.000000 0.000000 0.000000 0.942000
2.000 0.583950 1.000000 1.314670 0.459440 14.817000 10.670300 30.000000 2.124850 0.000000 0.000000 0.000000 0.000000
2.500 0.583950 1.000000 1.363360 0.408460 14.817000 10.670300 30.000000 14.381810 0.000000 0.000000 0.000000 0.000000
3.000 0.583950 1.000000 1.412050 0.364210 14.817000 10.670300 30.000000 14.381810 0.000000 0.000000 0.000000 0.000000
3.500 0.583950 1.000000 1.460750 0.329840 14.817000 10.670300 30.000000 14.381810 0.000000 0.000000 0.000000 0.000000
4.000 0.583950 1.000000 1.509440 0.309120 14.817000 10.670300 30.000000 14.381810 0.000000 0.000000 0.000000 0.000000
4.500 0.583950 1.000000 1.558130 0.292510 14.817000 10.670300 30.000000 14.381810 0.000000 0.000000 0.000000 0.000000
5.000 0.583950 1.000000 1.606820 0.547360 14.817000 10.670300 30.000000 14.381810 0.000000 0.000000 0.000000 0.000000
""")
CONSTANTS = {"m_c": 7.1,
"xcro": 2.0,
"xinto": 10.0, # Used in subduction Interface class
"alpha": 2.0,
"beta": 0.6,
"Imin": 1.0,
"Imax": 12.0,
"m_sc": 6.3} # Used in Subduction Slab class
# Impedence ratio factors for each site class
IMF = {1: (1. + 0.8 * 2.73) / 3.5, # IMF 0.91
2: 3.07 / 3.0,
3: (1. + 0.9 * 1.76) / 2.5,
4: (1. + 0.6 * 2.02) / 3.0}
[docs]class ZhaoEtAl2016AscSiteSigma(ZhaoEtAl2016Asc):
"""
Adaption of the Zhao et al (2016a) GMPE for active shallow crust
events for the case when within-event variability is dependent on site
class
"""
[docs] def get_stddevs(self, C, n_sites, idx, stddev_types):
"""
Returns site class specific standard deviation
"""
stddevs = []
tau = C["tau"] + np.zeros(n_sites)
phi = np.zeros(n_sites)
for i in range(1, 5):
phi[idx[i]] += C["sc{:g}_sigma_S".format(i)]
for stddev_type in stddev_types:
assert stddev_type in self.DEFINED_FOR_STANDARD_DEVIATION_TYPES
if stddev_type == const.StdDev.TOTAL:
stddevs.append(np.sqrt(phi ** 2. + tau ** 2.))
elif stddev_type == const.StdDev.INTRA_EVENT:
stddevs.append(phi)
elif stddev_type == const.StdDev.INTER_EVENT:
stddevs.append(tau)
return stddevs
[docs]class ZhaoEtAl2016UpperMantle(ZhaoEtAl2016Asc):
"""
Adaptation of the Zhao et al. (2016a) GMPE for the upper mantle events
"""
#: Supported tectonic region type is upper mantle
DEFINED_FOR_TECTONIC_REGION_TYPE = const.TRT.UPPER_MANTLE
[docs] def get_sof_term(self, C, rup):
"""
In the case of the upper mantle events separate coefficients
are considered for normal, reverse and strike-slip
"""
if rup.rake <= -45.0 and rup.rake >= -135.0:
# Normal faulting
return C["FN_UM"]
elif rup.rake > 45.0 and rup.rake < 135.0:
# Reverse faulting
return C["FRV_UM"]
else:
# No adjustment for strike-slip faulting
return 0.0
[docs] def get_depth_term(self, C, rup):
"""
No top of rupture depth is considered for upper mantle events
"""
return 0.0
[docs] def get_distance_term(self, C, dists, rup):
"""
Returns the distance attenuation term
"""
x_ij = dists.rrup
gn_exp = np.exp(C["c1"] + 6.5 * C["c2"])
g_n = C["gcrN"] * np.log(self.CONSTANTS["xcro"] + 30. + gn_exp) *\
np.ones_like(x_ij)
idx = x_ij <= 30.0
if np.any(idx):
g_n[idx] = C["gcrN"] * np.log(self.CONSTANTS["xcro"] +
x_ij[idx] + gn_exp)
c_m = min(rup.mag, self.CONSTANTS["m_c"])
r_ij = self.CONSTANTS["xcro"] + x_ij + np.exp(C["c1"] + C["c2"] * c_m)
return C["gUM"] * np.log(r_ij) +\
C["gcrL"] * np.log(x_ij + 200.0) +\
g_n + C["eum"] * x_ij + C["ecrV"] * dists.rvolc + C["gamma_S"]
# For Upper Mantle
SITE_COEFFS = CoeffsTable(sa_damping=5, table="""\
imt LnAmax1D1 LnAmax1D2 LnAmax1D3 LnAmax1D4 Src1D1 Src1D2 Src1D3 Src1D4 fsr1 fsr2 fsr3 fsr4
pga 0.650220 0.709730 0.644340 0.404280 8.429000 1.913680 1.117140 0.836440 1.000000 1.000000 1.000000 1.000000
0.005 0.650220 0.709730 0.644340 0.404280 8.429000 1.913680 1.117140 0.836440 1.000000 1.000000 1.000000 1.000000
0.010 0.651810 0.706790 0.646240 0.404280 8.090000 1.882560 1.114440 0.836440 1.000000 1.000000 1.000000 1.000000
0.020 0.653620 0.694650 0.638650 0.387890 6.992000 1.778610 1.124370 0.830000 1.000000 1.000000 1.000000 1.150000
0.030 0.654670 0.687550 0.634210 0.378300 6.350000 1.717810 1.130170 0.826240 1.000000 1.000000 1.000000 0.800000
0.040 0.652850 0.698920 0.606040 0.317370 4.883000 2.052340 1.150800 0.767580 1.000000 1.000000 1.000000 0.613000
0.050 0.672640 0.701370 0.617160 0.309340 5.043000 2.387130 1.239710 0.786320 1.000000 1.000000 1.000000 0.542000
0.060 0.699660 0.724450 0.637970 0.325300 6.271000 2.833990 1.348190 0.837750 1.000000 1.000000 1.000000 0.534000
0.070 0.717130 0.743430 0.654370 0.354120 7.667000 3.294470 1.451810 0.926160 1.000000 1.000000 1.000000 0.583000
0.080 0.716030 0.785980 0.680190 0.392820 9.034000 3.990910 1.583150 1.022280 1.000000 1.000000 0.830000 0.643000
0.090 0.725610 0.797210 0.708890 0.421840 11.251000 4.465760 1.732920 1.118020 1.000000 1.000000 0.706000 0.674000
0.100 0.742000 0.816680 0.718810 0.437360 14.817000 5.045610 1.841340 1.165780 1.000000 1.000000 0.800000 0.694000
0.120 0.762360 0.845230 0.725810 0.472080 14.817000 5.899600 2.030290 1.285510 0.000000 1.000000 0.759000 0.763000
0.140 0.752150 0.782960 0.745250 0.512780 14.817000 5.053530 2.281330 1.398080 0.000000 0.767000 0.715000 0.684000
0.150 0.738190 0.794800 0.761030 0.534320 14.817000 5.204900 2.444130 1.443270 0.000000 0.708000 0.686000 0.645000
0.160 0.719110 0.808610 0.768130 0.550220 14.817000 5.386940 2.580170 1.471770 0.000000 0.657000 0.644000 0.610000
0.180 0.654080 0.843310 0.756900 0.572790 14.817000 5.871650 2.741610 1.546940 0.000000 0.568000 0.549000 0.532000
0.200 0.583950 0.877700 0.717850 0.596740 14.817000 6.573910 2.825870 1.644010 0.000000 0.509000 0.434000 0.468000
0.250 0.583950 0.937670 0.654700 0.611360 14.817000 8.500000 2.718930 1.790130 0.000000 0.433000 0.292000 0.291000
0.300 0.583950 0.950000 0.696190 0.626380 14.817000 10.670300 2.417590 1.823450 0.000000 0.337000 0.275000 0.275000
0.350 0.583950 1.000000 0.779070 0.630120 14.817000 10.670300 2.303750 1.790370 0.000000 0.337000 0.295000 0.295000
0.400 0.583950 1.000000 0.827760 0.647730 14.817000 10.670300 2.236250 1.768440 0.000000 0.337000 0.293000 0.293000
0.450 0.583950 1.000000 0.876450 0.641520 14.817000 10.670300 2.216780 1.675390 0.000000 0.000000 0.000000 0.000000
0.500 0.583950 1.000000 0.925140 0.655820 14.817000 10.670300 2.243380 1.625390 0.000000 0.000000 0.000000 0.000000
0.600 0.583950 1.000000 0.973830 0.686680 14.817000 10.670300 2.805350 1.524530 0.000000 0.000000 0.000000 0.000000
0.700 0.583950 1.000000 1.022520 0.705600 14.817000 10.670300 6.658390 1.397240 0.000000 0.000000 0.000000 0.000000
0.800 0.583950 1.000000 1.071220 0.714290 14.817000 10.670300 30.000000 1.320290 0.000000 0.000000 0.000000 0.000000
0.900 0.583950 1.000000 1.119910 0.703880 14.817000 10.670300 30.000000 1.266370 0.000000 0.000000 0.000000 0.000000
1.000 0.583950 1.000000 1.168600 0.678130 14.817000 10.670300 30.000000 1.226800 0.000000 0.000000 0.000000 0.000000
1.250 0.583950 1.000000 1.217290 0.611190 14.817000 10.670300 30.000000 1.220650 0.000000 0.000000 0.000000 0.000000
1.500 0.583950 1.000000 1.265980 0.547360 14.817000 10.670300 30.000000 1.318050 0.000000 0.000000 0.000000 0.000000
2.000 0.583950 1.000000 1.314670 0.459440 14.817000 10.670300 30.000000 2.124850 0.000000 0.000000 0.000000 0.000000
2.500 0.583950 1.000000 1.363360 0.408460 14.817000 10.670300 30.000000 14.381810 0.000000 0.000000 0.000000 0.000000
3.000 0.583950 1.000000 1.412050 0.364210 14.817000 10.670300 30.000000 14.381810 0.000000 0.000000 0.000000 0.000000
3.500 0.583950 1.000000 1.460750 0.329840 14.817000 10.670300 30.000000 14.381810 0.000000 0.000000 0.000000 0.000000
4.000 0.583950 1.000000 1.509440 0.309120 14.817000 10.670300 30.000000 14.381810 0.000000 0.000000 0.000000 0.000000
4.500 0.583950 1.000000 1.558130 0.292510 14.817000 10.670300 30.000000 14.381810 0.000000 0.000000 0.000000 0.000000
5.000 0.583950 1.000000 1.606820 0.547360 14.817000 10.670300 30.000000 14.381810 0.000000 0.000000 0.000000 0.000000
""")
[docs]class ZhaoEtAl2016UpperMantleSiteSigma(ZhaoEtAl2016UpperMantle):
"""
Adaption of the Zhao et al (2016a) GMPE for upper mantle events for the
case when within-event variability is dependent on site class
"""
[docs] def get_stddevs(self, C, n_sites, idx, stddev_types):
"""
Returns site class specific standard deviation
"""
stddevs = []
tau = C["tau"] + np.zeros(n_sites)
phi = np.zeros(n_sites)
for i in range(1, 5):
phi[idx[i]] += C["sc{:g}_sigma_S".format(i)]
for stddev_type in stddev_types:
assert stddev_type in self.DEFINED_FOR_STANDARD_DEVIATION_TYPES
if stddev_type == const.StdDev.TOTAL:
stddevs.append(np.sqrt(phi ** 2. + tau ** 2.))
elif stddev_type == const.StdDev.INTRA_EVENT:
stddevs.append(phi)
elif stddev_type == const.StdDev.INTER_EVENT:
stddevs.append(tau)
return stddevs
[docs]class ZhaoEtAl2016SInter(ZhaoEtAl2016Asc):
"""
Implements the subduction interface GMPE of Zhao et al (2016b)
Zhao, J. X., Liang, X., Jiang, F., Xing, H., Zhu, M., Hou, R., Zhang, Y.,
Lan, X., Rhoades, D. A., Irikura, K., Fukushima, Y.,
Somerville, P. (2016b), "Ground-Motion Prediction Equations for
Subduction Interface Earthquakes in Japan Using Site Class and Simple
Geometric Attenuation Functions", Bulletin of the Seismological
Society of America, 106(4), 1518-1534
Main version with standard deviations independent of site term
"""
#: Supported tectonic region type is subduction interface
DEFINED_FOR_TECTONIC_REGION_TYPE = const.TRT.SUBDUCTION_INTERFACE
#: Required rupture parameters are magnitude and top-of-rupture depth
REQUIRES_RUPTURE_PARAMETERS = set(('mag', 'ztor'))
[docs] def get_magnitude_scaling_term(self, C, rup):
"""
Returns magnitude scaling term, which is dependent on top of rupture
depth - as described in equations 1 and 2
"""
if rup.ztor > 25.0:
# Deep interface events
c_int = C["cint"]
else:
c_int = C["cintS"]
if rup.mag <= self.CONSTANTS["m_c"]:
return c_int * rup.mag
else:
return (c_int * self.CONSTANTS["m_c"]) +\
(C["dint"] * (rup.mag - self.CONSTANTS["m_c"]))
[docs] def get_sof_term(self, C, rup):
"""
No style of faulting dependence here
"""
return 0.0
[docs] def get_depth_term(self, C, rup):
"""
Returns depth term (dependent on top of rupture depth) as given
in equations 1 and 2
"""
return (C["bint"] * rup.ztor)
[docs] def get_distance_term(self, C, dists, rup):
"""
Returns distance scaling term, dependent on top of rupture depth,
as described in equation 6
"""
x_ij = dists.rrup
# Get r_ij - distance for geometric spreading (equations 4 & 5)
c_m = min(rup.mag, self.CONSTANTS["m_c"])
r_ij = self.CONSTANTS["xinto"] + x_ij +\
np.exp(C["alpha"] + C["beta"] * c_m)
# Get factors common to both shallow and deep
dist_term = C["gint"] * np.log(r_ij) + C["eintV"] * dists.rvolc +\
C["gammaint"]
if rup.ztor < 25.:
# Shallow events have geometric and anelastic attenuation term
dist_term += (C["gintLS"] * np.log(x_ij + 200.0)
+ C["eintS"] * x_ij) + C["gamma_ints"]
else:
# Deep events do not have an anelastic attenuation term
dist_term += (C["gintLD"] * np.log(x_ij + 200.0))
return dist_term
def _get_ln_a_n_max(self, C, n_sites, idx, rup):
"""
Defines the rock site amplification defined in equations 7a and 7b
For events deeper than 25 km the rock-site factor is slightly different
for site classes SCII, SCIII, SCIV
"""
ln_a_n_max = C["lnSC1AM"] * np.ones(n_sites)
for i in [2, 3, 4]:
if np.any(idx[i]):
# For deep events site classes 5, 6 and 7 are used
if rup.ztor > 25.0:
loc = i + 3
else:
loc = i
ln_a_n_max[idx[i]] += C["S{:g}".format(loc)]
return ln_a_n_max
def _get_ln_sf(self, C, C_SITE, idx, n_sites, rup):
"""
Returns the log SF term required for equation 9
For events deeper than 25 km the rock-site factor is slightly different
for site classes SCII, SCIII, SCIV
"""
ln_sf = np.zeros(n_sites)
for i in range(1, 5):
ln_sf_i = (C["lnSC1AM"] - C_SITE["LnAmax1D{:g}".format(i)])
if i > 1:
if rup.ztor > 25.0:
# For deep events site classes 5, 6, and 7 are used
loc = i + 3
else:
# For shallow events the conventional approach applies
loc = i
ln_sf_i += C["S{:g}".format(loc)]
ln_sf[idx[i]] += ln_sf_i
return ln_sf
# Coefficients table taken from spreadsheet supplied by the author
COEFFS = CoeffsTable(sa_damping=5, table="""\
imt alpha beta cint cintS dint gamma_ints bint gint gintLD gintLS eintV eintS gammaint S2 S3 S4 S5 S6 S7 lnSC1AM sigma tau sigma_T sc1_sigma_S sc1_tau_S sc1_sigma_ST sc2_sigma_S sc2_tau_T sc2_sigma_ST sc3_sigma_S sc3_tau_S sc3_sigma_ST sc4_sigma_S sc4_tau_S sc4_sigma_ST
pga -5.3011890395444800 1.151 1.09973 1.31479 0.553 -3.89528 0.01999 -2.05587 0.545410000000000 1.133644136349760 -0.01123 -0.00628 -4.49858 0.31288 -0.0043099999999998 0.22838 0.31288 -0.0043099999999998 0.22838 0.3059054757248630 0.553 0.377 0.669 0.398 0.419 0.571 0.417 0.424 0.597 0.409 0.376 0.544 0.415 0.346 0.574
0.005 -5.3011890395444800 1.151 1.09973 1.31479 0.553 -3.89528 0.01999 -2.05587 0.545410000000000 1.133644136349760 -0.01123 -0.00628 -4.49858 0.31288 -0.0043099999999998 0.22838 0.31288 -0.0043099999999998 0.22838 0.3059054757248630 0.553 0.377 0.669 0.398 0.419 0.571 0.417 0.424 0.597 0.409 0.376 0.544 0.415 0.346 0.574
0.010 -5.2884351314221800 1.151 1.09848 1.31739 0.553 -3.89528 0.01999 -2.06565 0.549750000000000 1.133647532106160 -0.01125 -0.00625 -4.45894 0.30850 -0.0111500000000002 0.22305 0.30850 -0.0111500000000002 0.22305 0.2211189175034500 0.554 0.377 0.670 0.397 0.420 0.571 0.417 0.424 0.598 0.409 0.377 0.545 0.415 0.348 0.543
0.020 -5.2756812232998800 1.151 1.09227 1.31919 0.553 -3.89528 0.01999 -2.10231 0.561715000000000 1.133645539711160 -0.01127 -0.00616 -4.25807 0.29297 -0.0216599999999998 0.20892 0.29297 -0.0216599999999998 0.20892 0.1391818272825550 0.553 0.384 0.673 0.395 0.425 0.574 0.417 0.425 0.598 0.408 0.376 0.544 0.416 0.352 0.544
0.030 -5.2682206653106900 1.151 1.10690 1.34096 0.553 -3.89528 0.02066 -2.19232 0.578870000000000 1.133638051189290 -0.01158 -0.00572 -3.91800 0.22869 -0.1129091455512500 0.13314 0.22869 -0.0829091455512501 0.16314 0.0928148352701154 0.555 0.397 0.682 0.389 0.444 0.586 0.418 0.422 0.595 0.409 0.386 0.553 0.417 0.355 0.546
0.040 -5.2629273151775800 1.151 1.11578 1.38051 0.553 -3.89528 0.02308 -2.24640 0.493330000000000 0.988105377906839 -0.01203 -0.00532 -3.11423 0.16316 -0.1887005077743300 0.06959 0.16316 -0.1887005077743300 0.05959 0.0626979621246270 0.565 0.428 0.709 0.387 0.474 0.611 0.420 0.434 0.604 0.413 0.393 0.557 0.420 0.362 0.555
0.050 -5.2588214738333900 1.151 1.10234 1.43246 0.553 -3.89528 0.02709 -2.29341 0.491000000000000 0.904439734658832 -0.01256 -0.00503 -2.76035 0.12129 -0.2283118706921090 0.02845 0.12129 -0.2283118706921090 0.00845 0.0360112352938003 0.570 0.463 0.734 0.387 0.486 0.620 0.422 0.440 0.608 0.409 0.384 0.554 0.422 0.375 0.562
0.060 -5.2554667571883900 1.151 1.08611 1.46238 0.553 -3.89528 0.02972 -2.31172 0.508475000000000 0.887664211721920 -0.01312 -0.00528 -2.64088 0.12345 -0.2192024788530120 0.00070 0.12345 -0.2192024788530120 0.00070 0.0372315296271544 0.583 0.488 0.760 0.397 0.500 0.632 0.416 0.449 0.614 0.401 0.379 0.550 0.423 0.388 0.571
0.070 -5.2526303846795600 1.151 1.07291 1.47120 0.553 -3.89528 0.03207 -2.31101 0.527500000000000 0.904912893427476 -0.01359 -0.00569 -2.65622 0.13970 -0.1900877585085220 -0.00949 0.13970 -0.1900877585085220 0.00551 0.0491648511532895 0.602 0.501 0.784 0.403 0.531 0.658 0.412 0.466 0.628 0.394 0.377 0.558 0.420 0.410 0.587
0.080 -5.2501734070552700 1.151 1.06383 1.46432 0.553 -3.89461 0.03196 -2.28778 0.545953146186990 0.942062229717542 -0.01382 -0.00619 -2.75266 0.16385 -0.1549929496094120 -0.00485 0.16385 -0.1549929496094120 0.01415 0.0983337731804003 0.614 0.501 0.793 0.413 0.553 0.681 0.412 0.489 0.647 0.389 0.374 0.552 0.423 0.416 0.591
0.090 -5.2480061991992000 1.151 1.05863 1.44695 0.553 -3.90179 0.02972 -2.24680 0.563067703920420 0.986519531877540 -0.01393 -0.00673 -2.89920 0.20501 -0.1125034957053890 0.03047 0.20501 -0.1125034957053890 0.04047 0.1753764826899090 0.625 0.495 0.797 0.422 0.559 0.688 0.412 0.521 0.672 0.391 0.381 0.557 0.427 0.435 0.604
0.100 -5.2460675657110900 1.151 1.05673 1.42323 0.553 -3.90765 0.02789 -2.20410 0.576153228856880 1.035532649999650 -0.01395 -0.00718 -3.07698 0.24449 -0.0750832172171707 0.06080 0.24449 -0.0750832172171707 0.07080 0.2447184227049870 0.637 0.478 0.796 0.429 0.562 0.693 0.418 0.540 0.687 0.394 0.409 0.574 0.426 0.469 0.632
0.120 -5.2427128490660800 1.151 1.06051 1.36833 0.553 -3.91644 0.02470 -2.12011 0.592575859588849 1.135292002539040 -0.01381 -0.00793 -3.48283 0.32284 0.0150002886394358 0.14231 0.32284 0.0150002886394358 0.14031 0.3366973661131230 0.646 0.453 0.789 0.440 0.564 0.697 0.429 0.537 0.684 0.428 0.480 0.635 0.445 0.517 0.670
0.140 -5.2398764765572600 1.151 1.07135 1.31558 0.553 -3.92273 0.02117 -2.04337 0.609834200714791 1.234237180364630 -0.01351 -0.00853 -3.91612 0.40120 0.0969875152845021 0.22696 0.40120 0.0969875152845021 0.20196 0.4219918629306410 0.654 0.412 0.773 0.441 0.556 0.692 0.435 0.560 0.712 0.433 0.440 0.599 0.442 0.507 0.670
0.150 -5.2386070077219000 1.151 1.07856 1.29277 0.553 -3.92526 0.01951 -2.01088 0.619605756588544 1.281341268532020 -0.01333 -0.00879 -4.13481 0.43622 0.1458772263727260 0.25758 0.43622 0.1358772263727260 0.23258 0.4561999113437010 0.659 0.404 0.773 0.450 0.552 0.690 0.440 0.574 0.724 0.420 0.427 0.592 0.440 0.504 0.663
0.160 -5.2374194989329700 1.151 1.08664 1.27319 0.553 -3.92753 0.01793 -1.98301 0.630846000686806 1.326585227368230 -0.01312 -0.00902 -4.35243 0.46736 0.1879291503650610 0.30748 0.46736 0.1729291503650610 0.26248 0.4857829542311720 0.663 0.398 0.774 0.452 0.548 0.688 0.444 0.573 0.726 0.424 0.429 0.589 0.437 0.496 0.661
0.180 -5.2352522910768900 1.151 1.10470 1.24828 0.553 -3.93130 0.01505 -1.94607 0.661957417503245 1.411336369650280 -0.01269 -0.00927 -4.78030 0.51201 0.2515429577444550 0.35974 0.51201 0.2415429577444550 0.31974 0.5337861886971570 0.672 0.387 0.776 0.454 0.534 0.681 0.447 0.557 0.719 0.446 0.436 0.600 0.436 0.488 0.667
0.200 -5.2333136575887900 1.151 1.12443 1.23715 0.553 -3.93446 0.01255 -1.92702 0.699775273836698 1.488543999616950 -0.01223 -0.00942 -5.19439 0.53926 0.3029826435245730 0.40313 0.53926 0.3029826435245730 0.37313 0.5700442187746950 0.678 0.382 0.778 0.462 0.528 0.679 0.450 0.532 0.710 0.441 0.437 0.602 0.432 0.481 0.663
0.250 -5.2292078162446100 1.151 1.17689 1.22387 0.553 -3.94068 0.00769 -1.89876 0.784470723616457 1.652089063569910 -0.01108 -0.00959 -6.15803 0.58602 0.4269231521787280 0.50765 0.58602 0.4269231521787280 0.48765 0.6247924955873390 0.659 0.365 0.753 0.474 0.493 0.650 0.470 0.551 0.723 0.459 0.398 0.572 0.432 0.448 0.637
0.300 -5.2258530995996000 1.151 1.22970 1.22846 0.553 -3.94547 0.00438 -1.89141 0.859387831257168 1.781266639239090 -0.00998 -0.00952 -7.02003 0.60465 0.5161970850888910 0.57779 0.60465 0.5161970850888910 0.57779 0.6505259707420690 0.640 0.348 0.729 0.472 0.471 0.632 0.475 0.509 0.695 0.437 0.421 0.579 0.432 0.442 0.623
0.350 -5.2230167270907800 1.151 1.28058 1.24219 0.553 -3.94943 0.00215 -1.89299 0.923381083222217 1.884375114633360 -0.00898 -0.00933 -7.79153 0.60640 0.5695125616089470 0.63821 0.60640 0.5795125616089470 0.64821 0.6617360212939630 0.634 0.360 0.729 0.468 0.474 0.629 0.478 0.485 0.683 0.444 0.457 0.607 0.432 0.442 0.618
0.400 -5.2205597494664900 1.151 1.32873 1.26077 0.553 -3.95273 0.00000 -1.89531 0.980106570565885 1.967625670226430 -0.00808 -0.00911 -8.49551 0.60277 0.6236634605365610 0.70323 0.60277 0.6336634605365610 0.71323 0.6649532313027650 0.627 0.354 0.720 0.457 0.449 0.610 0.484 0.496 0.683 0.453 0.463 0.607 0.416 0.435 0.607
0.450 -5.2183925416104100 1.151 1.37394 1.28190 0.553 -3.95558 0.00000 -1.90578 1.022203929823320 2.035537053673680 -0.00727 -0.00888 -9.11351 0.58043 0.6581416274761160 0.75082 0.58043 0.6581416274761160 0.75082 0.6654266030879300 0.620 0.363 0.719 0.450 0.423 0.595 0.477 0.495 0.679 0.477 0.447 0.593 0.408 0.418 0.592
0.500 -5.2164539081223000 1.151 1.41630 1.30432 0.553 -3.95795 0.00000 -1.91467 1.058737908474810 2.091423339372400 -0.00656 -0.00866 -9.68515 0.55692 0.6867116556863830 0.79378 0.55692 0.6867116556863830 0.79378 0.6635290829305090 0.612 0.364 0.712 0.445 0.415 0.583 0.470 0.473 0.658 0.472 0.435 0.579 0.409 0.432 0.603
0.600 -5.2130991914773000 1.151 1.49310 1.35019 0.553 -3.96181 0.00000 -1.92743 1.118034334349720 2.176403198118720 -0.00534 -0.00824 -10.68946 0.50969 0.7122000042516460 0.84953 0.50969 0.7122000042516460 0.84953 0.6564608147533650 0.613 0.379 0.720 0.448 0.419 0.576 0.460 0.450 0.636 0.459 0.434 0.580 0.407 0.442 0.598
0.700 -5.2102628189684700 1.151 1.56069 1.39523 0.560 -3.96483 0.00000 -1.93452 1.163018824473690 2.235996913698260 -0.00437 -0.00787 -11.54596 0.46501 0.7123896083711760 0.87979 0.46501 0.7123896083711760 0.87979 0.6475225398028140 0.625 0.393 0.739 0.440 0.429 0.583 0.460 0.454 0.634 0.461 0.448 0.591 0.403 0.454 0.613
0.800 -5.2078058413441800 1.151 1.62055 1.43824 0.580 -3.96729 0.00000 -1.93739 1.197347867256100 2.278310911386920 -0.00359 -0.00755 -12.28722 0.42440 0.6993827558539090 0.89539 0.42440 0.6993827558539090 0.89539 0.6378151725112400 0.628 0.396 0.743 0.441 0.427 0.580 0.460 0.448 0.631 0.457 0.440 0.590 0.407 0.468 0.627
0.900 -5.2056386334881000 1.151 1.67390 1.47881 0.602 -3.96962 0.00000 -1.93725 1.223616172806850 2.308465259401480 -0.00296 -0.00726 -12.93630 0.38841 0.6799630664371980 0.90256 0.38841 0.6799630664371980 0.90256 0.6275299463931780 0.628 0.397 0.743 0.435 0.433 0.586 0.456 0.446 0.622 0.449 0.445 0.593 0.408 0.485 0.633
1.000 -5.2037000000000000 1.151 1.72171 1.51685 0.622 -3.97202 0.00000 -1.93505 1.243670000000000 2.329852618176160 -0.00244 -0.00700 -13.51000 0.35703 0.6582841883037910 0.90528 0.35703 0.6582841883037910 0.90528 0.6166948032035810 0.633 0.403 0.750 0.427 0.453 0.606 0.448 0.445 0.617 0.442 0.422 0.578 0.408 0.498 0.642
1.250 -5.1995941586558200 1.151 1.82188 1.60148 0.667 -3.97949 0.00000 -1.92471 1.272450000000000 2.358529248389570 -0.00153 -0.00644 -14.69034 0.29673 0.6192828473123780 0.91793 0.29673 0.6192828473123780 0.91793 0.5871136982629510 0.636 0.404 0.753 0.413 0.457 0.613 0.442 0.457 0.624 0.428 0.428 0.587 0.412 0.507 0.657
1.500 -5.1962394420108100 1.151 1.90081 1.67280 0.705 -3.99048 0.00000 -1.91189 1.285410000000000 2.366526445593360 -0.00097 -0.00597 -15.60298 0.25785 0.5828781633806380 0.92130 0.25785 0.5828781633806380 0.92130 0.5541350591719840 0.644 0.392 0.754 0.416 0.464 0.619 0.448 0.448 0.628 0.418 0.434 0.604 0.418 0.505 0.653
2.000 -5.1909460918777000 1.151 2.01482 1.78372 0.768 -4.02652 0.00000 -1.88859 1.288300000000000 2.355357695068850 -0.00043 -0.00518 -16.90010 0.22262 0.5262044723719100 0.91709 0.22262 0.5262044723719100 0.91709 0.4823188293988830 0.635 0.382 0.741 0.409 0.464 0.623 0.437 0.458 0.633 0.406 0.470 0.625 0.413 0.524 0.665
2.500 -5.1868402505335200 1.151 2.08892 1.86241 0.820 -4.08299 0.00000 -1.87252 1.277270000000000 2.331055616102490 -0.00023 -0.00451 -17.73658 0.21840 0.4871991042586300 0.90551 0.21840 0.4871991042586300 0.90551 0.4108485356913110 0.619 0.393 0.734 0.398 0.441 0.599 0.429 0.443 0.623 0.387 0.446 0.612 0.414 0.529 0.676
3.000 -5.1834855338885100 1.151 2.13574 1.91711 0.863 -4.15939 0.00000 -1.86345 1.260490000000000 2.304087350833550 -0.00016 -0.00393 -18.27135 0.21595 0.4569807448342530 0.88668 0.21595 0.4569807448342530 0.88668 0.3476970039189440 0.599 0.385 0.712 0.390 0.412 0.572 0.423 0.441 0.623 0.369 0.400 0.565 0.412 0.524 0.668
3.500 -5.1806491613796900 1.151 2.16254 1.95322 0.902 -4.25418 0.00000 -1.85969 1.241065000000000 2.277938604223720 0.00000 -0.00344 -18.59264 0.21595 0.4280893607115190 0.85880 0.21595 0.4280893607115190 0.85880 0.2981089167617070 0.581 0.376 0.692 0.386 0.394 0.553 0.410 0.424 0.608 0.376 0.353 0.536 0.402 0.535 0.672
4.000 -5.1781921837553900 1.151 2.17393 1.97450 0.935 -4.36584 0.00000 -1.85953 1.220255000000000 2.253713119173820 0.00000 -0.00302 -18.75474 0.21595 0.3952406635315720 0.82025 0.21595 0.3952406635315720 0.82025 0.2652398018598770 0.568 0.377 0.682 0.377 0.376 0.541 0.410 0.409 0.603 0.362 0.365 0.540 0.394 0.542 0.671
4.500 -5.1760249758993200 1.151 2.17301 1.98361 0.966 -4.49267 0.00000 -1.86151 1.198585000000000 2.231608900335380 0.00000 -0.00267 -18.79351 0.21595 0.3548841036361290 0.76990 0.21595 0.3548841036361290 0.76990 0.2508917433151500 0.551 0.376 0.667 0.360 0.363 0.528 0.412 0.393 0.591 0.368 0.349 0.519 0.385 0.530 0.652
5.000 -5.1740863424112100 1.151 2.16199 1.98255 0.994 -4.63313 0.00000 -1.86451 1.176285000000000 2.211490595470470 0.00000 -0.00240 -18.73393 0.21595 0.3046830024209600 0.70705 0.21595 0.3046830024209600 0.70705 0.2365436847704230 0.563 0.374 0.675 0.361 0.346 0.525 0.447 0.408 0.615 0.381 0.293 0.476 0.380 0.514 0.633
""")
# Coefficients specific to the site amplification
SITE_COEFFS = CoeffsTable(sa_damping=5, table="""\
imt LnAmax1D1 LnAmax1D2 LnAmax1D3 LnAmax1D4 Src1D1 Src1D2 Src1D3 Src1D4 fsr1 fsr2 fsr3 fsr4
pga 0.650220 0.709730 0.644340 0.404280 8.429000 1.913680 1.117140 0.836440 1.0000 1.0000 1.1650 1.0000
0.005 0.650220 0.709730 0.644340 0.404280 8.429000 1.913680 1.117140 0.836440 1.0000 1.0000 1.1650 1.0000
0.010 0.651810 0.706790 0.646240 0.404280 8.090000 1.882560 1.114440 0.836440 1.0000 1.0000 0.9440 1.0000
0.020 0.653620 0.694650 0.638650 0.387890 6.992000 1.778610 1.124370 0.830000 1.0000 1.0000 1.0120 1.0000
0.030 0.654670 0.687550 0.634210 0.378300 6.350000 1.717810 1.130170 0.826240 1.0000 0.9990 1.1000 1.0000
0.040 0.652850 0.698920 0.606040 0.317370 4.883000 2.052340 1.150800 0.767580 1.0000 0.8430 0.9590 0.5570
0.050 0.672640 0.701370 0.617160 0.309340 5.043000 2.387130 1.239710 0.786320 1.0000 0.6630 0.8890 0.5430
0.060 0.699660 0.724450 0.637970 0.325300 6.271000 2.833990 1.348190 0.837750 1.0000 0.8410 0.9460 0.5740
0.070 0.717130 0.743430 0.654370 0.354120 7.667000 3.294470 1.451810 0.926160 1.0000 1.0290 1.0060 0.6480
0.080 0.716030 0.785980 0.680190 0.392820 9.034000 3.990910 1.583150 1.022280 1.0000 1.2350 1.0650 0.7210
0.090 0.725610 0.797210 0.708890 0.421840 11.251000 4.465760 1.732920 1.118020 1.0000 1.1440 1.0930 0.8090
0.100 0.742000 0.816680 0.718810 0.437360 14.817000 5.045610 1.841340 1.165780 1.0000 1.0920 1.0770 1.0150
0.120 0.762360 0.845230 0.725810 0.472080 14.817000 5.899600 2.030290 1.285510 0.0000 0.9450 1.0360 0.9720
0.140 0.752150 0.782960 0.745250 0.512780 14.817000 5.053530 2.281330 1.398080 0.0000 0.6240 0.8950 0.9670
0.150 0.738190 0.794800 0.761030 0.534320 14.817000 5.204900 2.444130 1.443270 0.0000 0.5770 0.8600 0.9630
0.160 0.719110 0.808610 0.768130 0.550220 14.817000 5.386940 2.580170 1.471770 0.0000 0.5450 0.8220 0.9530
0.180 0.654080 0.843310 0.756900 0.572790 14.817000 5.871650 2.741610 1.546940 0.0000 0.5270 0.7430 0.9670
0.200 0.583950 0.877700 0.717850 0.596740 14.817000 6.573910 2.825870 1.644010 0.0000 0.5460 0.6630 1.0050
0.250 0.583950 0.937670 0.654700 0.611360 14.817000 8.500000 2.718930 1.790130 0.0000 0.5960 0.4870 1.0450
0.300 0.583950 0.950000 0.696190 0.626380 14.817000 10.670300 2.417590 1.823450 0.0000 0.6230 0.4470 1.0350
0.350 0.583950 1.000000 0.779070 0.630120 14.817000 10.670300 2.303750 1.790370 0.0000 0.7010 0.4730 1.0080
0.400 0.583950 1.000000 0.827760 0.647730 14.817000 10.670300 2.236250 1.768440 0.0000 0.7080 0.4870 1.0070
0.450 0.583950 1.000000 0.876450 0.641520 14.817000 10.670300 2.216780 1.675390 0.0000 0.7370 0.5110 0.9810
0.500 0.583950 1.000000 0.925140 0.655820 14.817000 10.670300 2.243380 1.625390 0.0000 0.7480 0.5360 0.9900
0.600 0.583950 1.000000 0.973830 0.686680 14.817000 10.670300 2.805350 1.524530 0.0000 0.7280 0.5400 1.0160
0.700 0.583950 1.000000 1.022520 0.705600 14.817000 10.670300 6.658390 1.397240 0.0000 0.6340 0.4770 1.0220
0.800 0.583950 1.000000 1.071220 0.714290 14.817000 10.670300 30.000000 1.320290 0.0000 0.0000 0.0000 1.0230
0.900 0.583950 1.000000 1.119910 0.703880 14.817000 10.670300 30.000000 1.266370 0.0000 0.0000 0.0000 0.9970
1.000 0.583950 1.000000 1.168600 0.678130 14.817000 10.670300 30.000000 1.226800 0.0000 0.0000 0.0000 0.9480
1.250 0.583950 1.000000 1.217290 0.611190 14.817000 10.670300 30.000000 1.220650 0.0000 0.0000 0.0000 0.8020
1.500 0.583950 1.000000 1.265980 0.547360 14.817000 10.670300 30.000000 1.318050 0.0000 0.0000 0.0000 0.0000
2.000 0.583950 1.000000 1.314670 0.459440 14.817000 10.670300 30.000000 2.124850 0.0000 0.0000 0.0000 0.0000
2.500 0.583950 1.000000 1.363360 0.408460 14.817000 10.670300 30.000000 14.381810 0.0000 0.0000 0.0000 0.0000
3.000 0.583950 1.000000 1.412050 0.364210 14.817000 10.670300 30.000000 14.381810 0.0000 0.0000 0.0000 0.0000
3.500 0.583950 1.000000 1.460750 0.329840 14.817000 10.670300 30.000000 14.381810 0.0000 0.0000 0.0000 0.0000
4.000 0.583950 1.000000 1.509440 0.309120 14.817000 10.670300 30.000000 14.381810 0.0000 0.0000 0.0000 0.0000
4.500 0.583950 1.000000 1.558130 0.292510 14.817000 10.670300 30.000000 14.381810 0.0000 0.0000 0.0000 0.0000
5.000 0.583950 1.000000 1.606820 0.547360 14.817000 10.670300 30.000000 14.381810 0.0000 0.0000 0.0000 0.0000
""")
[docs]class ZhaoEtAl2016SInterSiteSigma(ZhaoEtAl2016SInter):
"""
Subclass of the Zhao et al. (2016b) subduction interface GMPE for the
case of site-dependent within-event variability
"""
[docs] def get_stddevs(self, C, n_sites, idx, stddev_types):
"""
"""
stddevs = []
tau = C["tau"] + np.zeros(n_sites)
phi = np.zeros(n_sites)
for i in range(1, 5):
phi[idx[i]] += C["sc{:g}_sigma_S".format(i)]
for stddev_type in stddev_types:
assert stddev_type in self.DEFINED_FOR_STANDARD_DEVIATION_TYPES
if stddev_type == const.StdDev.TOTAL:
stddevs.append(np.sqrt(phi ** 2. + tau ** 2.))
elif stddev_type == const.StdDev.INTRA_EVENT:
stddevs.append(phi)
elif stddev_type == const.StdDev.INTER_EVENT:
stddevs.append(tau)
return stddevs
[docs]class ZhaoEtAl2016SSlab(ZhaoEtAl2016Asc):
"""
Implements the subduction slab GMPE of Zhao et al (2016c)
Zhao, J. X., Jiang, F., Shi, P., Xing, H., Huang, H., Hou, R.,
Zhang, Y., Yu, P., Lan, X., Rhoades, D. A., Somerville, P. G., Irikura, K.,
Fukushima, Y. (2016c), "Ground-Motion Prediction Equations for
Subduction Slab Earthquakes in Japan Using Site Class and Simple
Geometric Attenuation Functions", Bulletin of the Seismological
Society of America, 106(4), 1535-1551
Main version with standard deviations independent of site term
"""
#: Supported tectonic region type is subduction inslab
DEFINED_FOR_TECTONIC_REGION_TYPE = const.TRT.SUBDUCTION_INTRASLAB
#: Required rupture parameters are magnitude and top of rupture depth.
REQUIRES_RUPTURE_PARAMETERS = set(('mag', 'ztor'))
[docs] def get_magnitude_scaling_term(self, C, rup):
"""
Returns the magnitude scaling defined in equation 1
"""
m_c = self.CONSTANTS["m_c"]
if rup.mag <= m_c:
return C["cSL"] * rup.mag +\
C["cSL2"] * ((rup.mag - self.CONSTANTS["m_sc"]) ** 2.)
else:
return C["cSL"] * m_c +\
C["cSL2"] * ((m_c - self.CONSTANTS["m_sc"]) ** 2.) +\
C["dSL"] * (rup.mag - m_c)
[docs] def get_sof_term(self, C, rup):
"""
No style of faulting dependence here
"""
return 0.0
[docs] def get_depth_term(self, C, rup):
"""
Returns depth term (dependent on top of rupture depth) as given
in equations 1
Note that there is a ztor cap of 100 km that is introduced in the
Fortran code but not mentioned in the original paper!
"""
if rup.ztor > 100.0:
return C["bSLH"] * 100.0
else:
return C["bSLH"] * rup.ztor
[docs] def get_distance_term(self, C, dists, rup):
"""
Returns the distance scaling term in equation 2a
Note that the paper describes a lower and upper cap on Rvolc that
is not found in the Fortran code, and is thus neglected here.
"""
x_ij = dists.rrup
# Get anelastic scaling term in quation 5
if rup.ztor >= 50.:
qslh = C["eSLH"] * (0.02 * rup.ztor - 1.0)
else:
qslh = 0.0
# r_volc = np.copy(dists.rvolc)
# r_volc[np.logical_and(r_volc > 0.0, r_volc <= 12.0)] = 12.0
# r_volc[r_volc >= 80.0] = 80.0
# Get r_ij - distance for geometric spreading (equations 3 and 4)
c_m = min(rup.mag, self.CONSTANTS["m_c"])
r_ij = x_ij + np.exp(C["alpha"] + C["beta"] * c_m)
return C["gSL"] * np.log(r_ij) + \
C["gLL"] * np.log(x_ij + 200.) +\
C["eSL"] * x_ij + qslh * x_ij +\
C["eSLV"] * dists.rvolc + C["gamma"]
# Coefficients table taken from spreadsheet supplied by the author
COEFFS = CoeffsTable(sa_damping=5, table="""\
imt alpha beta cSL cSL2 dSL bSLH gSL gLL eSLV eSL eSLH gamma S2 S3 S4 lnSC1AM sigma tau sigma_T sc1_sigma_S sc1_tau_S sc1_sigma_ST sc2_sigma_S sc2_tau_S sc2_sigma_ST sc3_sigma_S sc3_tau_S sc3_sigma_ST sc4_sigma_S sc4_tau_S sc4_sigma_ST
pga -5.30118903954448 1.151 1.44758 0.37625 0.42646 0.01825668401347 -1.98471 1.12071 -0.01499 -0.00340 -0.000501 -9.879560 0.232000 0.143711 0.147037 0.3230000000000000 0.5870 0.4570 0.7440 0.3981 0.5107 0.6475 0.4174 0.4494 0.6133 0.4091 0.4306 0.5940 0.4152 0.4217 0.5918
0.005 -5.30118903954448 1.151 1.44758 0.37625 0.42646 0.01825668401347 -1.98471 1.12071 -0.01499 -0.00340 -0.000501 -9.879560 0.232000 0.143711 0.147037 0.3230000000000000 0.5870 0.4570 0.7440 0.3981 0.5107 0.6475 0.4174 0.4494 0.6133 0.4091 0.4306 0.5940 0.4152 0.4217 0.5918
0.010 -5.28843513142218 1.151 1.45400 0.38099 0.42075 0.01825668401347 -1.96360 1.03278 -0.01503 -0.00331 -0.000501 -9.512690 0.228860 0.139780 0.132840 0.2050000000000000 0.5870 0.4590 0.7460 0.3968 0.5166 0.6514 0.4173 0.4499 0.6136 0.4086 0.4314 0.5941 0.4149 0.4183 0.5892
0.020 -5.27568122329988 1.151 1.46625 0.39101 0.40055 0.01825668401347 -1.91839 0.94715 -0.01517 -0.00345 -0.000501 -9.266260 0.218250 0.126000 0.144310 0.0830000000000000 0.5870 0.4650 0.7490 0.3951 0.5180 0.6515 0.4172 0.4488 0.6128 0.4082 0.4313 0.5938 0.4158 0.4245 0.5942
0.030 -5.26822066531069 1.151 1.49246 0.41976 0.36433 0.01825668401347 -1.89271 0.93420 -0.01567 -0.00391 -0.000501 -9.331500 0.187370 0.061640 0.066000 0.0410000000000000 0.5880 0.4800 0.7590 0.3894 0.5369 0.6633 0.4179 0.4492 0.6135 0.4088 0.4300 0.5933 0.4166 0.4224 0.5933
0.040 -5.26292731517758 1.151 1.50129 0.45746 0.32072 0.01825668401347 -1.87260 0.97168 -0.01616 -0.00454 -0.000501 -9.507980 0.123320 -0.017050 -0.017140 0.0340000000000000 0.6000 0.5200 0.7940 0.3873 0.5716 0.6905 0.4201 0.4559 0.6200 0.4125 0.4281 0.5945 0.4201 0.4308 0.6017
0.050 -5.25882147383339 1.151 1.51051 0.48601 0.30000 0.01825668401347 -1.85351 1.01492 -0.01676 -0.00510 -0.000501 -9.728580 0.072070 -0.063310 -0.073120 0.0460000000000000 0.6070 0.5550 0.8230 0.3874 0.5857 0.7022 0.4219 0.4794 0.6386 0.4089 0.4286 0.5924 0.4217 0.4387 0.6085
0.060 -5.25546675718839 1.151 1.51380 0.50311 0.31147 0.01825668401347 -1.83395 1.06854 -0.01722 -0.00552 -0.000501 -9.966280 0.027010 -0.101030 -0.119550 0.0690000000000000 0.6230 0.5840 0.8540 0.3970 0.6131 0.7304 0.4160 0.5074 0.6561 0.4012 0.4485 0.6018 0.4227 0.4450 0.6138
0.070 -5.25263038467956 1.151 1.51111 0.50704 0.32673 0.01825668401347 -1.81345 1.13401 -0.01752 -0.00588 -0.000487 -10.225830 -0.006210 -0.146840 -0.160060 0.0980000000000000 0.6380 0.5980 0.8750 0.4027 0.6326 0.7499 0.4121 0.5317 0.6727 0.3942 0.4556 0.6025 0.4203 0.4732 0.6329
0.080 -5.25017340705527 1.151 1.50406 0.50004 0.34289 0.01825668401347 -1.79189 1.20364 -0.01768 -0.00615 -0.000479 -10.551110 0.015650 -0.144790 -0.124340 0.1320000000000000 0.6510 0.5980 0.8840 0.4130 0.6434 0.7646 0.4124 0.5554 0.6918 0.3887 0.4554 0.5987 0.4230 0.4905 0.6477
0.090 -5.24800619919920 1.151 1.49423 0.48071 0.35921 0.01825668401347 -1.76931 1.25808 -0.01772 -0.00635 -0.000476 -10.807210 0.050890 -0.126660 -0.072930 0.1690000000000000 0.6630 0.5850 0.8840 0.4221 0.6307 0.7589 0.4122 0.5677 0.7015 0.3905 0.4752 0.6151 0.4269 0.5177 0.6710
0.100 -5.24606756571109 1.151 1.48300 0.45759 0.37000 0.01825668401347 -1.74581 1.30112 -0.01768 -0.00652 -0.000478 -11.021900 0.095600 -0.093190 -0.014580 0.2080000000000000 0.6740 0.5670 0.8810 0.4286 0.6258 0.7585 0.4177 0.5662 0.7036 0.3936 0.5067 0.6416 0.4263 0.5593 0.7033
0.120 -5.24271284906608 1.151 1.45559 0.41355 0.40606 0.01825668401347 -1.73746 1.39137 -0.01742 -0.00660 -0.000489 -11.365300 0.200370 -0.008780 0.082520 0.2880000000000000 0.6900 0.5340 0.8720 0.4402 0.6135 0.7551 0.4285 0.5695 0.7127 0.4279 0.5479 0.6951 0.4445 0.5762 0.7277
0.140 -5.23987647655726 1.151 1.44277 0.37828 0.43450 0.01825668401347 -1.74463 1.47084 -0.01700 -0.00652 -0.000508 -11.730390 0.303720 0.089260 0.171510 0.3700000000000000 0.6920 0.5040 0.8560 0.4411 0.6126 0.7549 0.4347 0.5971 0.7386 0.4326 0.5072 0.6666 0.4424 0.5608 0.7143
0.150 -5.23860700772190 1.151 1.43314 0.36308 0.45000 0.01825668401347 -1.74972 1.50784 -0.01676 -0.00647 -0.000520 -11.880130 0.342840 0.136020 0.209320 0.4120000000000000 0.6960 0.4860 0.8500 0.4495 0.6023 0.7515 0.4399 0.5985 0.7428 0.4200 0.4896 0.6450 0.4401 0.5546 0.7080
0.160 -5.23741949893297 1.151 1.43253 0.34919 0.46055 0.01825668401347 -1.76259 1.54326 -0.01649 -0.00636 -0.000532 -12.056370 0.374000 0.177510 0.241170 0.4530000000000000 0.6970 0.4650 0.8380 0.4518 0.5987 0.7500 0.4443 0.5970 0.7442 0.4239 0.4939 0.6509 0.4369 0.5517 0.7037
0.180 -5.23525229107689 1.151 1.43710 0.32464 0.48439 0.01825668401347 -1.78989 1.60985 -0.01594 -0.00614 -0.000559 -12.420440 0.427000 0.253090 0.298960 0.5350000000000000 0.7040 0.4300 0.8250 0.4540 0.5990 0.7516 0.4474 0.6005 0.7489 0.4463 0.5003 0.6704 0.4361 0.5534 0.7046
0.200 -5.23331365758879 1.151 1.44781 0.30358 0.50900 0.01825668401347 -1.82110 1.67146 -0.01537 -0.00590 -0.000588 -12.785420 0.462970 0.320050 0.345910 0.6060000000000000 0.7130 0.4060 0.8210 0.4622 0.5897 0.7493 0.4499 0.5962 0.7469 0.4411 0.4907 0.6598 0.4322 0.5617 0.7088
0.250 -5.22920781624461 1.151 1.48260 0.26174 0.55500 0.01825668401347 -1.90412 1.80738 -0.01395 -0.00526 -0.000667 -13.635370 0.508560 0.453040 0.442310 0.6700000000000000 0.7110 0.3850 0.8080 0.4743 0.5545 0.7297 0.4702 0.6010 0.7631 0.4587 0.4500 0.6426 0.4319 0.5065 0.6657
0.300 -5.22585309959960 1.151 1.51881 0.23036 0.59300 0.01825668401347 -1.98439 1.92242 -0.01261 -0.00468 -0.000749 -14.380860 0.507760 0.548750 0.517820 0.7100000000000000 0.6830 0.3650 0.7750 0.4723 0.5320 0.7114 0.4749 0.5574 0.7322 0.4369 0.4989 0.6631 0.4323 0.4889 0.6526
0.350 -5.22301672709078 1.151 1.55291 0.20580 0.62500 0.01825668401347 -2.05756 2.02102 -0.01139 -0.00415 -0.000831 -15.035110 0.497100 0.617130 0.575960 0.7189367707551090 0.6650 0.3730 0.7620 0.4677 0.5050 0.6883 0.4780 0.5178 0.7047 0.4435 0.5315 0.6922 0.4318 0.4630 0.6331
0.400 -5.22055974946649 1.151 1.58443 0.18597 0.65200 0.01825668401347 -2.12282 2.10642 -0.01029 -0.00369 -0.000912 -15.615990 0.480650 0.666340 0.622390 0.7056101598849490 0.6570 0.3830 0.7610 0.4574 0.4792 0.6625 0.4836 0.4878 0.6868 0.4530 0.5458 0.7093 0.4155 0.4661 0.6244
0.450 -5.21839254161041 1.151 1.61360 0.16960 0.67500 0.01825668401347 -2.18047 2.18097 -0.00931 -0.00327 -0.000993 -16.138300 0.461590 0.701110 0.659760 0.6928932372909810 0.6470 0.3910 0.7560 0.4500 0.4608 0.6441 0.4772 0.4790 0.6761 0.4771 0.5294 0.7127 0.4079 0.4613 0.6158
0.500 -5.21645390812230 1.151 1.64075 0.15585 0.69500 0.01825668401347 -2.23118 2.24651 -0.00843 -0.00290 -0.001071 -16.613239 0.442240 0.725580 0.690660 0.6807534705053220 0.6400 0.4030 0.7560 0.4452 0.4506 0.6334 0.4696 0.4713 0.6653 0.4718 0.4947 0.6836 0.4087 0.4665 0.6202
0.600 -5.21309919147730 1.151 1.69020 0.13405 0.72900 0.01825668401347 -2.31475 2.35602 -0.00694 -0.00227 -0.001239 -17.452980 0.405370 0.752940 0.737960 0.6580415130288750 0.6330 0.4120 0.7550 0.4480 0.4355 0.6248 0.4597 0.4671 0.6553 0.4590 0.4701 0.6570 0.4065 0.4335 0.5943
0.700 -5.21026281896847 1.151 1.73450 0.11757 0.75600 0.01825668401347 -2.37885 2.44331 -0.00574 -0.00178 -0.001393 -18.180950 0.373420 0.762450 0.772260 0.6371531353146750 0.6330 0.4320 0.7660 0.4396 0.4299 0.6149 0.4595 0.4642 0.6532 0.4611 0.4727 0.6604 0.4031 0.4300 0.5894
0.800 -5.20780584134418 1.151 1.77474 0.10476 0.77800 0.01825668401347 -2.42769 2.51391 -0.00477 -0.00139 -0.001535 -18.824989 0.346230 0.761190 0.797360 0.6178103254041100 0.6360 0.4360 0.7710 0.4412 0.4269 0.6139 0.4596 0.4559 0.6474 0.4573 0.4565 0.6461 0.4067 0.4543 0.6097
0.900 -5.20563863348810 1.151 1.81162 0.09458 0.79600 0.01825668401347 -2.46450 2.57166 -0.00398 -0.00109 -0.001664 -19.403130 0.323640 0.753840 0.816160 0.5997867391697460 0.6360 0.4370 0.7720 0.4345 0.4314 0.6123 0.4559 0.4659 0.6519 0.4490 0.4391 0.6280 0.4082 0.4556 0.6117
1.000 -5.20370000000000 1.151 1.84561 0.08636 0.81200 0.01825668401347 -2.49170 2.61931 -0.00333 -0.00086 -0.001781 -19.927660 0.304790 0.742790 0.830090 0.5829000000000000 0.6370 0.4360 0.7720 0.4270 0.4400 0.6131 0.4479 0.4663 0.6466 0.4419 0.4415 0.6246 0.4080 0.4619 0.6163
1.250 -5.19959415865582 1.151 1.92015 0.07173 0.84100 0.01807627932941 -2.52758 2.70638 -0.00215 -0.00052 -0.001989 -21.058180 0.270260 0.708330 0.850360 0.5447531267038370 0.6350 0.4440 0.7750 0.4126 0.4401 0.6032 0.4422 0.4755 0.6493 0.4279 0.4321 0.6081 0.4123 0.4437 0.6057
1.500 -5.19623944201081 1.151 1.98274 0.06258 0.86100 0.01785908731221 -2.53359 2.76244 -0.00142 -0.00043 -0.002134 -21.996330 0.248310 0.672560 0.857320 0.5111822983011660 0.6450 0.4480 0.7850 0.4157 0.4485 0.6115 0.4477 0.4730 0.6513 0.4177 0.4605 0.6217 0.4178 0.4432 0.6090
2.000 -5.19094609187770 1.151 2.08214 0.05327 0.88400 0.01717728835933 -2.49565 2.82205 -0.00067 -0.00070 -0.002245 -23.488390 0.225290 0.610670 0.849910 0.4538170835899950 0.6330 0.4240 0.7620 0.4089 0.4409 0.6013 0.4371 0.4624 0.6362 0.4061 0.4720 0.6227 0.4133 0.4392 0.6031
2.500 -5.18684025053352 1.151 2.15841 0.05036 0.90000 0.01628389643418 -2.42623 2.84475 -0.00039 -0.00127 -0.002188 -24.647409 0.215400 0.564030 0.827570 0.4056165742693500 0.6080 0.4130 0.7350 0.3979 0.4247 0.5820 0.4293 0.4273 0.6057 0.3868 0.4841 0.6197 0.4139 0.4267 0.5945
3.000 -5.18348553388851 1.151 2.22046 0.04536 0.90000 0.01549254642056 -2.34726 2.84988 -0.00030 -0.00198 -0.002068 -25.597130 0.211540 0.526120 0.799110 0.3638313271186210 0.5820 0.4070 0.7110 0.3901 0.3963 0.5561 0.4227 0.4071 0.5869 0.3688 0.4470 0.5795 0.4117 0.4309 0.5959
3.500 -5.18064916137969 1.151 2.27406 0.04536 0.90000 0.01489171351523 -2.27002 2.84667 -0.00026 -0.00271 -0.001926 -26.409969 0.209756 0.497660 0.767820 0.3268167127622250 0.5620 0.3940 0.6870 0.3858 0.3856 0.5454 0.4104 0.4011 0.5738 0.3763 0.4339 0.5743 0.4016 0.4127 0.5758
4.000 -5.17819218375539 1.151 2.32307 0.04536 0.90000 0.01458020631717 -2.19947 2.83992 -0.00021 -0.00341 -0.001798 -27.131809 0.208748 0.476850 0.735940 0.2935047212753080 0.5400 0.3810 0.6600 0.3766 0.3728 0.5299 0.4104 0.3862 0.5636 0.3621 0.4064 0.5443 0.3943 0.3906 0.5550
4.500 -5.17602497589932 1.151 2.37009 0.04536 0.90000 0.01458710652043 -2.12528 2.82802 -0.00021 -0.00421 -0.001701 -27.792990 0.207737 0.462230 0.704070 0.2630000000000000 0.5250 0.3650 0.6400 0.3604 0.3635 0.5119 0.4116 0.3731 0.5555 0.3675 0.3831 0.5309 0.3853 0.3664 0.5317
5.000 -5.17408634241121 1.151 2.37009 0.04536 0.90000 0.01458710652043 -2.02646 2.82521 -0.00021 -0.00500 -0.001575 -28.313459 0.206721 0.452670 0.672200 0.2350000000000000 0.5220 0.3780 0.6450 0.3612 0.3588 0.5091 0.4469 0.3219 0.5508 0.3806 0.3163 0.4949 0.3799 0.3235 0.4990
""")
# Coefficients specific to the site amplification
SITE_COEFFS = CoeffsTable(sa_damping=5, table="""\
imt LnAmax1D1 LnAmax1D2 LnAmax1D3 LnAmax1D4 Src1D1 Src1D2 Src1D3 Src1D4 fsr1 fsr2 fsr3 fsr4
pga 0.650220 0.709730 0.644340 0.404280 8.429000 1.913680 1.117140 0.836440 1.0000 1.0000 1.0000 1.0000
0.005 0.650220 0.709730 0.644340 0.404280 8.429000 1.913680 1.117140 0.836440 1.0000 1.0000 1.0000 1.0000
0.010 0.651810 0.706790 0.646240 0.404280 8.090000 1.882560 1.114440 0.836440 1.0000 1.0000 1.0000 1.0000
0.020 0.653620 0.694650 0.638650 0.387890 6.992000 1.778610 1.124370 0.830000 1.0000 1.0000 1.0000 1.0500
0.030 0.654670 0.687550 0.634210 0.378300 6.350000 1.717810 1.130170 0.826240 1.0000 1.0000 1.0000 0.5800
0.040 0.652850 0.698920 0.606040 0.317370 4.883000 2.052340 1.150800 0.767580 1.0000 1.0060 1.0000 0.4820
0.050 0.672640 0.701370 0.617160 0.309340 5.043000 2.387130 1.239710 0.786320 1.0000 0.8510 1.0000 0.4720
0.060 0.699660 0.724450 0.637970 0.325300 6.271000 2.833990 1.348190 0.837750 1.0000 0.8030 1.0440 0.5060
0.070 0.717130 0.743430 0.654370 0.354120 7.667000 3.294470 1.451810 0.926160 1.0000 0.9180 0.9750 0.5870
0.080 0.716030 0.785980 0.680190 0.392820 9.034000 3.990910 1.583150 1.022280 1.0000 1.0620 0.9640 0.6830
0.090 0.725610 0.797210 0.708890 0.421840 11.251000 4.465760 1.732920 1.118020 1.0000 1.1060 0.9800 0.7820
0.100 0.742000 0.816680 0.718810 0.437360 14.817000 5.045610 1.841340 1.165780 1.0000 1.0710 0.9700 0.8230
0.120 0.762360 0.845230 0.725810 0.472080 14.817000 5.899600 2.030290 1.285510 0.0000 0.9515 1.0220 1.0290
0.140 0.752150 0.782960 0.745250 0.512780 14.817000 5.053530 2.281330 1.398080 0.0000 0.6720 0.8890 0.9910
0.150 0.738190 0.794800 0.761030 0.534320 14.817000 5.204900 2.444130 1.443270 0.0000 0.6310 0.8610 0.9830
0.160 0.719110 0.808610 0.768130 0.550220 14.817000 5.386940 2.580170 1.471770 0.0000 0.6000 0.8310 0.9730
0.180 0.654080 0.843310 0.756900 0.572790 14.817000 5.871650 2.741610 1.546940 0.0000 0.5710 0.7480 0.9790
0.200 0.583950 0.877700 0.717850 0.596740 14.817000 6.573910 2.825870 1.644010 0.0000 0.5650 0.6500 1.0060
0.250 0.583950 0.937670 0.654700 0.611360 14.817000 8.500000 2.718930 1.790130 0.0000 0.6010 0.4790 1.0270
0.300 0.583950 0.950000 0.696190 0.626380 14.817000 10.670300 2.417590 1.823450 0.0000 0.5790 0.4490 1.0210
0.350 0.583950 1.000000 0.779070 0.630120 14.817000 10.670300 2.303750 1.790370 0.0000 0.6790 0.4820 1.0030
0.400 0.583950 1.000000 0.827760 0.647730 14.817000 10.670300 2.236250 1.768440 0.0000 0.6550 0.4990 1.0100
0.450 0.583950 1.000000 0.876450 0.641520 14.817000 10.670300 2.216780 1.675390 0.0000 0.6150 0.5150 0.9850
0.500 0.583950 1.000000 0.925140 0.655820 14.817000 10.670300 2.243380 1.625390 0.0000 0.5500 0.5300 0.9900
0.600 0.583950 1.000000 0.973830 0.686680 14.817000 10.670300 2.805350 1.524530 0.0000 0.0000 0.5300 1.0060
0.700 0.583950 1.000000 1.022520 0.705600 14.817000 10.670300 6.658390 1.397240 0.0000 0.0000 0.4990 1.0000
0.800 0.583950 1.000000 1.071220 0.714290 14.817000 10.670300 30.000000 1.320290 0.0000 0.0000 0.3690 1.0000
0.900 0.583950 1.000000 1.119910 0.703880 14.817000 10.670300 30.000000 1.266370 0.0000 0.0000 0.3000 0.9600
1.000 0.583950 1.000000 1.168600 0.678130 14.817000 10.670300 30.000000 1.226800 0.0000 0.0000 0.2000 0.9040
1.250 0.583950 1.000000 1.217290 0.611190 14.817000 10.670300 30.000000 1.220650 0.0000 0.0000 0.0000 0.7380
1.500 0.583950 1.000000 1.265980 0.547360 14.817000 10.670300 30.000000 1.318050 0.0000 0.0000 0.0000 0.5350
2.000 0.583950 1.000000 1.314670 0.459440 14.817000 10.670300 30.000000 2.124850 0.0000 0.0000 0.0000 0.3580
2.500 0.583950 1.000000 1.363360 0.408460 14.817000 10.670300 30.000000 14.381810 0.0000 0.0000 0.0000 0.0000
3.000 0.583950 1.000000 1.412050 0.364210 14.817000 10.670300 30.000000 14.381810 0.0000 0.0000 0.0000 0.0000
3.500 0.583950 1.000000 1.460750 0.329840 14.817000 10.670300 30.000000 14.381810 0.0000 0.0000 0.0000 0.0000
4.000 0.583950 1.000000 1.509440 0.309120 14.817000 10.670300 30.000000 14.381810 0.0000 0.0000 0.0000 0.0000
4.500 0.583950 1.000000 1.558130 0.292510 14.817000 10.670300 30.000000 14.381810 0.0000 0.0000 0.0000 0.0000
5.000 0.583950 1.000000 1.606820 0.547360 14.817000 10.670300 30.000000 14.381810 0.0000 0.0000 0.0000 0.0000
""")
[docs]class ZhaoEtAl2016SSlabSiteSigma(ZhaoEtAl2016SSlab):
"""
Subclass of the Zhao et al. (2016c) subduction in-slab GMPE for the
case of site-dependent within-event variability
"""
[docs] def get_stddevs(self, C, n_sites, idx, stddev_types):
"""
Returns the intra-event standard deviation calibrated for the
specific site class
"""
stddevs = []
tau = C["tau"] + np.zeros(n_sites)
phi = np.zeros(n_sites)
for i in range(1, 5):
phi[idx[i]] += C["sc{:g}_sigma_S".format(i)]
for stddev_type in stddev_types:
assert stddev_type in self.DEFINED_FOR_STANDARD_DEVIATION_TYPES
if stddev_type == const.StdDev.TOTAL:
stddevs.append(np.sqrt(phi ** 2. + tau ** 2.))
elif stddev_type == const.StdDev.INTRA_EVENT:
stddevs.append(phi)
elif stddev_type == const.StdDev.INTER_EVENT:
stddevs.append(tau)
return stddevs