AUTOPIC is a tool to automatically pick phases on events registered into the database. AUTOPICK can also be used for many events in connection with program AUTO, see below.

If an event file (S-file) has any readings, the AUTOPIC program will not reread in order to not destroy old picks. The automatic readings in the file are marked with an A after the weight column to indicate automatic pick. Each pick is evaluated by using the signal to noise ratio and an indication of the quality is given with the weight. The program will run on all waveform files given in an S-file. Each time the program runs, there is a file called autopic.out containing information about the run. If there are any 3-component stations, an azimuth will also be calculated, and the S-phase will be more reliable. The AUTOPIC program can also be used from EEV by typing Z (will run program AUTOPIC). When it is used from EEV, there is always an output in the S-file, which will be grouped at the bottom of the file, making it possible to compare manual and automatic readings. THE S-FILE MUST THEN BE EDITED MANUALLY IN ORDER TO REMOVE DOUBLE READINGS. The program requires an input parameter file in the working directory or DAT with the name AUTOPIC.INP. The program will first look in the working directory. The parameters in that file are explained below. NOTE: The file is formatted, data must be in columns exactly as shwown and no tabs must be used. The program uses a 4-pole filter running one way. This might result in phases being picked a bit late. However, it seems more accurate than the earlier version where the filter run both ways and picks were often far too early. The program is made mainly by Bent Ruud. For more information about how it works, see Ruud and Husebye [1992]; Ruud et al. [1988]. Description of parameters


Example of input file AUTOPIC.INP



A Program for automatic processing: Phase picking, location, magnitude and fault plane solution determination.

This program run the picking program AUTOPHASE first (optionally AUTOPIC), locates events with or without outlier removal with HYP and then determines magnitudes with AUTOMAG. Automatic amplitudes for fault plane solutions can be done optionally with AUTORATIO. Fault plane solutions can also be done optionally with FPFIT or HASH, note velocty is the hardwired option. The default is to do AUTOPHASE, HYP with outlier removal and AUTOMAG using default parameters. Individual steps, like doing magnitudes, can optionally be deselected and some non default parameters for AUTOMAG can be selected. All changes are made using arguments, see list below. Thus the program can also be used to do any one of the operations. However to do only e.g. magnitudes, a more rational choice would be to use AUTOMAG only. Similarly for doing only locations, it is more logical to use HYP. The input can be all SEISAN types: a file, a data base or an index file. Outlier removal by HYP can be deselected. An alternative to outlier removal by HYP is to use the HYP option for residual weighting. For more detail of the individual programs, see program descriptions elsewhere in SEISAN manual. There are always two output files: auto.out: all events with the final results with updated values of location and magnitude. auto.log: a summary of what has been done to each event. If input is from a data base or index file, the results are also written back to the data base, overwriting what was there from before. The ID line is not updated by AUTO. When picking automatically, a few bad picks can throw the solution off so for a local event, the distance might wrongly be very large. It is then important that the range used for distance weighting in HYP is large so the initial wrong location can be made and the outlier rejection gets a chance to eliminate bad picks. Doing fault plane solution is only recommend with very good data. It is then important to use at least the default values set for both minimum number of polarities and maximum az gap. Note the gap is not the gap as in a hypocenter solutions but the gap when polarities are plotted on the projection on the focal sphere.

See list of argument below

 Examples of argument use:

auto          : pick phases, locate and do magnitudes
auto l m      : only pick phases with AUTOPHASE
auto j l m    : only pick phases with AUTOPIC
auto i l       : only do magnitudes
auto s 30     : all defaults except that spectral window for magnitude is 30 s
auto l m i f  : only do fault plane solutions with FPFIT
auto l m i ar h f: automatic amplitudes and fault plane solution with HASH

Writing auto help will give all the possible arguments:

c:\seismo\PRO>auto help
 Write auto help to get list of arguments. If no
 arguments, all defaults are used. This is autophase,
 location with outlier rejection and ML and spectral Mw

  s xxx : do spectrum for Mw, xxx is window length, default
          is 20 s. if 0, no spectrum. Default is spectrum

  w xxx : amplitude for Ml, xxx is window length
          if zero, no amplitude. Default is do amplitude
          default window length is 50 s

  n or e: use N or E component for magnitude, respectively.
          Default is Z

  p     : use P for Mw, default is S

  l     : do not locate, def. is to locate with rejection
          if magnitudes are done, location is also done
          after determination of amps and spectra in order
          to update magnitudes on header line

  r     : do not remove outliers when locating

  m     : no magnitude, default is to do magnitude

  i     : no phase pick, default is to do phase pick

  j     : if phase pick, use autopic, default is autophase

  ar    : do autoratio, default is not to do

  af  xx xx  : filter for autoratio, default is 2-4 Hz

  ad  xx   : max distance for autoratio, def 100km

  at  xx  : time window for autoratio, def 2 s

  ag  xx   : ground motion for autoratio, def 0 for none

  az    : if given, use z for s in autoratio

  f     : do fault plane solution with fpfit, def.  is not

  h     : do fault plane solution with hash, def.  is not

  n     : minimum number of polarities for fps, default 10

  g     : maximum gap for fps, default 180

When using amplitude ratios for fps, velecity is the harwired default.

Peter Voss 2022-08-29T10:10Z | Please contact pv@geus.dk regarding the webpage, it is hosted by The Geological Survey of Denmark and Greenland (GEUS) | xxxxxxxxxxxyyyyyyyyyyyy