How to Make .pz files for use in localmag, gmew, etc.
Artist: Dr. Stephen R. McNutt
Obtain Dataless Seed Volume
Sources:
- IRIS BUD: http://www.iris.edu/bud_stuff/bud/bud_start.pl?BUDDIR=/budnas/virtualnets/ALL
- Geofon WebDC: http://webdc.eu/arclink/query?sesskey=ea8effeb
Install rdseed
The following instructions are specific to Ubuntu.
Install Dependencies:
sudo apt-get install cmake sudo apt-get install g++-multilib sudo updatedb
Download rdseed: http://www.iris.edu/forms/rdseed_request.htm
Unpack tar ball:
tar -xvf rdseed...tar.gz cd rdseed...
Compile it:
make
troubles? see README
Run rdseed
./<rdseed_directory>/rdseed -pf <dataless_file>
"An important fact to note about the SAC PZ output is that acceleration and velocity responses are converted to displacement to conform to the SAC convention." -rdseed manual (included with distribution)
NOTE: Poles and Zeros outputted in radians. Constant in counts/meters. Earthworm requires radians and counts/nanometers, respectively.
Enter file(s) entitled e.g., SAC_PZs_CX_HMBCX_BHZ2007.320.00.00.00.0000_2599.365.23.59.59.99999
* **********************************
* NETWORK (KNETWK): CX
* STATION (KSTNM): HMBCX
* LOCATION (KHOLE):
* CHANNEL (KCMPNM): BHZ
* CREATED : 2012-07-17T19:30:36
* START : 2007-11-16T00:00:00
* END : 2599-12-31T23:59:59
* DESCRIPTION : Plate Boundary Station HMBCX, Chile
* LATITUDE : -20.278220
* LONGITUDE : -69.887910
* ELEVATION : 1152.0
* DEPTH : 0.0
* DIP : 0.0
* AZIMUTH : 0.0
* SAMPLE RATE : 20.0
* INPUT UNIT : M
* OUTPUT UNIT : COUNTS
* INSTTYPE : Streckeisen STS-2/N
* INSTGAIN : 1.500000e+03 (M/S)
* COMMENT : N/A
* SENSITIVITY : 6.291450e+08 (M/S)
* A0 : 6.007700e+07
* **********************************
ZEROS 3
+0.000000e+00 +0.000000e+00
+0.000000e+00 +0.000000e+00
+0.000000e+00 +0.000000e+00
POLES 5
-3.700400e-02 +3.701600e-02
-3.700400e-02 -3.701600e-02
-2.513300e+02 +0.000000e+00
-1.310400e+02 -4.672900e+02
-1.310400e+02 +4.672900e+02
CONSTANT +3.779715e+16
Make pz file
Confirm the INPUT UNIT is M (meters) [NOTE: CONSTANT given in counts/meter] and then cp from ZEROS to CONSTANT and paste in your Earthworm response file directory (e.g., /opt/earthworm/osopdata/params/responses/) a file entitled, e.g., HMBCX_BHZ_CX.pz:
Now convert the numbers to nm by subtracting 9 units to the exponent (to convert counts/meters to counts/nanometers, 1 m = 1E9 nm so we must DIVIDE by 1E9):
counts/meter = meter/1E9 nanometer = counts/nanometer
ZEROS 3
+0.000000e+00 +0.000000e+00
+0.000000e+00 +0.000000e+00
+0.000000e+00 +0.000000e+00
POLES 5
-3.700400e-02 +3.701600e-02
-3.700400e-02 -3.701600e-02
-2.513300e+02 +0.000000e+00
-1.310400e+02 -4.672900e+02
-1.310400e+02 +4.672900e+02
CONSTANT +3.779715e+07
Additional Resources
- http://vps.isti.com/trac/ew/wiki/localmag, see section entitled "Response File Example"
- Scherbaum, F. (2007) Of Poles and Zeros: Fundamentals of Digital Seismology (Modern Approaches in Geophysics)
- Youtube "Z-Domain- Pole and Zero plots": http://www.youtube.com/watch?v=m5TP2uG_O2M&feature=relmfu
- Youtube "Z-Domain- 3d View": http://www.youtube.com/watch?v=mELFJIZ9cEY&feature=relmfu
- http://www.ee.up.ac.za/main/_media/en/undergrad/subjects/eli220/polezero.pdf
- Poles and Zeros GUI: http://www.osop.com.pa/?page_id=121
Attachments
- 739px-Calibration.jpg (86.6 KB) - added by branden 14 months ago.

