pygmm.atkinson_boore_2006.AtkinsonBoore2006¶

class pygmm.atkinson_boore_2006.AtkinsonBoore2006(scenario)[source]¶

Atkinson and Boore (2006, [Atkinson and Boore, 2006]) model.

Developed for the Eastern North America with a reference velocity of 760 or 2000 m/s.

Parameters:: scenario (pygmm.model.Scenario) – earthquake scenario

NAME = 'Atkinson and Boore (2006)'¶: Long name of the model

ABBREV = 'AB06'¶: Short name of the model

COEFF = {'bc': rec.array([(-2. , -3.09 , 0.828 , -0.008541, -2.381, 0.1782, -1.643 , 0.1628 , 0.581 , -0.145 , -0.0002337), (-1. , -1.662 , 1.05 , -0.06035 , -2.496, 0.184 , -2.301 , 0.25 , 0.1268 , -0.08704, -0.0004266), ( 0. , 0.5233, 0.9686, -0.06196 , -2.439, 0.1465, -2.335 , 0.1912 , -0.08695, -0.08285, -0.0006304), ( 0.025, 1.052 , 0.903 , -0.05768 , -2.571, 0.1483, -2.652 , 0.2065 , -0.4084 , -0.05769, -0.0005122), ( 0.031, 1.191 , 0.8884, -0.05642 , -2.577, 0.1451, -2.84 , 0.2121 , -0.437 , -0.05866, -0.0004329), ( 0.04 , 1.261 , 0.8789, -0.05515 , -2.536, 0.1388, -2.994 , 0.2158 , -0.3908 , -0.06746, -0.0003881), ( 0.05 , 1.209 , 0.883 , -0.05441 , -2.44 , 0.1295, -3.035 , 0.2133 , -0.2098 , -0.08997, -0.0004145), ( 0.063, 1.109 , 0.8875, -0.05386 , -2.334, 0.1229, -2.881 , 0.2007 , -0.03189, -0.1069 , -0.0005483), ( 0.079, 0.9667, 0.9033, -0.05476 , -2.249, 0.1215, -2.53 , 0.1775 , 0.1001 , -0.1147 , -0.0007724), ( 0.1 , 0.7818, 0.9235, -0.05555 , -2.165, 0.1191, -2.097 , 0.1483 , 0.2847 , -0.1319 , -0.0009897), ( 0.125, 0.5356, 0.9647, -0.05835 , -2.11 , 0.1205, -1.672 , 0.1156 , 0.3433 , -0.1322 , -0.00113 ), ( 0.158, 0.1194, 1.057 , -0.06473 , -2.054, 0.119 , -1.355 , 0.0916 , 0.5164 , -0.1503 , -0.001178 ), ( 0.199, -0.3056, 1.156 , -0.07211 , -2.038, 0.122 , -1.147 , 0.07375, 0.5082 , -0.143 , -0.00114 ), ( 0.251, -0.8756, 1.293 , -0.08193 , -2.014, 0.1226, -1.027 , 0.06341, 0.5808 , -0.1491 , -0.001053 ), ( 0.315, -1.56 , 1.455 , -0.09312 , -1.977, 0.1209, -0.9466, 0.05576, 0.6499 , -0.1558 , -0.0009552), ( 0.397, -2.281 , 1.629 , -0.1054 , -1.967, 0.1227, -0.888 , 0.05033, 0.6839 , -0.1582 , -0.0008587), ( 0.5 , -3.007 , 1.803 , -0.1178 , -1.982, 0.1274, -0.8466, 0.04698, 0.667 , -0.1546 , -0.0007676), ( 0.629, -3.748 , 1.973 , -0.1294 , -1.997, 0.1313, -0.8417, 0.0482 , 0.6772 , -0.1557 , -0.0006763), ( 0.794, -4.446 , 2.119 , -0.1387 , -2.009, 0.1356, -0.8576, 0.04976, 0.7084 , -0.1589 , -0.0005751), ( 1. , -5.058 , 2.233 , -0.1454 , -2.03 , 0.1408, -0.8744, 0.05412, 0.7922 , -0.1697 , -0.0004886), ( 1.25 , -5.489 , 2.289 , -0.1476 , -2.081, 0.1501, -0.9 , 0.05794, 0.8208 , -0.1719 , -0.000407 ), ( 1.587, -5.754 , 2.287 , -0.145 , -2.131, 0.1582, -0.9568, 0.06762, 0.867 , -0.1789 , -0.0003429), ( 2. , -5.853 , 2.233 , -0.1385 , -2.195, 0.1688, -1.037 , 0.08002, 0.8666 , -0.179 , -0.000286 ), ( 2.5 , -5.8 , 2.126 , -0.1278 , -2.257, 0.179 , -1.123 , 0.09539, 0.8911 , -0.1797 , -0.0002601), ( 3.125, -5.59 , 1.972 , -0.1136 , -2.331, 0.1908, -1.204 , 0.1099 , 0.8449 , -0.1723 , -0.0002452), ( 4. , -5.256 , 1.787 , -0.09785 , -2.435, 0.2068, -1.307 , 0.121 , 0.734 , -0.156 , -0.0001959), ( 5. , -4.852 , 1.58 , -0.08066 , -2.53 , 0.2216, -1.426 , 0.1361 , 0.634 , -0.1413 , -0.0001608)], dtype=[('period', '<f8'), ('c_1', '<f8'), ('c_2', '<f8'), ('c_3', '<f8'), ('c_4', '<f8'), ('c_5', '<f8'), ('c_6', '<f8'), ('c_7', '<f8'), ('c_8', '<f8'), ('c_9', '<f8'), ('c_10', '<f8')]), 'rock': rec.array([(-2. , -2.82 , 0.7491, -0.004196, -2.611, 0.2157, -1.671 , 0.1643 , 0.7452 , -0.1547 , -1.917e-04), (-1. , -1.442 , 0.9909, -0.05848 , -2.701, 0.2155, -2.436 , 0.2659 , 0.08479, -0.06927, -3.734e-04), ( 0. , 0.9069, 0.983 , -0.06595 , -2.698, 0.1594, -2.795 , 0.212 , -0.3011 , -0.06532, -4.484e-04), ( 0.025, 1.522 , 0.9597, -0.06351 , -2.813, 0.1458, -3.654 , 0.2362 , -0.6544 , -0.055 , -4.848e-05), ( 0.031, 1.436 , 0.9592, -0.06276 , -2.714, 0.14 , -3.728 , 0.2343 , -0.543 , -0.06448, -3.230e-05), ( 0.04 , 1.264 , 0.968 , -0.06232 , -2.581, 0.1317, -3.644 , 0.2276 , -0.3506 , -0.08126, -1.225e-04), ( 0.05 , 1.105 , 0.9719, -0.06197 , -2.466, 0.1276, -3.39 , 0.2144 , -0.1391 , -0.09839, -3.167e-04), ( 0.063, 0.9109, 0.9802, -0.06208 , -2.36 , 0.1263, -2.972 , 0.191 , 0.1069 , -0.1173 , -5.786e-04), ( 0.079, 0.6906, 0.9974, -0.06276 , -2.262, 0.1246, -2.487 , 0.1636 , 0.2139 , -0.1207 , -8.469e-04), ( 0.1 , 0.4797, 1.017 , -0.06404 , -2.201, 0.127 , -2.007 , 0.1326 , 0.3371 , -0.1266 , -1.047e-03), ( 0.125, 0.2144, 1.054 , -0.06664 , -2.154, 0.1295, -1.608 , 0.1046 , 0.4273 , -0.1303 , -1.153e-03), ( 0.158, -0.1455, 1.123 , -0.07143 , -2.116, 0.1302, -1.303 , 0.08311, 0.5617 , -0.1438 , -1.182e-03), ( 0.199, -0.6153, 1.227 , -0.07886 , -2.087, 0.1312, -1.12 , 0.06788, 0.6055 , -0.1459 , -1.125e-03), ( 0.251, -1.121 , 1.342 , -0.08722 , -2.082, 0.1349, -0.9714, 0.05628, 0.614 , -0.1432 , -1.055e-03), ( 0.315, -1.721 , 1.483 , -0.09739 , -2.08 , 0.1382, -0.8893, 0.04869, 0.6101 , -0.1389 , -9.538e-04), ( 0.397, -2.437 , 1.649 , -0.1084 , -2.051, 0.1363, -0.8426, 0.04483, 0.7386 , -0.1557 , -8.509e-04), ( 0.5 , -3.216 , 1.826 , -0.1201 , -2.018, 0.1344, -0.8134, 0.04437, 0.8839 , -0.1751 , -7.704e-04), ( 0.629, -3.917 , 1.987 , -0.1314 , -2.045, 0.1419, -0.7818, 0.04297, 0.7878 , -0.159 , -6.948e-04), ( 0.794, -4.604 , 2.132 , -0.1406 , -2.062, 0.1468, -0.7974, 0.04345, 0.7748 , -0.1558 , -5.790e-04), ( 1. , -5.272 , 2.264 , -0.1483 , -2.069, 0.1497, -0.8132, 0.04666, 0.8262 , -0.1622 , -4.862e-04), ( 1.25 , -5.724 , 2.324 , -0.1505 , -2.104, 0.1565, -0.8202, 0.05186, 0.8563 , -0.1661 , -4.329e-04), ( 1.587, -6.043 , 2.342 , -0.1496 , -2.157, 0.1662, -0.8704, 0.06047, 0.9207 , -0.1734 , -3.748e-04), ( 2. , -6.183 , 2.302 , -0.1442 , -2.223, 0.177 , -0.937 , 0.07067, 0.9518 , -0.1768 , -3.220e-04), ( 2.5 , -6.169 , 2.211 , -0.1348 , -2.299, 0.1898, -0.986 , 0.0786 , 0.9683 , -0.1765 , -2.823e-04), ( 3.125, -6.038 , 2.08 , -0.1221 , -2.367, 0.2002, -1.073 , 0.0895 , 1.002 , -0.1803 , -2.306e-04), ( 4. , -5.791 , 1.916 , -0.1071 , -2.441, 0.2113, -1.162 , 0.1018 , 1.012 , -0.1824 , -2.010e-04), ( 5. , -5.408 , 1.714 , -0.09012 , -2.537, 0.2267, -1.268 , 0.1162 , 0.9792 , -0.1767 , -1.757e-04)], dtype=[('period', '<f8'), ('c_1', '<f8'), ('c_2', '<f8'), ('c_3', '<f8'), ('c_4', '<f8'), ('c_5', '<f8'), ('c_6', '<f8'), ('c_7', '<f8'), ('c_8', '<f8'), ('c_9', '<f8'), ('c_10', '<f8')])}¶

PERIODS = array([-2. , -1. , 0. , 0.025, 0.031, 0.04 , 0.05 , 0.063, 0.079, 0.1 , 0.125, 0.158, 0.199, 0.251, 0.315, 0.397, 0.5 , 0.629, 0.794, 1. , 1.25 , 1.587, 2. , 2.5 , 3.125, 4. , 5. ])¶: Indices of the periods

COEFF_SITE = rec.array([(-1. , -0.6 , -0.495, -0.06 ), ( 0. , -0.361, -0.641, -0.144), ( 0.025 , -0.33 , -0.624, -0.115), ( 0.0313, -0.322, -0.618, -0.108), ( 0.04 , -0.314, -0.609, -0.105), ( 0.05 , -0.286, -0.643, -0.105), ( 0.0629, -0.249, -0.642, -0.105), ( 0.0794, -0.232, -0.637, -0.117), ( 0.1 , -0.25 , -0.595, -0.132), ( 0.125 , -0.26 , -0.56 , -0.14 ), ( 0.1587, -0.28 , -0.528, -0.185), ( 0.2 , -0.306, -0.521, -0.185), ( 0.25 , -0.39 , -0.518, -0.16 ), ( 0.3125, -0.445, -0.513, -0.13 ), ( 0.4 , -0.5 , -0.508, -0.095), ( 0.5 , -0.6 , -0.495, -0.06 ), ( 0.625 , -0.67 , -0.48 , -0.031), ( 0.7692, -0.69 , -0.465, -0.002), ( 1. , -0.7 , -0.44 , 0. ), ( 1.5873, -0.726, -0.395, 0. ), ( 2. , -0.73 , -0.375, 0. ), ( 3.125 , -0.74 , -0.33 , 0. ), ( 4. , -0.745, -0.31 , 0. ), ( 5. , -0.752, -0.3 , 0. )], dtype=[('period', '<f8'), ('b_lin', '<f8'), ('b_1', '<f8'), ('b_2', '<f8')])¶

COEFF_SF = rec.array([(-1. , 0.11, 2. , 5.5 ), ( 0.0101, 0.15, 0.5 , 5.5 ), ( 0.025 , 0.15, 0. , 5. ), ( 0.0315, 0.15, 0. , 5. ), ( 0.0396, 0.15, 0. , 5. ), ( 0.0499, 0.15, 0. , 5. ), ( 0.0629, 0.15, 0.17, 5.17), ( 0.0791, 0.15, 0.34, 5.34), ( 0.0996, 0.15, 0.5 , 5.5 ), ( 0.1255, 0.15, 1.15, 5.67), ( 0.158 , 0.15, 1.85, 5.84), ( 0.1988, 0.15, 2.5 , 6. ), ( 0.2506, 0.15, 2.9 , 6.12), ( 0.3155, 0.15, 3.3 , 6.25), ( 0.3968, 0.15, 3.65, 6.37), ( 0.5 , 0.15, 4. , 6.5 ), ( 0.6289, 0.15, 4.17, 6.7 ), ( 0.7937, 0.15, 4.34, 6.95), ( 1. , 0.15, 4.5 , 7.2 ), ( 1.25 , 0.15, 4.67, 7.45), ( 1.5873, 0.15, 4.84, 7.7 ), ( 2. , 0.15, 5. , 8. ), ( 2.5 , 0.15, 5.25, 8.12), ( 3.125 , 0.15, 5.5 , 8.25), ( 4. , 0.15, 5.75, 8.37), ( 5. , 0.15, 6. , 8.5 )], dtype=[('period', '<f8'), ('delta', '<f8'), ('m_1', '<f8'), ('m_h', '<f8')])¶

INDEX_PGD = 0¶: Index of the peak ground displacement

INDEX_PGV = 1¶: Index of the peak ground velocity

INDEX_PGA = 2¶: Index of the peak ground acceleration

INDICES_PSA = array([ 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26])¶: Indices for the spectral accelerations

PARAMS = [<pygmm.model.NumericParameter object>, <pygmm.model.NumericParameter object>, <pygmm.model.NumericParameter object>]¶: Model parameters

__init__(scenario)[source]¶

Initialize the model.

Parameters:: scenario (Scenario)

LIMITS = {}¶: Limits of model applicability

PGD_SCALE = 1.0¶: Scale factor to apply to get PGD in cm

PGV_SCALE = 1.0¶: Scale factor to apply to get PGV in cm/sec

interp_ln_spec_accels(periods, kind='linear')¶

Interpolate the spectral acceleration.

Interpolation of the spectral acceleration is done in natural log space.

Parameters:

periods (array_like) – spectral periods to interpolate the response.
kind (str, optional) – see scipy.interpolate.interp1d() for description of kind. Options include: ‘linear’ (default), ‘nearest’, ‘zero’, ‘slinear’, ‘quadratic’, and ‘cubic’

Returns:

ln_spec_accels – interpolated spectral accelerations

Return type:

np.ndarray

interp_ln_stds(periods, kind='linear')¶

Interpolate the logarithmic standard deviation.

Interpolate the logarithmic standard deviation (\(\sigma_{\ln}\)) of spectral acceleration at the provided damping at specified periods.

Parameters:

periods (array_like) – spectral periods to interpolate the response.
kind (str, optional) – see scipy.interpolate.interp1d() for description of kind. Options include: ‘linear’ (default), ‘nearest’, ‘zero’, ‘slinear’, ‘quadratic’, and ‘cubic’

Returns:

ln_stds – interpolated logarithmic standard deviations

Return type:

np.ndarray

interp_spec_accels(periods, kind='linear')¶

Interpolate the spectral acceleration.

Interpolation of the spectral acceleration is done in natural log space.

Parameters:

periods (array_like) – spectral periods to interpolate the response.
kind (str, optional) – see scipy.interpolate.interp1d() for description of kind. Options include: ‘linear’ (default), ‘nearest’, ‘zero’, ‘slinear’, ‘quadratic’, and ‘cubic’

Returns:

spec_accels – interpolated spectral accelerations

Return type:

np.ndarray

property ln_std_pga: float¶: Peak ground accelaration log-standard deviation.

property ln_std_pgd: float¶: Peak ground displacement log-standard deviation.

property ln_std_pgv: float¶: Peak ground velocity log-standard deviation.

property ln_stds: ndarray¶: Pseudo-spectral accelerations log-standard deviation.

property periods: ndarray¶: Periods specified by the model.

property pga: float¶: Peak ground acceleration (PGA) computed by the model (g).

property pgd: float¶: Peak ground displacement (PGD) computed by the model (cm).

property pgv: float¶: Peak ground velocity (PGV) computed by the model (cm/sec).

property scenario¶

property spec_accels: ndarray¶: Pseudo-spectral accelerations computed by the model (g).