7.1.10 Preventing a local minimum

For some data, the iterations will stop in a local minimum. This can partly be checked with the DRMS output in the print.out file. The DRMS output gives the RMS in 6 points around the hypocenter at distance DRMS. DRMS is set with RESET TEST (13). In order to check for a local minimum, a more extensive grids search can be done. This is done by using RESET TEST 119-121. Setting these parameters, a grid search can be done in 3D around the hypocenter in a user selected grid size. The grid search is done after the last iteration. If a point is found with a lower RMS than given by the iterative process, a new iterative process is started with that point as start location. The grid search can be particularly useful for local earthquakes, however, in the current implementation it is also done for distant earthquakes using the same parameters. Note that using a large grid might significantly slow down the processing time.

When using the grid search option for start location, an extra ouput is given in the print.out file:

Grid rms smaller than solution rms 1.953 3.391

Hypocentral dist. (km) from epicenter to grid point min. 28.3

Depth change 10.00

Iterate again with new start location: 62.303 5.760 15.0

This information gives how far away the minimum is and what the RMS is.

The most likely local minimum is happening for the depth. The start depth (on control line in STATION0.HYP) is therefore important and should be what is expected for the area. If there is a large variation in depths, like in a subduction zone, and the network is small, it might be an advantage to use several start depths. This can be done in two ways:

1: Set a range of start depths on the control line. In this case a complete location is done with each start depth. Setting many depth will increase the processing time proportionally.

2: Set RESET TEST(96) to 1.0. With this parameter set, the program will, once a location is found, look for the depth with the smallest RMS using the loction found. A new location is now made with the a start depth location at the minimum rms depth. This procedure is fast since the mimumum RMS depth is found with a fixed epicenter and no new location is found, however, since a new location is not found, it will not always work as well as option 1. The procedure for this option is to start with a depth 3 times the current depth (set to at least 40km) and then calculate rms in 40 steps from that depth to the surface. In SEISAN before version 12.0, the procedure was to start with current depth and the go in steps of 10 to the surface but this did not work well with shallow earthquakes.

Peter Voss : Fri Nov 12 10:33:10 UTC 2021