Web Requests¶
Data requests use get_data!
as a wrapper to either FDSN or IRIS data services; for live streaming, see SeedLink.
-
get_data!(S, method, channels; KWs)
-
S = get_data(method, channels; KWs)
--src
using the server list (also available by typing ?seis_www
).?chanspec
at the Julia prompt for more info.?SeisIO.KW
.--s
: Start time--t
: Termination (end) timeExamples¶
get_data!(S, "FDSN", "UW.SEP..EHZ,UW.SHW..EHZ,UW.HSR..EHZ", "IRIS", t=(-600))
: using FDSNWS, get the last 10 minutes of data from three short-period vertical-component channels at Mt. St. Helens, USA.get_data!(S, "IRIS", "CC.PALM..EHN", "IRIS", t=(-120), f="sacbl")
: using IRISWS, fetch the last two minutes of data from component EHN, station PALM (Palmer Lift (Mt. Hood), OR, USA,), network CC (USGS Cascade Volcano Observatory, Vancouver, WA, USA), in bigendian SAC format, and merge into SeisData structure S.get_data!(S, "FDSN", "CC.TIMB..EHZ", "IRIS", t=(-600), w=true)
: using FDSNWS, get the last 10 minutes of data from channel EHZ, station TIMB (Timberline Lodge, OR, USA), save the data directly to disk, and add it to SeisData structure S.S = get_data("FDSN", "HV.MOKD..HHZ", "IRIS", s="2012-01-01T00:00:00", t=(-3600))
: using FDSNWS, fill a new SeisData structure S with an hour of data ending at 2012-01-01, 00:00:00 UTC, from HV.MOKD..HHZ (USGS Hawai’i Volcano Observatory).
FDSN Queries¶
The International Federation of Digital Seismograph Networks (FDSN) is a global organization that supports seismology research. The FDSN web protocol offers near-real-time access to data from thousands of instruments across the world.
FDSN queries in SeisIO are highly customizable; see data keywords list and channel id syntax.
Data Query¶
-
get_data!(S, "FDSN", channels; KWs)
-
S = get_data("FDSN", channels; KWs)
FDSN data query with get_data! wrapper.
--s
: Start time--t
: Termination (end) timexml_file
: Name of XML file to save station metadataStation Query¶
-
FDSNsta!(S, chans, KW)
-
S = FDSNsta(chans, KW)
Fill channels chans of SeisData structure S with information retrieved from remote station XML files by web query.
Event Header Query¶
-
H = FDSNevq(ot)
Shared keywords: evw, rad, reg, mag, nev, src, to, v, w
Multi-server query for the event(s) with origin time(s) closest to ot. Returns a SeisHdr.
Notes:
- Specify ot as a string formatted YYYY-MM-DDThh:mm:ss in UTC (e.g. “2001-02-08T18:54:32”). Returns a SeisHdr array.
- Incomplete string queries are read to the nearest fully-specified time constraint; thus, FDSNevq(“2001-02-08”) returns the nearest event to 2001-02-08T00:00:00.
- If no event is found in the specified search window, FDSNevq exits with an error.
Event Header and Data Query¶
-
Ev = FDSNevt(ot::String, chans::String)
Get trace data for the event closest to origin time ot on channels chans. Returns a SeisEvent.
--len
: desired record length in minutes.IRIS Queries¶
Incorporated Research Institutions for Seismology (IRIS) is a consortium of universities dedicated to the operation of science facilities for the acquisition, management, and distribution of seismological data.
Data Query Features¶
- Stage zero gains are removed from trace data; all IRIS data will appear to have a gain of 1.0.
- IRISWS disallows wildcards in channel IDs.
- Channel spec must include the net, sta, cha fields; thus, CHA = “CC.VALT..BHZ” is OK; CHA = “CC.VALT” is not.
Phase Onset Query¶
Command-line interface to IRIS online travel time calculator, which calls TauP [1-2]. Returns a matrix of strings.
Specify Δ in decimal degrees, z in km with + = down.
-model
: velocity model (defaults to “iasp91”)References
- Crotwell, H. P., Owens, T. J., & Ritsema, J. (1999). The TauP Toolkit: Flexible seismic travel-time and ray-path utilities, SRL 70(2), 154-160.
- TauP manual: http://www.seis.sc.edu/downloads/TauP/taup.pdf