AUTOSIG is a program to perform some automatic processing. The program includes routines for P-phase picking, determination of signal duration, amplitude determination, determination of spectral parameters (Ottemöller and Havskov, 2003) and determination of distance type (local, teleseismic, noise). The program can still do with improvement. The input to the program can be either a parametric Nordic file (with one or several events) or waveform files. In both these cases, the output is written to the autosig.out file. Additional output files are autosig.trace and autosig.err, which will help to find potential problems. Alternatively, the program can also be started from EEV; the output is then directly written to the S-file. The input parameters are defined in the file autosig.par, which is located either in the DAT or the working directory.
Following are descriptions of the automatic processing routines:
The phase picking is based on changes in the STA/LTA ratio. A band pass filter can be specified. The routine gives options to enhance the changes in the signal before computing the STA/LTA ratio. It is recommended to first remove the linear trend and then to compute the characteristic function which is given by y**2+k*(dy/dt)**2, which enhances changes in both amplitude and frequency content. Then the STA/LTA ratio is computed to detect changes in the signal. The routine can also compute the squared STA/LTA. When a change is detected (STA/LTA ratio above trigger level), it is tested whether the signal spectral amplitudes are significantly higher (factor of 2 in amplitude) than the pre-signal noise spectral amplitudes. This is done to avoid triggering on spikes.
The signal duration is determined by comparing the signal amplitudes with the amplitudes of the pre-signal noise. The duration is determined by the point from which the ratio of these amplitudes is lower than a given value. A filter is applied if specified in the parameter file.
Routine finds maximum amplitude between two peaks.
The routine computes the displacement amplitude spectrum for P or S waves (see section 8.12) and, using either a converging grid search or a genetic algorithm determines the seismic moment and the corner frequency by minimizing the difference between observed and synthetic source spectra. The frequency band is determined by comparison with the pre-signal spectrum. The grid search is generally more cost effective and produces better results. The method is described in Ottemöller and Havskov (2003). The displacement spectrum is corrected for geometrical spreading and attenuation (both along the travel path and near surface). Therefore, the hypocentral distance has to be known. The time domain window for extracting the data from the trace can be given by either a group velocity (Vg=distance/travel time) window or a fixed window in seconds around the phase pick.
The routine determines whether the signal is from a local or teleseismic event, or noise. If signal spectral amplitudes are not significantly higher than pre-signal noise amplitudes, it is assumed that the signal is noise. Otherwise the amplitudes at two selected frequencies given by `DIST FREQ SELECT' are compared, the rules are (f1 f2):
Spec signal amp(f1) - Spec noise amp(f1)
Spec signal amp(f2) - Spec noise amp(f2): teleseismic
Spec signal amp(f2) Spec noise amp(f2): local
There are a few command line options that can be used to run autosig in non-interactive mode, syntax is
autosig -infile filename [-spec on/off -phase on/off -clear on/off]
-spec on/off: determine spectral parameters if option given
-clear on/off: remove phases from input S-file before start if option given
-phase on/off: detect phases if option given
-infile file : give name of input file, either S-file or waveform file
-help: get help
Note: When running the program the first time and the hypocenter location is not known, the determination of spectral parameters is not done. To run the determination of spectral parameters, the hypocenter location has to be given in the S-file.
The meaning of most parameters in the parameter file is clear from the keyword. The spectral parameters are as described in the MULPLT section. Other parameters that need explanation are:
AUTO PHASE, AUTO SPECTRUM and AUTO AMPLITUDE: Logical flag to activate phase picking, spectral analysis and amplitude reading, respectively (1. for true)
GA POPULATION SIZE: Number of elements in the population, used only if SEARCH ALGORITHM is 1.
GA GENERATIONS: Number of generations in one run, used only if SEARCH ALGORITHM is 1.
Note: Increasing GA POPULATION SIZE and GA GENERATIONS will increase the computation time.
GRID NLOOP: Number of loops in converging grid search for spectral parameters, used only if SEARCH ALGORITHM is 2. Resolution increases with every loop.
NGRID FREQUENCY: Number of grid points in search for corner frequency, used only if SEARCH ALGORITHM is 2.
NGRID SPECTRAL AMP: Number of grid points in search for spectral amplitude, used only if SEARCH ALGORITHM is 2.
NORM: Norm for computation of residuals in spectral fitting can be set, however, tests show that 1 or 2 produce the same result, and generally default of 1 can be used.
SEARCH ALGORITHM: Defines whether genetic algorithm (1) or converging grid search (2) should be used. Converging grid search is recommended.
SELECT PHASE: Defines, which phase to use for spectral analysis, choices are: 0 for P by AUTOSIG, 1 for computed P arrival for given location, 2 for computed S arrival, 3 for P from s-file, 4 for S from s-file or 5 for S or P from s-file.
SEPCTRUM F LOW: Lower limit of frequency band to be used.
SPECDURATION CHOICE: The time window for computation of the spectrum can be given either as a time window starting from the phase onset (0.) or can be defined by a group velocity window (1.).
SPECTRUM P LENGTH: Duration in seconds of signal starting from P arrival.
SPECTRUM S LENGTH: Duration in seconds of signal starting from S arrival.
SPECTRUM PRE LENGTH: Duration in seconds of signal to be included prior to phase arrival.
GROUP VEL WINDOW P: Range of group velocities defining time window to be used for P spectrum. Time window is given by (distance/group velocity)
GROUP VEL WINDOW S: Range of group velocities defining time window to be used for S spectrum. Time window is given by (distance/group velocity)
STALTA NREC/REC: There are two STA/LTA algorithms, recursive (0.) and non-recursive (1.).
STATION LINE: One line with processing parameters for phase detection is given for each channel. The parameters are (also see example below):
Example of the parameter file autosig.par:
# KEYWORD............Comments.............Par 1.....Par 2 # # spectral parameters # SPECTRAL S-Q0 Q0 440.0 SPECTRAL P-Q0 Q0 85.0 SPECTRAL S-QALPHA Q = Q0**Qalpha 0.70 SPECTRAL P-QALPHA Q = Q0**Qalpha 0.70 SPECTRAL KAPPA 0.02 SPECTRAL P-VELOCITY P velocity 6.2 SPECTRAL S-VELOCITY S velocity 3.6 SPECTRAL DENSITY Density 2.8 # # auto signal processing parameters # REMOVE MEAN 1. for true 0. REMOVE TREND 1. for true 1. CHAR FUNCT 1. for true 1. K IN CHAR FUNCT K IN X=Y^2+K*Y'^2. 3. STALTA NREC/REC rec 0./ non-rec 1. 1. STALTA SQUARE/ABS square 1. 1. AUTOCODA SQUARE/ABS square 1. 0. AUTO PHASE 1. for true 0. only if no phase AUTO SPECTRUM 1. for true 1. AUTO AMPLITUDE 1. for true 1. AUTO LOCATE 1. for true 0. NORM 1. SEARCH ALGORITHM 1.=GA 2.=GRID 2. # # window selection # SPECTRUM P LENGTH in seconds 5. SPECTRUM S LENGTH in seconds 2. SPECTRUM PRE LENGTH in seconds .5 GROUP VEL WINDOW P 5.0 6.5 GROUP VEL WINDOW S 2.7 3.7 SPECDURATION CHOICE 0:SPEC. P/S LEN. 1. 1:GROUP VEL W. P/S # # select phase # SELECT PHASE 0:auto P 2. 1:synth P 2:synth S 3:P from file only 4:S from file only 5:S or P from file --- for 3-5, if no phase and AUTO PHASE is 1. use auto phase pick --- DIST FREQ SELECT 1. 10. # # parameters used in the genetic algorithm searching for spectral parameters # GA POPULATION SIZE 50. GA GENERATIONS 250. # # grid dimension in case of grid search # GRID NLOOP 5. NGRID FREQUENCY 100. NGRID SPECTRAL AMP 100. # # low filter limit to use for auto spectrum SEPCTRUM F LOW .05 # # parameters controlling output # CREATE WAVEOUT 1. SPEC OVERWRITE 0. PHASE OVERWRITE 0. # # station parameters # # STAT- COMP -sta-- -lta-- -ratio mincod -dtrle fill filh STATION MOL S Z 3.0 20.0 10.0 2.5 1.5 5.0 10.0 STATION MOL B Z 3.0 20.0 10.0 2.5 1.5 5.0 10.0 STATION HYA S Z 3.0 20.0 10.0 2.5 1.5 5.0 10.0 STATION LSA L Z 3.0 20.0 10.0 10.0 9999.0 00.1 10.1 STATION CHTO L Z 3.0 20.0 10.0 10.0 9999.0 00.1 10.1 STATION XAN L Z 3.0 20.0 10.0 10.0 9999.0 00.1 10.1