Catalog.estimate_mc_maxc#

Catalog.estimate_mc_maxc(delta_m: float | None = None, correction_factor: float = 0.2) → float#

Returns the completeness magnitude (mc) estimate using the maximum curvature method.

Source:

  • Wiemer, S. and Wyss, M., 2000. Minimum magnitude of completeness in earthquake catalogs: Examples from Alaska, the western United States, and Japan. Bulletin of the Seismological Society of America, 90(4), pp.859-869.

  • Woessner, J. and Wiemer, S., 2005. Assessing the quality of earthquake catalogues: Estimating the magnitude of completeness and its uncertainty. Bulletin of the Seismological Society of America, 95(2), pp.684-698.

Attention

Catalog.mc will be replaced by the mc return value.

Parameters:

  • magnitudes – Array of magnitudes to test.

  • delta_m – Bin size of discretized magnitudes. Catalog needs to be rounded to bins beforehand. Either given as parameter or taken from the object attribute.

  • correction_factor – Correction factor for the maximum curvature method (default value after Woessner & Wiemer 2005).

Returns:

mc – Estimated completeness magnitude.

Examples

>>> from seismostats import Catalog
>>> simple_catalog = Catalog.from_dict({
...     'magnitude': [2.3, 1.2, 1.5, 1.2, 1.7, 1.1, 1.2, 1.5,
...                   1.8, 1.6, 1.2, 1.5, 1.2, 1.7, 1.6, 1.1,
...                   1.1, 1.2, 2.0, 1.1, 1.2, 1.1, 1.2, 1.6,
...                   1.9, 1.3, 1.7, 1.3, 1.0, 1.2, 1.7, 1.3,
...                   1.3, 1.1, 1.5, 1.4]})
>>> simple_catalog.delta_m = 0.1
>>> simple_catalog.estimate_mc_maxc()
>>> simple_catalog.mc

1.4