# -*- coding: utf-8 -*-
# vim: tabstop=4 shiftwidth=4 softtabstop=4
#
# Copyright (C) 2017-2019 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:`MunsonThurber1997`
:class:`MunsonThurber1997Hawaii`
:class:`MunsonThurber1997Vector`.
"""
import numpy as np
from openquake.hazardlib.gsim.base import GMPE
from openquake.hazardlib import const
from openquake.hazardlib.imt import PGA, SA
[docs]class MunsonThurber1997(GMPE):
"""
Implements GMPE developed by Clifford G. Munson and Clifford H. Thurber
and published as "Analysis of the Attenuation of Strong Ground Motion
on the Island of Hawaii" (1997, Bulletin of the Seismological Society
of America, Vol. 87, No. 4, pp. 954-960).
"""
#: Supported tectonic region type is volcanic,
#: see paragraph 'Introduction', page 99.
DEFINED_FOR_TECTONIC_REGION_TYPE = const.TRT.VOLCANIC
#: Supported intensity measure types is spectral acceleration,
#: see table 3, pag. 110
DEFINED_FOR_INTENSITY_MEASURE_TYPES = set([PGA])
#: Supported intensity measure component is maximum horizontal
#: :attr:`~openquake.hazardlib.const.IMC.VECTORIAL`,
DEFINED_FOR_INTENSITY_MEASURE_COMPONENT = const.IMC.HORIZONTAL
#: Supported standard deviation type is total
DEFINED_FOR_STANDARD_DEVIATION_TYPES = set([
const.StdDev.TOTAL
])
#: Required site parameters is Vs30.
#: See paragraph 'Predictor Variables', pag 103
REQUIRES_SITES_PARAMETERS = set(('vs30', ))
#: Required rupture parameter is magnitude
REQUIRES_RUPTURE_PARAMETERS = set(('mag', ))
#: Required distance measure is hypocentral distance
#: see page 18 in Atkinson and Boore's manuscript
REQUIRES_DISTANCES = set(('rjb', ))
[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.
"""
# Distance term
R = np.sqrt(dists.rjb ** 2 + 11.29 ** 2)
# Magnitude term
M = rup.mag - 6
# Site term only distinguishes between lava and ash;
# since ash sites have Vs30 in the range 60-200m/s,
# we use this upper value as class separator
S = np.zeros(R.shape)
S[sites.vs30 <= 200] = 1
# Mean ground motion (log10)
mean = (0.518 + 0.387*M - np.log10(R) - 0.00256*R + 0.335*S)
# Converting to natural log
mean /= np.log10(np.e)
# Check for standard deviation type
assert all(stddev_type in self.DEFINED_FOR_STANDARD_DEVIATION_TYPES
for stddev_type in stddev_types)
# Constant (total) standard deviation
stddevs = [0.237/np.log10(np.e) + np.zeros(R.shape)]
return mean, stddevs
[docs]class MunsonThurber1997Hawaii(MunsonThurber1997):
"""
Modifies :class:`MunsonThurber1997` for use with the USGS Hawaii seismic
hazard map of Klein FW, Frankel AD,Mueller CS, Wesson RL, Okubo PG.
Seismic-hazard maps for Hawaii. US Geological Survey; 2000.
"""
DEFINED_FOR_INTENSITY_MEASURE_TYPES = set([PGA, SA])
[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.
"""
# assign constant
log10e = np.log10(np.e)
# Distance term
R = np.sqrt(dists.rjb ** 2 + 11.29 ** 2)
# Magnitude term
M = rup.mag - 6
# Site term only distinguishes between lava and ash;
# since ash sites have Vs30 in the range 60-200m/s,
# we use this upper value as class separator
S = np.zeros(R.shape)
S[sites.vs30 <= 200] = 1
# Mean ground motion (natural log)
# call super
mean, stddevs = super().get_mean_and_stddevs(sites, rup, dists,
imt, stddev_types)
if rup.mag > 7. and rup.mag <= 7.7:
mean = (0.171 * (1 - M)) / log10e + mean
elif rup.mag > 7.7:
mean = (0.1512 + 0.387 * (1 - M)) / log10e + mean
# define natural log of SA 0.3 sec and 0.2 sec
if isinstance(imt, SA):
if imt.period == 0.3:
mean = np.log(2.2) + mean
if imt.period == 0.2:
mean = np.log(2.5) + mean
return mean, stddevs
[docs]class MunsonThurber1997Vector(MunsonThurber1997):
"""
Modification of the original base class to correct mean ground motion
to geometric mean of horizontal components (Beyer and Bommer, 2006)
"""
#: Supported intensity measure component is geometric mean of horizontal
#: :attr:`~openquake.hazardlib.const.IMC.VECTORIAL`,
DEFINED_FOR_INTENSITY_MEASURE_COMPONENT = const.IMC.VECTORIAL
[docs] def get_mean_and_stddevs(self, sites, rup, dists, imt, stddev_types):
mean, stddevs = super().get_mean_and_stddevs(sites, rup, dists,
imt, stddev_types)
# Conversion to geometric mean of horizontal components
# using the coefficient in Beyer and Bommer, 2006
mean += np.log(1.1)
return mean, stddevs