obspy-seismology

# ObsPy Seismology A comprehensive skill for seismological data acquisition, processing, and analysis using the ObsPy framework. Covers FDSN waveform retrieval, instrument response removal, P/S phase arrival picking, basic moment tensor inversion, and seismicity visualisation. ## When to Use This Skill Use this skill when you need to: - Download seismogram waveforms from IRIS, GEOFON, ORFEUS, or any FDSN-compliant data centre - Remove instrument response and convert raw counts to physical units (m/s, m/s², Pa) - Filter, decimate, and taper seismic traces before analysis - Automatically pick P and S phase arrivals with STA/LTA or kurtosis detectors - Estimate earthquake source parameters (origin time, location, focal mechanism) - Build seismicity catalogues and map earthquake distributions - Compute and visualise spectrograms and particle motion diagrams This skill is **not** appropriate for: - Real-time continuous acquisition from hardware digitisers (use SeisComP or Earthworm instead) - Full waveform tomography (use SPECFEM or SALVUS skill) - Distributed seismic arrays with >100 stations in a single session ## Background & Key Concepts ### ObsPy Data Model ObsPy organises seismic data in three nested containers: | Class | Contains | Analogy | |---|---|---| | `Stream` | list of `Trace` objects | a multi-channel recording session | | `Trace` | 1-D NumPy array + `Stats` | a single channel | | `Stats` | dict-like metadata | network, station, channel, start time, sa...