3. Example of how to run automated multi-event manager

This example shows how the code can be applied to undertake shear-wave splitting analysis for many events, with the data from the SKS example.

Here, the S-picks are defined as 10 s into the sac data for each event.

Data is from:

J Asplet, J Wookey, M Kendall. (2020) “A potential post-perovskite province in D″ beneath the Eastern Pacific: evidence from new analysis of discrepant SKS–SKKS shear-wave splitting”. GJI.

%load_ext autoreload
%autoreload 2

# Import modules:
import swspy
import obspy
from obspy import UTCDateTime
import numpy as np
%matplotlib notebook
import matplotlib.pyplot as plt
import glob
import os, sys
import pandas as pd

3.1. Specify parameters for processing:

First one specifies all the relevent parameters for the splitting analysis. This changes various parameters from their default values to something more appropriate for a particular dataset.

Note: For a detailed description of each parameter/attribute, do help(proc_many_events()) or read the documentation.

# Specify data management parameters:
datadir = "data"
outdir = "outputs"

# Setup automation object and set key splitting data processing parameters:
sws_proc_obj = swspy.automate.proc_many_events()
sws_proc_obj.filter = True
sws_proc_obj.filter_freq_min_max = [0.01, 0.5]
sws_proc_obj.overall_win_start_pre_fast_S_pick = 4.0
sws_proc_obj.win_S_pick_tolerance = 1.0
sws_proc_obj.overall_win_start_post_fast_S_pick = 30.0
sws_proc_obj.rotate_step_deg = 1.0
sws_proc_obj.max_t_shift_s = 5.0
sws_proc_obj.n_win = 10
sws_proc_obj.downsample_factor = 1 #4
sws_proc_obj.coord_system = "ZNE"
sws_proc_obj.sws_method = "EV"
sws_proc_obj.output_plots = False

3.2. Run the processing for multiple events:

After specifying the various parameters, one can run the analysis over multiple events.

# And run events through automated processing:
S_pick_time_after_start_s = 10.0 # Time of S pick after start of SAC trace
sws_proc_obj.run_events_sws_fmt(datadir, outdir, S_pick_time_after_start_s=S_pick_time_after_start_s)

Processing for event UID: 116A_2006360_122621_SKKS (1/11)
Saved sws result to: outputs/data/20061226125151_sws_result.csv
Saved sws wfs to: outputs/data/20061226125151_wfs_*.mseed
Processing for event UID: COR_2008321_170232_SKS (2/11)
Saved sws result to: outputs/data/20081116172713_sws_result.csv
Saved sws wfs to: outputs/data/20081116172713_wfs_*.mseed
Processing for event UID: DAN_2003174_121231_ScS (3/11)
Saved sws result to: outputs/data/20030623123100_sws_result.csv
Saved sws wfs to: outputs/data/20030623123100_wfs_*.mseed
Processing for event UID: FACU_2009297_144044_SKS (4/11)
Saved sws result to: outputs/data/20091024150503_sws_result.csv
Saved sws wfs to: outputs/data/20091024150503_wfs_*.mseed
Processing for event UID: HUMO_2008321_170232_SKS (5/11)
Saved sws result to: outputs/data/20081116172718_sws_result.csv
Saved sws wfs to: outputs/data/20081116172718_wfs_*.mseed
Processing for event UID: IRON_2009297_144044_SKS (6/11)
Saved sws result to: outputs/data/20091024150513_sws_result.csv
Saved sws wfs to: outputs/data/20091024150513_wfs_*.mseed
Processing for event UID: K20A_2009003_223342_SKKS (7/11)
Saved sws result to: outputs/data/20090103225943_sws_result.csv
Saved sws wfs to: outputs/data/20090103225943_wfs_*.mseed
Processing for event UID: L07A_2007256_094844_SKS (8/11)
Saved sws result to: outputs/data/20070913101323_sws_result.csv
Saved sws wfs to: outputs/data/20070913101323_wfs_*.mseed
Processing for event UID: L24A_2009003_194355_SKKS (9/11)
Saved sws result to: outputs/data/20090103201022_sws_result.csv
Saved sws wfs to: outputs/data/20090103201022_wfs_*.mseed
Processing for event UID: NE81_2006360_122621_SKKS (10/11)
Saved sws result to: outputs/data/20061226125216_sws_result.csv
Saved sws wfs to: outputs/data/20061226125216_wfs_*.mseed
Processing for event UID: RDM_2003174_121231_ScS (11/11)
Saved sws result to: outputs/data/20030623123110_sws_result.csv
Saved sws wfs to: outputs/data/20030623123110_wfs_*.mseed
Finished processing shear-wave splitting for data in: data
Data saved to: outputs