FAILURE MAP
← Case archive

FA-71341 / Seismic magnitude estimation / Open access

Hutton-Boore local magnitude: zero-to-peak halving · case 01

Every station magnitude is 0.30 units too large.

Verified by executionVariant 1 · 8 checks per implementationDownload source bundle ↓JSON ↗

ROOT CAUSE

The peak-to-peak amplitude is used where the scale is calibrated on zero-to-peak.

VERIFIED REPAIR

Halve the peak-to-peak amplitude.

Unsuccessful approach: Dividing by sqrt(2) treats the reading as an RMS amplitude and leaves a 0.15 bias.

Case contract

Input the Wood-Anderson peak-to-peak amplitude in nm, epicentral distance and depth in km. Non-positive amplitude returns None. A (mm, zero-to-peak) = p2p/2 * 1e-6; r = hypocentral distance; ML = log10(A) + 1.11 log10(r/100) + 0.00189 (r - 100) + 3.0, rounded to 0.01.

Why this case matters

Regional networks publish ML for most small and moderate earthquakes.

1 / The failure

Exit 1
"""Failure Map reference implementation. Python standard library only."""
import json
import math
N = 1
observations = []
def solve(amp_p2p_nm, epi_km, depth_km):
    if amp_p2p_nm <= 0:
        return None
    a_mm = amp_p2p_nm * 1e-6
    r = math.hypot(epi_km, depth_km)
    ml = math.log10(a_mm) + 1.11 * math.log10(r / 100.0) + 0.00189 * (r - 100.0) + 3.0
    return round(ml, 2)
def check(label, actual, expected):
    observations.append({"check": label, "actual": actual, "expected": expected, "passed": actual == expected})
fixtures = [[['0 nm p2p at 0 km (depth 8)', [0, 0, 8], None], ['0 nm p2p at 12 km (depth 5)', [0, 12, 5], None], ['0 nm p2p at 40 km (depth 10)', [0, 40, 10], None], ['0 nm p2p at 100 km (depth 0)', [0, 100, 0], None], ['0 nm p2p at 100 km (depth 15)', [0, 100, 15], None], ['0 nm p2p at 250 km (depth 12)', [0, 250, 12], None], ['150 nm p2p at 0 km (depth 8)', [150, 0, 8], -2.52], ['150 nm p2p at 12 km (depth 5)', [150, 12, 5], -2.27]], [['150 nm p2p at 0 km (depth 8)', [150, 0, 8], -2.52], ['150 nm p2p at 12 km (depth 5)', [150, 12, 5], -2.27], ['150 nm p2p at 40 km (depth 10)', [150, 40, 10], -1.66], ['150 nm p2p at 100 km (depth 0)', [150, 100, 0], -1.12], ['150 nm p2p at 100 km (depth 15)', [150, 100, 15], -1.12], ['150 nm p2p at 250 km (depth 12)', [150, 250, 12], -0.4], ['150 nm p2p at 480 km (depth 30)', [150, 480, 30], 0.35], ['2000 nm p2p at 0 km (depth 8)', [2000, 0, 8], -1.39]], [['150 nm p2p at 480 km (depth 30)', [150, 480, 30], 0.35], ['2000 nm p2p at 0 km (depth 8)', [2000, 0, 8], -1.39], ['2000 nm p2p at 12 km (depth 5)', [2000, 12, 5], -1.15], ['2000 nm p2p at 40 km (depth 10)', [2000, 40, 10], -0.54], ['2000 nm p2p at 100 km (depth 0)', [2000, 100, 0], 0.0], ['2000 nm p2p at 100 km (depth 15)', [2000, 100, 15], 0.01], ['2000 nm p2p at 250 km (depth 12)', [2000, 250, 12], 0.73], ['2000 nm p2p at 480 km (depth 30)', [2000, 480, 30], 1.48]], [['2000 nm p2p at 40 km (depth 10)', [2000, 40, 10], -0.54], ['48000 nm p2p at 0 km (depth 8)', [48000, 0, 8], -0.01], ['48000 nm p2p at 12 km (depth 5)', [48000, 12, 5], 0.23], ['48000 nm p2p at 40 km (depth 10)', [48000, 40, 10], 0.84], ['48000 nm p2p at 100 km (depth 0)', [48000, 100, 0], 1.38], ['48000 nm p2p at 100 km (depth 15)', [48000, 100, 15], 1.39], ['48000 nm p2p at 250 km (depth 12)', [48000, 250, 12], 2.11], ['48000 nm p2p at 480 km (depth 30)', [48000, 480, 30], 2.86]], [['2000 nm p2p at 250 km (depth 12)', [2000, 250, 12], 0.73], ['900000 nm p2p at 0 km (depth 8)', [900000, 0, 8], 1.26], ['900000 nm p2p at 12 km (depth 5)', [900000, 12, 5], 1.51], ['900000 nm p2p at 40 km (depth 10)', [900000, 40, 10], 2.12], ['900000 nm p2p at 100 km (depth 0)', [900000, 100, 0], 2.65], ['900000 nm p2p at 100 km (depth 15)', [900000, 100, 15], 2.66], ['900000 nm p2p at 250 km (depth 12)', [900000, 250, 12], 3.38], ['900000 nm p2p at 480 km (depth 30)', [900000, 480, 30], 4.13]]]
for label, args, expected in fixtures[N - 1]:
    check(label, solve(*args), expected)
print(json.dumps({"observations": observations, "passed": all(x["passed"] for x in observations)}, ensure_ascii=False))
raise SystemExit(0 if all(x["passed"] for x in observations) else 1)
Boundary fixtureActualExpectedOutcome
0 nm p2p at 0 km (depth 8)NoneNonePassed
0 nm p2p at 12 km (depth 5)NoneNonePassed
0 nm p2p at 40 km (depth 10)NoneNonePassed
0 nm p2p at 100 km (depth 0)NoneNonePassed
0 nm p2p at 100 km (depth 15)NoneNonePassed
0 nm p2p at 250 km (depth 12)NoneNonePassed
150 nm p2p at 0 km (depth 8)-2.22-2.52Failed
150 nm p2p at 12 km (depth 5)-1.97-2.27Failed

SHA-256 / a60a1e21ff7a0299c28a4df06884e0c39725b3f1ffde65db0bf002de9e99e7ab

2 / The unsuccessful fix

Exit 1
"""Failure Map reference implementation. Python standard library only."""
import json
import math
N = 1
observations = []
def solve(amp_p2p_nm, epi_km, depth_km):
    if amp_p2p_nm <= 0:
        return None
    a_mm = amp_p2p_nm / 2.0 ** 0.5 * 1e-6
    r = math.hypot(epi_km, depth_km)
    ml = math.log10(a_mm) + 1.11 * math.log10(r / 100.0) + 0.00189 * (r - 100.0) + 3.0
    return round(ml, 2)
def check(label, actual, expected):
    observations.append({"check": label, "actual": actual, "expected": expected, "passed": actual == expected})
fixtures = [[['0 nm p2p at 0 km (depth 8)', [0, 0, 8], None], ['0 nm p2p at 12 km (depth 5)', [0, 12, 5], None], ['0 nm p2p at 40 km (depth 10)', [0, 40, 10], None], ['0 nm p2p at 100 km (depth 0)', [0, 100, 0], None], ['0 nm p2p at 100 km (depth 15)', [0, 100, 15], None], ['0 nm p2p at 250 km (depth 12)', [0, 250, 12], None], ['150 nm p2p at 0 km (depth 8)', [150, 0, 8], -2.52], ['150 nm p2p at 12 km (depth 5)', [150, 12, 5], -2.27]], [['150 nm p2p at 0 km (depth 8)', [150, 0, 8], -2.52], ['150 nm p2p at 12 km (depth 5)', [150, 12, 5], -2.27], ['150 nm p2p at 40 km (depth 10)', [150, 40, 10], -1.66], ['150 nm p2p at 100 km (depth 0)', [150, 100, 0], -1.12], ['150 nm p2p at 100 km (depth 15)', [150, 100, 15], -1.12], ['150 nm p2p at 250 km (depth 12)', [150, 250, 12], -0.4], ['150 nm p2p at 480 km (depth 30)', [150, 480, 30], 0.35], ['2000 nm p2p at 0 km (depth 8)', [2000, 0, 8], -1.39]], [['150 nm p2p at 480 km (depth 30)', [150, 480, 30], 0.35], ['2000 nm p2p at 0 km (depth 8)', [2000, 0, 8], -1.39], ['2000 nm p2p at 12 km (depth 5)', [2000, 12, 5], -1.15], ['2000 nm p2p at 40 km (depth 10)', [2000, 40, 10], -0.54], ['2000 nm p2p at 100 km (depth 0)', [2000, 100, 0], 0.0], ['2000 nm p2p at 100 km (depth 15)', [2000, 100, 15], 0.01], ['2000 nm p2p at 250 km (depth 12)', [2000, 250, 12], 0.73], ['2000 nm p2p at 480 km (depth 30)', [2000, 480, 30], 1.48]], [['2000 nm p2p at 40 km (depth 10)', [2000, 40, 10], -0.54], ['48000 nm p2p at 0 km (depth 8)', [48000, 0, 8], -0.01], ['48000 nm p2p at 12 km (depth 5)', [48000, 12, 5], 0.23], ['48000 nm p2p at 40 km (depth 10)', [48000, 40, 10], 0.84], ['48000 nm p2p at 100 km (depth 0)', [48000, 100, 0], 1.38], ['48000 nm p2p at 100 km (depth 15)', [48000, 100, 15], 1.39], ['48000 nm p2p at 250 km (depth 12)', [48000, 250, 12], 2.11], ['48000 nm p2p at 480 km (depth 30)', [48000, 480, 30], 2.86]], [['2000 nm p2p at 250 km (depth 12)', [2000, 250, 12], 0.73], ['900000 nm p2p at 0 km (depth 8)', [900000, 0, 8], 1.26], ['900000 nm p2p at 12 km (depth 5)', [900000, 12, 5], 1.51], ['900000 nm p2p at 40 km (depth 10)', [900000, 40, 10], 2.12], ['900000 nm p2p at 100 km (depth 0)', [900000, 100, 0], 2.65], ['900000 nm p2p at 100 km (depth 15)', [900000, 100, 15], 2.66], ['900000 nm p2p at 250 km (depth 12)', [900000, 250, 12], 3.38], ['900000 nm p2p at 480 km (depth 30)', [900000, 480, 30], 4.13]]]
for label, args, expected in fixtures[N - 1]:
    check(label, solve(*args), expected)
print(json.dumps({"observations": observations, "passed": all(x["passed"] for x in observations)}, ensure_ascii=False))
raise SystemExit(0 if all(x["passed"] for x in observations) else 1)
Boundary fixtureActualExpectedOutcome
0 nm p2p at 0 km (depth 8)NoneNonePassed
0 nm p2p at 12 km (depth 5)NoneNonePassed
0 nm p2p at 40 km (depth 10)NoneNonePassed
0 nm p2p at 100 km (depth 0)NoneNonePassed
0 nm p2p at 100 km (depth 15)NoneNonePassed
0 nm p2p at 250 km (depth 12)NoneNonePassed
150 nm p2p at 0 km (depth 8)-2.37-2.52Failed
150 nm p2p at 12 km (depth 5)-2.12-2.27Failed

SHA-256 / 0253128f4ec7e6e1b9353c884530413b0965c6160bce737fdd8298d782940e8e

3 / The verified repair

Exit 0
"""Failure Map reference implementation. Python standard library only."""
import json
import math
N = 1
observations = []
def solve(amp_p2p_nm, epi_km, depth_km):
    if amp_p2p_nm <= 0:
        return None
    a_mm = amp_p2p_nm / 2.0 * 1e-6
    r = math.hypot(epi_km, depth_km)
    ml = math.log10(a_mm) + 1.11 * math.log10(r / 100.0) + 0.00189 * (r - 100.0) + 3.0
    return round(ml, 2)
def check(label, actual, expected):
    observations.append({"check": label, "actual": actual, "expected": expected, "passed": actual == expected})
fixtures = [[['0 nm p2p at 0 km (depth 8)', [0, 0, 8], None], ['0 nm p2p at 12 km (depth 5)', [0, 12, 5], None], ['0 nm p2p at 40 km (depth 10)', [0, 40, 10], None], ['0 nm p2p at 100 km (depth 0)', [0, 100, 0], None], ['0 nm p2p at 100 km (depth 15)', [0, 100, 15], None], ['0 nm p2p at 250 km (depth 12)', [0, 250, 12], None], ['150 nm p2p at 0 km (depth 8)', [150, 0, 8], -2.52], ['150 nm p2p at 12 km (depth 5)', [150, 12, 5], -2.27]], [['150 nm p2p at 0 km (depth 8)', [150, 0, 8], -2.52], ['150 nm p2p at 12 km (depth 5)', [150, 12, 5], -2.27], ['150 nm p2p at 40 km (depth 10)', [150, 40, 10], -1.66], ['150 nm p2p at 100 km (depth 0)', [150, 100, 0], -1.12], ['150 nm p2p at 100 km (depth 15)', [150, 100, 15], -1.12], ['150 nm p2p at 250 km (depth 12)', [150, 250, 12], -0.4], ['150 nm p2p at 480 km (depth 30)', [150, 480, 30], 0.35], ['2000 nm p2p at 0 km (depth 8)', [2000, 0, 8], -1.39]], [['150 nm p2p at 480 km (depth 30)', [150, 480, 30], 0.35], ['2000 nm p2p at 0 km (depth 8)', [2000, 0, 8], -1.39], ['2000 nm p2p at 12 km (depth 5)', [2000, 12, 5], -1.15], ['2000 nm p2p at 40 km (depth 10)', [2000, 40, 10], -0.54], ['2000 nm p2p at 100 km (depth 0)', [2000, 100, 0], 0.0], ['2000 nm p2p at 100 km (depth 15)', [2000, 100, 15], 0.01], ['2000 nm p2p at 250 km (depth 12)', [2000, 250, 12], 0.73], ['2000 nm p2p at 480 km (depth 30)', [2000, 480, 30], 1.48]], [['2000 nm p2p at 40 km (depth 10)', [2000, 40, 10], -0.54], ['48000 nm p2p at 0 km (depth 8)', [48000, 0, 8], -0.01], ['48000 nm p2p at 12 km (depth 5)', [48000, 12, 5], 0.23], ['48000 nm p2p at 40 km (depth 10)', [48000, 40, 10], 0.84], ['48000 nm p2p at 100 km (depth 0)', [48000, 100, 0], 1.38], ['48000 nm p2p at 100 km (depth 15)', [48000, 100, 15], 1.39], ['48000 nm p2p at 250 km (depth 12)', [48000, 250, 12], 2.11], ['48000 nm p2p at 480 km (depth 30)', [48000, 480, 30], 2.86]], [['2000 nm p2p at 250 km (depth 12)', [2000, 250, 12], 0.73], ['900000 nm p2p at 0 km (depth 8)', [900000, 0, 8], 1.26], ['900000 nm p2p at 12 km (depth 5)', [900000, 12, 5], 1.51], ['900000 nm p2p at 40 km (depth 10)', [900000, 40, 10], 2.12], ['900000 nm p2p at 100 km (depth 0)', [900000, 100, 0], 2.65], ['900000 nm p2p at 100 km (depth 15)', [900000, 100, 15], 2.66], ['900000 nm p2p at 250 km (depth 12)', [900000, 250, 12], 3.38], ['900000 nm p2p at 480 km (depth 30)', [900000, 480, 30], 4.13]]]
for label, args, expected in fixtures[N - 1]:
    check(label, solve(*args), expected)
print(json.dumps({"observations": observations, "passed": all(x["passed"] for x in observations)}, ensure_ascii=False))
raise SystemExit(0 if all(x["passed"] for x in observations) else 1)
Boundary fixtureActualExpectedOutcome
0 nm p2p at 0 km (depth 8)NoneNonePassed
0 nm p2p at 12 km (depth 5)NoneNonePassed
0 nm p2p at 40 km (depth 10)NoneNonePassed
0 nm p2p at 100 km (depth 0)NoneNonePassed
0 nm p2p at 100 km (depth 15)NoneNonePassed
0 nm p2p at 250 km (depth 12)NoneNonePassed
150 nm p2p at 0 km (depth 8)-2.52-2.52Passed
150 nm p2p at 12 km (depth 5)-2.27-2.27Passed

SHA-256 / 1102858606c7ea41f4d5d387292062eb5bec55d66786048019918b7357636f38

Verification & scope

Stipulated deterministic teaching model of a seismological magnitude procedure; constants and tables are fixed by the contract and no claim of agency or standards conformance is made. This reproducer isolates one failure mechanism. Results cover the supplied fixtures. Variants within a family share a test contract and should remain grouped when constructing evaluation splits. Related mechanisms with a shared evaluation_group must also remain together; these controlled models are not independent production incidents.

Observations recorded using Python 3.12.14 at 2026-09-29T14:48:28.759155+00:00.

Case digest / edfb6002e5125deca84ba20591815454750d9a0204760ec7fb0d16266c05d4bc