Loading skill
Install any skill in seconds. Free to start, no credit card required.
Get Started Free →brycewang-stanford/seismology-data-guide
Earthquake data analysis, seismogram processing, and seismic research
Passed 19 of 24 test cases with the skill, and 15 of 24 without it.
and 1 more in the full report
4 fixed15 passed either way5 still failing0 broke
Tokens, turns and time describe resource use, not a dollar price.
One of 4 test cases the skill fixed. Same model, same prompt, with and without the skill.
Body-wave magnitude (mb) is derived from the amplitude of compression P-waves measured near what target frequency?
| Checked | Without skill | With skill |
|---|---|---|
| Specifies P-wave amplitude measurement at 1 Hz | Missed | Passed |
The instructions your model reads when the skill loads.
A skill for processing seismic data, analyzing earthquake catalogs, and working with seismograms using standard tools in observational seismology. Covers data retrieval from global networks, waveform processing with ObsPy, magnitude estimation, focal mechanism analysis, and seismic hazard assessment.
| Data Center | Abbreviation | Coverage | Access |
|---|---|---|---|
| IRIS Data Management Center | IRIS DMC | Global broadband | FDSN Web Services |
| European Integrated Data Archive | EIDA | European networks | FDSN Web Services |
| USGS Earthquake Hazards Program | USGS EHP | Global catalog | API + ComCat |
| International Seismological Centre | ISC | Global bulletin | ISC web services |
| NIED F-net | F-net | Japan broadband | NIED website |
from obspy.clients.fdsn import Client
from obspy import UTCDateTime
client = Client("IRIS")
# Fetch earthquake catalog for a region and time window
catalog = client.get_events(
starttime=UTCDateTime("2024-01-01"),
endtime=UTCDateTime("2024-12-31"),
minmagnitude=5.0,
maxmagnitude=9.0,
minlatitude=30.0, maxlatitude=45.0,
minlongitude=125.0, maxlongitude=150.0,
orderby="magnitude",
)
print(f"Found {len(catalog)} events")
for event in catalog[:5]:
origin = event.preferred_origin()
mag = event.preferred_magnitude()
print(f" M{mag.mag:.1f} {origin.time} "
test cases passed, vs 15 / 24 without
More skills from the same source and nearby in the registry, with their headline benchmark lift.