next up previous contents index
Next: 8.11 Data manipulation commands Up: 8. Trace plotting, phase Previous: 8.9 Instrument correction and   Contents   Index

8.10 Determine azimuth of arrival (3 comp or array) and component rotation

Azimuth of arrival from 3-component stations, h(Azim)

If a 3 component station is available, the azimuth of arrival can be determined using the method developed by Roberts et al. [1989]. Display any of the 3 components and press h (Azim). Then select a zoom window around the P-arrival of a few secs duration for the analysis. The 3 components will now be displayed below in order Z, N and E and the calculated azimuth, apparent velocity and correlation will be displayed at the bottom line. In order to check the stability of the estimate, try different windows and filters. Often, a filter must be used to get reliable results. The displayed azimuth and apparent velocity is only saved in the S-file when an associated phase is picked. THAT PHASE MUST BE PICKED ON THE SINGLE UPPER TRACE SEEN ON THE SAME SCREEN. If there is none, use I or E. The velocity estimate is not very reliable and is dependent on the local velocities. In order to calculate the apparent velocity, the P-velocity of the top layer must be given. The default value is 5.0 km/sec, but another value can be set in the MULPLT.DEF file. To get a good estimate, the correlation coefficient should be as high as possible and positive. The quality of the obtained azimuth can be tested by locating the event with the calculated azimuth weighted out and observe the azimuth residual. Figure 8.10 shows an example.

Azimuth and apparent velocity from array data, FK analysis f(FK)

Using this command, the traces seen on the screen will be put into the FK program and an FK plot will be displayed. The azimuth and apparent velocity with the highest correlation is selected by the program, however any other value can be manually selected. The values will ONLY enter the S-file if associated with a phase in the same way as amplitudes are picked. For more details, see section 40.

Rotated seismograms

Option U(Rotat) will rotate the horizontal components for the next plot if the two horizontal components are available. The rotation will display the radial component instead of the N-component and the transverse component instead of the E-component. The back-azimuth used is displayed above the trace. All channels will be displayed rotated until u(Rotat) is pressed again This means that phases can be picked and spectra made with the rotated channel. When picking phases on rotated signals, these will appear in the S-file with components R or T instead of N and E respectively. This also means that only if the rotated signals are shown, will the phases read on rotated channels appear on the plot. The station back-azimuth is obtained in the following way: If a hypocenter is given in the header line, the angles are calculated using the current STATIONx.HYP file. If no hypocenter is available, the angle will be read from the S-file under column observed azimuth (47-51) (if not blank) and the azimuth residual will be added. This option permits the user to first determine the azimuth with the 3-component option and then rotate the signals with the determined azimuth. Finally, if no observed azimuth is available, the event to station azimuth + 180 deg. will be used if available (column 77-79). If no back-azimuth can be found, no rotation is done and an angle of 999 deg. is displayed. If in single trace mode and choosing the 3-component option AND the rotate option, the user will be prompted for a rotation angle and the rotated channels will be shown in the usual 3-component plot, however, the azimuth determined is done with the unrotated channels.

PROBLEM: In general, the R-channel will use the response of the N-channel and the T-channel will use the response of the E-channel so for instrument response removal to be correct, the 2 channels must have the same response curve.


next up previous contents index
Next: 8.11 Data manipulation commands Up: 8. Trace plotting, phase Previous: 8.9 Instrument correction and   Contents   Index
Peter Voss : Mon Feb 27 10:16:12 UTC 2017