Hi,
I am working with data from a Kinemetrics FBA-11 accelerometer network. I have used detrend(‘demean’) to remove drift but it seems that there are still artifacts from the instrument. I want to use instrument correction in obspy and I am wondering if there is a central database that provides the needed PAZ etc?
Many thanks
John:
Happy New Year.
There is PDCC, IRIS’s Java-based metadata management tool:
http://www.iris.edu/dms/nodes/dmc/manuals/pdcc_3-0/
http://www.iris.edu/forms/pdcc_request.htm
This uses a central database called the Nominal Response Library.
or visit gempa GmbH’s SMP web-based tool:
https://seismo.gempa.de/smp
Kind Regards from New York,
You can also directly access the IRIS NRL (nominal response library) through this URL: http://www.iris.edu/NRL/
Thank you both very much and a happy new year!
I've used the resp file you mentioned, found at
and it looks much better: I'm interested in displacement (FBA-11
accelerometers on 6th and 8th floor of a building). If I detrend the
stream and integrate twice (without correction) I get
If I now use instrument correction:
seedresp = {'filename':
'/home/hd/Dropbox/projects/current/RESP.XX.NS235..BNZ.FBA11.DC_50.2_5V.1G',#RESP
filename
'date': UTCDateTime("2007-08-09T07:57:49Z"),
'units': 'DIS' # Units to return response in ('DIS', 'VEL' or ACC)}
st_orig.simulate(paz_remove=None,seedresp=seedresp)
the new plot I get is 6th and 8th floor with resp file correction - Imgur which looks much
better. However, this is still not correct. There is an animation
available for the this data stream
http://goldengate.ce.caltech.edu/~kohler/Factor_data/2007/221/2007_221.avi
and it moves much more before it returns to its starting position than
one sees in 6th and 8th floor with resp file correction - Imgur
Many thanks for your help!
Hi John,
looks like blown up low frequency noise to me. Try lowering the water
level (e.g. 60, 30 or even down to 10) and/or pre-filtering in frequency
domain (to the part of the spectrum that is of interest, use option
"pre_filt").
Set the parameters for your Stream.simulate() call, you can look up the
options at the corresponding low level function here:
http://docs.obspy.org/master/packages/autogen/obspy.signal.invsim.seisSim.html?highlight=seissim#obspy.signal.invsim.seisSim
best,
Tobias
Thx!!!
I am now prefiltering:
fl1 = 0
fl2 = 0.45
fl3 = 5.
fl4 = 500.
and use water_level=60.
What confuses is me is that calling Stream.simulate() with seedresp
set to units: 'VEL' and subsequently calling Stream.integrate() to get
the DIStance, gives something completely different than calling
Stream.simulate() with seedresp set to units: 'DIS'.
with units 'VEL': Corrected instrument for VEL - Imgur
with units 'VEL' after Stream.integrate(): Corrected instrument for VEL, then Stream.integrate() - Imgur
with units 'DIS': Corrected instrument for DIS - Imgur
I can't say what is going wrong for you (with just pictures) but i can
say that it works in principle when using correct response information
and proper options and pre/post-processing:
#!/usr/bin/env python
from obspy import read
import matplotlib.pyplot as plt
kwargs = dict(water_level=60, pre_filt=(1, 3, 20, 30))
tr = read()[0]
tr.simulate(seedresp={"filename": "RESP.BW.RJOB..EHZ", "units": "DIS"},
**kwargs)
tr.filter(type="highpass", freq=0.1)
plt.plot(tr.data[2:])
tr = read()[0]
tr.simulate(seedresp={"filename": "RESP.BW.RJOB..EHZ", "units": "VEL"},
**kwargs)
tr.filter(type="highpass", freq=0.1)
tr.integrate()
plt.plot(tr.data[1:])
tr = read()[0]
tr.simulate(seedresp={"filename": "RESP.BW.RJOB..EHZ", "units": "ACC"},
**kwargs)
tr.filter(type="highpass", freq=0.1)
tr.integrate()
tr.integrate()
plt.plot(tr.data)
plt.show()
The necessary response can be prepared as RESP like:
wget
https://raw.github.com/obspy/obspy/master/obspy/xseed/tests/data/dataless.seed.BW_RJOB
obspy-dataless2resp dataless.seed.BW_RJOB
I assume there is something wrong with your processing chain and the
response information you put together (e.g. digitizer?). You should end
up in m/s or m and your numbers are a bit high for that.
best,
Tobias
That was extremely helpful, thx! I now end up with (nearly) the same plot after correcting either directly for DIS or alternatively for ACC and then integrating twice (http://pastebin.com/esniGcGY). However, I’m still dubious about the actual numbers: the data is from a Kinemetrics FBA-11 acclerometer on the 6th floor of a building, hence they are way too high to make sense in meters?
My goal is to simply reproduce the animation
http://goldengate.ce.caltech.edu/~kohler/Factor_data/2007/221/2007_221.avi
with obspy. I’m using the resp file from http://www.iris.edu/NRL/sensors/kinemetrics/kinemetrics_fba11_sensors.htm and the
data http://goldengate.ce.caltech.edu/~kohler/Factor_data/2007/221/ which is from an earthquake with ML=4.6 that was 30km away so I don’t see where things are going wrong.
Hello John,
Just a wild thought that you might have already dealt with, but did you include the sensitivity of the digitizer in your resp file (Volts/count)? I think the NRL response library assumes a digitizer with a unity sensitivity(1 Volt/count). You could always find the digitizer you are using in the NRL and include that piece into stage 2 of your resp file.
Best,
Adam
Hi John,
You also need to know the gain in counts/V of the digitizer and account for that. Right now you’re only correcting for the sensor..
-Renate
Thx! That must be it, the NRL resp file was indeed using 1V/count. However, all http://factor.gps.caltech.edu/instrumentation says is that they are Quanterra 4128s digitizers and according to the NRL response library http://www.iris.edu/NRL/dataloggers/quanterra/q412x/quanterra_Q412x_dataloggers.htm these can be user configured. Is there any way I can find out what setting was used from the SAC files on http://goldengate.ce.caltech.edu/~kohler/Factor_data/2007/221/?
Another basic question: incorporating the digitizer gain would only affect the amplitude however the resp files about the digitzers at NRL seem to have much more information about response coefficients. Wouldn’t I have to include this information as well?
Hi John,
The only way to figure out digitizer settings without meta-data (such as a dataless SEED or stationXML file) is to contact the original deployer and ask, however, you are in luck. These data are also archived at the IRIS DMC so you can get the full meta-data there. Go to this link http://www.iris.edu/mda/FA, select any of the stations, then select “view RESP”…
-Renate
Thanks Renate, that’s brilliant! I’m now using the RESP file from http://www.iris.edu/mda/FA but I’m getting a warning from Obspy:
WARNING: FIR normalized: sum[coef]=9.781111E-01; FA FABD 6N HNE
While the shape looks good, I’m dubious about the amplitude numbers…the maximum change of position during the event on the 6th floor would have only been 0.6mm?
http://imgur.com/RFclca6
(The first row is uncorrected data, the 2nd instrument corrected for ACC with RESP and the third instrument corrected for DIS with RESP)
The (obs)python code I used is http://pastebin.com/QDPG7r0y