pygmm.macedo_abrahamson_liu_2021.MacedoAbrahamsonLiu2021

class pygmm.macedo_abrahamson_liu_2021.MacedoAbrahamsonLiu2021(scenario, pga_model=None, pga=None)[source]

Macedo, Abrahamson, and Liu (2021, [Macedo et al., 2021]) CAV model.

Conditional ground-motion model (CGMM) and scenario-based models for cumulative absolute velocity (CAV, m/s) for shallow crustal tectonic settings. Calibrated on the NGA-West2 database (~14,000 recordings, 287 events, M 3.0–7.9).

Two evaluation modes are supported:

Conditional mode (pga_model=None): CAV is estimated from a user-supplied PGA value (in g) using Eq. (12). The aleatory variability is \(\sigma = \sqrt{\tau^2 + \phi^2} \approx 0.31\).

Scenario mode (pga_model is one of 'ASK14', 'BSSA14', 'CB14', 'CY14', 'I14'): The median PGA and its lognormal standard deviation are taken from the chosen NGA-West2 GMM, then combined with the conditional model via propagation of errors (Eq. 17) to yield a full scenario-based CAV estimate.

Parameters:
  • scenario (pygmm.model.Scenario) – Earthquake scenario. When using scenario mode the scenario must also contain all parameters required by the chosen backbone PGA GMM (e.g. dist_x, dip, depth_tor for CY14).

  • pga_model (str or None, optional) – Backbone NGA-West2 GMM for scenario mode. Accepted values: 'ASK14', 'BSSA14', 'CB14', 'CY14', 'I14'. Default None selects conditional mode.

  • pga (float or None, optional) – Peak ground acceleration in g for conditional mode. Required when pga_model is None.

References

NAME = 'Macedo, Abrahamson, & Liu (2021)'

Long name of the model

ABBREV = 'MAL21'

Short name of the model

C1 = 1.79
C2 = 0.67
C3 = 0.57
C4 = -0.47
C5 = -0.0026
C6 = 0.17
TAU = 0.17
PHI = 0.26
SIGMA_COND = 0.31064449134018135
SUPPORTED_PGA_MODELS = ('ASK14', 'BSSA14', 'CB14', 'CY14', 'I14')
PARAMS = [<pygmm.model.NumericParameter object>, <pygmm.model.NumericParameter object>, <pygmm.model.NumericParameter object>, <pygmm.model.NumericParameter object>, <pygmm.model.NumericParameter object>, <pygmm.model.CategoricalParameter object>]

Model parameters

__init__(scenario, pga_model=None, pga=None)[source]

Initialize the model.

property ln_cav

Natural log of the median CAV (m/s).

property cav

Median CAV (m/s).

property ln_std

Total lognormal standard deviation of CAV.

property tau

Between-event standard deviation (conditional mode only).

property phi

Within-event standard deviation (conditional mode only).

property cav_plus_sigma

exp(ln_cav + ln_std).

Type:

84th-percentile CAV (m/s)

property cav_minus_sigma

exp(ln_cav - ln_std).

Type:

16th-percentile CAV (m/s)

LIMITS = {}

Limits of model applicability

property scenario