Instrument response files for seismic stations in South Korea

Citation

If you use it, please cite the following article in your work: 

Lim and Kim (2020), A dataset of seismic sensor responses of South Korea seismic stations. Journal of the Geological Society of Korea, 56(4), 515-524, http://dx.doi.org/10.14770/jgsk.2020.56.4.515 (Korean with English abstract)

Example 1 (Seismic Analysis Code)

Examples to deconvolve the instrument response from raw data are below, using the Seismic Analysis Code (SAC, version sac-101.6a). 

1. A unit of output is 'm/s'

r $input
rtr
rmean
taper
TRANS FROM POLEZERO S $pzfile TO none FREQ 0.05 0.1 1.5 3.0
w vel.sac
q

2. A unit of output is 'm'

r $input
rtr
rmean
taper
TRANS FROM POLEZERO S $pzfile TO none FREQ 0.05 0.1 1.5 3.0
TRANS FROM VEL TO NONE
w disp.sac
q

3. A unit of output is 'm/s2'

r $input
rtr
rmean
taper
TRANS FROM POLEZERO S $pzfile TO none FREQ 0.05 0.1 1.5 3.0
TRANS FROM VEL TO ACC
w acc.sac
q

Note that the unit of the output is incoherent with the SAC header 'IDEP'.

Example 2 (StationXML with obspy)

#!/usr/bin/env python3

from obspy import read_inventory
from obspy import read

st = read('/path/to/waveforms')

inv = read_inventory('ksgn.xml')

st.detrend(type="linear")
st.detrend(type="demean")

pref = [0.1, 0.2, 15, 20]

#output = m/s**2
st2 = st.copy()
st2.remove_response(inventory=inv,output='ACC',pre_filt=pref,taper=True,zero_mean=False)
st2[0].write('acc.sac',format="SAC")
#...

#output = m/s
st2 = st.copy()
st2.remove_response(inventory=inv,output='VEL',pre_filt=pref,taper=True,zero_mean=False)
st2[0].write('vel.sac',format="SAC")
#...

#output = m
st2 = st.copy()
st2.remove_response(inventory=inv,output='DISP',pre_filt=pref,taper=True,zero_mean=False)
st2[0].write('disp.sac',format="SAC")
#...

Update notes

• [1.02] Typo in station name of KS.HA2B is modified (KS.HA2B -> HA2B)
• [1.01]
  • Files for the sensor JC-V100 are provided from Won-Young Kim (2 Dec 2021)
  • Station KS.HA2B is added
  • StationXML file is added (ksgn.xml)
• [1.00] The files for the KG network are made based on the logs until 12 April 2019 (personal communication with the Korea Institute of Geoscience and Mineral Resources). 

Pending works

Reflecting a change in logger of KS.DGY2

Others

Velocity seismometer https://doi.org/10.5281/zenodo.3700312
Accelerometer https://doi.org/10.5281/zenodo.3872436

License

This distribution follows the Creative Commons Attribution 4.0 International. You are free to modify and redistribute the files.