Wolfram Function Repository

Instant-use add-on functions for the Wolfram Language

Function Repository Resource:

 SeismographData

Source Notebook

Obtain data from a seismographic station for a given date range

Contributed by: Arnoud Buzing

ResourceFunction["SeismographData"][]

lists all supported seismographic clients.

ResourceFunction["SeismographData"][client]

lists all supported networks and stations for the given client.

ResourceFunction["SeismographData"][client,network]

lists all supported stations for the given client and network.

ResourceFunction["SeismographData"][client,network,station,location,channel,{start,end}]

gives seismographic data for the given parameters.

ResourceFunction["SeismographData"][client,network,station,{start,end}]

gives seismographic data for the given parameters, for any location and channel.

Details

Supported client data services include "IRIS".
Supported network codes include "IU".
SeismographData encapsulates functionality from the ObsPy package.

Examples

Basic Examples

List all supported seismographic clients:

In[1]:=
ResourceFunction["SeismographData"][]
Out[1]=

Get a full list of networks and stations provided by the IRIS client (this may take a long time):

In[2]:=
AbsoluteTiming[ResourceFunction["SeismographData"]["IRIS"]]
Out[2]=

Get all stations for the "IU" network:

In[3]:=
ResourceFunction["SeismographData"]["IRIS", "IU"]
Out[3]=

Get information for a specific station:

In[4]:=
ResourceFunction["SeismographData"]["IRIS", "IU", "ANMO"]
Out[4]=

Get all seismographic information for the Albuquerque station for the given date range and any location and channel:

In[5]:=
timeseries = ResourceFunction["SeismographData"]["IRIS", "IU", "ANMO", "*", "*", {DateObject[{2010, 2, 27, 6, 45, 0}, "Instant", "Gregorian", "UTC"], DateObject[{2010, 2, 27, 7, 45, 0}, "Instant", "Gregorian", "UTC"]}];

Plot the first three of the seismographic timeseries:

In[6]:=
GraphicsGrid[ Map[{DateListPlot[#, Sequence[ PlotRange -> All, AspectRatio -> 1/4]]} &, Take[timeseries, 3]]]
Out[6]=

You can also omit the "*" wildcards for location and channel:

In[7]:=
timeseries = ResourceFunction["SeismographData"]["IRIS", "IU", "ANMO", {DateObject[{2010, 2, 27, 6, 45, 0}, "Instant", "Gregorian", "UTC"], DateObject[{2010, 2, 27, 7, 45, 0}, "Instant", "Gregorian", "UTC"]}];
In[8]:=
GraphicsGrid[ Map[{DateListPlot[#, Sequence[ PlotRange -> All, AspectRatio -> 1/4]]} &, Take[timeseries, 3]]]
Out[8]=

Publisher

Arnoud Buzing

Requirements

Wolfram Language 14.0 (January 2024) or above

Version History

  • 1.0.0 – 23 October 2024

Related Resources

License Information