GetDPZ: Retrieve Pre-Calculated Discrete Instrument Response

In TDD: Time-Domain Deconvolution of Seismometer Response

Description

Discrete instrument responses for common seismometers and sample rates have been pre-calculated and included in this package. GetDPZ retrieves them.

Usage

GetDPZ(w, dt)

Arguments

w	Vector of indices of seismometers used (see Details)
dt	Sample intervals corresponding to instruments in w

Details

Seismometers are numbered as follows:

Broadband Seismometers: 1. Streckeisen STS-1 (360 s) 2. Trillium 240 (generation 1) 3. Trillium 240 (generation 2) 4. Guralp CMG-3T 5. Streckeisen STS-2 (generation 1) 6. Streckeisen STS-2 (generation 2) 7. Streckeisen STS-2 (generation 3) 8. Trillium 120 9. Compact Trillium

Intermediate Seismometers: 10. Trillium 40 11. Guralp CMG-3ESP 12. Guralp CMG-40T (30 s) 13. Streckeisen STS-1 (20 s)

Short-period Seismometers: 14: Guralp CMG-40T (1 s)

Value

List of instrument responses corresponding to instruments and sample intervals given in w and dt.

Note

The STS-2 and Trillium 240 come in multiple generations, each with a slightly different response. For the Trillium 240, serial numbers less than 400 belong to generation 1, while serial numbers greater than or equal to 400 are in generation 2. To determine which generation an STS-2 is, see http://www.iris.edu/NRL/sensors/streckeisen/streckeisen_sts2_sensors.htm .

Many short-period instruments were intentionally omitted because their responses depend on installation-specific parameters and are therefore not completely standardized. Given that discrete responses of short-period instruments can be calculated quickly, I consider the convenience of having pre-calculated responses for these instruments to not be worth the risk of the user selecting the wrong response and getting inaccurate results.

Author(s)

Jake Anderson

References

Sources for all instrument responses are given in the comments of GetPZ; the reference list is too long to include here.

See Also

GetPZ MakeDPZ

Examples

# responses for 3T sampling at 1 Hz and 40T (30 s) sampling at 40 Hz
DPZLIST = GetDPZ(c(4, 12), c(1, 0.025))

Example output

Loading required package: signal

Attaching package: 'signal'

The following objects are masked from 'package:stats':

    filter, poly

Loading required package: RSEIS
Loading required package: pracma

Attaching package: 'pracma'

The following objects are masked from 'package:RSEIS':

    detrend, hypot, logspace, peaks, trapz

The following objects are masked from 'package:signal':

    conv, ifft, interp1, pchip, polyval, roots


Attaching package: 'TDD'

The following object is masked from 'package:RSEIS':

    ReadInstr

TDD documentation built on May 2, 2019, 4:51 a.m.