Seismogram deconvolution in SAC Good day! I have a problem replicating a result. Let us consider the top plot here:



On the plot, velocity is in mm/s. I think the problem is to remove instrument response correctly. I downloaded seismogram and pole-zero instrument response. I want to obtain velocity data. The instrument response is already given in velocity. I was able to get synthetic seismograms as shown in the figure, but I have a problem working with real data.

I need a signal to be filtered in 0.001–0.004 Hz

Using IRIS SAC software, I do:

SAC&gt; r KYK_BHZ
SAC&gt; rmean
SAC&gt; rtrend
SAC&gt; taper
SAC&gt; trans from polezero subtype KYK_LHZ.zp
SAC&gt; bp p 2 n 4 c 0.001 0.004



I do not set FREQlimits for transfer function at all here, but I tried to do it like 0.0009 0.001 0.004 0.005 or 0.0009 0.001 1 2, but it also doesn't work.
My synthetic seismogram is 1 Hz data, so I don't do any  resampling.
I tried to do trans from polezero subtype KYK_LHZ.zp to vel, but it also doesn't work.


Both the data and zero-pole files are here.

Thanks!
I typed in the commands as shown using your data in SAC. I've uploaded my result.  The plot window shows the result of the commands run in the highlighted red area. Is this what you are trying to do? The plot xlim is set from 1750 s to 2700 s.