1. Introduction

Unpack the section_03.tgz with the command

gunzip -c section_03.tgz | tar xvf - 

This will create the directory Section_03 with subdirectories Section_03/Forward and Section_03/Processing.

2. Data set in absolute time

Rather that viewing the data set in term of what is seen in real time, consider it as for all time, or essentially, previous acquired history that must be processed. A segment of that data stream is show in the next figure. For simplicity, it is assumed that the trial epicenter is the true epicenter.

3. Processing time windows using a group velocity window

In this set an origin time is assumed and then a time window is computed based on the travel time of the minimum and maximum group velocities to the station. In first example, the assumes origin time is offset -5, 0 or +5 minutes about the true origin time of the signal. The consideration of 3 windows is used to illustrate the procedure.

The next figure displays the sections of data used for the data compression.

-5 minutes 0 minutes +5 minutes

It is obvious that if the window diverges significantly from the range of the actualy data, then the signal compression and stacking will not yield significant values.

The next figure illustrates this for offsets of -10, -5, 0, 5 and 10 minutes about the true origin time. The figures are annotated with the assumed origin time and the maximum amplitude of the envelope. The assumed origin time is indicated by the red vertical bar on each plot.

Stacked compressed signals for different assumed origin times

The first thing to notice that one almost always sees a peak at the assumed origin time. this is because other windows have signal that is compressed, although not perfectly.

4. Processing time windows without a group velocity window

Here we do nto apply a group velocity window to the traces. Instead a moving window is slected, and the inverse dispersion oeprator is applied to the entire window. The windows for -1, 0 and +10 minutes about the true origin time are displayed in the next figure.

-10 minutes 0 minutes +10 minutes

Note that these plots start at 00:51:59.566, 01:01:59.566 and 01:11:59.566, respectively. This is illustrated in the next figure which show the synthetics trace for IUANTO plotted in real time for different start times of the cut window, as indicated.

The result of the compression at -10, -5, 0, +5 and +10 minute offsets are shown in the next figure.

For the windows considered, the compression and stacking results are virtually the same. There is some soouce of signal at roughly 01:02. The compression projects gthe observed signal back in absolute time, and since the cut signals ahd the same information, we are not surprised with this plot.

Signals commensurate with the assumed dispersion are mapped back to the origin time. Other signals are spread out and offset. In this figure the initial wavefield, think P waves, propagate faster than the assumed phase velocities. thus they are mapped and spread out later in time.

Final plot compares the two stacks. the red trace is the result of applying the group velocity window, which the blue trace uses the entire trace. Near the true origin time, the stack amplitudes are about the same. The group velocity windowed processing does not yield as sharp a pulse. The use of a wider window and a different taper will be investigaated.