Source code for openquake.hazardlib.scalerel.thingbaijam2017

# -*- coding: utf-8 -*-
# vim: tabstop=4 shiftwidth=4 softtabstop=4
#
# Copyright (C) 2012-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 :mod:`openquake.hazardlib.scalerel.thingbaijam2017` implements
:class:`Thingbaijam2017_Interface`
:class:`Thingbaijam2017_Crustal`
"""

from numpy import log10
from openquake.hazardlib.scalerel.base import BaseMSRSigma, BaseASRSigma


[docs]class ThingbaijamInterface(BaseMSRSigma, BaseASRSigma): """ Thingbaijam, K. K. S., P. M. Mai, and K. Goda (2017). New Empirical Earthquake Source-Scaling Laws. Bulletin of the Seismological Society of America, 107(5), pp 2225-2946, doi: 10.1785/0120170017. Implements both magnitude-area and area-magnitude scaling relationships. """
[docs] def get_median_area(self, mag, rake): """ Calculates median fault area from magnitude. """ return 10**(-3.292+0.949*mag)
[docs] def get_std_dev_area(self, mag, rake): """ Returns std """ return 10**0.150
[docs] def get_median_mag(self, area, rake): """ Returns magnitude for a given fault area """ return (log10(area)+3.292)/0.949
[docs] def get_std_dev_mag(self, area, rake): """ Returns std """ return 0.150
[docs]class ThingbaijamStrikeSlip(BaseMSRSigma, BaseASRSigma): """ Thingbaijam, K. K. S., P. M. Mai, and K. Goda (2017). New Empirical Earthquake Source-Scaling Laws. Bulletin of the Seismological Society of America, 107(5), pp 2225-2946, doi: 10.1785/0120170017. Implements both magnitude-area and area-magnitude scaling relationships. """
[docs] def get_median_area(self, mag, rake): """ Calculates median fault area from magnitude. """ return 10**(-3.486+0.942*mag)
[docs] def get_std_dev_area(self, mag, rake): """ Returns std """ return 10**0.184
[docs] def get_median_mag(self, area, rake): """ Returns magnitude for a given fault area """ return (log10(area)+3.486)/0.942
[docs] def get_std_dev_mag(self, area, rake): """ Returns std """ return 0.184