.. _lambda-sec-analysis: Interactive Array Analysis ========================== |vappname| implements :ref:`different methods ` to detect seismic or infrasound phases interactively by an operator. As a result slowness, backazimuth and arrival time of the detected phase or time and location of an origin are determined. If possible phase hints are added. The set of configuration parameters is available in the :ref:`lambda description `. Read the section :ref:`lambda-sec_configuration` for more details. .. hint :: Multiple processing and parameter widgets are available in |vappname|. Read the section on :ref:`lambda-sec_widget-arrangement` for learning how to select and arrange these widgets interactively. Usual array analysis in |vappname| consists of #. Selection of * :ref:`array `, * :ref:`data and event` if available, #. Data pre-processing including filtering as described in section :ref:`Interactive data processing `. #. Application of :ref:`array methods ` for signal detection to one or more arrays or 3C stations as described in the following sections: * :ref:`single time windows ` * :ref:`sliding time windows ` #. Joint analysis of the detections in the :ref:`cart ` of |vappname|: phase selection and relocation, #. Selection and committing the solutions in the :ref:`cart `. .. _lambda-sec-data-singletw: Single time-window detection ---------------------------- The data can be processed in s single time window which is to be selected specifically. Within the time window one of the methods described below can be applied typically results in a single detection which can be selected, analysed and provided to further processing. Besides classical beam methods, |vappname| also provides conventional phase picking similar to :program:`scolv` but also 3-component (3C) polarization analysis. .. note :: The methods for single time-window analysis are started by pressing buttons in the upper trace window of the Waveform widget. .. _lambda-sec-data-singletw-interactive: Interactive beamforming ~~~~~~~~~~~~~~~~~~~~~~~ Values for slowness and backazimuth can be provided manually. Then the beam is be formed from the loaded waveforms. The beam is shown in the upper zoom trace widget along with the traces of the stations. #. Open Beam parameters window and Timeshift plot: View -> Windows, #. Select frequency range in ATF Settings window, press *Apply*, #. Toggle data filtering: press *Filter* button, #. Adjust backazimuth and slowness in Beam parameters window. The corresponding time shifts can be viewed in the window named *Timeshift plot*. #. Optionally invalidate beam. .. _lambda-fig-interactivebeams: .. figure:: /base/media/lambda/lambda_array_interactivebeams.png :width: 12cm :align: center Selection of analysis methods and parameters. .. _lambda-sec-data-singletw-picking: Phase picking ~~~~~~~~~~~~~ You may interactively pick phase arrivals on single 1C or 3C stations. The arrival times may be added to the :ref:`cart ` together with previously determined backazimuth and slowness. #. Empty beam parameters. Otherwise they will be added to the picks! #. Select a trace, #. Select a phase type, #. Make a pick: click with mouse on trace, #. Optionally: set uncertainty by dragging the triangles, #. Add to :ref:`cart ` or remove. .. _lambda-fig-picker: .. figure:: /base/media/lambda/lambda-singletw-picking.png :width: 12cm :align: center .. _lambda-sec-data-singletw-3c: 3C-polarization analysis ~~~~~~~~~~~~~~~~~~~~~~~~ Instead of using multi-station arrays, the arrival time, backazimuth and incidence angle can be measured on 3-component (3C) station. Given an effective velocity :math:`v` of the uppermost Earth layer, the horizontal slowness :math:`p_H` is derived from the incidence angle :math:`i`: .. math:: :label: lambda-slowness p_H = \frac{\sin i}{v}. #. Select a 3C station - may be the same as for picking, #. Activate Detail view, #. Adjust the time window, #. Click on *POL*, #. Adjust values, #. Optional: Export ascii file for custom processing, #. Press *Apply*, #. Add the arrival time of the phase, #. Add to :ref:`cart ` or remove. .. _lambda-fig-3c: .. figure:: /base/media/lambda/lambda-singletw-3cpol.png :width: 12cm :align: center 3C polarization analysis. .. _lambda-sec-data-singletw-fk: Slowness - backazimuth analysis ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Slowness - backazimuth analysis applied to multiply stations is available for single time-window analysis. The methods :ref:`lambda-sec-methods-fk` and :ref:`lambda-sec-methods-beamp` are available. #. Select many or all traces: Press :kbd:`A`. #. Set the time window large enough / slowness small enough. #. Press OK to confirm or #. Press FK/BP and adjust values, then start again #. Make an additional pick on the beam (assigned to the reference station) #. Add to :ref:`cart ` or remove. .. _lambda-fig-single-tw-fk: .. figure:: /base/media/lambda/lambda-singletw-fk.png :width: 12cm :align: center Slowness - backazimuth analysis in single time window. PMCC ~~~~ Select PMCC for single time-window analysis. Currently not implemented. .. _lambda-sec-data-slidingtw: Sliding time-window detection ----------------------------- The data can be processed continuously by shifting time windows over a longer time series. Within a time window one of the methods described below can be applied. Sliding-time window processing typically results in multiple detections which can be individually selected, analysed and provided to further processing. #. Select the data, #. Adjust and apply the processing parameters, #. Select travel time parameters, e.g. for :ref:`backprojection ` or phase hints, #. Start a method. The methods for sliding time-window analysis are started with buttons in the lower trace window of the Waveform widget. .. _lambda-fig-slidingtw: .. figure:: /base/media/lambda/lambda-sliding-traces.png :width: 12cm :align: center Prepare data and select a method. .. _lambda-sec-analysis-manual: Inversion of manual picks ~~~~~~~~~~~~~~~~~~~~~~~~~ The arrival times of manually picked phases can be inverted for slowness, backazimuth and arrival time at the reference station. The :ref:`plane-wave fitting ` method is applied. .. _lambda-fig-manual-inversion: .. figure:: /base/media/lambda/lambda-manual.png :width: 12cm :align: center Phase picking and inversion in plane-wave fitting window. Derived values are shown in the upper right. For comparison the contour plot indicates the semblance prediction based on the derived parameters. The maximum should be close the derived value. #. Start the manual data processing by clicking on the *MPWF* button for *manual plane-wave fitting* in the Waveform widget, #. Prepare the data: * Zoom in to the desired time window. The upper pick window can be controlled separately from the lower overview window, * Apply filtering by clicking on the filter button, * Press :kbd:`m` to toggle the stacked view, * View spectrograms for enhancing visibility, #. Start picking by pressing 1 or by clicking on the pick icon: * Minima or maxima can be picked **automatically based on the amplitude extrema**. For picking the extrema, drag the time window on the upper pick trace with the mouse. This generally results in very accurate phase picks but the absolute arrival time is delayed. Therefore, later picking on the beam trace is available. * To set picks **manually** at the arrival time click on the upper trace in the pick window. #. Pick multiple traces: Scroll though the traces using the mouse wheel or the :kbd:`arrow up` and :kbd:`arrow down` keys, #. Invert the picks for beam parameters, #. Set pick on the beam trace if required, #. Optionally: Select *New set* in the pull down Result selection menu to start another pick set and repeat picking and inversion, #. After inversion chose and review a pick set can be chosen in the pull down selection menu, #. Integrate the solutions in any of the possible ways: * Press the Apply button to apply the solution in the main |vappname| window, * Press Commit and Add to cart to provide the solution to the :ref:`cart`, * Press Commit and Send to send picks to the |scname| system, * Press Commit and check the box *Send as artifical origin* to send picks with an artificial origin to the |scname| system. Open the new origin, e.g. in :program:`scolv` for further process the solution. .. _lambda-fig-commit: .. figure:: /base/media/lambda/lambda-manual-send.png :width: 10cm :align: center Add picks to cart or send an artificial origin. .. _lambda-sec-slidingtw-fk: Slowness - wavenumber analysis ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Once the data are loaded and pre-processed the phases can be detected based on slowness / wavenumber analysis. The methods :ref:`lambda-sec-methods-fk` and :ref:`lambda-sec-methods-beamp` are available. #. In the Waveforms window select the type in the Type window along with the slowness / wavenumber ranges. Click on Apply finalize the selection. #. For choosing the slowness / wavenumber analysis, double-click on the button for **phase detection** in the Waveforms window. An new window will open which can be use simultaneously. #. For detection set the following parameters in the new window: * **Method**: Choose :ref:`F-K analysis ` or :ref:`beampbacking ` for detection. * **Unified sampling rate:** The sampling interval of the data should not be larger than the differential time shifts between pairs of stations. Therefore the data can be resampled to to a unified sampling rate. Upsampling and downsampling may be useful for small-aperture arrays ad large-aperture arrays, respectively. The sampling rate depends on the array geometry as well as on the slowness and the frequency of the expected signals. Upsampling increases the accuracy of beam packing but slows down the analysis. * **Maximum channel delay:** Time delay allowing to wait for complete data. * **Noise window:** Length of noise window used before detection to calculate the :term:`SNR` * **Processing time window:** Length of data processing window used for stacking. * **Processing step:** Time increment between 2 subsequent processing time windows. * **Semblance trigger thresholds:** Threshold of semblance to declare a detection. * **Minium SNR:** Minimum :term:`SNR` which must be reached at the position of the :term:`AIC` detection for declaring a pick. #. To start the analysis, click the start button. When the compuation is finished, the semblance and the beam traces will be shown along with the beam detections. Beam detections are plotted on top of semblance, slowness, backazimuth representing maximum values for each time step. #. Select the detection result for reviewing. The selected detection is highlighted in the slowness plot by changing the color. The corresponding beam parameters are shown to the right. The pick and waveform parameters update automatically when selecting another detection. The beam waveform is computed based on the beam parameters. The phase pick is found on the beam trace using the :term:`AIC` method. It is associated to the reference station of the array. #. Compare the detection diagram with the array response function (ATF) for identifying the minor maxima which may be caused by the ATF. For showing the ATF on top of the detection diagram, right click on the diagram and choose "Show ATF contours". The color of the circles on the contour indicate the semblance level of the contour. Right-click on the detection diagram for saving it to a file. .. _lambda-fig-fkbp: .. figure:: /base/media/lambda/lambda-sliding-FK-5.png :width: 12cm :align: center Interactive beam packing window. #. To send the pick to the :ref:`cart` or to the |scname| messaging click on "Commit". This will open the inspector window in which the pick parameters and the pick XML code can be inspected. After inspection click on the "Send" button to eventually send the parameters or cancel the operation. .. figure:: /base/media/lambda/lambda_sl-az-bp-processed-commit.png :width: 12cm :align: center Inspector window showing some details of the solution to be committed. .. figure:: /base/media/lambda/lambda_beampacking-processed_inspector_2.png :width: 12cm :align: center Inspector window showing the full list of objects to be committed in :term:`scml`. #. After sending the pick is available in the |scname| system and can be used for further analysis. E.g. it can be used for relocating events in scolv. As the beam pick is associated to the reference station of the array it can also be viewed there. The beam pick can be identified by the method **B**. .. figure:: /base/media/lambda/scolv_beampacking-pick.png :width: 12cm :align: center Beam pick in the scolv picker window. .. _lambda-sec-slidingtw-pmcc: PMCC ~~~~ PMCC is available for sliding time-window analysis. The :ref:`PMCC ` method is applied. #. Select PMCC in the lower trace window of the Waveform widget for sliding time-window analysis. #. Adjust the processing parameters as for :ref:`lambda-sec-slidingtw-fk`. #. Adjust the family criteria as well to join solutions belonging to the same phase. #. Press *Start* to calculate the parameters. #. Select a family. #. Apply or commit the solution of the selected family and send to :ref:`cart ` or to the |scname| messaging. .. _lambda-fig-pmcc: .. figure:: /base/media/lambda/lambda-sliding-pmcc-5.png :width: 12cm :align: center PMCC: Sliding time-window analysis. .. _lambda-sec-slidingtw-backproj: Backprojection ~~~~~~~~~~~~~~ |vappname| can also be used to detect and to locate an earthquake at the same time. The :ref:`back projection ` method is applied based on waveform stacking along pre-defined traveltime curves for given source points. #. **Array widget or Map window:** Select the array. #. **Location grid or map widget:** Select the range of the location grid by manually entering the values or by :ref:`selection on the map ` on the map. The grid resolution is adapted from the ATF settings widget. Finally, press *Apply*. .. figure:: /base/media/lambda/lambda_backporj-location-grid.png :width: 12cm :align: center Map window and Location grid widget: Define the location grid. #. **Waveform window:** * Choose the array and the waveform data * Select the traveltime interface and the traveltime table for calculating the traveltime from each grid point to each array element. * For starting the backprojection analysis, double-click on the *BPRO* button for **backprojection**. An new window Backprojection widget will open which can be use simultaneously. #. **Backprojection widget:** * For detection and location set the following parameters: * **Method**: Choose :ref:`F-K analysis ` or :ref:`beampbacking ` * **Unified sampling rate:** The sampling interval of the data should not be larger than the differential time shifts between pairs of stations. Therefore the data can be resampled to to a unified sampling rate. Upsampling and downsampling may be useful for small-aperture arrays ad large-aperture arrays, respectively. The sampling rate depends on the array geometry as well as on the slowness and the frequency of the expected signals. Upsampling increases the accuracy of beam packing but slows down the analysis. * **Maximum channel delay:** Time delay allowing to wait for complete data. * **Noise window:** Length of noise window used before detection to calculate the :term:`SNR` * **Processing time window:** Length of data processing window used for stacking. * **Processing step:** Time increment between 2 subsequent processing time windows. * **Depth:** Depth in kilometer at which the event is search for within the location grid. * **Semblance trigger thresholds:** Threshold of semblance to declare a detection. * **Minium SNR:** Minimum :term:`SNR` which must be reached at the position of the :term:`AIC` detection for declaring a pick. * To start the analysis, click the *Start* button. When the computation is finished, the semblance and the beam traces will be shown along with the beam detections. Beam detections are plotted on top of semblance, slowness, backazimuth representing maximum values for each time step. * Select the detection result for reviewing. The selected detection is highlighted in the slowness plot by changing the color. The corresponding origin and beam parameters are shown to the right. The parameters and waveforms update automatically when selecting another result. The beam waveform is computed based on the computed beam parameters for the reference station. The phase pick is found on the beam trace using the :term:`AIC` method. It is associated to the reference station of the array. .. figure:: /base/media/lambda/lambda_backporj-results.png :width: 12cm :align: center Processing window for beampacking with results. * To send the origin with the pick to the :ref:`cart ` or the |scname| messaging system click on "Commit". Committing opens the inspector window in which the origin parameters and the pick XML code can be inspected. After inspection click on the "Send" button to eventually send the parameters or cancel the operation. .. figure:: /base/media/lambda/lambda_backproj_origin.png :width: 10cm :align: center Evaluate and send origin. #. **Map Widget:** View the location grid diagram for evaluating the computed location on the map. For showing the semblance values and the origin location on top of the map in the *Map* widget, right click on the diagram and choose "Show on map". You may also save the semblance map to a file. .. figure:: /base/media/lambda/lambda_backproj_map.png :width: 12cm :align: center Show semblance grid on map. #. **Cart widget:** After sending the origin parameters are either available in * the :ref:`cart ` for further evaluation, * the |scname| system for further analysis. E.g. adding other picks and for relocating events in scolv. As the beam pick is associated to the reference station of the array it can also be viewed there. The beam pick can be identified by the method **B** shown along with the pick on the traces. Phase type identification ------------------------- **not yet implemented** .. _lambda-sec-slidingtw-vespagram: Vespagrams ---------- Vespagrams can be computed based on a selected solution after the slowness-backazimuth, PMCC or backprojection analysis. #. Select a solution #. For better visibility increase the row height of the vespgram above the traces of the processing widget: :kbd:`Shift` + :kbd:`y`, #. Keep or adjust the backazimuth and the frequency range, #. Adjust the vespagram settings: * time window, * slowness range, #. Press *Update vespagram* for computing the vespagram. .. figure:: /base/media/lambda/lambda_vespagram.png :width: 12cm :align: center Vespagram based on backazimuth from solution.