rFDSN-package: Get seismic data from the International Federation of Digital...

In rFDSN: Get Seismic Data from the International Federation of Digital Seismograph Networks

Description

Gather information about seismic networks (stations, channels, locations, etc), download seismic time series in miniseed format, and find historical earthquake data (origins, magnitudes, etc).

Details

Package:	rFDSN
Type:	Package
Version:	0.0.0
Date:	2014-09-22
License:	GPL v3

Author(s)

Daniel C. Bowman [e-mail]

References

FDSN Web Services:
http://www.fdsn.org/webservices/
List of base URLS for data download:
http://www.fdsn.org/webservices/datacenters/

Examples

#Find all stations with operating vertical seismic channels
#within 2 degrees of Socorro, New Mexico,
#between 2010 and 2013
#available through the IRIS DMC server.

library(rFDSN)

## Not run: 

#IRIS station server
base.url <- "http://service.iris.edu/fdsnws/station/1/"

#Selection criteria
parameters <- list(
    name = c(
    "latitude",
    "longitude",
    "maxradius",
    "cha",
    "start",
    "end"),
    value = c(
    34.068110, #Latitude
    -106.901847, #Longitude
    2, #Degrees from center point
    "*Z",#Channel (with wildcards)
    "2010-01-01", #Start time
    "2014-01-01") #End time
    )
 
networks <- FDSNGetNetworks(base.url, parameters)

#Download pressure data at station Y22D on August 9, 2014

#IRIS data download server
base.url <- "http://service.iris.edu/fdsnws/dataselect/1/"
parameters <- list(
    name = c(
    "sta",
    "cha",
    "start",
    "end"),
    value = c(
    "Y22D", #Station name
    "BDF", #SEED channel name
    "2014-08-09T00:00:00.000", #Midnight GMT
    "2014-08-09T23:59:59.999") #Midnight GMT
   )
miniseed.file <- FDSNGetTimeSeries(base.url, parameters, save.file = "Y22D_BDF.mseed")

#Get all earthquakes within 2 degrees of Socorro, new Mexico, 
#between 2010 and 2013

#IRIS event server
base.url <- "http://service.iris.edu/fdsnws/event/1/"

#Selection criteria
parameters <- list(
    name = c(
    "latitude",
    "longitude",
    "maxradius",
    "start",
    "end"),
    value = c(
    34.068110, #Latitude
    -106.901847, #Longitude
    2, #Degrees from center point
    "2010-01-01", #Start time
    "2014-01-01") #End time
    )

events <- FDSNGetEvents(base.url, parameters)


## End(Not run)

rFDSN documentation built on May 1, 2019, 6:51 p.m.