This electronic supplement embodies a Matlab™ implementation of our algorithms and guidance as revised two years after publication. It also provides ANSS guidance, Earth-noise models, and earthquake examples for co-plot.
The original bundle published in 2010 contained serious typographical errors that cause the two- and three-instrument (Holcomb and Sleeman) methods to produce incorrect results (too-high noise levels in the presence of ambient Earth noise). The revised bundle (2012) corrects these errors, adds more test examples to verify and demonstrate a wider range of applications, and subsumes numerous minor corrections, clarifications, and general cleanup.
Both the old and new scripts require Matlab™ version R2007a or later; http://www.mathworks.com/.
Key filenames, function names, and terse descriptions are provided in Table S1 (Unix™ annotation style). The primary routine is GST_OpRngDiagrams(), for which users should generate driving routines similar to the examples provided for the test data, as files named Exec_*.m. These driving routines prepare the data, set up response information, and so forth, and then call routine GST_OpRngDiagrams(), finally saving the results as desired (we advise using the comprehensive saves shown in the examples).
In practice, users will modify at least the directory and file names in Exec_*.m as well as instrument clipping and response information for their own applications. Additional user guidance is given in file GST_OpRngDiagrams_Instructions.txt and in the extensive comments at the head of GST_OpRngDiagrams().
In operation, GST_OpRngDiagrams() generates two PostScript™ plot files, an ORD and an rPSD like those illustrated in Figures 2 and 3 of the paper. At the user's request, it also plots the input time series. Press “s” for “save” when the cursor appears in the Matlab™ plot (or its legend). GST_OpRngDiagrams() then returns the various vectors computed for those plots, which the user may save, use, or ignore.
Although the plots generated by our scripts are usable as created, we recommend modest graphical editing, for example with Adobe™ Illustrator™, to make them more presentable. In particular, the user may wish to move or modify plot legends, which commonly obscure some results, and correct the square-root symbol in rPSD plots, which can suffer change due to font incompatibilities among various software systems and platforms. Both such edits were performed in preparing Figures 2 and 3, for example (using Illustrator™ on a Macintosh™).
Our Matlab™-based system is very likely to evolve as any additional bugs are corrected and features are added, so the electronic supplement may become obsolete. We will maintain an updated version of that content at a USGS anonymous ftp() site, which should be used in preference to this e-Supplement except immediately after distribution in early 2012. The USGS-supported version of the package can be found at:
<ftp://anonymous:@aslftp.cr.usgs.gov/pub/users/jrevans/OperatingRange/>
Finally, if users discover missing components, bugs, out-of-date documentation, or other problems, or have difficulty using the scripts they are asked to contact the first author at jrevans@usgs.gov to obtain assistance and provide feedback.
| Filename(s) | Function Name | Comment |
|---|---|---|
| GST_OpRngDiagrams.m | GST_OpRngDiagrams() | Primary routine |
| Exec_GST_OpRngDiagrams_ … … Example_1inst_*.m | main() | User's routines for single- sensor (or DAU) method |
| Exec_GST_OpRngDiagrams_ … … Example_2inst_*.m | User's routines for two-sensor method | |
| Exec_GST_OpRngDiagrams_ … … Example_3inst_*.m | User's routines for three-sensor method | |
| ANSS_Guidance*.m | ANSS_Guidance() | ANSS (and author) instrument guidelines |
| GST_NormalizedWelch … … _[1-3]inst.m | GST_NormalizedWelch … … _[1-3]inst() | Welch-method PSD computation for 1-, 2-, or 3-instrument methods |
| Other ../*/*.m | GST_HalfOctavePower() through unepoch() | Subroutines used by GST_OpRngDiagrams() |
| [12]*.mat | input data files | Used by Exec*.m |
| NoiseSpectra/*.mat | Earth-noise vectors | PSD, rPSD, and half- octave rms values of various noise models |
| GST_OpRngDiagrams … … _Instructions.txt | text | Instructions for users (cf. GST_OpRngDiagrams.m) |
[ Back ]