FA-71446 / Seismic magnitude estimation / Open access

Teleseismic body-wave magnitude: short-period limit · case 01

3.0 s readings are discarded.

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

ROOT CAUSE

The period limit excludes exactly 3 s.

VERIFIED REPAIR

Accept periods up to and including 3.0 s.

Unsuccessful approach: Tightening to 2 s discards valid 2-3 s readings.

Case contract

P amplitude in micrometres, period in s, epicentral distance in degrees. Valid only for 20 <= delta <= 100 and 0 < period <= 3.0 (else None). Q(delta) is linearly interpolated from [(20,6.0),(30,6.5),(40,6.8),(60,6.9),(80,7.0),(100,7.4)]. mb = log10(A/T) + Q, rounded 0.01.

Why this case matters

mb is the first teleseismic magnitude available for global events.

1 / The failure

Exit 1

"""Failure Map reference implementation. Python standard library only."""
import json
import math
N = 1
observations = []
def solve(amp_um, period_s, delta_deg):
    table = [(20, 6.0), (30, 6.5), (40, 6.8), (60, 6.9), (80, 7.0), (100, 7.4)]
    if not (20 <= delta_deg <= 100):
        return None
    if period_s <= 0 or period_s >= 3.0:
        return None
    q = None
    for (d0, q0), (d1, q1) in zip(table, table[1:]):
        if d0 <= delta_deg <= d1:
            q = q0 + (q1 - q0) * (delta_deg - d0) / (d1 - d0)
            break
    return round(math.log10(amp_um / period_s) + q, 2)
def check(label, actual, expected):
    observations.append({"check": label, "actual": actual, "expected": expected, "passed": actual == expected})
fixtures = [[['0.05 um 0.5 s at 15 deg', [0.05, 0.5, 15], None], ['0.05 um 0.5 s at 20 deg', [0.05, 0.5, 20], 5.0], ['0.05 um 0.5 s at 25 deg', [0.05, 0.5, 25], 5.25], ['0.05 um 0.5 s at 45 deg', [0.05, 0.5, 45], 5.83], ['0.05 um 0.5 s at 70 deg', [0.05, 0.5, 70], 5.95], ['0.05 um 0.5 s at 99 deg', [0.05, 0.5, 99], 6.38], ['0.05 um 3.0 s at 20 deg', [0.05, 3.0, 20], 4.22], ['0.05 um 3.0 s at 25 deg', [0.05, 3.0, 25], 4.47]], [['0.05 um 0.5 s at 101 deg', [0.05, 0.5, 101], None], ['0.05 um 1.0 s at 15 deg', [0.05, 1.0, 15], None], ['0.05 um 1.0 s at 20 deg', [0.05, 1.0, 20], 4.7], ['0.05 um 1.0 s at 25 deg', [0.05, 1.0, 25], 4.95], ['0.05 um 1.0 s at 45 deg', [0.05, 1.0, 45], 5.52], ['0.05 um 1.0 s at 70 deg', [0.05, 1.0, 70], 5.65], ['0.05 um 3.0 s at 70 deg', [0.05, 3.0, 70], 5.17], ['0.4 um 3.0 s at 20 deg', [0.4, 3.0, 20], 5.12]], [['0.05 um 1.0 s at 100 deg', [0.05, 1.0, 100], 6.1], ['0.05 um 1.0 s at 101 deg', [0.05, 1.0, 101], None], ['0.05 um 2.0 s at 15 deg', [0.05, 2.0, 15], None], ['0.05 um 2.0 s at 20 deg', [0.05, 2.0, 20], 4.4], ['0.05 um 2.0 s at 25 deg', [0.05, 2.0, 25], 4.65], ['0.05 um 2.0 s at 45 deg', [0.05, 2.0, 45], 5.22], ['0.4 um 3.0 s at 20 deg', [0.4, 3.0, 20], 5.12], ['0.4 um 3.0 s at 100 deg', [0.4, 3.0, 100], 6.52]], [['0.05 um 2.0 s at 99 deg', [0.05, 2.0, 99], 5.78], ['0.05 um 2.0 s at 100 deg', [0.05, 2.0, 100], 5.8], ['0.05 um 2.0 s at 101 deg', [0.05, 2.0, 101], None], ['0.05 um 3.0 s at 15 deg', [0.05, 3.0, 15], None], ['0.05 um 3.0 s at 20 deg', [0.05, 3.0, 20], 4.22], ['0.05 um 3.0 s at 25 deg', [0.05, 3.0, 25], 4.47], ['0.4 um 3.0 s at 70 deg', [0.4, 3.0, 70], 6.07], ['2.5 um 3.0 s at 99 deg', [2.5, 3.0, 99], 7.3]], [['0.05 um 3.0 s at 70 deg', [0.05, 3.0, 70], 5.17], ['0.05 um 3.0 s at 99 deg', [0.05, 3.0, 99], 5.6], ['0.05 um 3.0 s at 100 deg', [0.05, 3.0, 100], 5.62], ['0.05 um 3.0 s at 101 deg', [0.05, 3.0, 101], None], ['0.05 um 3.5 s at 15 deg', [0.05, 3.5, 15], None], ['0.05 um 3.5 s at 20 deg', [0.05, 3.5, 20], None], ['0.05 um 3.5 s at 25 deg', [0.05, 3.5, 25], None], ['2.5 um 3.0 s at 20 deg', [2.5, 3.0, 20], 5.92]]]
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 fixture	Actual	Expected	Outcome
0.05 um 0.5 s at 15 deg	`None`	`None`	Passed
0.05 um 0.5 s at 20 deg	`5.0`	`5.0`	Passed
0.05 um 0.5 s at 25 deg	`5.25`	`5.25`	Passed
0.05 um 0.5 s at 45 deg	`5.83`	`5.83`	Passed
0.05 um 0.5 s at 70 deg	`5.95`	`5.95`	Passed
0.05 um 0.5 s at 99 deg	`6.38`	`6.38`	Passed
0.05 um 3.0 s at 20 deg	`None`	`4.22`	Failed
0.05 um 3.0 s at 25 deg	`None`	`4.47`	Failed

SHA-256 / c4c33bdb534bd669a9017787a36dc73509c2a782e661fccb083f26e627e1bdfc

2 / The unsuccessful fix

Exit 1

"""Failure Map reference implementation. Python standard library only."""
import json
import math
N = 1
observations = []
def solve(amp_um, period_s, delta_deg):
    table = [(20, 6.0), (30, 6.5), (40, 6.8), (60, 6.9), (80, 7.0), (100, 7.4)]
    if not (20 <= delta_deg <= 100):
        return None
    if period_s <= 0 or period_s > 2.0:
        return None
    q = None
    for (d0, q0), (d1, q1) in zip(table, table[1:]):
        if d0 <= delta_deg <= d1:
            q = q0 + (q1 - q0) * (delta_deg - d0) / (d1 - d0)
            break
    return round(math.log10(amp_um / period_s) + q, 2)
def check(label, actual, expected):
    observations.append({"check": label, "actual": actual, "expected": expected, "passed": actual == expected})
fixtures = [[['0.05 um 0.5 s at 15 deg', [0.05, 0.5, 15], None], ['0.05 um 0.5 s at 20 deg', [0.05, 0.5, 20], 5.0], ['0.05 um 0.5 s at 25 deg', [0.05, 0.5, 25], 5.25], ['0.05 um 0.5 s at 45 deg', [0.05, 0.5, 45], 5.83], ['0.05 um 0.5 s at 70 deg', [0.05, 0.5, 70], 5.95], ['0.05 um 0.5 s at 99 deg', [0.05, 0.5, 99], 6.38], ['0.05 um 3.0 s at 20 deg', [0.05, 3.0, 20], 4.22], ['0.05 um 3.0 s at 25 deg', [0.05, 3.0, 25], 4.47]], [['0.05 um 0.5 s at 101 deg', [0.05, 0.5, 101], None], ['0.05 um 1.0 s at 15 deg', [0.05, 1.0, 15], None], ['0.05 um 1.0 s at 20 deg', [0.05, 1.0, 20], 4.7], ['0.05 um 1.0 s at 25 deg', [0.05, 1.0, 25], 4.95], ['0.05 um 1.0 s at 45 deg', [0.05, 1.0, 45], 5.52], ['0.05 um 1.0 s at 70 deg', [0.05, 1.0, 70], 5.65], ['0.05 um 3.0 s at 70 deg', [0.05, 3.0, 70], 5.17], ['0.4 um 3.0 s at 20 deg', [0.4, 3.0, 20], 5.12]], [['0.05 um 1.0 s at 100 deg', [0.05, 1.0, 100], 6.1], ['0.05 um 1.0 s at 101 deg', [0.05, 1.0, 101], None], ['0.05 um 2.0 s at 15 deg', [0.05, 2.0, 15], None], ['0.05 um 2.0 s at 20 deg', [0.05, 2.0, 20], 4.4], ['0.05 um 2.0 s at 25 deg', [0.05, 2.0, 25], 4.65], ['0.05 um 2.0 s at 45 deg', [0.05, 2.0, 45], 5.22], ['0.4 um 3.0 s at 20 deg', [0.4, 3.0, 20], 5.12], ['0.4 um 3.0 s at 100 deg', [0.4, 3.0, 100], 6.52]], [['0.05 um 2.0 s at 99 deg', [0.05, 2.0, 99], 5.78], ['0.05 um 2.0 s at 100 deg', [0.05, 2.0, 100], 5.8], ['0.05 um 2.0 s at 101 deg', [0.05, 2.0, 101], None], ['0.05 um 3.0 s at 15 deg', [0.05, 3.0, 15], None], ['0.05 um 3.0 s at 20 deg', [0.05, 3.0, 20], 4.22], ['0.05 um 3.0 s at 25 deg', [0.05, 3.0, 25], 4.47], ['0.4 um 3.0 s at 70 deg', [0.4, 3.0, 70], 6.07], ['2.5 um 3.0 s at 99 deg', [2.5, 3.0, 99], 7.3]], [['0.05 um 3.0 s at 70 deg', [0.05, 3.0, 70], 5.17], ['0.05 um 3.0 s at 99 deg', [0.05, 3.0, 99], 5.6], ['0.05 um 3.0 s at 100 deg', [0.05, 3.0, 100], 5.62], ['0.05 um 3.0 s at 101 deg', [0.05, 3.0, 101], None], ['0.05 um 3.5 s at 15 deg', [0.05, 3.5, 15], None], ['0.05 um 3.5 s at 20 deg', [0.05, 3.5, 20], None], ['0.05 um 3.5 s at 25 deg', [0.05, 3.5, 25], None], ['2.5 um 3.0 s at 20 deg', [2.5, 3.0, 20], 5.92]]]
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 fixture	Actual	Expected	Outcome
0.05 um 0.5 s at 15 deg	`None`	`None`	Passed
0.05 um 0.5 s at 20 deg	`5.0`	`5.0`	Passed
0.05 um 0.5 s at 25 deg	`5.25`	`5.25`	Passed
0.05 um 0.5 s at 45 deg	`5.83`	`5.83`	Passed
0.05 um 0.5 s at 70 deg	`5.95`	`5.95`	Passed
0.05 um 0.5 s at 99 deg	`6.38`	`6.38`	Passed
0.05 um 3.0 s at 20 deg	`None`	`4.22`	Failed
0.05 um 3.0 s at 25 deg	`None`	`4.47`	Failed

SHA-256 / fe8023364a9837f5945d5c40ea6af15756213b64479bf1465ff04dd4ed6985ed

3 / The verified repair

Exit 0

"""Failure Map reference implementation. Python standard library only."""
import json
import math
N = 1
observations = []
def solve(amp_um, period_s, delta_deg):
    table = [(20, 6.0), (30, 6.5), (40, 6.8), (60, 6.9), (80, 7.0), (100, 7.4)]
    if not (20 <= delta_deg <= 100):
        return None
    if period_s <= 0 or period_s > 3.0:
        return None
    q = None
    for (d0, q0), (d1, q1) in zip(table, table[1:]):
        if d0 <= delta_deg <= d1:
            q = q0 + (q1 - q0) * (delta_deg - d0) / (d1 - d0)
            break
    return round(math.log10(amp_um / period_s) + q, 2)
def check(label, actual, expected):
    observations.append({"check": label, "actual": actual, "expected": expected, "passed": actual == expected})
fixtures = [[['0.05 um 0.5 s at 15 deg', [0.05, 0.5, 15], None], ['0.05 um 0.5 s at 20 deg', [0.05, 0.5, 20], 5.0], ['0.05 um 0.5 s at 25 deg', [0.05, 0.5, 25], 5.25], ['0.05 um 0.5 s at 45 deg', [0.05, 0.5, 45], 5.83], ['0.05 um 0.5 s at 70 deg', [0.05, 0.5, 70], 5.95], ['0.05 um 0.5 s at 99 deg', [0.05, 0.5, 99], 6.38], ['0.05 um 3.0 s at 20 deg', [0.05, 3.0, 20], 4.22], ['0.05 um 3.0 s at 25 deg', [0.05, 3.0, 25], 4.47]], [['0.05 um 0.5 s at 101 deg', [0.05, 0.5, 101], None], ['0.05 um 1.0 s at 15 deg', [0.05, 1.0, 15], None], ['0.05 um 1.0 s at 20 deg', [0.05, 1.0, 20], 4.7], ['0.05 um 1.0 s at 25 deg', [0.05, 1.0, 25], 4.95], ['0.05 um 1.0 s at 45 deg', [0.05, 1.0, 45], 5.52], ['0.05 um 1.0 s at 70 deg', [0.05, 1.0, 70], 5.65], ['0.05 um 3.0 s at 70 deg', [0.05, 3.0, 70], 5.17], ['0.4 um 3.0 s at 20 deg', [0.4, 3.0, 20], 5.12]], [['0.05 um 1.0 s at 100 deg', [0.05, 1.0, 100], 6.1], ['0.05 um 1.0 s at 101 deg', [0.05, 1.0, 101], None], ['0.05 um 2.0 s at 15 deg', [0.05, 2.0, 15], None], ['0.05 um 2.0 s at 20 deg', [0.05, 2.0, 20], 4.4], ['0.05 um 2.0 s at 25 deg', [0.05, 2.0, 25], 4.65], ['0.05 um 2.0 s at 45 deg', [0.05, 2.0, 45], 5.22], ['0.4 um 3.0 s at 20 deg', [0.4, 3.0, 20], 5.12], ['0.4 um 3.0 s at 100 deg', [0.4, 3.0, 100], 6.52]], [['0.05 um 2.0 s at 99 deg', [0.05, 2.0, 99], 5.78], ['0.05 um 2.0 s at 100 deg', [0.05, 2.0, 100], 5.8], ['0.05 um 2.0 s at 101 deg', [0.05, 2.0, 101], None], ['0.05 um 3.0 s at 15 deg', [0.05, 3.0, 15], None], ['0.05 um 3.0 s at 20 deg', [0.05, 3.0, 20], 4.22], ['0.05 um 3.0 s at 25 deg', [0.05, 3.0, 25], 4.47], ['0.4 um 3.0 s at 70 deg', [0.4, 3.0, 70], 6.07], ['2.5 um 3.0 s at 99 deg', [2.5, 3.0, 99], 7.3]], [['0.05 um 3.0 s at 70 deg', [0.05, 3.0, 70], 5.17], ['0.05 um 3.0 s at 99 deg', [0.05, 3.0, 99], 5.6], ['0.05 um 3.0 s at 100 deg', [0.05, 3.0, 100], 5.62], ['0.05 um 3.0 s at 101 deg', [0.05, 3.0, 101], None], ['0.05 um 3.5 s at 15 deg', [0.05, 3.5, 15], None], ['0.05 um 3.5 s at 20 deg', [0.05, 3.5, 20], None], ['0.05 um 3.5 s at 25 deg', [0.05, 3.5, 25], None], ['2.5 um 3.0 s at 20 deg', [2.5, 3.0, 20], 5.92]]]
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 fixture	Actual	Expected	Outcome
0.05 um 0.5 s at 15 deg	`None`	`None`	Passed
0.05 um 0.5 s at 20 deg	`5.0`	`5.0`	Passed
0.05 um 0.5 s at 25 deg	`5.25`	`5.25`	Passed
0.05 um 0.5 s at 45 deg	`5.83`	`5.83`	Passed
0.05 um 0.5 s at 70 deg	`5.95`	`5.95`	Passed
0.05 um 0.5 s at 99 deg	`6.38`	`6.38`	Passed
0.05 um 3.0 s at 20 deg	`4.22`	`4.22`	Passed
0.05 um 3.0 s at 25 deg	`4.47`	`4.47`	Passed

SHA-256 / 8d3055c9a2d022599829b801b58865f1ae241beb2647abcf37779ab75e3cdbb8

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:29.755377+00:00.

Case digest / ef24767fd5731034ef33f71f7d18fa8c66b1c2492ba4323113c77b6adfc014cf