Using Geoscience Skills

Using Geoscience Skills

Meta-skill for discovering, routing, and composing the geoscience skills library. This skill maps user intent to domain skills, workflow skills, slash commands, and agents.

Domain Routing Table

Match user intent keywords to the appropriate domain skill.

| Keywords / Triggers | Skill | Domain | |---------------------|-------|--------| | SEG-Y, seismic traces, trace headers, inline, crossline | segyio | Seismic I/O | | waveform, earthquake, FDSN, seismogram, miniSEED | obspy | Seismology | | surface wave, dispersion, Rayleigh, Love wave | disba | Seismology | | LAS, well logs, wireline, borehole curves | lasio | Well Logs | | DLIS, RP66, array logs, modern well data | dlisio | Well Logs | | well analysis, curve QC, multi-well, despike | welly | Well Logs | | petrophysics, Sw, porosity, formation evaluation | petropy | Petrophysics | | lithology, stratigraphy, striplog, facies log | striplog | Stratigraphy | | 3D model, geology, implicit surface, faults | gempy | 3D Modelling | | fold modelling, structural frame, Loop3D | loopstructural | 3D Modelling | | GIS, spatial data prep, borehole to GemPy | gemgis | GIS Preprocessing | | inversion, DC resistivity, magnetics, gravity, EM | simpeg | Inversion | | ERT, SRT, IP, near-surface inversion | pygimli | Inversion | | PDE, wave equation, finite differences, stencil | devito | Simulation | | linear operator, inverse problem, sparsity | pylops | Inverse Problems | | gravity, magnetic, Bouguer, upward continuation | harmonica | Potential Fields | | AVO, Zoeppritz, Gassmann, fluid substitution, wavelet | bruges | Rock Physics | | gridding, interpolation, spatial, Verde | verde | Spatial Analysis | | variogram, kriging, GSLIB, geostatistics | geostatspy | Geostatistics | | variogram fitting, scikit-learn style geostat | scikit-gstat | Geostatistics | | spatial regression, GWR, GNNWR, non-stationarity, coefficient mapping | gnnwr | Spatial Regression | | groundwater, time series, pumping test | pastas | Hydrology | | landscape, erosion, surface processes, DEM | landlab | Surface Processes | | stereonet, strike, dip, poles, structural | mplstereonet | Structural Geology | | geochemistry, REE, spider diagram, ternary | pyrolite | Geochemistry | | GPR, ground-penetrating radar, radargram | gprpy | Near-Surface | | magnetotellurics, MT, impedance tensor | mtpy | Near-Surface | | NetCDF, xarray, multi-dimensional, climate | xarray | Data Formats | | 3D visualization, mesh, VTK, point cloud | pyvista | Visualization | | data download, sample data, cache, fetch | pooch | Utilities |

Workflow Skills

Workflow skills chain multiple domain skills into end-to-end pipelines.

| Workflow | Slash Command | Skill Chain | |----------|---------------|-------------| | Seismic Interpretation | /seismic-workflow | segyio -> obspy -> bruges -> disba -> pyvista | | Well Log Evaluation | /well-analysis | lasio/dlisio -> welly -> petropy -> striplog -> pyvista | | Geological Modelling | /model-3d | gemgis -> gempy/loopstructural -> pyvista | | Geophysical Inversion | /inversion-workflow | simpeg/pygimli -> verde -> pyvista | | Rock Physics & AVO | /rock-physics | lasio/welly -> bruges -> segyio |

Available Agents

| Agent | Purpose | Typical Trigger | |-------|---------|-----------------| | data-qc-reviewer | Automated data quality checks across formats | "QC my data", "check data quality" | | geoscience-mentor | Guided explanations of geoscience concepts and methods | "explain", "teach me", "what is" |

All 30 Domain Skills by Category

Seismic and Seismology

• segyio -- SEG-Y file I/O, trace and header access
• obspy -- seismological waveform processing, FDSN services
• disba -- surface wave dispersion (Rayleigh, Love)

Well Log Analysis

• lasio -- LAS file reading and writing
• dlisio -- DLIS/RP66 binary well log parsing
• welly -- well data analysis, curve QC, multi-well projects
• petropy -- petrophysical analysis, formation evaluation
• striplog -- lithological and stratigraphic log display

3D Geological Modelling

• gempy -- implicit 3D geological modelling
• loopstructural -- 3D modelling with fold and fault support
• gemgis -- spatial data preprocessing for GemPy

Geophysical Inversion

• simpeg -- multi-method geophysical inversion framework
• pygimli -- ERT, SRT, IP inversion with simple API
• devito -- symbolic PDE solver for wave propagation
• pylops -- linear operators for inverse problems

Potential Fields and Rock Physics

• harmonica -- gravity and magnetic data processing
• bruges -- AVO, Gassmann, wavelets, elastic moduli

Spatial Analysis and Geostatistics

• verde -- spatial gridding and interpolation
• geostatspy -- variograms, kriging (GSLIB-style)
• scikit-gstat -- geostatistics with scikit-learn API
• gnnwr -- geographically weighted neural network regression

Hydrology and Surface Processes

• pastas -- groundwater time series modelling
• landlab -- landscape evolution modelling

Structural Geology and Geochemistry

• mplstereonet -- stereonet plots for orientation data
• pyrolite -- geochemical analysis and diagrams

Near-Surface Geophysics

• gprpy -- GPR data processing
• mtpy -- magnetotelluric data analysis

Data Formats and Visualization

• xarray -- NetCDF, multi-dimensional labeled arrays
• pyvista -- 3D mesh visualization and analysis
• pooch -- data file fetching and caching

Skill Composition Rules

Chain skills when a task spans multiple stages of a geoscience workflow.

Composition Patterns

Data Loading -> Processing -> Modelling -> Visualization

1. Always start with a data I/O skill (segyio, lasio, dlisio, xarray)
2. Use processing skills for QC and transformation (welly, obspy, verde)
3. Apply domain modelling (bruges, gempy, simpeg, pygimli)
4. Finish with visualization (pyvista, matplotlib via domain skill)

When to Chain vs Use Standalone

| Scenario | Approach | |----------|----------| | Single file format question | Standalone domain skill | | End-to-end analysis pipeline | Workflow skill to orchestrate | | Data QC across formats | data-qc-reviewer agent | | Concept explanation | geoscience-mentor agent | | Multi-library code generation | Chain domain skills in order |

Dependency Awareness

When composing skills, respect data flow:

# Correct: segyio loads, obspy processes, bruges models
import segyio
import obspy
from bruges.reflection import zoeppritz

# Load with segyio
with segyio.open('seismic.sgy') as f:
    data = f.trace[:]

# Process with obspy (convert to Stream if needed)
# Model with bruges
Rpp = zoeppritz(vp1, vs1, rho1, vp2, vs2, rho2, theta)

When to Use This Skill

This is the discovery and routing skill. Use it when:

• Starting a new geoscience coding session and unsure which library to use
• A user request spans multiple geoscience domains
• You need to find the right slash command or workflow for a task
• Composing multiple domain skills into a pipeline
• Looking up which skill handles a specific file format or analysis type

This skill does not perform any analysis itself. It directs to the appropriate domain skill, workflow skill, or agent for execution.

Quick Decision Tree

User wants to...
  |
  +-- Load/write a file? --> Check format:
  |     SEG-Y -> segyio    LAS -> lasio    DLIS -> dlisio
  |     NetCDF -> xarray    VTK -> pyvista
  |
  +-- Process signals? --> obspy (seismology), welly (well logs)
  |
  +-- Build a model?
  |     3D geology -> gempy or loopstructural
  |     Rock physics -> bruges
  |     Inversion -> simpeg or pygimli
  |
  +-- Visualize results? --> pyvista (3D), matplotlib (2D via domain skill)
  |
  +-- Run a full pipeline? --> Use workflow skills above

Source & license

This open-source skill is cataloged on AgentStack and links to its original source — we do not rehost the code.

• Author: SteadfastAsArt
• Source: SteadfastAsArt/geoscience-skills
• License: MIT

Install and usage instructions live in the source repository linked above.