Data access

<THIS PAGE IS A WORK IN PROGRESS - PLEASE REPORT ISSUES  / CHECK BACK OFTEN! >

Datasets across the AlpArray Region - the AlpArray Seismic Network

The AlpArray region is defined as the region within 250 km distance from a smoothed 800m contour line around the Alps. The polygon describing this is shown below

New image

Within this region, there are stations operated by existing permanent networks that are complimented by the dense AlpArray Seismic Network temporary experiment. Together, the permanent and temporary networks are the AlpArray Seismic Network (AASN).

With very few exceptions, all broadband stations operating within the AlpArray limits belonging to permanent seismic networks are available to the community (and hence also the AASN) without restriction, and the highest available sampling rates - in general this means HH? channels with sampling rates of at least 100sps. All these permanent stations have realtime communications and data is available in the archive typically within minutes of realtime. Data can be downloaded from these stations using EIDA tools ( EIDA data access portal, arclink_fetch, fdsnws). These stations are operated under their usual network codes, eg. CH for the National Swiss Network, FR for RESIF etc.

The AlpArray temporary network has network code Z3, all stations in the temporary stations use this network code independent of where or whom has installed and operates the station. All Z3 stations are broadband - at least 30s corner period, though the majority are at least 120s, and are sampled at at least 100sps. Around half of Z3 stations stream data in realtime or with small delays, the rest have no communications. For stations with no communications, its is targeted to make data available within 3 months maximum. All Z3 data are permanently archived at EIDA nodes, and the waveforms are available using arclink_fetch and the EIDA data access portal.

The official period of the AASN started on 1.1.2016.

Data Restriction for Z3 network

Z3 station data is available as soon as it is included in EIDA with password protection via EIDA arclink services to registered members of institutions part of the Core Group in order to carry out research following approved topics.

At latest one year after operation of the AlpArray Seismic Network has ended (by decision of the Core Group) the waveform data will be freely shared among the entire AlpArray Working Group.

Three years after operation of the AlpArray Seismic Network has ended (by decision of the Core Group) the waveform data will be freely available to the public.

Virtual Network Code _ALPARRAY 

In order to identify and facilitate easy and rapid access to all AlpArray Seismic Network stations, we have created a virtual network for AlpArray - _ALPARRAY. This network code can be used using both the EIDA data access portal and arclink_fetch.

The rules defining which stations are included in _ALPARRAY are:

- include all permanent stations that are within the 250km distance polygon AND providing HH? channel data AND  operational at any time between 1.1.2016 and 1.1.2020

- include all Z3 network code stations (temporary AlpArray stations)

There are a number of exceptions that are also included (whitelisted):
- G.ECH, RD.LOR, RD.ORIF, RD.PGF (stations only recording at 20sps)
- CR.BRJN, CR.KALN, CR.MOSL, CR.NVLJ, CR.OZLJ, CR.PTJ, CR.RIY, CR.UDBI, CR.ZAG (stations only recording at 50sps)
- CR.HVAR, CR.KIJV, CR.MORI, CR.STON (50sps AND outside polygon)
- HS.ERNH (station just outside polygon)

Data Access - General Information

  • Data from the permanent networks in the region is openly available from EIDA and can be accessed using all tools below as well as fdsnws without password
  • AlpArray Z3 data is restricted to participants in the AASN Core Group of AlpArray
  • Z3 Access Credentials are in the form of a unique keyword that is provided to each AlpArray Seismic Network (AASN) Core Member – for Current list of Core Members, see http://www.alparray.ethz.ch/seismic_network/backbone/management/
  • John Clinton ( jclinton AT sed.ethz.ch ) provides credentials to the contact person for each AASN Core Member – see Science Council at http://www.alparray.ethz.ch/organisation/management/ for the list of contact persons. If you are from an AASN Core Member and would like access to the Z3 data, get in touch with your contact person.

Data Access - Summary of Data Archival 

  • all temporary stations are collected using network code Z3 (the EASI Complementary Experiment in contrast is XT)
  • Z3 waveform data is archived across various EIDA Nodes. The EIDA node archiving a particular station depends on the agency that collects the data rather than the country the stations are placed in.
    • SED / ETHZ EIDA node archives SED/SEG ETHZ data
    • INGV node archives INGV data
    • RESIF node archives data collected by French institutions
    • GFZ / BGR / LMU nodes archives data collected by German institutions
    • ODC node archives data collected by institutions from Austria, Hungary, Czech Republic.
  • waveform data is available (only using password) via the EIDA arclink portal or arclink_fetch
  • station metadata can be explored using EIDA fdsnws (see http://www.seismo.ethz.ch/prod/archiv_seisdata/index_EN for details on fdsnws:station). Eg.:

http://eida.ethz.ch/fdsnws/station/1/query?&net=Z3&station=*&format=text&level=sta&startafter=2015-01-01&nodata=404

Data Access - Access Details 

  • instead of your e-mail address, use keyword uniquely provided to each AASN Core Member
  • this keyword allows you access all closed Z3 data, and can also be used also to download unrestricted (open) data from permanent networks.
  • Waveform data is retrieved in either mseed of full seed format. Metadata can be retrieved as Dataless Seed or inventoryXML.
  1. EIDA Portal (target: exploring datasets, irregular download)
  1. arclink_fetch (target: advanced users - for scripting and batch jobs)
  • an example for Z3, downloading 5 minutes of full seed HHZ data from the M4.6 Bosnia Earthquake on 2016/2/14:
$ arclink_fetch -k fseed -u XXXXXXX -o reqZ3_quake_20160214.seed -v reqZ3_quake_20160214.txt -g -p

where XXXXXX is the AASN keyword, txt is a text file in the same directory and contains the following lines:

2016,02,14,14,51,30 2016,02,14,14,56,30 Z3 * HHZ *
2016,02,14,14,51,30 2016,02,14,14,56,30 SL BOJS HHZ *
2016,02,14,14,51,30 2016,02,14,14,56,30 MN BLY HHZ *
  • the first line only can also be used to retrieve only the Z3 stations. Z3 can also be replaced by _ALPARRAY

 

  1. Obspy Arclink client (target: Obspy users - for scripting and interactive use)

Below a quick example on how to query data and visualize data with Obspy. A more complete documentation of Obspy refer to http://docs.obspy.org/ or more details specifically on the Arclink client and a comprehensive list of parameters, see http://docs.obspy.org/packages/obspy.clients.arclink.html#module-obspy.clients.arclink. The example below is for ObsPy 1.0.0.

  • Import Obspy and the client
import obspy
from obspy.clients.arclink import Client
# set to the secret user provided
user_Z3 = '<secret user>'
  • Example data
t1=obspy.UTCDateTime('2016-02-21')t2=obspy.UTCDateTime('2016-02-21T01:00')
net='IV'sta='ACER'loc=''cha='HH*'
# Or an example of Z3 datanet='Z3'sta='A313A'loc='*'cha='*'
  • Setup client for in routing mode
client = Client(user=user_Z3)
# or use one specific server for the first contact (routing)client = Client(user=user_Z3, host='eida.rm.ingv.it') 
  • Query data
st=client.get_waveforms(net,sta,loc,cha,t1,t2,route=True)
# print as summary (list of traces) of the retried Stream objectprint(st) # Quick plot of the tracesst.plot()  
  • Setup client for in proxy mode (use proxy server at port 18001)
client_proxy = Client(user=user_Z3, port=18001, institution='INGV')
  • Query data (proxy mode)
st=client_proxy.get_waveforms(net,sta,loc,cha,t1,t2,route=False)
print(st)st.plot()